Blurring of High-Resolution Data Shows that the Effect of Intrinsic Nucleosome Occupancy on Transcription Factor Binding is Mostly Regional, Not Local
[['Genomic DNA is largely covered in proteins , mostly in the form of nucleosomes .', 'Much of the remainder consists of chromatin-associated proteins , including enzymes that modify histones or DNA , or catalyze the rearrangement of nucleosomes , and sequence specific DNA binding proteins ( transcription factors ) that mediate the activation or repression of genes .', 'A deep understanding of gene regulation requires an understanding of how each of these cooperate and compete for access to genomic DNA .', 'That nucleosomes and TFs do , in fact , compete on a genomic-wide scale was substantiated several years ago by chromatin immunoprecipitation microarray experiments ( ChIP-chip ) , which determined the distribution of histones along the yeast genome .', '[1]–[3] These studies revealed an under-representation of nucleosomes in promoter regions , relative to transcribed regions .', 'In contrast , TFs are under-represented in transcribed regions and enriched in promoter regions .', 'Furthermore , nucleosome occupancy differs among promoters and these differences correlate with the probability of TF binding .', '[4] Competition between nucleosomes and transcription factors is a simple consequence of each having an inherent probability of binding to the same site .', 'Some transcription factors may be able to bind DNA on the outside surface of a nucleosome , but steric occlusion of sites on the inside and the sharp bending of DNA around the nucleosome preclude most transcription factors from binding to nucleosomal DNA .', 'As a consequence of this competition , nucleosomes can mediate interactions among transcription factors ( TFs ) in an entirely passive way .', 'For example , a binding motif that is close to a second , TF-occupied , site will tend to have a lower nucleosome occupancy than it would in the absence of the occupied site because nucleosome configurations that span both the motif and a nearby occupied site are disallowed .', 'This lower nucleosome occupancy translates into a higher effective binding affinity of the site .', 'In this scenario cooperative binding of factors is mediated not by direct interactions between the factors but by the passive effects of nucleosomes due to mutual competition .', 'This effect has been demonstrated experimentally .', '[5] Passive mediation of TF-TF interactions is one way in which nucleosomes affect the occupancy of TFs in the genome .', 'A second is through the intrinsic sequence specificity of nucleosomes themselves .', 'Nucleosomes lack highly specific amino acid side chain-to-base pair contacts that are characteristic of sequence-specific transcription factors , but they do have sequence preferences that are determined by the capacity of the DNA to be wrapped around the nucleosome .', 'One manifestation of this preference is a subtle but significant tendency towards a 10bp periodicity of certain dinucleotide steps .', '[2] , [6]–[9] This periodicity reflects the helical periodicity of DNA in nucleosomes and differences in the propensity for structural perturbations among different basepair-basepair steps .', '[10] More recently , attention has been drawn to longer A-rich regions , which are inhibitory for nucleosome binding [9] , [11] , [12] .', 'The near-exclusion of these sequences from the central regions of nucleosomes was first noted around the same time that dinucleotide preferences were discovered [7] but it has become clear only recently that these sequences contribute substantially to the ability to predict nucleosome binding .', '[9] , [11] , [12] Kaplan et al recently provided a simple and elegant demonstration that the intrinsic sequence specificity for nucleosome binding is a partial determinant of nucleosome positioning in vivo .', '[9] In their work , nucleosome positions were determined from deep sequencing of nuclease-digested chromatin prepared from yeast cells , and also from chromatin reconstituted in vitro from chicken histones and yeast genomic DNA .', 'Sequence features associated with nucleosome occupancy were found to be similar in real chromatin and in reconstituted chromatin .', 'Furthermore , genomic loci that are bound by transcription factors have , on average , lower nucleosome occupancies than unbound sites even when nucleosome occupancies are determined in vitro .', 'This suggests that intrinsic nucleosome binding preferences have some effect on the selection of binding sites by transcription factors in vivo .', 'Kaplan et al further reported that there are no significant differences in mapped nucleosome positions when chromatin was crosslinked compared to the more conventional procedure of mapping without crosslinking .', '[9] Mapping without crosslinking assumes that sites occupied by nucleosomes in vivo remain bound to those sites on the time scale of the assay .', 'The reported similarity of the data using the two procedures suggests that that is largely the case .', 'Here , we revisit the data of Kaplan et al , analyzing it in the context of existing ChIP-chip data [13] as well as new , more quantitative ChIP-qPCR experiments reported here .', 'We confirm a role for intrinsic nucleosome binding preferences in the binding of transcription factors , though our analyses suggest a rather low bound on how much information may be provided .', 'More interestingly , we show that the role played by sequence specific nucleosome binding lies not so much in the precise positioning of nucleosomes as in the determination of the average nucleosome density over a promoter-sized region .', 'We also show that there is a difference between crosslinked and uncrosslinked chromatin that is at least as informative with regard to the binding of transcription factors as either data set is alone .'], ['The analyses in this paper make extensive use of the nucleosome mapping data reported recently by Kaplan et al , so it was important to first establish that we can replicate and extend an analysis that they performed .', 'To that end , we took Abf1 motifs in the yeast genome and calculated a set of averaged profiles of nucleosome sequence tags around those motifs ( Figure 1A ) .', 'The most prominent of these nucleosome profiles , showing a very deep minimum for occupancy at the Abf1 motifs , was calculated for Abf1 motifs that are bound by the protein with high confidence ( ChIP-chip enrichment p-value≤1e-3 ) .', 'This particular profile likely corresponds to the single profile shown in Figure 4c of Kaplan et al , [9] though it has been calculated here over a wider genomic interval to show the enrichment of phased nucleosomes up to three units away from the Abf1 motifs .', 'We extended this analysis by testing other Abf1 sites for which the evidence for Abf1 binding is weaker .', 'Remarkably , as the stringency for defining Abf1 binding is relaxed , from the most stringent ( p≤1e-3 ) to the least ( p>0 . 5 ) , both the depletion of nucleosomes at the Abf1 motif and the enrichment of nucleosomes at adjacent flanking positions decrease , but not to the point where they disappear altogether ( Figure 1A ) .', 'This is reminiscent of the analyses of Tanay , who provided evidence for TF binding at ChIP enrichment p-values far worse than what would ordinarily be considered a meaningful threshold for binding .', '[14] To assess quantitatively the correlation between low nucleosome occupancy and TF binding we asked how well nucleosome tag counts correctly distinguish TF-bound sites from random sites selected from yeast promoters .', 'We use the area under the ROC curve ( ROC AUC: receiver operator characteristic area under the curve ) as a measure of this association .', '[15] As shown in Figure 1B , for even the lowest confidence Abf1 binding sites ( p>0 . 5 ) , the ROC AUC exceeds 0 . 5 , the value that is expected by chance .', 'This analysis was extended to the 41 yeast TFs for which there are at least 50 binding motifs bound according to the ChIP-chip data of Harbison et al ( p≤1e-3 ) and the subsequent motif analysis of MacIsaac et al . [13] , [16] Bootstrap analysis of ROC curve areas shows a significant association between TF occupancy and nucleosome depletion for most of the 41 TFs ( Figure 1C ) .', 'Kaplan et al used two different methods in their nucleosome mapping experiments , one involving formaldehyde crosslinking ( two replicates ) and the other a more traditional non-crosslinking protocol ( four replicates ) .', '[9] Crosslinking should be unnecessary if nucleosomes are sufficiently stable that nucleosomal locations are the same at the end of the assay as they are in vivo .', 'Having performed the experiment both ways , Kaplan et al deemed the two data sets to be sufficiently similar that they averaged all six replicates and used this single set of averaged tag counts for their analyses .', '[9] However , they make available the two separate datasets .', 'For Figure 1 we used data from uncrosslinked chromatin only , but we performed the analysis with data from crosslinked chromatin as well .', 'Our expectation was that nucleosome occupancies obtained from crosslinked chromatin would show , if anything , a stronger association with transcription factor binding than uncrosslinked chromatin because crosslinking would prevent the movement of nucleosomes into regions that are occupied by TFs in vivo but which might rearrange in the time course of the experiment .', 'Surprisingly , we found that crosslinking generally weakened the association between TF binding and nucleosome depletion at TF sites , rather than strengthening it .', '( Figure S1 )', 'To investigate this result further , we examined the crosslinked tag count distribution around Abf1 sites , as was done for Figure 1A using the data from uncrosslinked chromatin .', 'Inspection of the tag count distribution around bound Abf1 motifs reveals a remarkable concordance between the two data sets in the enrichment of nucleosome tag sequences flanking the binding site .', 'These correspond to a series of ( averaged ) phased nucleosomes adjacent to Abf1 binding sites ( Figure 2A ) .', 'However , at the Abf1 site itself there is a reduction in the amount of nucleosome depletion inferred from the crosslinked chromatin .', 'It is this slightly less dramatic effect that reduces the predictive value of nucleosome tags in predicting Abf1 binding sites when crosslinked data is used rather than uncrosslinked .', 'The bottom panel in Figure 2A shows the difference between the averaged tag counts in the crosslinked experiment and the averaged ( and normalized ) tag counts in the uncrosslinked experiment .', 'A clear peak in excess tag counts can be seen at Abf1 sites , suggesting that excess tag counts found in the crosslinked experiment are associated with TF binding .', 'To assess this more broadly , we normalized and subtracted , genome wide , the tag counts obtained with uncrosslinked chromatin from those obtained with crosslinked chromatin to produce a “nucleosome difference map” .', 'We then asked whether high tag counts in the difference map could predict TF binding in a manner analogous to how low tag counts in the raw nucleosome maps predict TF binding .', 'Remarkably , excess tags in the crosslinked chromatin are as strongly correlated with TF binding as nucleosome depletion is in the uncrosslinked sample ( Figure 2B ) .', 'This is reflected in the strong correlation between the two sets of ROC AUC values ( R\u200a=\u200a0 . 80 ) and in the slope of the line relating the values ( ∼1 . 02 ) .', 'If anything , the excess tags found in the crosslinked sample are more strongly correlated with TF binding than is nucleosome depletion in the uncrosslinked library alone .', '26 out of 41 TFs have higher ROC AUC values based on excess crosslinked tags , and the absolute value of the differences for those 26 are about 50% higher on average than for the remaining 15 .', 'We suspect this effect is due to the crosslinking of nucleosomes that span the binding site .', 'Crosslinking will trap transiently bound nucleosomes , and will likely do so more efficiently than for TFs because of the large number of amines ( lysines and arginines ) that lie in close proximity to DNA .', 'That crosslinking appears to be trapping nucleosomes over TF binding sites is illustrated by the appearance of a very strong nucleosome peak in the difference map that lies right on top of a set of Gal4 sites in the GAL1–GAL10 promoter ( Figure S2 ) .', 'The nucleosome is present when crosslinked , and nearly absent when not; a set of six neighboring nucleosome positions are scarcely affected by crosslinking .', 'The crosslinked nucleosome is much less prominent when cells are grown in galactose , presumably because Gal4 occupancy is higher under these conditions , resulting in less opportunity for a nucleosome to be crosslinked at that location .', 'Regardless of the mechanism , the association between TF bound sites and excess tag counts in crosslinked chromatin suggests that difference maps based on crosslinked and uncrosslinked chromatin might be used to identify non-histone DNA-binding sites without ChIP enrichment for particular proteins .', 'How such sites would compare to DNase hypersensitive sites or nucleosome poor regions defined by FAIRE [17] remains to be seen .', 'Kaplan et al . made the important observation that genomic loci that are bound by TFs in vivo tend to be also depleted for nucleosomes in reconstituted chromatin .', '[9] This shows that at least some of the low nucleosome occupancy observed at TF sites is intrinsic to the DNA binding specificity of nucleosomes and is not simply a consequence of competition by TF binding .', 'As a prelude to our analysis of resolution-sensitivity , we validated the observations made by Kaplan et al using the same TF binding data but with a different analytical measure ( ROC AUC vs . average tag counts ) .', 'We also obtained additional , higher accuracy TF binding data using ChIP-qPCR at selected binding motifs in order to establish more quantitatively the correlation between binding and nucleosome occupancy .', 'ROC AUC values are generally similar whether nucleosome occupancies are obtained in vivo or in vitro ( Figure 3A ) .', 'Not surprisingly , the values are somewhat higher with in vivo chromatin than with in vitro reconstituted chromatin for most TFs ( 33/41 , notably for Abf1 and Reb1 ) .', 'This indicates that some of the nucleosome depletion at binding sites is a consequence of TF binding and not really an effector of it .', 'Nevertheless , the correlation between in vivo and in vitro values is remarkably good ( R\u200a=\u200a0 . 79; R\u200a=\u200a0 . 90 if Abf1 and Reb1 are removed as outliers ) .', 'Neither the ROC AUC values nor the raw differences in tag counts used by Kaplan et al lend themselves to a simple interpretation in terms of the amount of TF binding information that lies in the intrinsic binding specificity of nucleosomes .', 'To assess more directly how much of an effect on TF binding is encoded by the intrinsic DNA binding specificity of nucleosomes , we determined the apparent binding occupancies of 107 perfect consensus binding sites in the genome using ChIP-qPCR ( Table S1 ) .', 'Between 16 and 33 consensus sites were assayed for each of four TFs ( Dig1 , Bas1 , Gcn4 and Nrg1 ) .', 'The four TFs are typical of those evaluated in Figure 3A , lying close to the best fit through the data , but have associations with nucleosome occupancy that are skewed to lower-than-average values .', 'This makes them particularly stringent targets for independent evaluation .', 'Figure 3B shows that ChIP enrichment values are correlated with nucleosome occupancy in the expected direction ( i . e . higher ChIP enrichment is correlated with lower nucleosome occupancy ) .', 'Correlation coefficients to the in vivo nucleosome data average 0 . 34 and range from 0 . 05 for Bas1 to 0 . 60 for Nrg1 .', 'All except Bas1 are significantly different than 0 ( i . e . , show a significant inverse correlation between ChIP enrichment values and nucleosome occupancy ) .', 'ChIP-qPCR enrichment values appear to be correlated with in vitro nucleosome data as well , but poorly .', '( Figure 3B; right hand set of panels ) .', 'For Gcn4 , for example , only ∼1 . 7% ( R2\u200a=\u200a0 . 132 ) of the variance in ChIP enrichment values is explained by intrinsic nucleosome binding , as defined by reconstituted chromatin , while 16% ( 0 . 402 ) of the variance is explained by nucleosome positions in vivo .', 'Overall , for the four TFs we assayed , we estimate that only about 5% of the variance associated with nucleosome occupancy differences in vivo is due to intrinsic nucleosome positioning; the rest is a consequence of other effects that determine chromatin structure in vivo .', 'The fact that all four TFs show correlations in the expected direction between ChIP-qPCR enrichment and nucleosome occupancy attests to the sensitivity of this analysis because the ChIP-chip based ROC AUC values for these same TFs are only marginally different than the value expected by chance .', 'It is possible that the ROC AUC values are underestimated due to the definition of unbound sites that we chose to use .', 'We chose to use random sites selected from promoter regions ( 600bp 5′ to ORFs ) thinking they would be more appropriate controls for the TF-bound sites , but the ROC AUC values obtained using this background are lower than what is obtained when sites randomly selected from throughout the genome are used instead ( data not shown ) .', 'A systematic underestimation of the true ROC AUC value would also explain why Nrg1 has an ROC AUC value below 0 . 5 , implying a direct association between nucleosome occupancy and binding , even though our ChIP-qPCR analysis unambiguously shows the expected inverse correlation .', 'While absolute ROC AUC values should be interpreted with caution , comparisons of ROC AUC values are valid because each calculation was performed using the same set of unbound sequences as background .', 'Since there are 41 TFs for which we have performed analyses using ChIP-chip data , we use those data and the ROC AUC metric for all subsequent analyses , rather than the correlation to ChIP-qPCR enrichment values , for which we have data for only four of the 41 TFs .', 'As discussed above , most of the 41 TFs are modestly associated with nucleosome depletion in vivo ( Figure 1C , Figure S3 ) .', 'This is consistent with the conclusion we reached previously using much lower resolution nucleosome occupancy data .', '[4] Those data were obtained from histone ChIP-chip experiments using DNA microarrays whose probes mostly corresponded to entire ORFs or intergenic regions .', '[1] The TF-bound regions were also low-resolution , having been obtained using the same microarray probes .', '[13] In the current analysis we use the same low resolution ChIP-chip data but the precise binding sites within those regions have now been inferred from motif analysis , as described by MacIsaac et al . [16] Thus , the TF binding data has effectively been made higher resolution through bioinformatics methods .', 'The nucleosome occupancy data available now is truly higher resolution .', 'If transcription factor binding depends sensitively on the positioning of nucleosomes , we would expect high resolution data to produce a stronger association between nucleosome depletion and TF binding .', 'To test this , we started with high-resolution data and simulated the effects of lower resolution data by averaging the high-resolution nucleosome occupancy data over windows of various sizes .', 'Figure 4 ( panels A and B ) shows the effect of averaging these data on the association between nucleosome occupancy in vivo and Abf1 binding .', 'Up to this point , we have used 15bp windows spanning TF binding sites ( and randomly selected promoter sites ) to calculate ROC curves and their areas .', 'If a substantially larger window is used instead ( 150bp ) , the ROC AUC is noticeably lower; as the window is expanded to 300bp , the ROC AUC drops even more precipitously , and by 600bp , there remains only a very small association between low nucleosome occupancy and Abf1 binding sites .', 'This dramatic drop in ROC AUC scores is expected because there are high occupancy nucleosomes flanking the Abf1 binding sites; a 300bp window includes the peaks of these nucleosomes , while a 600bp window includes all of those nucleosomes and a bit more of others .', 'Window sizes of 40bp or 75bp actually show somewhat higher ROC AUC scores than the 15bp window , reflecting the fact that the average bound Abf1 has a nucleosome-depleted window that is about 50–75bp ( Figure 4A ) .', 'We repeated this analysis for all 41 TFs , comparing the ROC AUC values obtained with 600bp windows to those obtained with 15bp windows ( Figure 4C; Figure S4 ) .', 'Abf1 and Reb1 , the two outliers in what is otherwise a good correlation between the effects of vivo chromatin vs . in vitro chromatin ( Figure 3A ) , are exceptionally strongly affected by the averaging of nucleosome position data .', 'This is because they have such a strong effect on local nucleosome density: displacement of nucleosomes at the binding site results in high nucleosome occupancy immediately adjacent to the site , and therefore there is a rapid regression to the mean nucleosome occupancy as the window size is expanded .', 'Although there are exceptions , and the effect is less dramatic , this is a trend that is seen for the set of TFs as a whole .', 'Altogether , a majority of TFs ( ∼60% ) show poorer association with nucleosome occupancy when those occupancy values are averaged over 600bp .', 'In addition , of the TFs that are most dramatically affected by averaging nucleosome occupancy over 600bp ( i . e . those with changes greater than two standard deviations larger than the mean ) , all are adversely affected by the averaging of nucleosome occupancy .', 'Remarkably , the opposite effect is observed when in vitro reconstituted chromatin is used in the calculations rather than in vivo chromatin .', 'An improvement in the association between binding and nucleosome occupancy is found for about three quarters of the TFs when in vitro nucleosome occupancies are averaged over 600bp .', 'Among TFs that show the most dramatic effects ( i . e . those exceeding the mean by 2σ ) , twice as many are improved as are made worse .', 'At a cutoff of 1σ , three times as many are improved as are made worse .', 'The number of TFs that are adversely affected by blurring of the in vivo data , and the number of TFs that are positively affected by blurring of the in vitro data , are each significantly different than the numbers expected by chance ( p∼0 . 02 ) .', 'The results with reconstituted chromatin are important because it is those data that are most relevant to an understanding of intrinsic nucleosome binding specificity and its effect on TF binding .', 'The analyses of in vivo nucleosome data serve as a kind of computational control , showing that simulation of low resolution data does indeed weaken the association with TF binding , as would be expected if the precise nucleosome location , as defined by high-resolution sequencing experiments , were relevant to TF binding .', 'The question then , is why is there an improvement in the correlation between TF binding and nucleosome occupancy when high-resolution data for the in vitro nucleosome date are averaged so as to simulate lower resolution data ?', 'To investigate this question further , we examined more closely the patterns of nucleosome occupancy around TF binding sites in vivo and in vitro .', 'To that end , we clustered the 41 TFs into five groups based on their in vivo and in vitro nucleosome profiles ( Methods; Figure 5A ) , and confirmed that members of these groups share similarities in their sensitivity to data blurring in vivo and in vitro ( Figure 5B ) .', 'Most of the TFs for which the blurring of in vitro nucleosome data improves the association with binding have at least one of two properties in their nucleosome occupancy profiles that can explain this result:', '( i ) the nucleosome poor region around the binding site is broad , such that averaging around the binding site provides greater contrast with random control sites or', '( ii ) the binding site is actually higher in nucleosome density than the surrounding regions , so blurring the data encompasses flanking regions that are lower in nucleosome density .', 'This latter set , in particular , suggest that the precise genomic position favored by nucleosomes is less relevant to TF binding than is overall nucleosome density in the region .', 'This is perhaps the case because nucleosomes that occlude binding sites can be displaced to nearby regions at little energetic cost .', 'In contrast , binding of Rap1/Fhl1 is correlated best with local nucleosome occupancy , and is adversely affected by blurring of the data .', 'These TFs , unlike all others , tend to have well occupied nucleosomes that immediately flank their binding sites in vitro .', 'As expected , TFs with similar patterns of nucleosome occupancy around their binding sites are also affected in similar ways by the averaging of nucleosome occupancy data ( Figure S4 ) .', 'As a further test of the effect of blurring high resolution nucleosome position data , we incorporated the data into a computational model that predicts TF binding to genomic regions .', '[18] The model uses position weight matrices ( PWMs ) to estimate Kd values for all sequence windows in the genome , and from those Kd values and an assumed protein concentration , it calculates the probability the protein is bound to at least one location within a genomic interval .', 'Previously , we showed that low resolution nucleosome occupancy data , obtained from histone ChIP-chip experiments , could be used as weights in the calculation of Kd values in this model , and that these weights improved the prediction of binding as verified via a ChIP-chip experiment .', '[4] Others have also shown the utility of incorporating nucleosome occupancy data in this way .', '[19] Here , we used the same weighting function and parameter values developed previously , but instead of applying weights based on large genomic regions , we used the base-pair resolution , in vitro nucleosome position data of Kaplan et al . [9] The ChIP-chip data of Harbison were used to evaluate the predictive efficacy of these weights .', '[13] Specifically , we first scored yeast promoters for the probability of binding based on our computational model and the PWMs for transcription factors .', '[18] We then evaluated how well this predicted binding identified genes whose promoters are bound experimentally .', 'The calculations were then repeated twice , once using weights based on in vitro nucleosome occupancy data at 15bp resolution , the other using weights based on averaging this nucleosome occupancy data over sliding 600bp windows .', 'Figure 6A illustrates the effect of these weighting schemes on two genomic regions , each containing a perfect consensus binding site for Gcn4 .', 'For each TF , we obtained three ROC AUC values that express how well binding is predicted: one based on the PWM alone; the second based on the PWM , but with genomic position weights determined by high resolution nucleosome position data; and the third based on the PWM and weights determined by simulated low resolution nucleosome position data ( i . e . high resolution data averaged over 600bp windows ) .', 'For TFs whose binding is well predicted by genomic sequence and the PWM alone , the inclusion of weights based on nucleosome occupancy evidently adds noise to the calculation , worsening the predictions .', 'However , for TFs whose binding is poorly predicted by sequence alone , the inclusion of binding affinity weights can substantially improve the prediction of binding ( Figure 6B ) .', 'Strikingly , the effect of intrinsic nucleosome position data on binding predictions is accentuated with the simulated low resolution data .', 'This is the opposite of what we would expect if precise nucleosome positioning were typically of great relevance to the binding of transcription factors , and it is the opposite of what we observed in most cases with the in vivo data .', 'Of course , the improvement in binding predictions with blurred data is for the set of bound promoters as a whole; within this set , some of the promoters bound by the TF fall in rank even if , overall , weighting improves the ROC AUC value ( Figure S5 ) .', 'For example , even though Yap5 is the most responsive TF to nucleosome averaging overall , 22 of the 88 promoters bound by Yap5 drop in rank when in vitro nucleosome data is averaged over a 600bp window .', 'Nevertheless , the overall effect of blurred nucleosome data on binding predictions provides additional support for our contention that it is the regional nucleosome occupancy that matters most in the localization of TF binding sites , not the precise position of the nucleosome .'], ['The yeast genome sequence ( Aug 2008 build ) and gene feature files were obtained from the Saccharomyces Genome Database ( SGD ) .', '[20] Promoter sequences were defined as the 600bp 5′ to the start of transcription of protein coding gene .', 'The genomic positions of ChIP microarray probes and their ChIP enrichment p-values using different epitope tagged transcription factors under normal growth conditions ( YPD media , 30°C ) were obtained from Harbison et al . [13] The nucleosome sequence tag maps of Kaplan et al . [9] were obtained from GEO ( accession number GSE13622 ) .', 'Sequence coordinates for each data set were converted to match the SGD Aug 2008 version of the genome based on a list of coordinate differences maintained at SGD .', 'Genomic loci deemed to be potential transcription factor ( TF ) binding sites on the basis of sequence analyses were obtained from MacIsaac et al . [16] .', 'For each potential TF binding site , the ChIP-chip probe spanning that site was identified and the ChIP enrichment p-value for that that probe was then assigned to the binding motif .', 'Binding motifs with p-values<0 . 001 , were classified as being bound by their respective TFs .', '41 TFs had ≥50 bound motifs by this criterion and were used for all subsequent analyses .', 'For the analysis of Abf1 binding , Abf1 motifs were further binned into p-value intervals 0 . 001–0 . 01 , 0 . 01–0 . 1 , 0 . 1–0 . 5 , >0 . 5 .', 'The nucleosome sequence tag data provided by Kaplan et al consists of a 5′ end , determined by sequencing , and a 3′ end 146bp away that is based on knowledge about the size of nucleosomes and on the preparation in the experiment of ∼150bp sized DNA by nuclease treatment and size-selection .', '[9] For simplicity , we refer throughout the paper to the inferred 146bp sequence as a ‘nucleosome tag’ , or simply a ‘tag’ .', 'The number of tags spanning a particular genomic basepair can be enumerated and is taken to be a measure of the nucleosome occupancy at that basepair .', 'For most analyses , tag counts were averaged over windows that were centered on a binding motif or on randomly selected basepairs from within promoters .', 'We refer to the windows as being of size 15 , 40 , 75 , 150 , 300 and 600 bp , though technically some of the windows are 1bp longer depending on whether the motif is of even length or odd .', 'Receiver operating characteristic ( ROC ) curves and the area under those curves ( AUC ) were used to quantify the ability of a predictor ( nucleosome tag counts ) to correctly classify sequences as TF bound ( defined by ChIP-chip ) or unbound ( randomly selected from yeast gene promoters ) .', 'Where error bars are shown for ROC AUC values , these were estimated from 1000-fold bootstrap re-sampling .', 'ROC AUC values were also used to quantify the predictive value of a TF binding prediction algorithm , with and without weights based on nucleosome occupancies .', 'In this case , the predictor is the estimated TF binding occupancy and the question is how well that value classifies promoters as TF bound or unbound .', 'Unbound sites were selected from 1000 randomly picked yeast promoters , defined as the 600bp region 5′ to the start of a gene .', 'For each of the 41 TFs , and for each of the two nucleosome datasets , we enumerated the number of nucleosomal sequences spanning each basepair in a 1200bp window .', 'The tag counts were averaged across the center of the profiles , and normalized to the mean value in that profile .', 'These 1200bp normalized windows were used to visualize the profiles for each TF ( Fig
['Introduction', 'Results', 'Methods', 'Discussion']
['Genome wide maps of nucleosome occupancy in yeast have recently been produced through deep sequencing of nuclease-protected DNA .', 'These maps have been obtained from both crosslinked and uncrosslinked chromatin in vivo , and from chromatin assembled from genomic DNA and nucleosomes in vitro .', 'Here , we analyze these maps in combination with existing ChIP-chip data , and with new ChIP-qPCR experiments reported here .', 'We show that the apparent nucleosome density in crosslinked chromatin , when compared to uncrosslinked chromatin , is preferentially increased at transcription factor ( TF ) binding sites , suggesting a strategy for mapping generic transcription factor binding sites that would not require immunoprecipitation of a particular factor .', 'We also confirm previous conclusions that the intrinsic , sequence dependent binding of nucleosomes helps determine the localization of TF binding sites .', 'However , we find that the association between low nucleosome occupancy and TF binding is typically greater if occupancy at a site is averaged over a 600bp window , rather than using the occupancy at the binding site itself .', 'We have also incorporated intrinsic nucleosome binding occupancies as weights in a computational model for TF binding , and by this measure as well we find better prediction if the high resolution nucleosome occupancy data is averaged over 600bp .', 'We suggest that the intrinsic DNA binding specificity of nucleosomes plays a role in TF binding site selection not so much through the specification of precise nucleosome positions that permit or occlude binding , but rather through the creation of low occupancy regions that can accommodate competition from TFs through rearrangement of nucleosomes .']
['Genomic DNA is largely covered by proteins that compete with one another for binding to regulatory sequences .', 'Most of these proteins are in the form of nucleosomes .', 'How nucleosomes come to occupy particular sites and thereby compete with sequence specific transcription factors is a central problem in developing a systems-level understanding of gene regulation .', 'Here , we performed a series of computational analyses using high-resolution nucleosome position data that has recently become available in yeast , thanks to advances in DNA sequencing technology .', 'Analysis of these data , combined with data on the location and occupancy of transcription factors genome-wide , shows that the precise location of nucleosomes as determined by nucleosome sequence specificity is often less important to transcription factor binding than the broader , regional occupancy of nucleosomes that is encoded in genomic DNA .', 'This result has implications for the evolution of DNA regulatory elements .']
['biochemistry/transcription and translation', 'genetics and genomics/gene expression', 'genetics and genomics/chromosome biology', 'computational biology/transcriptional regulation']