Genome-wide analysis of chromatin features identifies histone modification sensitive and insensitive yeast transcription factors
© Cheng et al.; licensee BioMed Central Ltd. 2011
Received: 24 May 2011
Accepted: 7 November 2011
Published: 7 November 2011
We propose a method to predict yeast transcription factor targets by integrating histone modification profiles with transcription factor binding motif information. It shows improved predictive power compared to a binding motif-only method. We find that transcription factors cluster into histone-sensitive and -insensitive classes. The target genes of histone-sensitive transcription factors have stronger histone modification signals than those of histone-insensitive ones. The two classes also differ in tendency to interact with histone modifiers, degree of connectivity in protein-protein interaction networks, position in the transcriptional regulation hierarchy, and in a number of additional features, indicating possible differences in their transcriptional regulation mechanisms.
Transcription factors regulate target gene expression through binding to specific genomic regions. In Saccharomyces cerevisiae, transcription factor (TF) binding sites (TFBSs) are often adjacent to and upstream of target loci due to the compact nature of the yeast genome [1, 2]. Upon binding, TFs interact with RNA polymerase II to activate or repress transcription. TFs also recruit chromatin modification enzymes to induce chromatin structure changes, which in turn affect the accessibility of factors to genomic DNA regions [3, 4]. The target genes of a TF change according to developmental, physiological and extracellular environmental conditions . In addition, TFs interact with each other through combinatorial binding . Uncovering TF target genes and inter-relationships between TFs for all different conditions is thus important for understanding gene expression regulation, but it is also a difficult task due to the scale of the problem.
Several different experimental methods have been developed to identify TFBSs. Chromatin immunoprecipitation followed by tiling array (ChIP-chip) has been widely used to identify TFBSs at the genomic scale [7–9]. More recently, high-throughput sequencing after chromatin immunoprecipitation (ChIP-Seq) has shown promise in identifying TFBSs at higher resolutions [10, 11]. With these methods, increasing amounts of TFBS data have been accumulated for different TFs in different species, cell types, conditions, and so on, which have started to unravel the global picture of gene expression regulation. In yeast, the TFBSs and target genes for an almost complete set of TFs have been mapped in common YPD medium using ChIP-chip . However, resources are still too limited to support a complete exploration of TF binding for all the combinations of cell types and conditions.
Many computational methods have also been proposed to predict TFBSs [12–19]. These methods are mostly based on the idea that the binding of a TF is mediated by the recognition of its binding motif represented as a position-specific scoring matrix (PSSM). PSSMs are usually discovered as enriched motifs from TFBSs in ChIP-chip or ChIP-seq experiments, or de novo from non-coding genomic sequences [5, 20]. Scanning and matching PSSMs in the genome constitute the core of these methods, which are then improved by incorporating information on motif conservation and TFBS co-localization. Nevertheless, these methods often lead to considerably high rates of false positives. Furthermore, most of these methods are not condition specific and thus do not reflect the dynamic nature of TF binding under different conditions.
Chromatin modifications can modulate the accessibility of DNA regions and affect the recruitment of TFs [3, 4, 21]. Both functions directly relate to transcription regulation by TFs. Genomic mapping of chromatin modifications in yeast using ChIP-chip has provided the opportunity to investigate their underlying relationships with TFBSs [22, 23]. Many chromatin modifications have been shown to be associated with transcription activation and repression [3, 4]. Recent studies have shown that incorporating histone modification data improves prediction of TFBSs in mouse and human [24, 25]. In these models, chromatin modifications generally provide non-TF-specific chromatin accessibility, while PSSMs determine TF-specific bindings.
Here we propose a method that integrates PSSMs and chromatin modifications to improve TF target gene predictions in yeast. Specifically, we trained individual support vector machine (SVM) models  for 203 yeast TFs using 2 types of features: the existence of PSSMs upstream of genes and chromatin modifications adjacent to the ATG start codons. The models were trained and tested using TF target genes from ChIP-chip experiments. Interestingly, we found that some yeast TFs were more sensitive to histone modifications than others: their target genes had relatively higher histone modification signals, and were better predicted when these signals were included in the models. We termed these two TF classes as histone-sensitive and histone-insensitive TFs for simplicity. We show that histone-sensitive and histone-insensitive TFs differ in many biological characteristics found in large-scale and small-scale experiments. Furthermore, we used the model to investigate condition specificity and TF-specificity of chromatin modifications as well as TF-TF co-operation. Our analysis helps understand the mechanism of gene expression regulation by TFs and chromatin modifications.
Differential histone modifications between ChIP-chip verified and non-verified TFBSs
Improving target gene prediction by combining histone modifications and PSSMs
For chromatin modifications, we used 11 histone modifications that covered most yeast ORF regions from ChIP-chip experiments . Since TFs bind to the upstream sequence of ORFs, we focused on histone modification signals in the 1 kb of sequence flanking transcription start sites (ATGs), because TFBSs were enriched in these regions. For TF PSSMs, two independent sets were obtained from previous studies. One set comprised PSSMs discovered using a sequence analysis-based method, basically looking for enriched motifs in the DNA regions upstream of all yeast ORFs . From approximately 5, 600 upstream sequences a total of 666 motifs have been discovered, among which 48 could be associated with known TFs according to the literature. The other set of PSSMs was based on ChIP-chip data . For each TF a target gene set was determined and then the binding motif for the TF was identified from the DNA region upstream of these genes.
Our results indicate that, for almost all TFs, the model using both histone modification and PSSM data achieves the best or nearly best performance (measured using the area under the receiver operator characteristic (ROC) curve (AUC)) compared to the other two models using histone modification data or PSSM data alone (Figure 2; Additional files 1 and 2). For example, we obtained an AUC of 0.89 for predicting target genes of the factor SUM1 when both the histone modification and PSSM data were used. If only the PSSM or histone modification information was used, however, the models resulted in lower AUCs (0.77 for PSSM only and 0.85 for histone modification only; Figure 2a; Additional file 1). The improved performance of the combined model was observed for both of the TF PSSM sets (Figure 2b, c; Additional files 1 and 2), indicating that the improvement does not rely on a particular source or quality of PSSMs. Interestingly, for some TFs, the histone-only model performed better than the PSSM-only model, while for some other TFs, the opposite was observed.
Transcription factors with improved prediction by including multiple histone modification datasets
We also investigated the positional effect of histone modification signals for target gene prediction. First, we observed that signals of different types of histone modifications showed different patterns at DNA regions around the ATG, suggesting that they might affect TF binding in different ways. Second, histone modification signals upstream of the ATG are generally more predictive than those downstream of it, as we have observed for both 500-bp and 1, 000-bp flanking region sizes (Figure 3c, d). This is somewhat expected because TFBSs are more enriched in the upstream regions of ORFs for the purpose of transcriptional regulation. It is also interesting to see that ATG flanking regions of 500 bp and 1, 000 bp result in almost the same performance (Figure 3c, d). Given the compact nature of the yeast genome, transcription start sites for most ORFs are located within the region 1, 000 bp upstream of the ATG . The comparable performance when using 500-bp ATG flanking region indicates that most discriminative histone modification signals for TF binding are embedded in this region. We also find that histone modifications in intergenic regions are more predictive than those in ORF coding regions (Figure 3e).
The number of PSSM-containing genes is another factor that might also affect prediction using histone modifications. We counted the number of genes with PSSMs in their promoters for each TF, and found no obvious relationship between this and the prediction power (AUC) when using histone modifications (Additional files 3 and 4). In fact, the number of PSSM-containing genes varies extensively from 79 to 5, 810 depending on the information content of the PSSMs, but its correlation with the prediction power is negligible (R = 0.1 for the Frankel et al. data set ).
Both histone modifications and TF binding are involved in regulation of gene expression. To rule out the possibility that the capability of histone modifications for predicting TF target genes is actually mediated by the expression level of genes, we examined the marginal effect of gene expression level on TF target prediction. Including expression levels as an additional predictor in the histone modification-based models does not significantly change the prediction accuracy for any TF. Furthermore, SVM models based on expression level alone can hardly predict targets for any of the TFs (Additional file 5). The only exception is FHL1 (prediction accuracy AUC = 0.79), which predominantly regulates the transcription of ribosomal protein genes; it is the extreme abundance of these genes that enables accurate prediction for this TF. Thus, the prediction power of histone modification models is not likely to be mediated by gene expression levels overall.
Condition specificity of the chromatin model
TFs bind to and regulate target gene expression in a complex and dynamic manner to coordinate biological processes [5, 33]. Chromatin modifications also change rapidly in response to stimuli from the extracellular environment . Therefore, chromatin modifications in one condition should match TF target binding in that condition but not other conditions.
We investigated the condition specificity of our chromatin model in two conditions, YPD medium and H2O2. We tested a total of 12 TFs, for which we had the PSSM, histone modification and TF target binding data under both conditions. For each of the TFs, we constructed two separate chromatin models: one model (model A) used PSSM and histone modifications under the YPD condition as features, while the other (model B) used PSSMs from the YPD condition but histone modifications under the H2O2 condition. The two models were then used to predict TF target binding under H2O2 conditions. It is generally believed that TFs keep their binding specificity PSSMs in different conditions and even over large evolutionary distances. Therefore, we used PSSMs from the YPD condition as a close approximation in model B, where no PSSM information is available under the H2O2 condition.
Relative importance of different histone modifications for target prediction
TFs bind to the upstream regions of target genes by recognizing their specific binding motif PSSMs. We then asked an analogous question: do the binding sites of TFs have specific histone modification profiles? To address this question, we calculated a histone modification profile for each TF by averaging upstream ATG histone modification signals over all its target genes. In our analysis, we included 25 different histone modifications from the two studies mentioned above [22, 23].
Correlations between pairs of histone modifications are shown in Figure 5c, based on their signals in the histone modification profiles for all TFs. Only pairs with strong correlation (r > 0.5 or < -0.5) were connected in the form of a correlation network. The dense connectivity in this network reveals strong pairwise redundancy of histone modification signals, which is also indicative of redundancy for predicting target genes.
We next examined the relative importance of each histone modification for predicting target genes of all TFs. Given a histone modification, we compared its signal difference between target and non-target genes of a TF. The signal difference is represented as t-statistics (see Materials and methods for details), which indicate the relative importance of that histone modification for predicting the target genes of a TF. A larger absolute value for a t-statistic indicates more importance. The t-statistics for all histone modifications form a TF-specific profile, denoted as differential modification profiles of the TF. Interestingly, histone-sensitive TFs and histone-insensitive TFs defined based on target histone modification profiles are also distinct according to their differential modification profiles (Figure 5b). This suggests that histone-sensitive and -insensitive TFs are actually robust clusters with different patterns of histone modifications in their target genes.
Histone modification sensitivity of transcription factors
Transcription factor histone sensitivity relates to hierarchical level in a regulatory network
Hierarchical levels in regulatory network
Number of histone-sensitive TFs
Number of histone-insensitive TFs
Interestingly, the histone sensitivity of TFs also indicates distinct co-regulation relationships. Two TFs are said to be co-regulatory if their sets of targets significantly overlap. Among the approximately 14, 000 possible TF pairs, we found 1, 440 significant co-regulatory relationships (P < 0.05, Fisher's exact test). Among the TFs involved in co-regulatory relationships, 64 are histone sensitive and 95 are histone insensitive. Of the 1, 440 significant co-regulatory pairs, 447 are between two histone-sensitive TFs, 437 between two histone-insensitive TFs, and 556 between one histone-sensitive TF and one histone-insensitive TF. Fisher's exact test showed that histone-sensitive TFs are more likely to be involved in a co-regulatory relationship than histone-insensitive TFs (P < 10-16). In summary, the histone-sensitive TFs reside mostly in the upper layers of the regulatory network, and tend to work and communicate with other TFs during transcriptional regulation.
Transcription factor histone sensitivity relates to cellular functions
ACE2, ASH1, CIN5, FKH1, FKH2, MBP1, MCM1, NDD1, RLM1, STB1, STE12, STP1, SWI4, SWI5, SWI6, TEC1
Heat shock or stress conditions
GAT1, MSN4, SKN7, YAP1
GLN3, HAL9, HMS2, HSF1, MGA1, MSN2, WAR1, USV1
PSSM predictability and cooperativity of transcription factors
Top 10 PSSM well-predictable transcription factors
PSSM sensitive TFs
For TFs that are weakly predictable using PSSM information, we hypothesized that these TFs may bind to their targets indirectly by cooperating with other TFs. If this is the case, we would expect to predict the target genes of such a TF accurately by using the PSSM of its cooperative TF. We tested this by using each PSSM to predict targets of all TFs (Additional files 1 and 2). The targets of most TFs were best predicted by their own PSSMs, but some TFs have their targets better predicted using other TFs' PSSMs. For example, YAP1's target genes were better predicted using CAD1's PSSM - the AUC increased from 0.71 to 0.75. Similarly, the target genes of INO4 were better predicted using INO2's PSSM than it's own, with the AUC increasing from 0.78 to 0.81. In fact, YAP1 and CAD1 work together in stress-induced transcriptional responses, and INO2 and INO4 form heteromeric complexes involved in phospholipid biosynthesis [37, 38]. We also found that the PSSMs of the cooperative TFs are actually quite similar, measured using a similarity score range from 0 to 1. The PSSMs of CAD1 and YAP1 have a similarity score of 0.72 (top 1% among all pairs), and those of INO2 and INO4 have a similarity score of 0.55 (top 5%). This further indicates the cooperation between the two TF pairs through indirect binding. Therefore, TF target gene prediction using 'cross-PSSMs' could help identify co-operative interactions between TFs. On the other hand, this suggests that using a TF's own PSSM may not always be best for predicting its target genes, especially when there is evidence it co-operates with another TF.
Chromatin model improves prediction of TF binding sites
We have examined our chromatin model using the TF target data from large-scale ChIP-chip experiments , and shown its effectiveness for predicting target genes. The study by Harbison et al.  investigated TF binding within yeast promoter regions only. However, technical advancement of ChIP-chip and ChIP-seq has enabled us to obtain the binding sites of a TF across the whole genome. Given these more high-resolution data, our chromatin model can also be used to predict TFBSs. For demonstration, we use the ChIP-seq data for STE12 as an example .
Comparison with previous methods
Comparison of several computational methods for target gene prediction
Number of target genes
PSSM alone (FIMO)
Among the previously published methods, EEL and Stubb take advantage of conservation of TF binding motifs between related species, and as shown they achieve relatively more accurate prediction results than FIMO, Cluster-Buster and MCAST. We also tried the 'Chromia' method proposed by Won et al. . Similar to our method, Chromia integrates histone modification and PSSM data sets but using a hidden Markov model. The method has shown impressive performance when applied to genome-wide ChIP-seq data in mouse for predicting TFBSs. However, when applied to the yeast data in our case, it does not result in good prediction due to the low coverage of the histone modification and TF binding data sets [5, 23]. For example, the arrays used for the Pokholok et al. ChIP-chip data contain approximately 42, 000 probes (60-mers), representing only about 20% of the yeast genome . The arrays used for identifying yeast TFBSs are essentially promoter arrays, covering only DNA regions around the transcription start site of yeast ORFs . In practice, our method requires only data for interested regions (for example, promoter regions), and thereby is more flexible and can be applied to a wide range of data sets.
We examine the effectiveness of these methods for predicting STE12 binding sites (see 'Chromatin model improves prediction of TFBSs' section). Chromia predicts STE12 binding sites with an AUC of 0.66 and a PPV of 5.6%, presumably also due to the low resolution of the histone modification data. Those motif-based methods perform similarly to our PSSM-based methods (Figure 7), and taking into account conservation does not lead to significant improvement for the STE12 case.
Histone-sensitive and -insensitive transcription factors
We classified the 203 yeast TFs used in our study into 68 histone-sensitive and 135 histone-insensitive TFs based on the upstream histone modification signals of their target genes. The two classes have generally opposite characteristics with regard to histone modification signals, expression levels, topology in regulatory networks and other biological features.
Steinfeld et al.  have discovered a list of TF-chromatin modifier interactions in yeast from genome-wide analysis of high-throughput experiments. Among the 35 TFs that interact with chromatin modifiers, 20 are classified as histone-sensitive and 15 as histone-insensitive according to our analysis. Namely, there is a significant enrichment of histone-sensitive TFs (P = 0.001, Fisher's exact test) in the chromatin modifier interacting TFs. This suggests that histone-sensitive TFs are more likely to interact with chromatin modifiers, consistent with the observation that their target genes tend to have stronger histone modification signals.
Hitone-sensitive TFs might target highly regulated genes. It is known that gene expression is regulated by specific TFs and their orchestrating chromatin modification enzymes. Thus, stronger histone modification signals upstream of the target genes of the histone-sensitive TFs might imply more intensive transcriptional regulation. Our results showed that cell cycle TFs were mostly histone-sensitive TFs, consistent with the previous knowledge that the cell cycle is highly regulated to achieve cyclical expression of genes.
Most histone modification data used in this study are derived from yeast grown on YPD medium. We found that histone-sensitive TFs tend to be active under the YPD medium condition, as indicated by a larger number of target genes and higher expression levels with respect to those insensitive TFs. It is possible that histones upstream of the target genes of the histone-sensitive TFs have more chance to be modified by histone modification enzymes because they are recruited by these more active TFs.
We found that histone-sensitive TFs were enriched in higher layers of the hierarchical regulatory network. This suggests that histone-sensitive TFs tend to be 'managers' that regulate other TFs and, for such a reason, their binding to target genes is highly regulated through histone modifications. Consistent with this hypothesis, we have observed stronger histone modification signals in the upstream regions of histone-sensitive TFs' target genes. We caution here that the classification of the histone-sensitive and -insensitive TFs is based only on the histone modification signals of their target genes and the contribution of these signals to target predictions. Further experiments might be worthwhile to investigate the relationship between the target selection of the TF classes and histone modifications in more detail.
A more recent paper in human showed the capability of using unique chromatin signatures to identify two distinct classes of genomic elements, active and poised enhancers . Consistently, here we find in yeast that the histone modifications are also informative for distinguishing ChIP-chip verified and non-verified binding sites. As we demonstrated, the method could be applied to identify target genes and the genome-wide binding sites of TFs.
Contribution of histone modifications to condition-specific target prediction
It is widely known that transcriptional regulation is condition specific in that TFs change their binding sites under different conditions. We show here that histone modification data are most predictive of TF target binding under the same condition. This is true especially for those TFs that are mostly active in specific stress conditions.
Because of limited resources, it is impossible to perform exhaustive experiments for every TF, cell type, and species and all possible conditions. Meanwhile, the TF binding recognition motifs, PSSMs, are generally thought to be non-condition-specific, not changing under different conditions . Thus, prediction based on motif information alone cannot provide condition-specific target genes for a TF. As an alternative method, we propose the feasibility of predicting target genes of a TF under a condition of interest by combining histone modification data under that condition with its PSSM. In this way, we can achieve much higher results than using PSSM alone. More importantly, since histone modifications reflect the chromatin states in a specific condition or tissue, the predicted targets are also condition and tissue specific.
Combinatorial interaction of TFs: direct and indirect binding
When a TF binds directly to the promoter regions of its target genes, the enriched motifs identified from its binding sites can be regarded as its own PSSM. However, TFs do not always act individually; sometimes they cooperate with (physically bind to) each other to form regulatory functional units, such as the yeast cell cycle complexes SBF (SWI4-SWI6) and MBF (MBP1-SWI6) . In these indirect binding cases, it is important to distinguish the TFs that are motif-recognizing and those that are not.
By examining the PSSM sensitivity of TFs, we were able to infer some possible combinatorial interactions between TFs. If TF A's targets are better predicted by using TF B's PSSM instead of its own PSSM, then this is an indication of potential cooperation between the two TFs. In particular, TF B directly binds to promoter regions through its PSSM, and TF A indirectly binds to promoter regions through physical binding to TF B . This is also referred to as indirect piggy-back binding .
PSSM sensitivity under indirect TF binding is important for our target gene prediction model. Instead of using a TF's own PSSM, the PSSM of another TF through which the TF binds should be used for more accurate predictions. Therefore, identifying those cases before using our model will be necessary to achieve better results.
Implications on gene expression regulation
We show in this study that incorporating chromatin modification information could substantially improve the prediction of TF target genes. In fact, chromatin modifications relate to gene expression regulation on two levels . First, chromatin is modified to form euchromatin, within which genes can be turned on and off, or heterochromatin, within which genes are silenced. Second, euchromatin is further modified by enzymes recruited by specific TFs to mark the 'on and off' status of transcription. We examined target versus non-target differential histone modification profiles for each individual TF, and observed TF-specific chromatin modifications marked in the target genes. Therefore, we suggest that chromatin modifications might function as both non-specific euchromatin marks and TF-specific regulatory marks. Our model takes advantage of the chromatin information from both of these levels.
However, the sequential order of chromatin modification and TF binding, in terms of time and causality, is still not quite clear. It is possible that one of them happens first and then drives the occurrence of the other. The other possibility is that the two events might interact in a feedback manner to regulate gene expression. More fine-tuned experiments in the future would be helpful for unraveling the time-dependent interaction between chromatin modification and TF binding.
Materials and methods
Chromatin modification data
The yeast histone modification data sets used in this study are basically from two sources. The first data set is from Pokholok et al.  and contains the profiles of 14 chromatin features under YPD or H2O2 conditions. These chromatin features include histone H3 and H4 occupation, H3K9ac, H3K14ac, H4K5ac8ac12ac16ac, H3K4me1, H3K4me2, H3K4me3, H3K36me3, and H3K79me3. The profiles of these features were measured by ChIP-chip experiments using over 40, 000 probes, which cover 85% of the yeast genome. We calculated the signal of each chromatin feature in the 1-kb upstream region of each ORF by averaging signals of all the probes within this region. Similarly, for each ORF the average signal of each feature in the 1-kb region downstream of the start codon was also calculated. We named these the upstream chromatin signal and downstream chromatin signal for ORFs, respectively.
The second data set is from Kurdistani et al. . These data contain levels of acetylation of 11 lysines in intergenic regions (IRs) as well as ORF regions. These profiles were also measured using ChIP-chip experiments. These 11 histone acetylations are H2AK7ac, H2BK11ac, H2bK16ac, H3K9ac, H3K14ac, H3K18ac, H3K23ac, H3K27ac, H4K8ac, H4K12ac and H4K16ac. We named the signals in IRs and ORFs as the IR chromatin signal and coding region chromatin signal for ORFs, respectively.
Target genes of yeast transcription factors
Target genes for 203 yeast TFs under various conditions (including YPD and H2O2) were identified using the ChIP-chip experiments by Harbison et al. . For each binding interaction, a probability score (P-value) was calculated, measuring the binding potential of a TF with the promoter region of a gene.
When TF target genes are determined according to ChIP-chip data, one needs to set a cutoff for P-values, which indicates the confidence of regulation of genes by TFs. A small (strict) P-value cutoff would result in fewer but more confident target genes, while a large (loose) P-value cutoff would do the opposite. For instance, there are 159 target genes for RAP1 using a cutoff value of 0.001, while the target gene number increases to 581 when a cutoff value of 0.05 is used. We therefore tested the influence of the P-value cutoff on our model performance. As shown in Figure 2a, our results indicate that a more stringent P-value cutoff (that is, a smaller target gene set) improves the prediction accuracy of our model. Moreover, at all cutoff values the models combining histone modification and PSSM data outperform the models using either of them alone. On the other hand, a more stringent cutoff results in less target genes. To ensure enough positive target genes for model training, we decided to use 0.01 as the P-value cutoff in our analysis.
Position-specific scoring matrices of transcription factors
Two sets of PSSMs for yeast TFs have been identified previously using different strategies [10, 18]. The first set was downloaded from , which was based on de novo motif finding in all yeast promoter sequences. The promoter DNA sequences (from the start codon of an ORF to 800 bp upstream) of all yeast ORFs were analyzed to identify enriched motifs using the AlignACE program . A total of 666 motifs were found, among which 51 can be associated with known yeast TFs. The occurrences and matching scores of these motifs in the promoter regions of all yeast genes were also provided by Beer and Tavazoie .
The second set of PSSMs was from , which is based on motif analysis of target promoters identified by the ChIP-chip experiment. Details on the motif discovery procedure can be found in . In brief, motifs for a TF were discovered by applying a suite of motif discovery programs to the intergenic sequences identified by the binding data for this factor. The resulting motifs were subsequently clustered, filtered and selected to give rise to a single PSSM that can best represent the motif of a factor. For some factors the above procedure failed to identify their motifs and in such cases motifs were derived from the literature or databases.
where i and j represent positions in PSSMs and four nucleotides, respectively. p i, j is then the weight at each PSSM position of each nucleotide, and p b is the background nucleotide frequency of the S. cerevisiae genome. Specifically, we use 37% as the GC content to calculate p b for each nucleotide.
where n is the length of PSSMs. If two PSSMs are of different lengths, we compare each possible alignment of the two PSSMs with no gap, and keep the maximum similarity from each alignment. The similarity score is in the range 0 to 1.
Searching promoters for known motifs
Given the list of PSSMs for TFs, we searched the promoters of all yeast genes for occurrences of these motifs using FIMO of the MEME suite . The promoter region was defined as the DNA region from the start codon to 800 bp upstream of an ORF. The cumulative matching score of the occurrences of a motif in the promoter region of a gene was calculated, which was subsequently used as a feature for predicting TF target genes.
Comparison of chromatin modifications between ChIP-chip-verified TFBSs and non-verified motif matching sites
We performed comparative analysis of chromatin modification differences between TFBSs verified by ChIP-chip experiments and those non-verified motif matching sites for TFs with available PSSMs. SWI4 is discussed here as an example. A list of binding sites (TFBSs) of the factor SWI4 was downloaded from the Saccharomyces Genome Database . This list contains 99 binding sites that are targeted by SWI4 in YPD medium according to the ChIP-chip data. We also collected a list of non-TFBSs by selecting DNA regions that are consistent with the SWI4 motif but not targeted by SWI4 in YPD medium as indicated by the ChIP-chip results (P > 0.4). These non-TFBSs were further filtered to ensure that there is no SWI4 TFBSs within the nearby 2-kb region, which ultimately resulted in 485 non-TFBSs for SWI4. All of the TFBSs and non-TFBSs are less than 20 bp in size. The levels of chromatin features in these TFBSs and non-TFBSs were calculated based on the intensities of probes covering them. Finally, the signal of the 14 chromatin features was compared between the TFBS and non-TFBS groups using the t-test.
Support vector machine model for transcription factor target prediction
For a TF in each gene we obtained the following features: cumulative matching score from motif searching, the upstream and downstream signal of 14 histone or histone modification profiles, and the IR and coding region signal of 11 histone acetylation profiles. Based on the ChIP-chip data from , genes were classified into target and non-target sets for a TF. All or subsets of these features were integrated using a SVM model  for predicting target genes of a TF. The radial kernel was used for training and predicting in the SVM classification model.
We evaluated the performances of the models using two-fold cross-validation. Specifically, we randomly split the data into two sets of equal sizes, a training set and a testing set. The model was then trained using the training set and applied to the testing set to predict target genes. The prediction power of the model was estimated based on the testing set. In general, the SVM model outputs a probability indicating how likely a gene is to be the target of a TF. By setting different cutoff values, we can balance the sensitivity (true positive rate) and specificity (true negative rate) of predictions of the model. The plot of the sensitivity versus '1 - specificity' is called the ROC, which can be used to show the classification accuracy of the SVM model. The AUC can be used to summarize the prediction power of the model. For each TF, we repeated this process 50 times and the average of the AUC values was calculated to represent the prediction accuracy.
Based on the above-described SVM classification method, we constructed and compared several different models for TF target prediction, each taking advantage of different features (histone modifications, PSSM information, gene expression levels) or a combination of these. First, we prepared several groups of histone modification features. Based on the Pokholok et al. data , we calculated the signals of 14 types of histone modification in the upstream (500 bp, 1, 000 bp) and the downstream (500 bp, 1, 000 bp) regions of the transcription start sites of all yeast genes. From the Kurdistani et al. data , we collected the signals of 11 types of histone modification in the intergenic and ORF regions. Second, we calculated the cumulative matching score of the occurrences of the PSSM of each TF in the promoter region of genes, resulting in the PSSM features. Third, these features together with the expression levels of genes were selected and used as the predictors to classify targets versus non-targets of TFs. For instance, in the histone + PSSM model for factor RAP1, the histone modification features and the PSSM feature of the RAP1 motif were selected as the predictors of the SVM. In 'Results', we mostly used the 1, 000-bp upstream and the 1, 000-bp downstream histone modification features as predictors, but the other groups of histone modification features were also investigated for comparison purposes. The SVM function in the R package 'e1071' was utilized to implement the models with default parameter settings. For the SVM the radial basis function kernel was used since it achieved the best performance according to the cross-validation results.
Clustering of TFs using target chromatin modification profile
For each TF, a target histone modification profile was calculated by averaging histone modification signals among all its targets. The 1-kb upstream chromatin modifications from Pokholok et al.  and intergenic chromatin modifications from Kurdistani et al.  were used. An unsupervised k-means clustering algorithm was performed to generate two TF clusters, histone-sensitive and -insensitive TFs, using their target histone modification profiles.
To understand the relative importance of each chromatin modification to target prediction, target-non-target differential histone modification profiles for TFs were calculated based on t-statistics. For each chromatin modification in a differential modification profile for a TF, modification signals for target genes and non-target genes were collected and a t-statistic calculated. The t-statistics in differential modification profiles indicated the directional significance of chromatin modifications to distinguish target genes.
Inferring interactions between transcription factors
The target genes identified by ChIP-chip experiments could be either direct or indirect targets of a TF. For example, if two TFs, A and B, interacte with each other, the ChIP-chip for A can potentially identify target genes of B as well. Conversely, the existence of TF B's motif would be informative for predicting target genes of TF A. We used the TF target prediction model with chromatin modifications and TFs' own PSSMs, and then compared the model's AUC values to those derived using models with chromatin modifications and other TFs' PSSMs. Models with improved AUC performances suggest better predictive power with PSSMs other than a TF's own PSSM. These cases might indicate interactions between TFs.
Application of previously reported methods
We ran all methods with their default parameter settings. Internal thresholding was turned off in all cases to report a full list of predictions with scores. PSSMs of ten TFs and upstream 1-kb DNA sequences for all annotated yeast S. cerevisiae ORFs were used as inputs to MCAST  and Cluster-Buster . Pairwise pre-aligned upstream 1-kb DNA sequences of all annotated S. cerevisiae and Saccharomyces paradoxus orthologous ORFs were used instead for running EEL  and Stubb . Predicted binding targets with respective scoring systems from the programs were collected for all ten TFs. ROC curves and AUCs were calculated based on the same thresholding scheme for all methods.
area under ROC curve
chromatin immunoprecipitation with microarray hybridization
chromatin immunoprecipitation with massively parallel DNA sequencing
open reading frame
positive predictive value
position-specific scoring matrix
receiver operator characteristic
support vector machine
transcription factor binding site.
We thank Dr Won, Kyoung-Jae and Dr Wang, Wei for providing 'Chromia' source code and suggestions about running the program. We also thank Robert Bjornson and Yale High Performance Computing for technology support. We acknowledge support from the NIH and from the AL Williams Professorship funds.
- Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, et al: Life with 6000 genes. Science. 1996, 274 (546): 563-547.Google Scholar
- Goffeau A: Four years of post-genomic life with 6, 000 yeast genes. FEBS Lett. 2000, 480: 37-41. 10.1016/S0014-5793(00)01775-0.PubMedView ArticleGoogle Scholar
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.PubMedView ArticleGoogle Scholar
- Li B, Carey M, Workman JL: The role of chromatin during transcription. Cell. 2007, 128: 707-719. 10.1016/j.cell.2007.01.015.PubMedView ArticleGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, et al: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104. 10.1038/nature02800.PubMedPubMed CentralView ArticleGoogle Scholar
- Kato M, Hata N, Banerjee N, Futcher B, Zhang MQ: Identifying combinatorial regulation of transcription factors and binding motifs. Genome Biol. 2004, 5: R56-10.1186/gb-2004-5-8-r56.PubMedPubMed CentralView ArticleGoogle Scholar
- Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E, et al: Genome-wide location and function of DNA binding proteins. Science. 2000, 290: 2306-2309. 10.1126/science.290.5500.2306.PubMedView ArticleGoogle Scholar
- Iyer VR, Horak CE, Scafe CS, Botstein D, Snyder M, Brown PO: Genomic binding sites of the yeast cell-cycle transcription factors SBF and MBF. Nature. 2001, 409: 533-538. 10.1038/35054095.PubMedView ArticleGoogle Scholar
- Horak CE, Snyder M: ChIP-chip: a genomic approach for identifying transcription factor binding sites. Methods Enzymol. 2002, 350: 469-483.PubMedView ArticleGoogle Scholar
- Robertson G, Hirst M, Bainbridge M, Bilenky M, Zhao Y, Zeng T, Euskirchen G, Bernier B, Varhol R, Delaney A, et al: Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing. Nat Methods. 2007, 4: 651-657. 10.1038/nmeth1068.PubMedView ArticleGoogle Scholar
- Johnson DS, Mortazavi A, Myers RM, Wold B: Genome-wide mapping of in vivo protein-DNA interactions. Science. 2007, 316: 1497-1502. 10.1126/science.1141319.PubMedView ArticleGoogle Scholar
- Stormo GD: DNA binding sites: representation and discovery. Bioinformatics. 2000, 16: 16-23. 10.1093/bioinformatics/16.1.16.PubMedView ArticleGoogle Scholar
- Frith MC, Li MC, Weng Z: Cluster-Buster: Finding dense clusters of motifs in DNA sequences. Nucleic Acids Res. 2003, 31: 3666-3668. 10.1093/nar/gkg540.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou Q, Wong WH: CisModule: de novo discovery of cis-regulatory modules by hierarchical mixture modeling. Proc Natl Acad Sci USA. 2004, 101: 12114-12119. 10.1073/pnas.0402858101.PubMedPubMed CentralView ArticleGoogle Scholar
- Gupta M, Liu JS: De novo cis-regulatory module elicitation for eukaryotic genomes. Proc Natl Acad Sci USA. 2005, 102: 7079-7084. 10.1073/pnas.0408743102.PubMedPubMed CentralView ArticleGoogle Scholar
- Sinha S, Liang Y, Siggia E: Stubb: a program for discovery and analysis of cis-regulatory modules. Nucleic Acids Res. 2006, 34: W555-559. 10.1093/nar/gkl224.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Loo P, Aerts S, Thienpont B, De Moor B, Moreau Y, Marynen P: ModuleMiner - improved computational detection of cis-regulatory modules: are there different modes of gene regulation in embryonic development and adult tissues?. Genome Biol. 2008, 9: R66-10.1186/gb-2008-9-4-r66.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie D, Cai J, Chia NY, Ng HH, Zhong S: Cross-species de novo identification of cis-regulatory modules with GibbsModule: application to gene regulation in embryonic stem cells. Genome Res. 2008, 18: 1325-1335. 10.1101/gr.072769.107.PubMedPubMed CentralView ArticleGoogle Scholar
- Rada-Iglesias A, Bajpai R, Swigut T, Brugmann SA, Flynn RA, Wysocka J: A unique chromatin signature uncovers early developmental enhancers in humans. Nature. 2011, 470: 279-283. 10.1038/nature09692.PubMedPubMed CentralView ArticleGoogle Scholar
- Beer MA, Tavazoie S: Predicting gene expression from sequence. Cell. 2004, 117: 185-198. 10.1016/S0092-8674(04)00304-6.PubMedView ArticleGoogle Scholar
- Berger SL: The complex language of chromatin regulation during transcription. Nature. 2007, 447: 407-412. 10.1038/nature05915.PubMedView ArticleGoogle Scholar
- Kurdistani SK, Tavazoie S, Grunstein M: Mapping global histone acetylation patterns to gene expression. Cell. 2004, 117: 721-733. 10.1016/j.cell.2004.05.023.PubMedView ArticleGoogle 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-527. 10.1016/j.cell.2005.06.026.PubMedView ArticleGoogle Scholar
- Won KJ, Ren B, Wang W: Genome-wide prediction of transcription factor binding sites using an integrated model. Genome Biol. 2010, 11: R7-10.1186/gb-2010-11-1-r7.PubMedPubMed CentralView ArticleGoogle Scholar
- Ernst J, Plasterer HL, Simon I, Bar-Joseph Z: Integrating multiple evidence sources to predict transcription factor binding in the human genome. Genome Res. 2010, 20: 526-536. 10.1101/gr.096305.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Cristianini N, Shawe-Taylor J: An Introduction to Support Vector Machines and other kernel-based learning methods. 2000View ArticleGoogle Scholar
- Pina B, Bruggemeier U, Beato M: Nucleosome positioning modulates accessibility of regulatory proteins to the mouse mammary tumor virus promoter. Cell. 1990, 60: 719-731. 10.1016/0092-8674(90)90087-U.PubMedView ArticleGoogle Scholar
- Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat Genet. 2004, 36: 900-905. 10.1038/ng1400.PubMedView ArticleGoogle Scholar
- Sekinger EA, Moqtaderi Z, Struhl K: Intrinsic histone-DNA interactions and low nucleosome density are important for preferential accessibility of promoter regions in yeast. Mol Cell. 2005, 18: 735-748. 10.1016/j.molcel.2005.05.003.PubMedView ArticleGoogle Scholar
- Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJ: Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005, 309: 626-630. 10.1126/science.1112178.PubMedView ArticleGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320: 1344-1349. 10.1126/science.1158441.PubMedPubMed CentralView ArticleGoogle Scholar
- The Fraenkel Motifs. [http://fraenkel.mit.edu/Harbison/release_v24]
- Ni L, Bruce C, Hart C, Leigh-Bell J, Gelperin D, Umansky L, Gerstein MB, Snyder M: Dynamic and complex transcription factor binding during an inducible response in yeast. Genes Dev. 2009, 23: 1351-1363. 10.1101/gad.1781909.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu H, Gerstein M: Genomic analysis of the hierarchical structure of regulatory networks. Proc Natl Acad Sci USA. 2006, 103: 14724-14731. 10.1073/pnas.0508637103.PubMedPubMed CentralView ArticleGoogle Scholar
- Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, 34: D535-539. 10.1093/nar/gkj109.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu WS, Li WH: Systematic identification of yeast cell cycle transcription factors using multiple data sources. BMC Bioinformatics. 2008, 9: 522-10.1186/1471-2105-9-522.PubMedPubMed CentralView ArticleGoogle Scholar
- Hirata D, Yano K, Miyakawa T: Stress-induced transcriptional activation mediated by YAP1 and YAP2 genes that encode the Jun family of transcriptional activators in Saccharomyces cerevisiae. Mol Gen Genet. 1994, 242: 250-256. 10.1007/BF00280413.PubMedView ArticleGoogle Scholar
- Schuller HJ, Richter K, Hoffmann B, Ebbert R, Schweizer E: DNA binding site of the yeast heteromeric Ino2p/Ino4p basic helix-loop-helix transcription factor: structural requirements as defined by saturation mutagenesis. FEBS Lett. 1995, 370: 149-152. 10.1016/0014-5793(95)00818-T.PubMedView ArticleGoogle Scholar
- Lefrancois P, Euskirchen GM, Auerbach RK, Rozowsky J, Gibson T, Yellman CM, Gerstein M, Snyder M: Efficient yeast ChIP-Seq using multiplex short-read DNA sequencing. BMC Genomics. 2009, 10: 37-10.1186/1471-2164-10-37.PubMedPubMed CentralView ArticleGoogle Scholar
- Bailey TL, Noble WS: Searching for statistically significant regulatory modules. Bioinformatics. 2003, 19 (Suppl 2): ii16-25. 10.1093/bioinformatics/btg1054.PubMedView ArticleGoogle Scholar
- Palin K, Taipale J, Ukkonen E: Locating potential enhancer elements by comparative genomics using the EEL software. Nat Protoc. 2006, 1: 368-374. 10.1038/nprot.2006.56.PubMedView ArticleGoogle Scholar
- Steinfeld I, Shamir R, Kupiec M: A genome-wide analysis in Saccharomyces cerevisiae demonstrates the influence of chromatin modifiers on transcription. Nat Genet. 2007, 39: 303-309. 10.1038/ng1965.PubMedView ArticleGoogle Scholar
- Gordan R, Hartemink AJ, Bulyk ML: Distinguishing direct versus indirect transcription factor-DNA interactions. Genome Res. 2009, 19: 2090-2100. 10.1101/gr.094144.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- The MEME Suite. [http://meme.sdsc.edu/meme/intro.html]
- Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe T, Schroeder M, et al: SGD: Saccharomyces Genome Database. Nucleic Acids Res. 1998, 26: 73-79. 10.1093/nar/26.1.73.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.