- Proceedings
- Open Access
A nucleosomal approach to inferring causal relationships of histone modifications
- Ngoc Tu Le^{1},
- Tu Bao Ho^{1},
- Bich Hai Ho^{2} and
- Dang Hung Tran^{3}
https://doi.org/10.1186/1471-2164-15-S1-S7
© Le et al.; licensee BioMed Central Ltd. 2014
- Published: 24 January 2014
Abstract
Motivation
Histone proteins are subject to various posttranslational modifications (PTMs). Elucidating their functional relationships is crucial toward understanding many biological processes. Bayesian network (BN)-based approaches have shown the advantage of revealing causal relationships, rather than simple cooccurrences, of PTMs. Previous works employing BNs to infer causal relationships of PTMs require that all confounders should be included. This assumption, however, is unavoidably violated given the fact that several modifications are often regulated by a common but unobserved factor. An existing non-parametric method can be applied to tackle the problem but the complexity and inflexibility make it impractical.
Results
We propose a novel BN-based method to infer causal relationships of histone modifications. First, from the evidence that nucleosome organization in vivo significantly affects the activities of PTM regulators working on chromatin substrate, hidden confounders of PTMs are selectively introduced by an information-theoretic criterion. Causal relationships are then inferred from a network model of both PTMs and the derived confounders. Application on human epigenomic data shows the advantage of the proposed method, in terms of computational performance and support from literature. Requiring less strict data assumptions also makes it more practical. Interestingly, analysis of the most significant relationships suggests that the proposed method can recover biologically relevant causal effects between histone modifications, which should be important for future investigation of histone crosstalk.
Keywords
- Bayesian Network
- Histone Modification
- Hide Variable
- Literature Support
- Conditional Mutual Information
Background
Genomes of higher organisms are organized into chromatin, a condensed structure of nucleosome units. Each of these units comprises a short piece of DNA wrapping around an octamer histone, containing two proteins of each type: H2A, H2B, H3, and H4 [1]. The histone protein is subject to various biochemical modifications, a.k.a. posttranslational modifications (PTMs), which have been shown to play crucial roles in many cellular processes, such as transcription and replication [2]. Defects of PTMs have also been implicated in determining cell fate and oncogenesis [3, 4]. The facts that PTMs may cause combinatorial effects on downstream events, and, by forming stable chromatin domains, properly pass modified states to the next generation [5, 6] suggest the existence of "histone codes" [7]. Therefore, revealing genome-wide PTM patterns and related functional implications would help increase our understanding of different DNA-mediated processes. For example, [8] discovered a common modification module of 17 modifications in human, suggesting their critical roles in gene regulation.
Advances in profiling techniques, such as ChIP-Chip and ChIP-Seq, have enabled the availability of genome-scale PTM data [8, 9], thus providing an unprecedented opportunity to decipher histone codes and their associated cis-regulatory elements. However, it also poses a great requirement for methods to understand such data. Many methods, ranging from clustering- to Hidden Markov Model (HMM)- to Bayesian network (BN)-based, have been developed to identify histone modifications patterns from ChIP-Chip and ChIP-Seq data [10–14]. Among them, BN-based approaches may help discover not only the cooccurrence but also the causal relationships of histone modifications [15]. This is especially important to understand histone crosstalk, a phenomenon that often occurs among different PTM events [16–18].
Bayesian network is a family of graphical models representing conditional independence of multiple variables [19]. First introduced to model gene regulatory networks (GRNs) from expression data [20], it has been widely used in reconstructing various biological networks, such as protein-protein interactions, protein signaling networks [21–23]. Likewise, there have been attempts to employ BNs to analyze histone modification data, in which compelled edges of the resulting models were considered causal relationships between PTMs [14, 24, 25]. Though useful, these works have a significant drawback: they require causal sufficiency assumption, i.e., all confounders of PTMs should be observed [26, 27]. This assumption, however, is unavoidably violated given the fact that some modifications can be regulated by enzymatic activity of a common but unobserved modifier [2].
In this work, we propose a novel BN-based method to infer causal relationships of PTMs that accounts for the existence of hidden confounders. First, an information-theoretic criterion is proposed to selectively introduce a pairwise hidden confounder (PHC) for each pair of PTMs. General hidden confounders (GHCs) are then derived from PHCs. The idea of deriving GHCs from PHCs has been presented in [31] to learn two-layer BNs with hidden variables. Differently, we based our approach on the evidence that chromatin in vivo imposes regulatory effects on the activities of PTM regulators. Thus, the criterion is proposed exploiting information about chromatin structure, i.e., nucleosome positioning. Matrix X is separately learned by a BN structure learning method. Compelled edges, i.e., causal relationships, are then derived from a network model of both PTMs and GHCs. Application on human epigenomic data of 38 histone modifications and histone variant H2A.Z, shows that the proposed method outperformed the non-parametric (Np) and the traditional one, which does not account for hidden confounders (noHidden), in terms of computational performance and literature support. Moreover, analysis of the most significant relationships shows that the proposed method can recover biologically relevant causal effects between histone modifications, such as H3K27Me3 → H3K9Me3, H3K4Me3 → H2AK5Ac, H4K8Ac → H2AZ. This is important for future investigation of histone crosstalk.
Methods
Information theory
where p(x, y), p(x), and p(y) are joint density function and marginal density functions of x and y, respectively.
where z = x − x_{ i }, d is the dimension of x, and Σ is the covariance matrix of z. When d = 1, equation (3) returns the estimated marginal density. When d = 2, it can be used to estimate the joint density function of bivariate variable (x, y). In our work, MI and CMI values were computed using a software package provided by [36].
Bayesian networks
Definition
in which p (x_{ i }| Pa_{ i }) corresponds to the local probability distribution of variable x_{ i }, and Pa_{ i } are x_{ i }'s parents.
D-separation property
BN structure learning
BN structure can be learned by score-based methods, aiming to identify the structure(s) that "best" describe the data. In this work, BDe score [37, 38] with uniform prior was used to measure the fitness of a candidate network. Because it is infeasible to search though all possible structures [39], greedy hill-climbing search combined with simulated annealing algorithm to avoid local maxima was employed.
Criterion for introducing PHCs
Proposition 1 Consider two PTMs, if each has its own (hidden) regulator, they will be d-separated given evidence on nucleosome positioning.
Proposition 2 Consider two PTMs, if they share a hidden regulator (in other words, a confounder), their d-separation property does not change upon the availability of nucleosome positioning evidence.
The results suggest that, given evidence on NucPos, the dependency level between two PTMs would not change if they share a hidden confounder, and would change (becoming "less" dependent) if each has its own (hidden) regulator. Using MI and CMI as the measures of dependency levels between two PTMs, we derive the following criterion for introducing a PHC for a pair of modifications, ptm1 and ptm2:
where α, β >0 are significant thresholds. These criteria will be used to derive PHCs for all pairs of PTMs.
Derivation of GHCs
From a set of PHCs derived in previous step, we define hidden confounder graph, an undirected graph whose nodes correspond to PTMs and edges to PHCs, implying that two nodes are connected if they share a PHC. Maximal clique algorithm is then applied on this graph, resulting in a set of maximal cliques, each corresponding to a GHC.
Causal relationship inference
To derive causal relationships of PTMs, we first combine BN received from structure learning step with GHCs, forming a network of PTMs and their hidden confounders. The edges among PTMs that share a GHC are then removed. Finally, the algorithm for finding compelled edges [26] is applied to the resulting structure, producing a set of compelled edges representing causal relationships of PTMs.
Data
Chromatin modification. CD4+ T cell data containing 20 methylations, 18 acetylations, and histone variant H2A.Z were retrieved from [9] and [8].
Nucleosome positioning data of resting CD4+ T cell was obtained from [45].
Gene set. UCSC Known Genes were retrieved from UCSC Genome Browser [46]. After removing genes with duplicated or without U133P2 probe IDs, 12456 genes were kept for analysis.
Results
Derivation of hidden confounders
Inference of PTM causal relationships
General scheme
BN structures were learned by Banjo (http://www.cs.duke.edu/~amink/software/banjo/), limited to 1, 300, 000 iterations because no significant improvement was achieved in further iteration (data not shown). The resulted structures were combined with 50 GHCs derived in previous step to produce a set of causal relationships. Significance scores were evaluated by bootstrapping method [20]. By which, original data was randomly bootstrapped N times, generating N bootstrapped datasets, and a set of causal relationships was derived for each. Significance score of each relationship was defined as the frequency of its appearance in N bootstrapped sets. In our experiment, N was set to 100.
For comparison, the implementation of Np by [30] was run on the same data. Because it only works with binary variables, the data were discretized into binary values by three schemes, using 70 (Scheme 1), 80 (Scheme 2), and 90 (Scheme 3) percentiles as thresholds. After receiving hidden confounders, the above procedure was employed to generate three sets of causal relationships, corresponding to each scheme.
Comparison
Performance
Performances of Np and our method.
Np(200 iterations) | hidden | |||
---|---|---|---|---|
Scheme1 | Scheme2 | Scheme3 | ||
Running time (sec.) | 5.0e+02 | 3.8e+02 | 2.8e+02 | 8.41 |
#Confounders | 22 | 30 | 17 | 50 |
Literature-based comparison
in which freq(ptm1, ptm2) is the frequency that both PTM terms appear together in PubMed abstracts, and freq(ptm_{ i }) is the frequency of each individually.
in which N is the number of the most significant shared concepts between two PTMs, sig_{1i}and sig_{2i}are the significant levels, assigned as point-wise mutual information values, of the associations between the ith shared concept and the two PTMs. All of these were retrieved through FACTA+ search with the list of the search terms given in supplementary information.
Comparison on the significance scores of three highly correlated PTM pairs reported in [11].
PTM pairs | hidden | noHidden | p − value |
---|---|---|---|
H3K27Ac-H3K4Me3 | 0.866 | 0.724 | 2.1e-10 |
H2AZ-H2BK120Ac | 0.002 | 0.002 | Nd |
H3K9Ac-H3K36Ac | 0.195 | 0.155 | Nd |
Comparison on the significance scores of 10 most correlated PTM pairs reported in [8].
PTM pairs | hidden | noHidden | p − value |
---|---|---|---|
H2BK5ac-H3K27ac | 0.677 | 0.481 | 6.26e-10 |
H2BK120ac-H2BK5ac | 0.594 | 0.301 | 6.11e-10 |
H2BK120ac-H4K91ac | 0.843 | 0.336 | 1.81e-15 |
H2BK5ac-H3K9ac | 0.524 | 0.416 | 2.36e-08 |
H3K79me2-H3K79me3 | 0.794 | 0.793 | Nd |
H2BK120ac-H3K27ac | 0.623 | 0.207 | 3.24e-13 |
H2BK120ac-H3K18ac | 0.61 | 0.196 | 1.55e-14 |
H3K18ac-H3K27ac | 0.453 | 0.19 | 1.28e-09 |
H2BK5ac-H3K18ac | 0.047 | 0.004 | 4.31e-08 |
H2BK5ac-H4K91ac | 0.294 | 0.283 | Nd |
Analysis and discussions
The top dominant histone modifications and significant Markov relations with corresponding dominance and significance scores (dScoreand C_{ 0 }, respectively) given by our method.
Modifications | dScore | Markov relations | C _{0} |
---|---|---|---|
H3K4Me3 | 4.473 | H3K23Ac → H3K14Ac | 1 |
H4K8Ac | 2.8836 | H4K8Ac → H2AZ | 1 |
H3K27Me1 | 2.7993 | H4K8Ac → H4K12Ac | 1 |
H3K27Ac | 2.6593 | H4K91Ac → H4K16Ac | 1 |
H4K5Ac | 2.3603 | H3K4Me2 → H3K23Ac | 1 |
H2BK120Ac | 2.2473 | H3K4Me3 → H2AZ | 1 |
H4K91Ac | 1.7744 | H3K4Me3 → H3K9Me1 | 1 |
H3K4Me1 | 1.6105 | H3R2Me1 → H3R2Me2 | 0.99 |
H3K27Me3 | 1.5533 | H3K27Me3 → H3K9Me3 | 0.98 |
H3K9Ac | 1.3325 | H3K27Me1 → H3K27Me3 | 0.96 |
Analyzing the top dominant modifications, we found that 8 out of 10 PTMs, {H3K4Me3, H3K27Ac, H2BK120Ac, H4K8Ac, H4K5Ac, H4K91Ac, H3K4Me1, H3K9Ac}, have been reported in the original research as important marks that appeared in the modification back-bone at promoters [8]. For the other two, H3K27Me3 is known as an important repressive mark, and H3K27Me1 as an active mark at promoters [9]. Interestingly, the result suggested the significant role of H2BK120Ac and its regulatory effect on H3K4Me3, an important modification mark of active promoters, through the chain H2BK120Ac → H3K18Ac → H3K4Me3. For a long time, the functions of H2B modifications, particularly H2BK120Ac, have remained obscure compared to other modifications [56]. Just recently there has been an indication that H2BK120Ac appears as an early modification mark in TSS regions and affects H2BK120Ub [57], a modification that regulates H3K4Me3 [58, 59], providing support for our finding. Investigation of the most significant Markov relations revealed that well-characterized modifications are mostly functionally related. For example, the N-terminal tail of histone H4 has four acetylated lysines: K5, K8, K12, K16, of which H4 K5Ac/K8Ac/K12Ac play a non-specific, cumulative regulatory role different from that of H4K16Ac [60]. In consistence with this observation, these modifications were predicted to be closely linked and separated from H4K16Ac in the resulting model: H4K5Ac → H4K8Ac, H4K5Ac → H4K12Ac, and H4K8Ac → H4K12Ac (one of the top 10 Markov relations). For other less well-known modifications, such as H3R2 methylations or H3K27 mono-methylation, the links might suggest novel biological understanding. While the relationship between H3R2Me1 → H3R2Me2 might reflect a directional equilibrium between mono- and di-methyl H3R2, the one between H3K27Me1 → H3K27Me3 might reflect their functional association through G9a methyltransferase, as recently reported by [61]. More interestingly, 4 out of 10 most significant Markov relations have already been reported to be causal in literature. [14] have shown evidences for causal relationships of H3K27Me3 → H3K9Me3 and H3K4Me3 → H2AZ. In [62], H3K9Me1/2 was shown to be demethylated by P HD finger protein 8 (PHF8), whose catalytic activity is in turn stimulated by H3K4Me3, suggesting the causal effect of H3K4Me3 on H3K9Me1, represented by the link H3K4Me3 → H3K9Me1. Also, the deposition of histone variant H2A.Z by SWR1 complex is known to be triggered by NuA4-mediated acetylation of histone H4 [63, 64]. Our model supported this observation with the relationship H4K8Ac → H2AZ. Additionally, causal effects have also been observed to support other relationships of the resulting model. For example, [14] have given evidence for the relationship H3K4Me3 → H3K36Me3. [65] have reported that the recruitment of MLL1, a histone methyltransferase responsible for H3K4 methylation, is required for the binding of TIP60 histone acetyltransferase, which catalytically acetylates H2AK5. In agreement, our model predicted the relationship H3K4Me3 → H2AK5Ac, suggesting causal effect of H3K4Me3 on H2AK5Ac.
Conclusion
Elucidation of functional relationships among histone modifications is crucial to understanding important chromatin-mediated processes. Previous BN-based approaches, however, have not taken into account the existence of hidden regulators when inferring causal relationships of PTMs. We tackled the problem by proposing a novel approach that exploits chromatin organizational information to capture the effect of PTM hidden regulators. Application on human epigenomic data showed the advantage of the proposed method over the previous ones. Moreover, it could recover biologically relevant causal relationships between histone modifications, which may be useful for future investigation of histone crosstalk.
Declarations
Acknowledgements
This work was supported by Vietnam's National Foundation for Science and Technology Development (NAFOSTED Project No.102.01−2011.05).
Declarations
The publication costs for this article were funded by Vietnam's National Foundation for Science and Technology Development (NAFOSTED Project No.102.01−2011.05).
This article has been published as part of BMC Genomics Volume 15 Supplement 1, 2014: Selected articles from the Twelfth Asia Pacific Bioinformatics Conference (APBC 2014): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/15/S1.
Authors’ Affiliations
References
- Luger K, Ku M, Jaffe DB, Issac B: Crystal structure of the nucleosome core particle at 2.8 A resolution. Nature. 1997, 389: 251-60. 10.1038/38444.PubMedView ArticleGoogle Scholar
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128 (4): 693-705. 10.1016/j.cell.2007.02.005.PubMedView ArticleGoogle Scholar
- Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, Lee W, Mendenhall E, O'Donovan A, Presser A, Russ C, Xie A, Meissner A, Wernig M, Jaenisch R, Nusbaum C, Lander ES, Bernstein BE: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448 (7151): 553-60.PubMedPubMed CentralView ArticleGoogle Scholar
- Fraga MF, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Loss of acetylation at Lys16 and trimethylation at Lys20 of histone H4 is a common hallmark of human cancer. Nature Genet. 2005, 37: 391-400. 10.1038/ng1531.PubMedView ArticleGoogle Scholar
- Bernstein BE, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : The mammalian epigenome. Cell. 2007, 128 (4): 669-81. 10.1016/j.cell.2007.01.033.PubMedView ArticleGoogle Scholar
- Hawkins RD, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Distinct epigenomic landscapes of pluripotent and lineage-committed human cells. Cell. 2010, 6 (5): 479-91.Google Scholar
- Strahl BD, Allis CD: The language of covalent histone modifications. Nature. 2000, 403 (6765): 41-5. 10.1038/47412.PubMedView ArticleGoogle Scholar
- Wang Z, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Combinatorial patterns of histone acetylations and methylations in the human genome. Nature Genet. 2008, 40 (7): 897-903. 10.1038/ng.154.PubMedPubMed CentralView ArticleGoogle Scholar
- Barski A, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-37. 10.1016/j.cell.2007.05.009.PubMedView ArticleGoogle Scholar
- Hon G, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : ChromaSig: a probabilistic approach to finding common chromatin signatures in the human genome. PLoS Comput Biol. 2008, 4 (10): e1000201-10.1371/journal.pcbi.1000201.PubMedPubMed CentralView ArticleGoogle Scholar
- Ucar D, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Combinatorial chromatin modification patterns in the human genome revealed by subspace clustering. Nucleic Acids Res. 2011, 39 (10): 4063-75. 10.1093/nar/gkr016.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaskchek R, Tanay A: Spatial clustering of multivariate genomic and epigenomic information. Proceedings of the 13th Annual International Conference on Research in Computational Molecular Biology, Springer-Verlag. 2009, 170-183.Google Scholar
- Ernst J, Kellis M: Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat Biotechnol. 2010, 28 (8): 817-25. 10.1038/nbt.1662.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu H, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Inferring causal relationships among different histone modifications and gene expression. Genome Res. 2008, 18 (8): 1314-1324. 10.1101/gr.073080.107.PubMedPubMed CentralView ArticleGoogle Scholar
- Hawkins RD, Hon GC, Ren B: Next-generation genomics: an integrative approach. Nat Rev Genet. 2010, 11 (7): 476-86.PubMedPubMed CentralGoogle Scholar
- Berger SL: The complex language of chromatin regulation during transcription. Nature. 2007, 447: 407-12. 10.1038/nature05915.PubMedView ArticleGoogle Scholar
- Latham JA, Dent SYR: Cross-regulation of histone modifications. Nat Struct Mol Biol. 2007, 14 (11): 1017-24. 10.1038/nsmb1307.PubMedView ArticleGoogle Scholar
- Suganuma T, Workman JL: Crosstalk among Histone Modifications. Cell. 2008, 135 (4): 604-7. 10.1016/j.cell.2008.10.036.PubMedView ArticleGoogle Scholar
- Jensen FV, Nielsen TD: Bayesian Networks and Decision Graphs (second edition). 2001, New York: Springer-VerlagView ArticleGoogle Scholar
- Friedman N, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Using Bayesian networks to analyze expression data. J Comput Biol. 2000, 7 (3-4): 601-20. 10.1089/106652700750050961.PubMedView ArticleGoogle Scholar
- Jansen R, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : A Bayesian Networks Approach for Predicting Protein-Protein Interactions from Genomic Data. Science. 2003, 302 (5644): 449-53. 10.1126/science.1087361.PubMedView ArticleGoogle Scholar
- Steensel BV, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Bayesian network analysis of targeting interactions in chromatin. Genome Res. 2009, 20 (2): 190-200.PubMedView ArticleGoogle Scholar
- Sachs K, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Causal protein-signaling networks derived from multiparameter single-cell data. Science. 2005, 308 (5127): 523-9.PubMedView ArticleGoogle Scholar
- Lv J, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Discovering Cooperative Relationships of Chromatin Modifications in Human T Cells Based on a Proposed Closeness Measure. PLoS One. 2010, 5 (12): e14219-10.1371/journal.pone.0014219.PubMedPubMed CentralView ArticleGoogle Scholar
- Le NT, Ho TB: Reconstruction of histone modification network from next-generation sequencing data. IEEE International Conference on Bioinformatics and Bioengineering. 2011Google Scholar
- Chickering DM: A transformational characterization of equivalent Bayesian network structures. Proceedings of 11th Conference on Uncertainty in Artificial Intelligence, Morgan Kaufmann Publishers. 1995, 87-89.Google Scholar
- Zhang J, Spirtes P: A Transformational Characterization of Markov Equivalence between DAGs with Latent Variables. Proceedings of 21th Conference on Uncertainty in Artificial Intelligence, UAUI Press. 2005, 667-674.Google Scholar
- Koller D, Friedman N: Probabilistic Graphical Models: Principles and Techniques. MIT Press. 2009Google Scholar
- Bannister AJ, Kouzarides T: Regulation of chromatin by histone modifications. Cell Res. 2011, 21 (3): 381-95. 10.1038/cr.2011.22.PubMedPubMed CentralView ArticleGoogle Scholar
- Wood F, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : A Non-Parametric Bayesian Method for Inferring Hidden Causes. Proceedings of the 22nd Conference on Uncertainty in Artificial Intelligence. 2006, 536-543.Google Scholar
- Martin JD, VanLehn K: Discrete factor analysis: Learning hidden variables in Bayesian networks. Technical report, Department of Computer Science, University of Pittsburgh. 1995Google Scholar
- Basso K, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Reverse engineering of regulatory networks in human B cells. Nat Genet. 2005, 37 (4): 382-90. 10.1038/ng1532.PubMedView ArticleGoogle Scholar
- Margolin AA, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics. 2006, 7 (S7):Google Scholar
- Meyer PE, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : minet: A R/Bioconductor package for inferring large transcriptional networks using mutual information. BMC Bioinformatics. 2008, 9 (461):Google Scholar
- Zhang X, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Inferring gene regulatory networks from gene expression data by path consistency algorithm based on conditional mutual information. Bioinformatics. 2012, 28: 98-104. 10.1093/bioinformatics/btr626.PubMedView ArticleGoogle Scholar
- Peng H, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Feature selection based on mutual information criteria of max-dependency, max-relevance, and min-redundancy. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2005, 27 (8): 1226-38.PubMedView ArticleGoogle Scholar
- Heckerman D, Geiger D, Chickering DM: Learning Bayesian networks: The combination of knowledge and statistical data. Machine Learning. 1995, 20: 197-243.Google Scholar
- Cooper G, Herskovits E: A Bayesian method for the induction of probabilistic networks from data. Machine Learning. 1992, 9: 309-47.Google Scholar
- Chickering DM: Learning Bayesian Networks is NP hard. Technical report, Redmond, WA: Microsoft Research. 1994Google Scholar
- Brown CE, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : The many HATs of transcription coactivators. Trends Biochem Sci. 2000, 25: 15-9. 10.1016/S0968-0004(99)01516-9.PubMedView ArticleGoogle Scholar
- Hager GL, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Influence of chromatin structure on the binding of transcription factors to DNA. Cold Spring Harb Symp Quant Biol. 1993, 58: 63-71. 10.1101/SQB.1993.058.01.010.PubMedView ArticleGoogle Scholar
- Kwon H, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Nucleosome disruption and enhancement of activator binding by a human SW1/SNF complex. Nature. 1994, 370 (6489): 477-81. 10.1038/370477a0.PubMedView ArticleGoogle Scholar
- Floer M, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : A RSC/Nucleosome Complex Determines Chromatin Architecture and Facilitates Activator Binding. Cell. 2010, 141 (3): 407-18. 10.1016/j.cell.2010.03.048.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang X, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Nucleosomes and the accessibility problem. Trends Genet. 2011, 27 (12): 487-92. 10.1016/j.tig.2011.09.001.PubMedView ArticleGoogle Scholar
- Schones DE, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132 (5): 887-98. 10.1016/j.cell.2008.02.022.PubMedView ArticleGoogle Scholar
- Karolchik D, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004, 132: D493-6.View ArticleGoogle Scholar
- Butte AJ, Kohane IS: Mutual Information Relevance Networks: Functional Genomic Clustering Using Pairwise Entropy Measurements. Pacific Symposium on Biocomputing. 2000, 5: 415-26.Google Scholar
- Tie F, A , B : CBP-mediated acetylation of histone H3 lysine 27 antagonizes Drosophila Polycomb silencing. Development. 2009, 136 (18): 3131-41. 10.1242/dev.037127.PubMedPubMed CentralView ArticleGoogle Scholar
- Morris SA, Banerjee R, Stratton CA: Identification of histone H3 lysine 36 acetylation as a highly conserved histone modification. J Biol Chem. 2007, 282 (10): 7632-40.PubMedPubMed CentralView ArticleGoogle Scholar
- Guillemette B, Banerjee R, Stratton CA: H3 lysine 4 is acetylated at active gene promoters and is regulated by H3 lysine 4 methylation. PLoS Genet. 2011, 7 (3): e1001354-10.1371/journal.pgen.1001354.PubMedPubMed CentralView ArticleGoogle Scholar
- Jenssen TK, Higuchi T, Goto T, Jaakkola TS: A literature network of human genes for high-throughput analysis of gene expression. Nat Genet. 2001, 28: 21-28.PubMedGoogle Scholar
- Djebbari A, Quackenbush J: Seeded Bayesian Networks: Constructing genetic networks from microarray data. BMC Systems Biology. 2008, 2 (57):Google Scholar
- Steele E, Tucker A, Schuemie PHM: Literature-based priors for gene regulatory networks. Bioinformatics. 2009, 25 (14): 1768-17748. 10.1093/bioinformatics/btp277.PubMedView ArticleGoogle Scholar
- Izzo A, Schneider R: Chatting histone modifications in mammals. Brief Funct Genomics. 2010, 9 (5-6): 429-43. 10.1093/bfgp/elq024.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsuruoka Y, Higuchi T, Goto T, Jaakkola TS: Discovering and visualizing indirect associations between biomedical concepts. Bioinformatics. 2011, 27 (13): i111-i119. 10.1093/bioinformatics/btr214.PubMedPubMed CentralView ArticleGoogle Scholar
- Wyrick JJ, Parra MA: The role of histone H2A and H2B post-translational modifications in transcription: a genomic perspective. Biochim Biophys Acta. 2009, 1789: 37-44. 10.1016/j.bbagrm.2008.07.001.PubMedView ArticleGoogle Scholar
- Gatta R, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : An acetylation-mono-ubiquitination switch on lysine 120 of H2B. Epigenetics. 2011, 6 (5): 630-7. 10.4161/epi.6.5.15623.PubMedView ArticleGoogle Scholar
- Lee JS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Histone crosstalk between H2B monoubiquitination and H3 methylation mediated by COMPASS. Cell. 2007, 131 (6): 1084-96. 10.1016/j.cell.2007.09.046.PubMedView ArticleGoogle Scholar
- Shukla A, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Histone methylation and ubiquitination with their cross-talk and roles in gene expression and stability. Cell Mol Life Sci. 2009, 66 (8): 1419-33. 10.1007/s00018-008-8605-1.PubMedView ArticleGoogle Scholar
- Dion MF, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Genomic characterization reveals a simple histone h4 acetylation code. Proc Natl Acad Sci USA. 2005, 102 (15): 5501-6. 10.1073/pnas.0500136102.PubMedPubMed CentralView ArticleGoogle Scholar
- Yoo K, Hennighausen L: EZH2 methyltransferase and H3K27 methylation in breast cancer. Int J Biol Sci. 2012, 8: 59-65.PubMedPubMed CentralView ArticleGoogle Scholar
- Feng W, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : PHF8 activates transcription of rRNA genes through H3K4me3 binding and H3K9me1/2 demethylation. Nat Struct Mol Biol. 2010, 17 (4): 445-50. 10.1038/nsmb.1778.PubMedView ArticleGoogle Scholar
- Altaf M, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : NuA4-dependent Acetylation of Nucleosomal Histones H4 and H2A Directly Stimulates Incorporation of H2A.Z by the SWR1 Complex. J Biol Chem. 2010, 285 (21): 15966-77. 10.1074/jbc.M110.117069.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou BO, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : SWR1 complex poises heterochromatin boundaries for antisilencing activity propagation. Mol Cell Biol. 2010, 30 (10): 2391-400. 10.1128/MCB.01106-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeong KW, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, P : Recognition of enhancer element specific histone methylation by TIP60 in transcriptional activation. Nat Struct Mol Biol. 2011, 18 (12): 1358-65. 10.1038/nsmb.2153.PubMedPubMed CentralView ArticleGoogle Scholar
Copyright
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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.