Nucleic Acids Research Advance Access originally published online on December 7, 2006
Nucleic Acids Research 2007 35(1):193-202; doi:10.1093/nar/gkl1059
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Nucleic Acids Research, 2007, Vol. 35, No. 1 193-202
© 2006 The Author(s)
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.
Genomics |
Genome-wide analysis of transcriptional dependence and probable target sites for Abf1 and Rap1 in Saccharomyces cerevisiae
1 Laboratory of Developmental Genetics, Wadsworth Center, New York State Department of Health Albany, NY 12201-2002, USA 2 Department of Biomedical Sciences, State University of New York at Albany School of Public Health Albany, NY 12201-2002, USA
*To whom correspondence should be addressed. Tel: +1 518 486 3116; Fax: +1 518 474 3181; Email: Randall.Morse{at}wadsworth.org
Received September 25, 2006. Revised October 19, 2006. Accepted November 7, 2006.
| ABSTRACT |
|---|
|
|
|---|
Abf1 and Rap1 are general regulatory factors (GRFs) that contribute to transcriptional activation of a large number of genes, as well as to replication, silencing and telomere structure in yeast. In spite of their widespread roles in transcription, the scope of their functional targets genome-wide has not been previously determined. Here, we use microarrays to examine the contribution of these essential GRFs to transcription genome-wide, by using ts mutants that dissociate from their binding sites at 37°C. We then combine this data with published ChIP-chip studies and motif analysis to identify probable direct targets for Abf1 and Rap1. We also identify a substantial number of genes likely to bind Rap1 or Abf1, but not affected by loss of GRF binding. Interestingly, the results strongly suggest that Rap1 can contribute to gene activation from farther upstream than can Abf1. Also, consistent with previous work, more genes that bind Abf1 are unaffected by loss of binding than those that bind Rap1. Finally, we show for several such genes that the Abf1 C-terminal region, which contains the putative activation domain, is not needed to confer this peculiar memory effect that allows continued transcription after loss of Abf1 binding.
| INTRODUCTION |
|---|
|
|
|---|
Abf1 and Rap1 are multipurpose, abundant, site-specific DNA-binding proteins that function in transcriptional activation, silencing and replication in the budding yeast Saccharomyces cerevisiae (13). Rap1 and Abf1 are essential for cell viability, and their multifunctional character has caused them to be termed General Regulatory Factors (GRFs). Binding sites for Abf1 and Rap1 are found in a large number of promoters, and have been shown to be important for activation of many of the corresponding genes. Genome-wide localization, or Chromatin immunoprecipitation (ChIP)-on-chip experiments indicate that Rap1 and Abf1 bind
200300 promoters each, and Rap1 has been found to bind with 122 out of 137 ribosomal protein (RP) genes (46). Genes containing Rap1 or Abf1 sites principally comprise RP genes and genes encoding proteins involved in amino acid biosynthesis, regulation of carbon source and sporulation. Many of these genes are among those most highly expressed in the yeast genome. The broad scope of promoters bounded by Abf1 and Rap1, and the importance of such promoters in basic cellular metabolism, underscore the key role that these GRFs play in gene control in yeast. Consistent with this notion, genome-wide interaction studies have found both Abf1 and Rap1 to be network hubs, with the implication that they play central and coordinating roles in cell function (7,8). Thus, identifying functional targets of Abf1 and Rap1 is important in constructing a circuit diagram of gene control in yeast. To this end, we report here the effect of loss of Abf1 and loss of Rap1 binding on genome-wide transcription. By combining our microarray results with genome-wide location results and sequence analyses, we identify probable direct targets of Rap1 and Abf1. Our results also indicate that many promoters that are controlled by Abf1 do not require continued binding by Abf1 for ongoing transcription, consistent with investigations of individual promoters (9,10). In contrast, loss of Rap1 at Rap1 binding promoters more commonly leads to significantly decreased transcription. Furthermore, we find that Rap1 activates genes from more distant sites than does Abf1, suggesting a previously unrecognized mechanistic difference between these two GRFs. Taken together, the results presented here provide new insight into the roles of Abf1 and Rap1 in genome-wide transcriptional regulation.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Plasmids
The expression vector for abf1(1592) ts was constructed by replacing the small NheIXhoI fragment from pRS315abf1-1 with the corresponding fragment from pRS315Abf1(1592), and the resulting plasmid shuttled into TMY86 as described previously (2).
Yeast strains and growth
The S.cerevisiae strains used in this study are listed in Table 1; all are derived from W303a. Transformations were done using a standard lithium acetate protocol (11), and yeast cells were grown at 30°C in complete synthetic medium (CSM) [6.7 g/l yeast nitrogen base without amino acids, 2% glucose and CSM dropout mixture (Bio101)] when selection was required, or in rich medium (YPD) (1% Bactoyeast extract, 2% Bactopeptone and 2% glucose). For temperature shift experiments, yeast were grown to mid-log phase in YPD medium at 2325°C and an equal volume of medium at 50°C was added to bring the culture rapidly to 37°C; growth was then continued for 12 h (as indicated in the text) at 37°C.
|
Strain RMY28
a, in which the Abf1 binding site in the promoter of RPS28A is mutated (replacing CGTCTAGAGTGAC with CATGG), was constructed by two-step replacement in GA208 using the URA3 gene (12). The modified RPS28A promoter from the resulting strain was verified by PCR amplification and sequencing of the promoter region.
ChIP and RNA analysis
ChIP and northern analysis were performed as described previously (9). For FLAG IP, 1.2 µl (6 µg) of antibody (Sigma) were used per IP reaction. Probes were prepared by PCR. For analysis of RPS28A mRNA, purified RNA was reverse transcribed with primers for RPS28A (GACGAGCTTCACGTTCAGATTCCATTAG) and PYK1 (TTTCGTGGTTTGGTGGGATT), using Bioscript reverse transcriptase (Bioline USA Inc., Randolph, MA). The resulting cDNA was amplified by PCR for 1520 cycles, after pilot experiments to determine the range over which exponential amplification was observed. For quantification, a low number of cycles (usually
15) was used, and the products separated by gel electrophoresis and analyzed by Southern blotting and Phosphorimager analysis, as described previously (9). A control experiment using an rps28a
yeast strain yielded no product using this protocol.
Microarray analysis
RNA was prepared from exponentially growing yeast cells (A600 = 0.8 1.2) grown in YPD using the Masterpure yeast RNA purification kit (Epicenter Technology, Madison, WI). RNA was further purified using the RNeasy purification kit (Qiagen). Processing and hybridization using Affymetrix S98 microarrays (Affymetrix, Santa Clara, CA) were performed as per manufacturer's instructions. Changes in gene expression were calculated by averaging log2 expression changes, and false discovery rates (FDRs) (13) were calculated in Genespring, after normalizing within each experiment (i.e. wild-type and ts mutant analyzed in parallel). Comparative analyses were done using Excel (Microsoft) and sequence analysis was performed using MEME (http://meme.sdsc.edu/meme/) and regulatory site analysis tools) (RSAT) (http://rsat.ccb.sickkids.ca/) (14,15). Sequence logos (16) were generated from position specific weight matrices (derived using RSAT) using the webserver at http://weblogo.berkeley.edu/ (17).
Microarray accession number
Microarray gene expression data are available at the Gene Expression Omnibus under accession no. GSE6073
[NCBI GEO]
, and are also available at www.wadsworth.org/resnres/bios/morse.htm and as Supplementary Data.
| RESULTS |
|---|
|
|
|---|
Genome-wide effects of loss of Abf1 or Rap1 binding
To assess the genome-wide effect of loss of Abf1 binding on transcription, we performed microarray analysis of mRNA from yeast harboring the abf1-1 ts mutation as well as from a congenic wild-type yeast strain after growth in rich medium (YPD) for 1 h at 37°C. RNA samples from three replicates were analyzed using Affymetrix microarrays. Similarly, genome-wide expression of wild-type and congenic rap1-2ts yeast incubated for 1 h at 37°C was also compared by microarray analysis. Neither abf1-1 nor rap1-2 yeast exhibit phenotypes that are characteristic of specific cell cycle arrest at the restrictive temperature (18,19). Consistent with this, no enrichment of 807 genes that exhibit cell cycle dependent regulation (20) was observed among genes showing altered transcription in abf1-1 yeast compared to wild-type at 37°C (data not shown). A modest enrichment was observed among genes showing decreased transcription in rap1 ts yeast compared to wild-type yeast at 37°C, with 76 out of 807 cell cycle regulated genes being among the 100 genes showing most strongly decreased transcription. However, this enrichment evidently did not reflect an altered cell cycle, as these genes showed dispersed peak expression patterns in the cell cycle, with peaks in S, G1/S, G1 and G2/M (20) (data not shown). Furthermore, the large majority of genes regulated by the cell cycle were not enriched among those showing altered expression in rap1 ts yeast compared to wild-type. Thus, cell cycle effects do not represent a major contribution to the observed changes in gene expression in abf1-1 or rap1 ts yeast.
Analysis of the microarray data indicates that global gene expression is affected much more by the loss of Rap1 binding than by loss of Abf1 binding (Figure 1). We found 235 genes (215 after removing cell cycle regulated genes) having changed expression upon loss of Abf1 binding with FDR <0.05, compared to 947 genes (833 after removing cell cycle regulated genes) with FDR <0.05 upon loss of Rap1 binding. Similarly, >10 times as many genes within these sets showed changed expression by 2-fold upon loss of Rap1 binding than upon loss of Abf1 binding. These results are not due to more genes having Rap1 than Abf1 binding sites, as genome-wide location analysis (ChIP-chip analysis) indicates a similar number of promoters binding Rap1 and Abf1 under our growth conditions (46) as well as in vitro (21). Another possible explanation is that the quality of the data for the Rap1 experiments is better than for the Abf1 experiments, resulting in a higher fraction of true positives having low FDRs for the former. However, even ignoring FDR values, we still observe considerably more genes showing changed expression upon loss of Rap1 than upon loss of Abf1 (Figure 1B). Analysis of this data to identify likely direct targets of Abf1 and Rap1 indicates that the larger effect genome-wide of loss of Rap1 compared to loss of Abf1 binding arises in part from Rap1 exerting a larger indirect effect on gene expression, and in part from fewer Abf1 binding than Rap1 binding targets being affected by loss of factor binding (see below).
|
We used the Funspec and T-profiler web tools to examine functional categories, promoter motifs and transcription factor binding sites derived from ChIP-chip data that were enriched in the gene sets responsive to loss of Abf1 or Rap1 binding (Figure 1C) (22,23). Funspec assigns genes in a specified list (in this case, genes most down-regulated upon loss of Abf1 or Rap1 binding) to functional categories and determines P-values for enriched categories, using Fisher's exact test (i.e. based on a hypergeometric distribution). T-profiler models gene activity as arising from a linear combination of transcription factor binding sites (consensus motifs) or transcription factors (ChIP-chip data) at each promoter, and models the entire set of microarray data on this basis, yielding an output of motifs or factors showing greatest contribution to gene expression changes. In both cases, the Abf1 data were enriched in mitochondrial ribosome biogenesis and the Rap1 data in ribosome biogenesis and function, and amino acid biosynthesis and metabolism. These represent known functions of Abf1 and Rap1 [(9) and references therein] and, therefore, this result indicates that gene expression changes measured by microarray analysis upon loss of binding accurately reflects the contribution of these GRFs to gene regulation.
In addition, the only consensus motif predicted to contribute to the change in expression in the Abf1 data set was that for Abf1, and Abf1 was one of only a few factors identified as contributing on the basis of ChIP-chip data (Figure 1C). In contrast, Rap1 was one of several factors identified as contributing to expression in the Rap1 dataset, with other identified factors (Met4, Fkh1, Met31, Gcn4, Fkh2 and Fhl1) or motifs (Pho4 and Gcr1) reflecting the known role of Rap1 in ribosomal biogenesis and amino acid biosynthesis (3,24). Several of these factors are indeed known to operate in conjunction with Rap1 in gene expression (2529), and it is possible that expression of some genes that bind Rap1 but do not suffer significant loss in expression upon loss of Rap1 binding is maintained by these or other auxiliary transcription factors. However, the finding of few such factors or motifs contributing to expression of the Abf1 data set argues against the idea that a small set of auxiliary factors that bind together with Abf1 at gene promoters is responsible for the continued expression of Abf1 binding genes upon loss of Abf1 binding. Interestingly, no Rap1 consensus motif emerged from T-profiler analysis, possibly reflecting the wide variety of sequences capable of binding to Rap1 (30). Alternatively, the Rap1 consensus motif may have been obscured by noise if a large proportion of affected genes were not direct targets of Rap1 (see below).
|
Finally, we note that we did not observe any significant telomeric clustering of genes showing increased expression upon loss of Rap1 or Abf1 binding using Pyxis, which employs nearest neighbor probability calculations to determine significant gene clustering (31) (data not shown). Although such an effect might have been predicted, based on the well known role of Rap1 in telomeric silencing (32), it has also been shown that native yeast telomeres do not exhibit gene silencing properties identical with those of the truncated telomeres used in many telomeric silencing studies (33). It is also possible that subtelomeric heterochromatin remains stable following the loss of Rap1 binding and therefore, exerts a continued repressive effect.
Identification of direct targets of Rap1 and Abf1
To identify genes likely to be directly controlled by Rap1, we first selected those genes identified as showing decreased expression by >1.5x (log2 change less than 0.59) in rap1 ts yeast compared to wild-type at 37°C (averages of three replicate microarray experiments) and binding Rap1 with P < 0.001 in asynchronous yeast grown in rich medium (YPD), the same conditions used in our microarray analyses [(6) and Figure 2]. Thirty-four such genes were found, 3 of which shared common promoters, leaving 31 unique promoters. These 31 promoters were analyzed for common sequence motifs using MEME (15) and two motifs identified. One of these was a T-rich motif, and the other was a clear match to previously reported Rap1 binding sites [(6,21,30) and Figure 2B]. We next used the position specific weight matrix (PSWM) representing the identified Rap1 binding site to search for Rap1 binding sites in these 31 promoters, and found 26 with significant matches to the Rap1 binding sequence. Using a protein binding microarray (PBM) protocol, 23 of these promoters were found to bind Rap1 in vitro (21), further supporting their assignment as highly probable direct targets of Rap1. Impressively, 18 of these 23 promoters contain Rap1 binding sites that are conserved among five closely related S.cerevisiae species, including 11 of 14 RP gene promoters (21). It seems likely in this case that the criterion of orthologous conservation employed is too stringent and excludes valid targets of Rap1. In addition, 2 of our original 31 promoters, CPA1 and ERG5, lacking a Rap1 binding site or having a marginal sequence match, bind Rap1 in the PBM assay but again are not conserved among related Saccharomyces species.
We similarly analyzed an additional 21 promoters showing decreased transcription by 2x or greater in rap1 ts compared to wild-type yeast at 37°C, and identified as binding Rap1 by Lieb et al. but not by Harbison et al. (5,6), using the PSWM defined above. (We used a slightly more stringent criterion for decreased transcription2x instead of 1.5xbecause P-values were not available as a quantitative criterion in the data of Lieb et al.) Twelve were found to have significant matches to the Rap1 binding sequence, and eight of these bound Rap1 in the in vitro binding PBM assay. In addition, 3 of the original 21 promoters bound Rap1 in vitro but did not have strong matches to the consensus site. Taking all these data together, we assigned 43 genes as probable direct targets of Rap1. These genes, together with their expression levels (34) and the location and sequence of the identified binding site(s), are listed in Supplementary Table 1. As expected, many of these genes are RP genes, and many are highly expressed. Interestingly, 27 of these 43 genes contain Rap1 sites >300 bp from the starting ATG, indicating that Rap1 exerts its activating effect on transcription over a relatively long distance compared to many yeast activators (Figure 3A). It is also notable that the large majority of highly expressed Rap1 targets identified here have their Rap1 binding site in the same orientation (as shown in Figure 2B).
|
We next identified additional putative targets of Rap1. These were defined as genes that (i) show decreased transcription by >2x in rap1 ts compared to wild-type yeast, and (ii) have a signficant match to the Rap1 binding sequence, as defined by its PSWM from above, or bind to Rap1 in vitro (PBM assay), but (iii) do not bind Rap1 in ChIP-chip studies. Of 391 unique promoters that showed at least 2-fold decreased transcription in rap1 ts compared to wild-type yeast but did not bind Rap1 in ChIP-chip studies (5,6), 30 contained sequences matching the Rap1 binding sequence and 16 of these bind Rap1 in vitro (PBM assay). An additional 15 promoters do not contain significant matches to the Rap1 binding site, but do bind Rap1 in vitro. Together, we categorize these 45 promoters as putative direct functional targets of Rap1 (Supplementary Table 2). Since these promoters derived from a set of 391 promoters showing 2-fold decreased transcription upon loss of Rap1 binding, this result suggests that a substantial fraction of genes affected upon loss of Rap1 binding are affected indirectly.
We also identified genes whose promoters are likely to bind Rap1 but that are unaffected upon loss of Rap1 binding. We selected 37 unique promoters for which the log2 change fell between 0.2 and 0.2, and which were identified as binding Rap1 with P < 0.001 (6). We eliminated from consideration four promoters for which a divergent transcript showed a decrease of at least 1.5-fold in rap1 ts yeast at 37°C. Significant matches to the Rap1 binding motif (Figure 2B) were found in 29 of the remaining 33 promoters, and 27 of these 29 also bind Rap1 in vitro (21). In addition, two sequences apparently lacking significant Rap1 binding sites, RPS26A and RPL33B, were found to bind Rap1 in vitro (21), yielding 31 genes identified as binding Rap1 but not affected by loss of Rap1 binding (Supplementary Table 3).
A similar analysis for Abf1 revealed that of 38 unique promoters showing decreased transcription (>1.5x) and binding Abf1 in ChIP-chip analysis (P < 0.001) (6), 37 had strong matches to the Abf1 motif identified from this set (Figure 2C and Supplementary Table 4). As with Rap1, the identified Abf1 motif is consistent with earlier work (6,21,35). Surprisingly, only 10 of these promoters were found to bind Abf1 in vitro using the PBM assay; eight of these promoters contained sites conserved in related Saccharomyces species (21). Furthermore, in contrast to promoters targeted by Rap1, the large majority (32 out of 37) of these promoters contained Abf1 binding sites within 300 bp, and most typically (29 out of 37) between 100 and 200 bp from the starting ATG, indicating that Abf1 may not activate transcription from distant locations as efficiently as Rap1 (Figure 3B). In addition, we found strong Abf1 binding sites in 37 of 66 unique promoters showing decreased transcription by 2x or greater but not identified as binding Abf1 in ChIP-chip studies (6); an additional two of these promoters lacking strong Abf1 binding sites were found to bind Abf1 in vitro (PBM) (21). These 39 genes are assigned as putative targets of Abf1 (Supplementary Table 5). Lastly, we identified genes whose promoters are likely to bind Abf1 but which are unaffected upon loss of Abf1 binding. We examined 92 genes for which the log2 change fell between 0.2 and 0.2, and which were identified as binding Abf1 with P < 0.001 (6). After removing from consideration five promoters associated with divergent genes that showed decreased (by >1.5-fold) transcription in abf1-1 ts yeast at 37°C, significant matches to the Abf1 PSWM (Figure 2C) were found in 79 of the remaining promoters, of which 47 bind Abf1 in vitro (21). In comparison, searching a set of 49 promoters binding Rap1 for matches to the Abf1 binding sequence yielded only three signficant matches. One additional gene that did not yield a significant match to the consensus binding site was found to bind Abf1 in vitro (21), giving a total of 80 genes that are likely to bind Abf1 but not be affected upon loss of Abf1 binding (Supplementary Table 6).
Expression of the RPS28A gene depends on an Abf1 binding site but is not affected by loss of Abf1 binding
Our microarray analysis indicated that only four genes showed >2-fold increased expression in abf1-1 ts yeast at 37°C compared to wild-type, and one of these is the ABF1 gene itself, consistent with results of Miyake et al. (36). This suggests that a possible explanation for the modest effect on global gene expression in abf1-1 ts yeast at 37°C is that the decreased affinity of the mutant Abf1 protein is offset by its increased abundance. To examine this possibility, we used ChIP to test Abf1-FLAG binding in the ts mutant compared to wild-type at 37°C. We observed substantially reduced Abf1 binding in the abf1-1 ts mutant compared to wild-type at 37°C at the RPS28A and SPT15 loci (Figure 4A and data not shown), consistent with our earlier results examining the HIS4-Abf1-MEL1 reporter gene, as well as with results of Miyake et al. (9,36). Expression of both RPS28A and SPT15 were found to be unaffected in abf1-1 compared to wild-type yeast at 37°C in our microarray experiments, and we confirmed this for the individual transcripts (Figure 4B and data not shown). Expression of SPT15, which encodes TATA-binding protein (TBP), was also shown by Schroeder and Weil to be unaffected by loss of Abf1 binding, although SPT15 promoter activity depends strongly on the Abf1 binding site (10,37). Deletion of the Abf1 binding site from the RPS28A promoter driving a GUS reporter gene resulted in about a 10-fold decrease in expression (38). We found that replacement of the Abf1 binding site with a mutated site in the chromosomal RPS28A gene resulted in 2- to 4-fold decreased expression (Figure 4B and data not shown). Thus, the RPS28A gene, like SPT15 and the modified HIS4 promoter (9), depends on the Abf1 binding site for normal levels of transcription, and yet loss of Abf1 binding from this site does not reduce its expression.
|
The C-terminal region of Abf1 is not required for continued transcription of Abf1-dependent genes after loss of Abf1 binding
One way in which an Abf1 binding site might contribute to gene activation without requiring continued occupancy by Abf1 would be if Abf1 recruited another protein that modified chromatin, or that remained bound to DNA even after loss of Abf1 binding. Since the C-terminal region of Abf1 can activate transcription and remodel chromatin in the context of protein fusions, this region is a logical candidate for effecting such hypothesized recruitment. We therefore tested whether the C-terminal region of Abf1 was necessary for ongoing transcription after loss of Abf1 binding.
We first introduced the abf1-1 point mutation into the abf1(1592) expression plasmid. When this plasmid was shuttled into yeast to replace the wild-type ABF1 gene, the resulting strain grew slowly and was temperature sensitive (Figure 5A and data not shown). We next examined transcripts of several genes that are known to depend on Abf1 binding sites for full expression in yeast expressing full-length Abf1, the ts mutant abf1-1, Abf1(1592), which lacks the C-terminal region, or the corresponding ts mutant from yeast grown at 2325°C or 1 h following a rapid shift to 37°C (Figure 5B). PRO3 showed a marked drop in transcript level in abf1-1 yeast at the restrictive temperature in the work of Miyake et al. (36) as well as in our microarray experiments, and so was examined as a positive control. Northern analysis showed that PRO3 mRNA levels decreased considerably after 1 h at 37°C in the ts mutants both in the background of full-length and the C-terminal deletion of Abf1, but not in the corresponding wild-type strains (Figure 5B). In contrast, neither QCR8 nor TCM1/RPL3, both of which depend on Abf1 binding sites for transcription (39,40), showed significantly decreased transcript levels at 37°C in either abf1-1 mutant, consistent with previous work (10), nor in the abf1(1592) ts mutant (Figure 5B). Similar results were also seen with SPT15 and RPS28A (data not shown), both of which rely on Abf1 binding sites for their expression [(10) and Figure 3B]. Furthermore, only PRO3 showed significantly reduced expression in yeast expressing Abf1(1592) compared to full-length Abf1. We conclude that whatever property Abf1 possesses that allows continued gene expression following loss of Abf1 binding does not reside exclusively in its C-terminal region.
|
| DISCUSSION |
|---|
|
|
|---|
We report here the effects of loss of Abf1 binding and loss of Rap1 binding on global gene expression. As both of these factors are essential, gene expression was compared in wild-type and rap1 ts or abf1-1 ts mutants after 1 h growth at 37°C following a rapid temperature shift. Previous work has shown by DMS footprinting and chromatin IP (for Abf1) that both mutant proteins rapidly dissociate from their DNA targets at this temperature in vitro and in vivo (9,10,18,19,36,41). Genes showing decreased expression upon loss of Abf1 or Rap1 binding were enriched in the expected functional categories, and were also enriched in Abf1 or Rap1 binding genes (although only slightly for the latter; see below) (Figure 1), attesting to the quality of the microarray data. In addition, our microarray results on gene expression changes in the abf1-1 mutant compared to wild-type yeast agreed well with previous work; 39 of the 51 (P < 1040) genes identified previously as showing decreased expression in abf1-1 yeast at the restrictive temperature also showed decreased expression by at least 1.5x in our study. However, our study differed from the previous work in that it included three biological replicates of both wild-type and mutant yeast, whereas the earlier work used three technical replicates from single mRNA isolations. The data in the current study therefore provides a better basis for statistical analysis of the microarray results. Regarding Rap1, although a previous study used microarrays to examine gene expression of yeast lacking the C-terminal silencing domain (42), to our knowledge this is the first work to examine the effect of loss of Rap1 binding on global gene expression.
To identify direct functional targets of Abf1 and Rap1, we used our microarray results together with motif analysis and published data on genome-wide in vivo and in vitro binding (5,6,21). This analysis led us to identify highly probable and putative direct targets for Abf1 and Rap1. Similar numbers of direct targets were found for both of these GRFs (Figure 2D), in agreement with the finding that roughly similar numbers of promoters bind both Abf1 and Rap1 in vitro and in vivo (6,21). However, a considerably smaller fraction of genes showing decreased transcription upon loss of Rap1 binding were identified as direct targets than of those showing decreased transcription upon loss of Abf1 binding (
9% compared to 29%). The reason for Rap1 apparently exerting a larger indirect effect than Abf1 on genome-wide transcription is not clear, but it does provide a possible explanation for the Rap1 consensus not emerging in the T-profiler analysis (Figure 1), as it may have been obscured by the greater amount of sequence noise caused by indirect effects on transcription. Interestingly, this analysis also yielded evidence that Rap1 binding sites typically function from a considerably greater distance from the transcription start site than do Abf1 sites; the mechanistic basis for this is completely unknown.
We also identified promoters likely to bind Abf1 or Rap1 that are not affected upon loss of binding (Supplementary Tables 3 and 6). Some of the Abf1 binding promoters in this class are probably not under control of Abf1 because the sites are too distant from the proximal promoter (Figure 3B), but this is probably not the case for Rap1 binding promoters, as Rap1 appears capable of activating promoters from sites far upstream (Figure 3A). It is also possible that some promoters bind Abf1 or Rap1 but are nevertheless not controlled by these GRFs, perhaps because other factors binding these promoters dominate their activation. However, such hypothetical promoters have not to our knowledge yet been verified (by mutagenesis of a validated binding site resulting in no change in transcription levels). Furthermore, there is no evident correlation between gene expression level and the effects of loss of Abf1 or Rap1 binding (Supplementary Tables 16). In contrast, several promoters that depend on Abf1 binding sites for fully activated levels of transcription have been shown not to yield decreased mRNA levels upon loss of Abf1 binding, including RPS28A (Figure 4) (10). What the mechanism is behind this memory effect [see references (9,10) for discussion] remains a topic for future study.
Several Abf1 binding sites that have been shown to function in transcriptional control and to bind Abf1 in vivo by DMS footprinting (SPT15, QCR8 and TCM1/RPL) (10) do not show statistically signficant Abf1 binding in ChIP-on-chip experiments, nor were their promoters found to bind Abf1 in vitro (6,21). Thus, these genes were not identified as having Abf1 binding promoters in our analysis, indicating that the criteria we used were perhaps overly stringent and likely to lead to a conservative estimate of Abf1 and Rap1 binding promoters. Indeed, transcription factor binding and contribution to gene expression may better be represented as a continuum than as a discrete On/Off function (43). Much additional work will be needed to achieve a quantitative model of genome-wide transcription factor binding and contribution to gene expression for individual transcription factors.
The persistent transcription induced by Abf1 does not require the protein's C-terminal region, which contains the putative activation domain, as we found the effect undiminished in yeast expressing Abf1 that lacks this region [Figure 5 and (2)]. Interestingly, expression of some but not all Abf1 binding genes is affected by loss of this region, suggesting that Abf1 may activate transcription by more than one mechanism [Figure 5 and (9,36)]. Our analysis (using MEME) of Abf1 binding promoters did not yield any motifs that distinguished those affected by loss of Abf1 binding from those that were not.
At present we do not know whether any Rap1-dependent promoters may also show persistent transcription following loss of Rap1 binding. Fewer promoters that are likely to bind Rap1 but that are not affected by loss of Rap1 binding were identified than for Abf1, and none have thus far been verified by mutagenesis of Rap1 binding sites, as has been done for Abf1 (3740). Promoters known to depend on Rap1 binding sites, such as HIS4 and many RP gene promoters (3,28) show decreased transcription in our microarray experiments, suggesting that at least most promoters that depend on Rap1 binding sites for full expression also require continued Rap1 binding. A more detailed investigation will be required to address this issue.
In summary, we have examined the effect on global gene expression of loss of binding of two GRFs, Abf1 and Rap1. Both of these GRFs bind to on the order of 200300 promoters; of these, we have identified 76 and 88 promoters likely to be direct functional targets of Abf1 and Rap1, respectively. Since both Abf1 and Rap1 have been identified as hubs in genome-wide interaction networks (7,8), these results should be useful in construction of circuit diagrams of the yeast transcriptome. We have also identified 80 promoters that are likely to bind Abf1 but are unaffected by loss of binding, and 31 such promoters for Rap1. At least some, and perhaps many, of the promoters that bind Abf1 but whose transcription is unaffected by loss of binding nonetheless depend on an Abf1 binding site (3740). The mechanism by which Abf1 marks such promoters remains enigmatic, but the phenomenon is reminiscent of regulation of many of the genes responsible for pattern formation in Drosophila, for which repressed or active states are maintained by PcG or Trx complexes even after the initial transcriptional repressors or activators have dissipated (44). It seems likely that further investigation into the mechanism by which Rap1 and Abf1 contribute to gene activation will provide new insights into fundamental aspects of transcriptional regulation.
| SUPPLEMENTARY DATA |
|---|
|
|
|---|
Supplementary Data are available at NAR online.
| ACKNOWLEDGEMENTS |
|---|
The authors gratefully acknowledge assistance from the Wadsworth Center Microarray and Molecular Genetics Cores. The authors also thank Tsuyoshi Miyake, Rong Li (University of Virginia Health Sciences Center) for generously providing yeast strains. This work was supported by grant RO1 GMS51993 from the National Institutes of Health. Funding to pay the Open Access publication charges for this article was provided by state funds appropriated by the New York State Department of Health.
Conflict of interest statement. None declared.
| REFERENCES |
|---|
|
|
|---|
- Morse, R.H. (2000) RAP, RAP, open up! new wrinkles for RAP1 in yeast Trends Genet, . 16, 5153[CrossRef][Web of Science][Medline] .
- Miyake, T., Loch, C.M., Li, R. (2002) Identification of a multifunctional domain in autonomously replicating sequence-binding factor 1 required for transcriptional activation, DNA replication, and gene silencing Mol. Cell. Biol, . 22, 505516
[Abstract/Free Full Text] . - Shore, D. (1994) RAP1: a protean regulator in yeast Trends Genet, . 10, 408412[CrossRef][Web of Science][Medline] .
- Lee, T.I., Rinaldi, N.J., Robert, F., Odom, D.T., Bar-Joseph, Z., Gerber, G.K., Hannett, N.M., Harbison, C.T., Thompson, C.M., Simon, I., et al. (2002) Transcriptional regulatory networks in Saccharomyces cerevisiae Science, 298, 799804
[Abstract/Free Full Text] . - Lieb, J.D., Liu, X., Botstein, D., Brown, P.O. (2001) Promoter-specific binding of Rap1 revealed by genome-wide maps of proteinDNA association Nature Genet, . 28, 327334[CrossRef][Web of Science][Medline] .
- Harbison, C.T., Gordon, D.B., Lee, T.I., Rinaldi, N.J., Macisaac, K.D., Danford, T.W., Hannett, N.M., Tagne, J.B., Reynolds, D.B., Yoo, J., et al. (2004) Transcriptional regulatory code of a eukaryotic genome Nature, 431, 99104[CrossRef][Medline] .
- Zhang, Z., Liu, C., Skogerbo, G., Zhu, X., Lu, H., Chen, L., Shi, B., Zhang, Y., Wang, J., Wu, T., et al. (2006) Dynamic changes in subgraph preference profiles of crucial transcription factors PLoS Comput. Biol, . 2, e47[CrossRef][Medline] .
- Luscombe, N.M., Babu, M.M., Yu, H., Snyder, M., Teichmann, S.A., Gerstein, M. (2004) Genomic analysis of regulatory network dynamics reveals large topological changes Nature, 431, 308312[CrossRef][Medline] .
- Yarragudi, A., Miyake, T., Li, R., Morse, R.H. (2004) Comparison of ABF1 and RAP1 in chromatin opening and transactivator potentiation in the budding yeast Saccharomyces cerevisiae Mol. Cell. Biol, . 24, 91529164
[Abstract/Free Full Text] . - Schroeder, S.C. and Weil, P.A. (1998) Genetic tests of the role of Abf1p in driving transcription of the yeast TATA box binding protein-encoding gene, SPT15 J. Biol. Chem, . 273, 1988419891
[Abstract/Free Full Text] . - Hill, J., Donald, K.A., Griffiths, D.E., Donald, G. (1991) DMSO-enhanced whole cell yeast transformation Nucleic Acids Res, . 19, 5791 [published erratum appears in Nucleic Acids Res 1991 Dec 11;19(23):6688]
[Free Full Text] . - Johnston, M., Riles, L., Hegemann, J.H. (2002) Gene disruption Meth. Enzymol, . 350, 290315[Web of Science][Medline] .
- Benjamini, Y. and Hochberg, Y. (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing J. Roy. Statist. Soc. B, 57, 289300 .
- van Helden, J. (2003) Regulatory sequence analysis tools Nucleic Acids Res, . 31, 35933596
[Abstract/Free Full Text] . - Bailey, T.L. and Elkan, C. (1994) Fitting a mixture model by expectation maximization to discover motifs in biopolymers Proc. Int. Conf. Intell. Syst. Mol. Biol, . 2, 2836[Medline] .
- Schneider, T.D. and Stephens, R.M. (1990) Sequence logos: a new way to display consensus sequences Nucleic Acids Res, . 18, 60976100
[Abstract/Free Full Text] . - Crooks, G.E., Hon, G., Chandonia, J.M., Brenner, S.E. (2004) WebLogo: a sequence logo generator Genome Res, . 14, 11881190
[Abstract/Free Full Text] . - Kurtz, S. and Shore, D. (1991) RAP1 protein activates and silences transcription of mating-type genes in yeast Genes Dev, . 5, 616628
[Abstract/Free Full Text] . - Rhode, P.R., Elsasser, S., Campbell, J.L. (1992) Role of multifunctional autonomously replicating sequence binding factor 1 in the initiation of DNA replication and transcriptional control in Saccharomyces cerevisiae Mol. Cell. Biol, . 12, 10641077
[Abstract/Free Full Text] . - Spellman, P.T., Sherlock, G., Zhang, M.Q., Iyer, V.R., Anders, K., Eisen, M.B., Brown, P.O., Botstein, D., Futcher, B. (1998) Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization Mol. Biol. Cell, . 9, 32733297
[Abstract/Free Full Text] . - Mukherjee, S., Berger, M.F., Jona, G., Wang, X.S., Muzzey, D., Snyder, M., Young, R.A., Bulyk, M.L. (2004) Rapid analysis of the DNA-binding specificities of transcription factors with DNA microarrays Nature Genet, . 36, 13311339[CrossRef][Web of Science][Medline] .
- Boorsma, A., Foat, B.C., Vis, D., Klis, F., Bussemaker, H.J. (2005) T-profiler: scoring the activity of predefined groups of genes using gene expression data Nucleic Acids Res, . 33, W592W595
[Abstract/Free Full Text] . - Robinson, M.D., Grigull, J., Mohammad, N., Hughes, T.R. (2002) FunSpec: a web-based cluster interpreter for yeast BMC Bioinformatics, 3, 35[CrossRef][Medline] .
- Warner, J.R. (1999) The economics of ribosome biosynthesis in yeast Trends Biochem. Sci, . 24, 437440[CrossRef][Web of Science][Medline] .
- Rudra, D., Zhao, Y., Warner, J.R. (2005) Central role of Ifh1p-Fhl1p interaction in the synthesis of yeast ribosomal proteins EMBO J, . 24, 533542[CrossRef][Web of Science][Medline] .
- Schawalder, S.B., Kabani, M., Howald, I., Choudhury, U., Werner, M., Shore, D. (2004) Growth-regulated recruitment of the essential yeast ribosomal protein gene activator Ifh1 Nature, 432, 10581061[CrossRef][Medline] .
- Wade, J.T., Hall, D.B., Struhl, K. (2004) The transcription factor Ifh1 is a key regulator of yeast ribosomal protein genes Nature, 432, 10541058[CrossRef][Medline] .
- Devlin, C., Tice-Baldwin, K., Shore, D., Arndt, K.T. (1991) RAP1 is required for BAS1/BAS2- and GCN4-dependent transcription of the yeast HIS4 gene Mol. Cell. Biol, . 11, 36423651
[Abstract/Free Full Text] . - Tornow, J., Zeng, X., Gao, W., Santangelo, G.M. (1993) GCR1, a transcriptional activator in Saccharomyces cerevisiae, complexes with RAP1 and can function without its DNA binding domain EMBO J, . 12, 24312437[Web of Science][Medline] .
- Pina, B., Fernandez-Larrea, J., Garcia-Reyero, N., Idrissi, F.Z. (2003) The different (sur)faces of Rap1p Mol. Genet. Genomics, 268, 791798[Web of Science][Medline] .
- Chang, C.F., Wai, K.M., Patterton, H.G. (2004) Calculating the statistical significance of physical clusters of co-regulated genes in the genome: the role of chromatin in domain-wide gene regulation Nucleic Acids Res, . 32, 17981807
[Abstract/Free Full Text] . - Grunstein, M. (1997) Molecular model for telomeric heterochromatin in yeast Curr. Opin. Cell. Biol, . 9, 383387[CrossRef][Web of Science][Medline] .
- Pryde, F.E. and Louis, E.J. (1999) Limitations of silencing at native yeast telomeres EMBO J, . 18, 25382550[CrossRef][Web of Science][Medline] .
- Holstege, F.C., Jennings, E.G., Wyrick, J.J., Lee, T.I., Hengartner, C.J., Green, M.R., Golub, T.R., Lander, E.S., Young, R.A. (1998) Dissecting the regulatory circuitry of a eukaryotic genome Cell, 95, 717728[CrossRef][Web of Science][Medline] .
- Beinoraviciute-Kellner, R., Lipps, G., Krauss, G. (2005) In vitro selection of DNA binding sites for ABF1 protein from Saccharomyces cerevisiae FEBS Lett, . 579, 45354540[CrossRef][Web of Science][Medline] .
- Miyake, T., Reese, J., Loch, C.M., Auble, D.T., Li, R. (2004) Genome-wide analysis of ARS (autonomously replicating sequence) binding factor 1 (Abf1p)-mediated transcriptional regulation in Saccharomyces cerevisiae J. Biol. Chem, . 279, 3486534872
[Abstract/Free Full Text] . - Schroeder, S.C. and Weil, P.A. (1998) Biochemical and genetic characterization of the dominant positive element driving transcription of the yeast TBP-encoding gene, SPT15 Nucleic Acids Res, . 26, 41864195
[Abstract/Free Full Text] . - Lascaris, R.F., Groot, E., Hoen, P.B., Mager, W.H., Planta, R.J. (2000) Different roles for Abf1p and a T-rich promoter element in nucleosome organization of the yeast RPS28A gene Nucleic Acids Res, . 28, 13901396
[Abstract/Free Full Text] . - de Winde, J.H. and Grivell, L.A. (1992) Global regulation of mitochondrial biogenesis in Saccharomyces cerevisiae: ABF1 and CPF1 play opposite roles in regulating expression of the QCR8 gene, which encodes subunit VIII of the mitochondrial ubiquinol-cytochrome c oxidoreductase Mol. Cell. Biol, . 12, 28722883
[Abstract/Free Full Text] . - Della Seta, F., Ciafre, S.A., Marck, C., Santoro, B., Presutti, C., Sentenac, A., Bozzoni, I. (1990) The ABF1 factor is the transcriptional activator of the L2 ribosomal protein genes in Saccharomyces cerevisiae Mol. Cell. Biol, . 10, 24372441
[Abstract/Free Full Text] . - Drazinic, C.M., Smerage, J.B., Lopez, M.C., Baker, H.V. (1996) Activation mechanism of the multifunctional transcription factor repressor-activator protein 1 (Rap1p) Mol. Cell. Biol, . 16, 31873196[Abstract] .
- Wyrick, J.J., Holstege, F.C., Jennings, E.G., Causton, H.C., Shore, D., Grunstein, M., Lander, E.S., Young, R.A. (1999) Chromosomal landscape of nucleosome-dependent gene expression and silencing in yeast Nature, 402, 418421[CrossRef][Medline] .
- Tanay, A. (2006) Extensive low-affinity transcriptional interactions in the yeast genome Genome Res, . 16, 962972
[Abstract/Free Full Text] . - Cernilogar, F.M. and Orlando, V. (2005) Epigenome programming by Polycomb and Trithorax proteins Biochem. Cell. Biol, . 83, 322331[CrossRef][Web of Science][Medline]
.
This article has been cited by other articles:
![]() |
S. A. Ansari, Q. He, and R. H. Morse Mediator complex association with constitutively transcribed genes in yeast PNAS, September 29, 2009; 106(39): 16734 - 16739. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Krushevskaya, H. Peterson, J. Reimand, M. Kull, and J. Vilo VisHiC--hierarchical functional enrichment analysis of microarray data Nucleic Acids Res., July 1, 2009; 37(suppl_2): W587 - W592. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Ezkurdia, L. Bartoli, P. Fariselli, R. Casadio, A. Valencia, and M. L. Tress Progress and challenges in predicting protein-protein interaction sites Brief Bioinform, May 1, 2009; 10(3): 233 - 246. [Abstract] [Full Text] [PDF] |
||||
![]() |
F.E. Frenkel and E.V. Korotkov Using Triplet Periodicity of Nucleotide Sequences for Finding Potential Reading Frame Shifts in Genes DNA Res, April 1, 2009; 16(2): 105 - 114. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Yu, J. B. Smirnova, E. C. Friedberg, B. Stillman, M. Akiyama, T. Owen-Hughes, R. Waters, and S. H. Reed ABF1-binding Sites Promote Efficient Global Genome Nucleotide Excision Repair J. Biol. Chem., January 9, 2009; 284(2): 966 - 973. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Sato, Y. Nakamura, T. Kaneko, E. Asamizu, T. Kato, M. Nakao, S. Sasamoto, A. Watanabe, A. Ono, K. Kawashima, et al. Genome Structure of the Legume, Lotus japonicus DNA Res, August 1, 2008; 15(4): 227 - 239. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. Schlecht, I. Erb, P. Demougin, N. Robine, V. Borde, E. v. Nimwegen, A. Nicolas, and M. Primig Genome-wide Expression Profiling, In Vivo DNA Binding Analysis, and Probabilistic Motif Prediction Reveal Novel Abf1 Target Genes during Fermentation, Respiration, and Sporulation in Yeast Mol. Biol. Cell, May 1, 2008; 19(5): 2193 - 2207. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Orchard and H. Hermjakob The HUPO proteomics standards initiative--easing communication and minimizing data loss in a changing world Brief Bioinform, March 1, 2008; 9(2): 166 - 173. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Yuan and K.-C. Li Context-dependent clustering for dynamic cellular state modeling of microarray gene expression Bioinformatics, November 15, 2007; 23(22): 3039 - 3047. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Kasahara, K. Ohtsuki, S. Ki, K. Aoyama, H. Takahashi, T. Kobayashi, K. Shirahige, and T. Kokubo Assembly of Regulatory Factors on rRNA and Ribosomal Protein Genes in Saccharomyces cerevisiae Mol. Cell. Biol., October 1, 2007; 27(19): 6686 - 6705. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. N. Koehler, N. Rachfall, and R. J. Rolfes Activation of the ADE Genes Requires the Chromatin Remodeling Complexes SAGA and SWI/SNF Eukaryot. Cell, August 1, 2007; 6(8): 1474 - 1485. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Labarga, F. Valentin, M. Anderson, and R. Lopez Web Services at the European Bioinformatics Institute Nucleic Acids Res., July 13, 2007; 35(suppl_2): W6 - W11. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||













