Skip to main content

Characterizing nucleosome dynamics from genomic and epigenetic information using rule induction learning

Abstract

Background

Eukaryotic genomes are packaged into chromatin, a compact structure containing fundamental repeating units, the nucleosomes. The mobility of nucleosomes plays important roles in many DNA-related processes by regulating the accessibility of regulatory elements to biological machineries. Although it has been known that various factors, such as DNA sequences, histone modifications, and chromatin remodelling complexes, could affect nucleosome stability, the mechanisms of how they regulate this stability are still unclear.

Results

In this paper, we propose a novel computational method based on rule induction learning to characterize nucleosome dynamics using both genomic and histone modification information. When applied on S. cerevisiae data, our method produced totally 98 rules characterizing nucleosome dynamics on chromosome III and promoter regions. Analyzing these rules we discovered that, some DNA motifs and post-translational modifications of histone proteins play significant roles in regulating nucleosome stability. Notably, these DNA motifs are strong determinants for nucleosome forming and inhibiting potential; and these histone modifications have strong relation with transcriptional activities, i.e. activation and repression. We also found some new patterns which may reflect the cooperation between these two factors in regulating the stability of nucleosomes.

Conclusion

DNA motifs and histone modifications can individually and, in some cases, cooperatively regulate nucleosome stability. This suggests additional insights into mechanisms by which cells control important biological processes, such as transcription, replication, and DNA repair.

Background

Genetic materials of eukaryotic organisms are packaged into chromatin inside cell nucleus. This compact structure has the form like a bead-on-string fiber containing fundamental repeating units, the nucleosomes. Each nucleosome is composed of 147 bp of DNA wrapped 1.65 turns around an octamer of histone proteins consisting of a central (H 3 - H 4)2 tetramer flanked on both side by two H 2A - H 2B dimers [1]. Since it was first recognized [2], there have been increasing evidences showing that chromatin plays a much more important role far beyond DNA compaction. By burying cis-regulatory elements under histone proteins and/or modifying related epigenetic information, chromatin imposes ubiquitous and profound effects on many DNA-based processes, including transcription, DNA repair and replication. To ensure faithfully copy both genetic and epigenetic information during replication or to facilitate the binding of Transcription Factors (TFs) to regulatory elements during transcription in the context of chromatin, cells have developed complicated biological pathways [3]. In these pathways, by regulating nucleosome stability cells can control the accessibility of underlying DNA sequences to biological machineries. For example, in replication, during the process known as parental histone segregation, pre-existing nucleosomes located ahead of replication forks are transiently disrupted from parental DNA strands and later transferred onto nascent DNA [3, 5]. In transcription, moving nucleosomes to different translational positions is known as one way to change the accessibility of nucleosomal DNA to TFs [4]. Also, promoter regions of actively transcribed genes are usually free of nucleosomes [7, 8]. So, understanding how cells regulate nucleosome stability will bring us additional insights into mechanisms of many important biological processes.

Nucleosome stability can be regulated by many factors, such as DNA sequences, histone modifications and histone variants, and chromatin remodelling complexes [9]. For example, DNA sequence is known as a reliable determinant for nucleosome preference, which can be used to predict nearly 50% of nucleosome positions [10], so it is likely to be an important factor in favouring or disfavouring nucleosome eviction. Histone variant H2A.Z (Htz1) is found to be preferentially enriched at promoters where some nucleosomes have to be quickly removed upon transcriptional activation [4]. Also, acetylated histones are shown to be easily dissociated from DNA [11, 12]. Chromatin remodelling complexes, such as Swi/Snf, act in concert with histone chaperones (e.g Asf1, Nap1) to displace histones from their original positions [4]. Although the complete list of factors has been fairly known, the mechanisms of how they act to mobilize nucleosome are still unclear.

Owing to recent advanced profiling techniques, such as ChIP-on-Chip and ChIP-Seq, we now have increasing amount of information about how nucleosomes and various kinds of histone modifications are distributed over the genomes of many organisms, including yeast, drosophila, and human [13, 8]. This opens up a chance for thorough investigation of nucleosome organization, its regulatory mechanisms and functions. Until now, there have been many works, both experimental and computational, concentrating on revealing the effects of factors stated above on nuclesome distribution [10, 13, 19, 20] but most of them have some common drawbacks. First, they mainly considered the effect of each factor separately while bypassing their combinatorial effects on nucleosome distribution. Second, although the distribution of destabilized nucleosomes is usually inhomogeneous throughout the genome and is known to have strong relation with transcriptional activities [13], it is still not well-characterized compared with that of stable nucleosomes.

There are several efforts trying to overcome these limitations. For example, Rippe et al. [21] and Schnitzler [22] investigated co-effects of DNA sequences and chromatin remodelling complexes; Widlund et al. [23] and Yang et al. [24] investigated co-effects of histone tails and DNA sequences on nucleosome distribution. Most of them, however, were based on experimental methods. More recently, Dai et al. [25] used both transcriptional interaction and genomic sequence information to computationally identify dynamic nucleosome distribution, but the number of works like this is still limit.

Enthused by these facts, in this paper, we propose a novel method for computationally characterizing nucleosome dynamics from both genomic sequences and histone modification profiles. Our method is based on induction rule learning adapted for subgroup discovery, which can discover sufficiently large and statistically meaningful subsets of population as shown in [26], so it is well suited for characterizing inhomogeneous distribution of destabilized nucleosomes. Moreover, by combining both genetic sequence and histone modification information, our method can discover the combinatorial nature of these two factors in regulating nucleosome stability. Our results on S. cerevisiae show that, some DNA motifs, which are reliable determinants for nucleosome forming/inhibiting potential, and post-translational modifications of histone proteins, which have strong relation with transcriptional activities, are likely to be more significant to nucleosome dynamics. We also found some patterns of cooperation between these DNA motifs and histone modifications in regulating nucleosome stability. Our results give additional insights into mechanisms of how cells regulate important biological processes, such as transcription, DNA repair and replication.

Results and discussion

Potentially significant motifs to nucleosome dynamics

DNA sequence has long been known to be a strong determinant for nucleosome formation potential, which can be used to identify nearly 50% of positioned nucleosomes in vivo, so it is likey to be an important factor affecting nuclesome stability. To determine DNA motifs which may be importantly related to nucleosome stability, two different approaches were applied (Section Method overview). In the first one, we used WordSpy [27] with the word length set to 6 to identify statistically significant motifs related to nucleosome states. The length of 6 was chosen because, as shown in some previous research [10, 19], nucleosome forming ability of DNA sequences may be decided mostly by short motifs, with length from 2 to 6. WordSpy uses dictionary-based approach so it is suitable to find short motifs among a group of DNA sequences [28]. Tables 1 and 2 show the 15 most significant motifs related to 2 states of nucleosomes found by WordSpy when run on chromosome III and promoter region data, respectively (complete lists are given in Additional File 1 and 2). The results show no big difference between important motifs of genetic regions and those of promoter regions. For example, both of them are enriched of dinucleotides TG/CA and this coincides with previous research [19], showing that TG/CA are highly flexible dinucleotides so they have large impact in imparting nucleosome forming ability. From the results given by WordSpy, it is difficult to identify motifs that may be important in discriminating nucleosome states. So, we used the second approach based on feature selection with Fisher criterion (Section Feature selection with Fisher criterion) to overcome this limitation. Table 3 shows 20 strongest discriminative motifs corresponding to chromosome III and promoter regions ranked by their F-score values (complete list is given in Additional File 3). Among them, dinucleotides are likely the most important motifs compared with the others in deciding nucleosome stability: 14 and 15 over 20 in chromosome III and promoter sequences, respectively. Moreover, among 10 strongest discriminative signals are AA/TT/AT/TA/CA/TG (for chromosome III) and AT/TT/CA/TG (for promoter regions), which are related with nucleosome forming (e.g. CA/TG) and inhibiting (e.g. AA/TT/AT/TA) potential of DNA sequences.

Table 1 Significant DNA motifs on chromosome III given by WordSpy
Table 2 Significant DNA motifs on promoter regions given by WordSpy
Table 3 Discriminative motifs ranked by F-scores

Significant histone modifications to nucleosome dynamics

Histone modification is one of the most important non-sequence regulatory factors of many chromatin-based processes and has also been known to affect nucleosome stability. To identify histone modifications potentially significant to nucleosome stability, we applied feature selection procedure, the same as what was done with DNA sequences, on the data of 12 different histone modifications corresponding to chromosome III and promoter regions (Section Data preparation). The result was ranked by F-score and given in Table 4. This result shows that, the first 9 modifications of chromosome III, including H3K14Ac/H4K5Ac/H3K4Me3/H4K12Ac/H3K4Me1/H3K9Ac/H2AK7Ac/H4K16Ac/H2BK16Ac, and the first 6 ones of promoter regions, including H3K4Me3/H3K9Ac/H3K18Ac/H4K16Ac/H4K12Ac/H4K8Ac, seem to be more important to nucleosome stability. Notably, all significant modifications in promoter regions are strongly related to transcriptional activation (e.g. H3K4Me3/H3K9Ac/H3K18Ac) and repression (e.g. H4K16Ac/H4K12Ac/H4K8Ac) [17, 18, 29]. That is also true with some significant modifications in chromosome III, where H3K4Me3/H3K9Ac and H4K12Ac/H4K16Ac/H2BK16Ac are known to have strong relation with transcriptional activation and repression, correspondingly.

Table 4 Histone modifications ranked by F-scores

Effects of DNA sequences and histone modifications on nucleosome dynamics

In order to see how DNA sequences and histone modifications affect nucleosome stability, we applied our method to the data containing significant DNA motifs and histone modifications identified above (Section Method overview). After filtering out uninteresting rules (Section Rule filtering), we received two sets of 60 rules (given in Additional File 4) and 38 rules (given in Additional File 5) characterizing nucleosome dynamics on chromosome III and promoter regions, correspondingly. Table 5 shows some selected rules from these rule sets. Analyzing these rules, we discovered that the enrichness of some specific DNA motifs has special impact on nucleosome stability. For example, nucleosomes bound by sequences enriched with AT/ATT/CTT are more stable (rules 1, 2, 6, 9, 10). This agrees with the result from [19], which said that sequences enriched with dinucleotides AT/TT have potential to inhibit nucleosome forming and deforming them on nucleosomes is more costly, so nucleosomes bound by these sequences may be more stable. Also, H3K9Ac/H3K18Ac/H3K4Me3 are known to have positive relation with transcriptional activation [17, 18, 29], so nucleosomes which are hyper-acetylated at H3K9/H3K18 and hyper-trimethylated at H3K4 seem to be more dynamic (rules 7, 8). In contrast, H4K12Ac is known to have positive relation with transcriptional repression [29], so H4K12 hyper-acetylated nucleosomes are more stable (rule 5) while H4K12 hypo-acetylated nucleosomes are more dynamic (rules 11, 12). However, there is no DNA pattern or post-translational modification showing dominant effect on nucleosome stability. Instead, there exist combinatorial effects, by DNA motifs themselves (rules 3, 4, 9) or by both DNA motifs and histone modifications (rules 2, 5, 7, 8, 10, 11, 12), on nucleosome stability. For example, if H3K4Me3 or H3K9Ac nucleosomes are located in regions enriched with ATT tri-nucleotide, they will become more stable (rules 2, 10); and even being located in regions enriched with AT dinucleotide, H4K12 hypo-acetylated nucleosomes still have potential of becoming unstable (rule 12). This agrees with the results from previous and recent works showing that the effects of histone acetylations depend on which lysines are acetylated and the locations of modified nuclesomes [30, 2]; and nucleosome positioning effect of DNA sequences is decided by the combination of nucleosome favouring and disfavouring motifs [19, 33].

Table 5 Selected rules characterizing nucleosome dynamics

Conclusion

Nucleosome dynamics plays important roles in many DNA-based processes and is regulated by many factors, such as DNA sequences, post-translational modifications of histone proteins, and chromatin remodelling complexes. However, most of the previous works only investigated the effect of individual factor while bypassing their combinatorial effects on the distribution of stable nucleosomes. In this paper, we proposed a novel method based on induction rule learning to computationally characterize nucleosome dynamics from both genomic and histone modification information. Our method is shown to be suitable for characterizing inhomogeneous distributions like that of destabilized nucleosomes; and by combining both genomic and histone modification information, it can discover potential co-effects of these two factors on nucleosome dynamics.

Our results on S. cerevisiae show that, some DNA motifs and histone modifications are more important in stabilizing and destabilizing nucleosomes. These DNA motifs and histone modifications are known to have strong relations with nucleosome forming/inhibiting potential and transcriptional activities, correspondingly. They not only act individually but also cooperate with each other by some specific patterns to combinatorially affect nucleosome stability.

Although our method is efficient in characterizing nucleosome dynamics, it produces a larger number of rules, of which many may be irrelevant. In the future, we need to develop a better method for filtering these uninteresting rules.

Methods

Data preparation

We used experimental data from Yuan et al. [13] and Liu et al. [17], which covered nearly 4% of yeast genome including chromosome III and 223 additional promoter regions, for our experiments. Data from Yuan contained 50-base DNA fragments tiled every 20 base pairs, and for each fragment we extracted its genomic sequence and HMM inferred state showing that it is nucleosomal sequence or not. Data extracted from Liu contained 12 different histone modification levels corresponding to DNA fragments above, including acetylations of H3K9, H3K14, H3K18, H4K5, H4K8, H4K12, H4K16, H2AK7, H2BK16 and mono-, di- and tri-methylations of H3K4. To investigate whether there exists any difference in characteristics of nucleosome dynamics between regulatory regions and genomic regions, we separated the data above into two datasets, corresponding to chromosome III and promoter regions. For each dataset, we filtered out data of linker regions to keep only nucleosomal data. Each nucleosome was assigned either as Well-positioned if it stretched from 6 to 8 fragments or as Delocalized if it stretched more than 9 fragments. Nucleosomes which had no histone modification values or delocalized nucleosomes whose lengths were longer than 350 base pairs were also treated as noise and removed. After these preprocessing steps, the dataset of chromosome III contained 997 well-positioned nucleosomes and 154 delocalized nucleosomes, the dataset of promoter regions contained 995 well-positioned nucleosomes and 69 delocalized nucleosomes. These two datasets were used for further analysis.

Method overview

In this work we aim at characterizing how DNA sequences and histone modifications affect nucleosome dynamics. To this end, we propose a novel method that takes significant DNA motifs and histone modifications along with nucleosome states as the input for the rule induction system to infer patterns which may represent the dependence of nucleosome stability on these two factors. Figure 1 depicts the overview of our method. At first, DNA motifs, which might be signigicantly related to nucleosome stability, were extracted from nucleosomal sequences by applying two different approaches. The first one was to find potentially conserved motifs related to nucleosome states using WordSpy, the software that has been shown to outperform other competing motif finding methods on benchmark datasets. The second one was to find motifs which could serve as discriminative information for two states of nucleosomes using feature selection function of Gist software package [34]. Motifs were ranked based on their important levels identified by Fisher criterion. Significant histone modifications were also extracted by applying the same feature selection procedure as the second approach above. We then constructed a decision table from these significant DNA motifs and histone modifications (see Figure 1) and used it as the input for CN2-SD rule induction system (Section Rule learning) to produce a set of rules. Some filtering procedures were applied to remove uninteresting rules and keep rules which may meaningfully characterize nucleosome dynamics.

Figure 1
figure 1

Method overview.

Feature selection with Fisher criterion

Feature selection is a process of selecting a subset of relevant features available from the data that most contribute to distinguishing instances from different classes. In our method, significant sequence and histone modification features related to two states of nucleosomes, Well-positioned and Delocalized, were identified and ranked by their Fisher scores (or F-score in short). This is one of statistical criteria that is simple, effective and independent of the choice of classification method. Because our method only concentrated on identifying features with highly discriminative strength instead of building any concrete classifiers so we chose F-score as the selection criterion. The discriminative strength of each feature is defined as following:

Given a dataset X with two classes, denote instances in class 1 as X1, and those in class 2 as X2. Assume is the average of the jth feature in Xk, the F-score of the jth feature is:

(1)

Where

(2)

The numerator indicates the discrimination between two classes, and the denominator indicates the scatter within each class. The larger the F-score is, the more likely this feature is more discriminative.

Rule learning

We consider this problem as a subgroup discovery problem and use a rule-based learning method for inducing rules. The problem of subgroup discovery can be defined as follows: given a population of individuals and a property of them, we are interested in finding population subgroups that are interesting with respect to the property of interest [26]. The induced rules usually have the form Cond → Class, where Class is a value of the property of interest, and Cond is a conjunction of attribute-value pairs selected from the features describing the training instances. In our work, Class has two values, Delocalized and Well-positioned. Attributes are significant histone modifications and DNA motifs as described above (Section Method overview).

Among several available rule induction systems, CN2 is a rule induction system implementing the separate-and-conquer strategy [35]. It learns a rule set by iteratively adding rules one at a time. Examples covered by a rule are removed from the search space before learning the next rule to add to the rule set. This is repeated until all examples are covered by at least one rule in the rule set or some stopping criteria is satisfied. Finally, CN2 can induce a set of independent rules, where each rule describes a specific subgroup of instances. However, CN2 only induces the first few rules discovered are usually interesting. Subsequently induced rules are obtained from biased example subsets, i.e., subsets including only positive examples that are not covered by previously induced rules. In 2004, Lavrac and her colleagues developed an improvement of CN2 for subgroup discovery, so-called CN2-SD [26]. The CN2-SD generalizes the covering algorithm by introducing example weights. Initially, all examples have a weight of 1.0. However, the weights of examples covered by a rule will not be set to 0 (they are not removed as in CN2), but instead will be reduced by a certain factor. The resulting number of rules is typically higher than with CN2, since most examples will be covered by more than one rule. CN2-SD is, therefore, better in learning local patterns, since the influence of previously covered patterns is reduced, but not completely ignored. In order to evaluate the rules with higher generality, CN2-SD also uses a weighted relative accuracy heuristic as presented in Equation 3. The weighted covering strategy tends to find rules that explain overlapped subgroups of instances in the search space, so the weighted relative accuracy heuristic produces highly general rules that express the knowledge contained in one specific subgroup. For these reasons, we utilize the CN2-SD in the rest of this paper for finding rules.

(3)

Rule filtering

Though the CN2-SD rule induction system uses a weighted covering strategy to restrict the redundancy of learned rules and guarantee the scanning of the whole search space, uninteresting rules are still produced [26, 36]. Let us assume that our rule r has a form: IF [Cond] THEN [ClassDistribution]. Where Cond = [motif1 = motifV al1 ∧ ... ∧ motif m = motifV al m ∧ histoneMod1 = hisV al1 ∧ ... ∧ histoneMod n = hisV al n ] with motif i is a DNA motif, motifV al i is enriched or low, histoneMod j is one kind of histone modification and hisV al j is hyper or neutral or hypo; ClassDistribution = [p, q] with p and q are the number of Well-positioned and Delocalized nucleosomes covered by r, respectively. We used several heuristics to filter out unexpected rules: rules that cover less than 2 positive examples or p/(p + q) < 0.8 if positive class is Delocalized and rules that cover less than 10 positive examples or q/(p + q) < 0.8 if positive class is Well-positioned (Positive class is the class characterized by the rule).

Note

Other papers from the meeting have been published as part of BMC Bioinformatics Volume 10 Supplement 15, 2009: Eighth International Conference on Bioinformatics (InCoB2009): Bioinformatics, available online at http://www.biomedcentral.com/1471-2105/10?issue=S15.

References

  1. Luger K, Mader AW, Richmond AK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particle at 2.8 A resolution. Nature. 1997, 389: 251-260. 10.1038/38444.

    Article  CAS  PubMed  Google Scholar 

  2. Kornberg RD, Thomas JO: Chromatin structure; oligomers of the histones. Science. 1974, 184 (139): 865-868. 10.1126/science.184.4139.865.

    Article  CAS  PubMed  Google Scholar 

  3. Groth A, Rocha W, Verreault A, Almouzni G: Chromatin Challenges during DNA Replication and Repair. Cell. 2007, 128 (4): 721-733. 10.1016/j.cell.2007.01.030.

    Article  CAS  PubMed  Google Scholar 

  4. Li B, Carey M, Workman JL: The Role of Chromatin during Transcription. Cell. 2007, 128 (4): 707-719. 10.1016/j.cell.2007.01.015.

    Article  CAS  PubMed  Google Scholar 

  5. Corpet A, Almouzni G: Making copies of chromatin: the challenge of nucleosomal organization and epigenetic information. Trends in Cell Biology. 2008, 19: 29-41. 10.1016/j.tcb.2008.10.002.

    Article  PubMed  Google Scholar 

  6. Probst AV, Dunleavy E, Almouzni G: Epigenetic inheritance during the cell cycle. Nature Reviews Molecular Cell Biology. 2009, 10: 192-206. 10.1038/nrm2640.

    Article  CAS  PubMed  Google Scholar 

  7. Lee CK, Shibata Y, Rao B, Strah BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nature Genetics. 2004, 36: 900-905. 10.1038/ng1400.

    Article  CAS  PubMed  Google Scholar 

  8. Henikoff S: Nucleosomes at active promoters: unforgettable loss. Cancer cell. 2007, 12 (5): 407-409. 10.1016/j.ccr.2007.10.024.

    Article  CAS  PubMed  Google Scholar 

  9. Henikoff S: Nucleosome destabilization in the epigenetic regulation of gene expression. Nature Reviews Genetics. 2008, 9: 15-26. 10.1038/nrg2206.

    Article  CAS  PubMed  Google Scholar 

  10. Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, Wang JPZ, Widom J: A genomic code for nucleosome positioning. Nature. 2006, 442 (7104): 772-778. 10.1038/nature04979.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Reinke H, Horz W: Histones Are First Hyperacetylated and Then Lose Contact with the Activated PHO5 Promoter. Molecular Cell. 2003, 11 (6): 1599-1607. 10.1016/S1097-2765(03)00186-2.

    Article  CAS  PubMed  Google Scholar 

  12. Zhao J, Diaz JH, Gross DS: Domain-Wide Displacement of Histones by Activated Heat Shock Factor Occurs Independently of Swi/Snf and Is Not Correlated with RNA Polymerase II Density. Molecular and Cellular Biology. 2005, 25 (20): 8985-8999. 10.1128/MCB.25.20.8985-8999.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. 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 (5734): 626-630. 10.1126/science.1112178.

    Article  CAS  PubMed  Google Scholar 

  14. Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: A high-resolution atlas of nucleosome occupancy in yeast. Nature Genetics. 2007, 39 (10): 1235-1244. 10.1038/ng2117.

    Article  CAS  PubMed  Google Scholar 

  15. Mavrich TN, Jiang C, Ioshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, Qi J, Glaser RL, Schuster SC, Gilmour DS, IstvanAlbert , Pugh BF: Nucleosome organization in the Drosophila genome. Nature. 2008, 453: 358-362. 10.1038/nature06929.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, Zeitlinger J, Lewitter F, Gifford DK, Young RA: Genome-wide map of nucleosome acetylation and methylation in yeast. Cell. 2005, 122 (4): 517-527. 10.1016/j.cell.2005.06.026.

    Article  CAS  PubMed  Google Scholar 

  17. Liu CL, Kaplan T, Kim M, Buratowski S, Schreiber SL, Friedman N, Rando OJ: Single-nucleosome mapping of histone modifications in S. cerevisiae. PLoS Biology. 2005, 3 (10): 10.1371/journal.pbio.0030328.

  18. Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-Resolution Profiling of Histone Methylations in the Human Genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.

    Article  CAS  PubMed  Google Scholar 

  19. Peckham HE, Thurman RE, Fu Y, Stamatoyannopoulos JA, Noble WS, Struhl K, Weng Z: Nucleosome positioning signals in genomic DNA. Genome Research. 2007, 17 (8): 1170-1177. 10.1101/gr.6101007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Zhang Y, Shin H, Song JS, Lei Y, Liu XS: Identifying Positioned Nucleosomes with Epigenetic Marks in Human from ChIP-Seq. BMC Genomics. 2008, 9: 537-10.1186/1471-2164-9-537.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Rippe K, Schrader A, Riede P, Strohner R, Lehmann E, Langst G: DNA sequence- and conformation-directed positioning of nucleosomes by chromatin-remodeling complexes. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (40): 15635-15640. 10.1073/pnas.0702430104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Schnitzler GR: Control of Nucleosome Positions by DNA Sequence and Remodeling Machines. Cell Biochemistry and Biophysics. 2008, 51 (2-3): 67-80. 10.1007/s12013-008-9015-6.

    Article  CAS  PubMed  Google Scholar 

  23. Widlund HR, Vitolo M, Thiriet C, Hayes JJ: DNA sequence-dependent contributions of core histone tails to nucleosome stability: differential effects of acetylation and proteolytic tail removal. Biochemistry. 2000, 39 (13): 3835-3841. 10.1021/bi991957l.

    Article  CAS  PubMed  Google Scholar 

  24. Yang Z, Zheng C, Hayes JJ: The core histone tail domains contribute to sequence-dependent nucleosome positioning. Journal of Biological Chemistry. 2007, 282 (11): 7930-7938. 10.1074/jbc.M610584200.

    Article  CAS  PubMed  Google Scholar 

  25. Dai Z, Dai X, Xiang Q, Feng J, Deng Y, Wang J, He C: Transcriptional interaction-assisted identification of dynamic nucleosome positioning. BMC Bioinformatics. 2009, 10 (Suppl l): S31-10.1186/1471-2105-10-S1-S31.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Lavrac N, Kavsek B, Flach P, Todorovski L: Subgroup discovery with CN2-SD. Journal of Machine Learning Research. 2004, 5: 153-188.

    Google Scholar 

  27. Wang G, Zhang W: A steganalysis-based approach to comprehensive identification and characterization of functional regulatory elements. Genome Biology. 2006, 7 (6):

  28. Das MK, Dai HK: A survey of DNA motif finding algorithms. BMC Bioinformatics. 2007, 8 (Suppl 7): S21-10.1186/1471-2105-8-S7-S21.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Kurdistani SK, Tavazoie S, Grunstein M: Mapping global histone acetylation patterns to gene expression. Cell. 2004, 117 (6): 721-733. 10.1016/j.cell.2004.05.023.

    Article  CAS  PubMed  Google Scholar 

  30. Hebbes TR, Thorne AW, Crane-Robinson C: A direct link between core histone acetylation and transcriptionally active chromatin. The EMBO Journal. 1988, 7 (5): 1395-1402.

    PubMed Central  CAS  PubMed  Google Scholar 

  31. Wang A, Kurdistani SK, Grunstein M: Requirement of Hos2 histone deacetylase for gene activity in yeast. Science. 2002, 298 (5597): 1412-1414. 10.1126/science.1077790.

    Article  CAS  PubMed  Google Scholar 

  32. de Nadal E, Zapater M, Alepuz PM, Sumoy L, Mas G, Posas F: The MAPK Hog1 recruits Rpd3 histone deacetylase to activate osmoresponsive genes. Nature. 2004, 427 (6972): 370-374. 10.1038/nature02258.

    Article  CAS  PubMed  Google Scholar 

  33. Jiang C, Pugh BF: Nucleosome positioning and gene regulation: advances through genomics. Nature Reviews Genetics. 2009, 10: 161-172. 10.1038/nrg2522.

    Article  CAS  PubMed  Google Scholar 

  34. Pavlidis P, Wapinski I, Noble WS: Support vector machine classification on the web. Bioinformatics. 2004, 20: 586-587. 10.1093/bioinformatics/btg461.

    Article  CAS  PubMed  Google Scholar 

  35. Clark P, Nibblet T: The CN2 induction algorithm. Machine Learning. 1989, 3: 261-283.

    Google Scholar 

  36. Pham TH, Clemente JC, Satou K, Ho TB: Computational discovery of transcriptional regulatory rules. Bioinformatics. 2005, 21: ii101-ii107. 10.1093/bioinformatics/bti1117.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to gratefully thank Prof. Nada Lavrac and Dr. Branko Kavsek for sharing CN2-SD software. The first and the third authors have been supported by Japanese Government Scholarship (Monbukagakusho) to study in Japan.

This article has been published as part of BMC Genomics Volume 10 Supplement 3, 2009: Eighth International Conference on Bioinformatics (InCoB2009): Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/10?issue=S3.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ngoc Tu Le.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NTL and TBH defined the research problem. NTL and DHT designed the experiment. NTL, TBH and DHT drafted the manuscript. All authors contributed to and approved the final version of the manuscript.

Electronic supplementary material

12864_2009_2585_MOESM1_ESM.txt

Additional file 1: The complete list of significant motifs on chromosome III given by WordSpy. This file contains statistically significant motifs, with lengths from 2 bp to 6 bp, given by WordSpy while run on chromosome III. (TXT 18 KB)

12864_2009_2585_MOESM2_ESM.txt

Additional file 2: The complete list of significant motifs on promoter regions given by WordSpy. This file contains statistically significant motifs, with lengths from 2 bp to 6 bp, given by WordSpy while run on promoter regions. (TXT 17 KB)

12864_2009_2585_MOESM3_ESM.txt

Additional file 3: The complete list of significantly discriminative motifs ranked by F-scores. This file contains motifs of chromosome III and promoter regions ranked by their discriminative powers based on F-score. (TXT 204 KB)

12864_2009_2585_MOESM4_ESM.txt

Additional file 4: Rules characterizing nucleosome dynamics on chromosome III. This file contains 60 rules characterizing nucleosome dynamics on chromosome III. (TXT 7 KB)

12864_2009_2585_MOESM5_ESM.txt

Additional file 5: Rules characterizing nucleosome dynamics on promoter regions. This file contains 38 rules characterizing nucleosome dynamics on promoter regions. (TXT 5 KB)

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Le, N.T., Ho, T.B. & Tran, D.H. Characterizing nucleosome dynamics from genomic and epigenetic information using rule induction learning. BMC Genomics 10 (Suppl 3), S27 (2009). https://doi.org/10.1186/1471-2164-10-S3-S27

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-S3-S27

Keywords