Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A. 2010;107:2926-2931. [PMID: 20133639 DOI: 10.1073/pnas.0909344107]

Down-regulation of the histone methyltransferase EZH2 contributes to the epigenetic programming of decidualizing human endometrial stromal cells

Differentiation of human endometrial stromal cells (HESC) into decidual cells represents a highly coordinated process essential for embryo implantation. We show that decidualizing HESC down-regulate the histone methyltransferase enhancer of Zeste homolog 2 (EZH2), resulting in declining levels of trimethylation of histone 3 on lysine 27 (H3K27me3) at the proximal promoters of key decidual marker genes PRL and IGFBP1. Loss of H3K27me3 was associated with a reciprocal enrichment in acetylation of the same lysine residue, indicating active remodeling from repressive to transcriptionally permissive chromatin. Chromatin immunoprecipitation coupled with DNA microarray analysis demonstrated that decidualization triggers genome-wide changes in H3K27me3 distribution that only partly overlap those observed upon EZH2 knockdown in undifferentiated HESC. Gene ontology revealed that gain of the repressive H3K27me3 mark in response to decidualization and upon EZH2 knockdown in undifferentiated cells was enriched at the promoter regions of genes involved in transcriptional regulation and growth/cell proliferation, respectively. However, loss of the H3K27me3 mark (indicating increased chromatin accessibility) in decidualizing cells and upon EZH2 knockdown occurred at selective loci enriched for genes functionally implicated in responses to stimulus. In agreement, EZH2 knockdown in undifferentiated HESC was sufficient to augment the induction of decidual marker genes in response to cyclic AMP and progesterone signaling. Thus, loss of EZH2-dependent methyltransferase activity in the endometrium is integral to the process of chromatin remodeling that enables the transition from a proliferative to a decidual phenotype in response to differentiation cues.

The dynamics of the nucleosome: thermal effects, external forces and ATP

With nucleosomes being tightly associated with the majority of eukaryotic DNA, it is essential that mechanisms are in place that can move nucleosomes 'out of the way'. A focus of current research comprises chromatin remodeling complexes, which are ATP-consuming protein complexes that, for example, pull or push nucleosomes along DNA. The precise mechanisms used by those complexes are not yet understood. Hints for possible mechanisms might be found among the various spontaneous fluctuations that nucleosomes show in the absence of remodelers. Thermal fluctuations induce the partial unwrapping of DNA from the nucleosomes and introduce twist or loop defects in the wrapped DNA, leading to nucleosome sliding along DNA. In this minireview, we discuss nucleosome dynamics from two angles. First, we describe the dynamical modes of nucleosomes in the absence of remodelers that are experimentally fairly well characterized and theoretically understood. Then, we discuss remodelers and describe recent insights about the possible schemes that they might use.

Epigenetics and the developmental origins of lung disease

The developmental origins of disease hypothesis have recently been expanded to include the early origins of lung disease, particularly early events that alter lung development. Intrauterine growth restriction (IUGR), preterm birth with the need for prolonged mechanical ventilation, and maternal tobacco smoke (MTS) or nicotine exposure produce neonatal and adult lung disease. These perinatal insults are characterized by alterations in alveolar formation and changes in the expression of genes that regulate alveolarization, including IGF1 and PPARγ. A potential mechanism for such changes in gene expression is epigenetics. IGF1 and PPARγ have altered epigenetic states in response to these perinatal insults. Identification of the specific epigenetic mechanisms involved in the developmental origin of lung disease may facilitate identification of molecular biomarkers with the potential to personalize respiratory disease risk assessment and treatment. The purpose of this review is to summarize what is known about the developmental origins of lung disease, the epigenetic contributions to lung disease, and areas that need further investigation.

Tight associations between transcription promoter type and epigenetic variation in histone positioning and modification

Background

Transcription promoters are fundamental genomic cis-elements controlling gene expression. They can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. Eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. For instance, it has been recently shown that histone with H3K9 acetylation (H3K9ac) is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone H3K9 and promoter architecture.

Results

Here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. We found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. Specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. This feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. More strikingly, the overall expression levels of genes associated with broad promoters (but not peak promoters) with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones.

Conclusion

These results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture.

CBP is required for environmental enrichment-induced neurogenesis and cognitive enhancement

The epigenetic changes of the chromatin represent an attractive molecular substrate for adaptation to the environment. We examined here the role of CREB-binding protein (CBP), a histone acetyltransferase involved in mental retardation, in the genesis and maintenance of long-lasting systemic and behavioural adaptations to environmental enrichment (EE). Morphological and behavioural analyses demonstrated that EE ameliorates deficits associated to CBP deficiency. However, CBP-deficient mice also showed a strong defect in environment-induced neurogenesis and impaired EE-mediated enhancement of spatial navigation and pattern separation ability. These defects correlated with an attenuation of the transcriptional programme induced in response to EE and with deficits in histone acetylation at the promoters of EE-regulated, neurogenesis-related genes. Additional experiments in CBP restricted and inducible knockout mice indicated that environment-induced adult neurogenesis is extrinsically regulated by CBP function in mature granule cells. Overall, our experiments demonstrate that the environment alters gene expression by impinging on activities involved in modifying the epigenome and identify CBP-dependent transcriptional neuroadaptation as an important mediator of EE-induced benefits, a finding with important implications for mental retardation therapeutics.

Quantification of histone modification ChIP-seq enrichment for data mining and machine learning applications

Background

The advent of ChIP-seq technology has made the investigation of epigenetic regulatory networks a computationally tractable problem. Several groups have applied statistical computing methods to ChIP-seq datasets to gain insight into the epigenetic regulation of transcription. However, methods for estimating enrichment levels in ChIP-seq data for these computational studies are understudied and variable. Since the conclusions drawn from these data mining and machine learning applications strongly depend on the enrichment level inputs, a comparison of estimation methods with respect to the performance of statistical models should be made.

Results

Various methods were used to estimate the gene-wise ChIP-seq enrichment levels for 20 histone methylations and the histone variant H2A.Z. The Multivariate Adaptive Regression Splines (MARS) algorithm was applied for each estimation method using the estimation of enrichment levels as predictors and gene expression levels as responses. The methods used to estimate enrichment levels included tag counting and model-based methods that were applied to whole genes and specific gene regions. These methods were also applied to various sizes of estimation windows. The MARS model performance was assessed with the Generalized Cross-Validation Score (GCV). We determined that model-based methods of enrichment estimation that spatially weight enrichment based on average patterns provided an improvement over tag counting methods. Also, methods that included information across the entire gene body provided improvement over methods that focus on a specific sub-region of the gene (e.g., the 5' or 3' region).

Conclusion

The performance of data mining and machine learning methods when applied to histone modification ChIP-seq data can be improved by using data across the entire gene body, and incorporating the spatial distribution of enrichment. Refinement of enrichment estimation ultimately improved accuracy of model predictions.

Tissue-specific prediction of directly regulated genes

Direct binding by a transcription factor (TF) to the proximal promoter of a gene is a strong evidence that the TF regulates the gene. Assaying the genome-wide binding of every TF in every cell type and condition is currently impractical. Histone modifications correlate with tissue/cell/condition-specific ('tissue specific') TF binding, so histone ChIP-seq data can be combined with traditional position weight matrix (PWM) methods to make tissue-specific predictions of TF-promoter interactions.

RESULTS

We use supervised learning to train a naïve Bayes predictor of TF-promoter binding. The predictor's features are the histone modification levels and a PWM-based score for the promoter. Training and testing uses sets of promoters labeled using TF ChIP-seq data, and we use cross-validation on 23 such datasets to measure the accuracy. A PWM+histone naïve Bayes predictor using a single histone modification (H3K4me3) is substantially more accurate than a PWM score or a conservation-based score (phylogenetic motif model). The naïve Bayes predictor is more accurate (on average) at all sensitivity levels, and makes only half as many false positive predictions at sensitivity levels from 10% to 80%. On average, it correctly predicts 80% of bound promoters at a false positive rate of 20%. Accuracy does not diminish when we test the predictor in a different cell type (and species) from training. Accuracy is barely diminished even when we train the predictor without using TF ChIP-seq data.

AVAILABILITY

Our tissue-specific predictor of promoters bound by a TF is called Dr Gene and is available at http://bioinformatics.org.au/drgene.

CONTACT

t.bailey@imb.uq.edu.au

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcription intermediary factor 1γ is a tumor suppressor in mouse and human chronic myelomonocytic leukemia

Transcription intermediary factor 1γ (TIF1γ) was suggested to play a role in erythropoiesis. However, how TIF1γ regulates the development of different blood cell lineages and whether TIF1γ is involved in human hematological malignancies remain to be determined. Here we have shown that TIF1γ was a tumor suppressor in mouse and human chronic myelomonocytic leukemia (CMML). Loss of Tif1g in mouse HSCs favored the expansion of the granulo-monocytic progenitor compartment. Furthermore, Tif1g deletion induced the age-dependent appearance of a cell-autonomous myeloproliferative disorder in mice that recapitulated essential characteristics of human CMML. TIF1γ was almost undetectable in leukemic cells of 35% of CMML patients. This downregulation was related to the hypermethylation of CpG sequences and specific histone modifications in the gene promoter. A demethylating agent restored the normal epigenetic status of the TIF1G promoter in human cells, which correlated with a reestablishment of TIF1γ expression. Together, these results demonstrate that TIF1G is an epigenetically regulated tumor suppressor gene in hematopoietic cells and suggest that changes in TIF1γ expression may be a biomarker of response to demethylating agents in CMML.

Revealing epigenetic patterns in gene regulation through integrative analysis of epigenetic interaction network

Chromatin components and DNA methylation play important roles in regulation of gene expression in mammalian genomes. However, the mechanism underlying how they regulate gene transcription, independently or synergistically, remains largely unknown. We constructed an epigenetic interaction network (EIN) of chromatin components, DNA methylation and gene expression by combining partial correlation coefficient with Pearson correlation coefficient. In EIN, we identified nine direct factors for gene expression. They constitute three interaction modules which synergistically affect gene expression. We introduced a new combination strategy to test how these direct factors in each module regulate gene expression synergistically. We found two inter-attracted patterns and one inter-repulsed patterns among the three modules. Furthermore, we identified 22 indirect factors for gene expression which have effect on gene expression via direct factors. DNA methylation, for example, could regulate gene expression through H3K4me3 and Pol II. Our approach has the potential to help in uncovering inherent relationships between epigenetic factors and gene transcription and guiding experiment.

Histone modification profiles are predictive for tissue/cell-type specific expression of both protein-coding and microRNA genes

Background

Gene expression is regulated at both the DNA sequence level and through modification of chromatin. However, the effect of chromatin on tissue/cell-type specific gene regulation (TCSR) is largely unknown. In this paper, we present a method to elucidate the relationship between histone modification/variation (HMV) and TCSR.

Results

A classifier for differentiating CD4+ T cell-specific genes from housekeeping genes using HMV data was built. We found HMV in both promoter and gene body regions to be predictive of genes which are targets of TCSR. For example, the histone modification types H3K4me3 and H3K27ac were identified as the most predictive for CpG-related promoters, whereas H3K4me3 and H3K79me3 were the most predictive for nonCpG-related promoters. However, genes targeted by TCSR can be predicted using other type of HMVs as well. Such redundancy implies that multiple type of underlying regulatory elements, such as enhancers or intragenic alternative promoters, which can regulate gene expression in a tissue/cell-type specific fashion, may be marked by the HMVs. Finally, we show that the predictive power of HMV for TCSR is not limited to protein-coding genes in CD4+ T cells, as we successfully predicted TCSR targeted genes in muscle cells, as well as microRNA genes with expression specific to CD4+ T cells, by the same classifier which was trained on HMV data of protein-coding genes in CD4+ T cells.

Conclusion

We have begun to understand the HMV patterns that guide gene expression in both tissue/cell-type specific and ubiquitous manner.

Epigenomic and RNA structural correlates of polyadenylation

Polyadenylation (poly(A)) of mRNA plays a critical role in regulating gene expression. Identifying the sequence, structural, and epigenomic determinants of poly(A) site usage is an important long term goal. Several cis elements that mediate poly(A) regulation have been identified. Highly used poly(A) sites are also known to have a greater nucleosome occupancy in the immediate downstream. However, a detailed exploration of additional epigenomic and mRNA structural correlates of poly(A) site usage has not been reported. Importantly, functional interaction between sequence, structure, and the epigenome in determining the poly(A) site usage is not known. We show that highly used poly(A) sites are positively associated with an mRNA structure that is energetically more favorable and one that better exposes a critical polyadenylation cis element. In exploring potential interplay between RNA and chromatin structure, we found that a stronger nucleosome occupancy downstream of poly(A) site strongly correlated with (1) a more favorable mRNA structure, and (2) a greater accumulation of RNA Polymerase II (PolII) at the poly(A) site. Further analysis suggested a causal relationship pointing from PolII accumulation to a stable RNA structure. Additionally, we found that distinct patterns of histone modifications characterize poly(A) sites and these epigenetic patterns alone can distinguish true poly(A) sites with ~76% accuracy and also discriminate between high and low usage poly(A) sites with ~74% accuracy. Our results suggest a causative link between chromatin structure and mRNA structure whereby a compacted chromatin downstream of the poly(A) site slows down the elongating transcript, thus facilitating the folding of nascent mRNA in a favorable structure at poly(A) site during transcription. Additionally we report hitherto unknown epigenomic correlates for poly(A) site usage.

The histone modifications governing TFF1 transcription mediated by estrogen receptor

Transcription regulation by histone modifications is a major contributing factor to the structural and functional diversity in biology. These modifications are encrypted as histone codes or histone languages and function to establish and maintain heritable epigenetic codes that define the identity and the fate of the cell. Despite recent advances revealing numerous histone modifications associated with transcription regulation, how such modifications dictate the process of transcription is not fully understood. Here we describe spatial and temporal analyses of the histone modifications that are introduced during estrogen receptor α (ERα)-activated transcription. We demonstrated that aborting RNA polymerase II caused a disruption of the histone modifications that are associated with transcription elongation but had a minimal effect on modifications deposited during transcription initiation. We also found that the histone H3S10 phosphorylation mark is catalyzed by mitogen- and stress-activated protein kinase 1 (MSK1) and is recognized by a 14-3-3ζ/14-3-3ε heterodimer through its interaction with H3K4 trimethyltransferase SMYD3 and the p52 subunit of TFIIH. We showed that H3S10 phosphorylation is a prerequisite for H3K4 trimethylation. In addition, we demonstrated that SET8/PR-Set7/KMT5A is required for ERα-regulated transcription and its catalyzed H4K20 monomethylation is implicated in both transcription initiation and elongation. Our experiments provide a relatively comprehensive analysis of histone modifications associated with ERα-regulated transcription and define the biological meaning of several key components of the histone code that governs ERα-regulated transcription.

Epigenetic regulation of key vascular genes and growth factors

The role of small RNAs in epigenetic regulation is an emerging field. This research may also open novel treatment strategies based on manipulation of the epigenetic status of the target tissues. Our objective is to review epigenetic regulation of key vascular genes and growth factors. Vascular endothelial growth factor A (VEGF-A) is one of the key players in regulating and maintaining cardiovascular functions and pathology. Although its epigenetic regulation is still not completely understood, expression of the VEGF gene can be manipulated by epigenetic mechanisms using small RNAs that are targeted to the gene promoter which results in the alteration of histone code. VEGF exerts its effects mostly through two receptors, VEGFR1 and VEGFR2, and their expression is also regulated by promoter DNA methylation in various cancer cells. These findings suggest the importance of epigenetic mechanisms in the regulation of vascular functions.

Computational study of associations between histone modification and protein-DNA binding in yeast genome by integrating diverse information

Background

In parallel with the quick development of high-throughput technologies, in vivo (vitro) experiments for genome-wide identification of protein-DNA interactions have been developed. Nevertheless, a few questions remain in the field, such as how to distinguish true protein-DNA binding (functional binding) from non-specific protein-DNA binding (non-functional binding). Previous researches tackled the problem by integrated analysis of multiple available sources. However, few systematic studies have been carried out to examine the possible relationships between histone modification and protein-DNA binding. Here this issue was investigated by using publicly available histone modification data in yeast.

Results

Two separate histone modification datasets were studied, at both the open reading frame (ORF) and the promoter region of binding targets for 37 yeast transcription factors. Both results revealed a distinct histone modification pattern between the functional protein-DNA binding sites and non-functional ones for almost half of all TFs tested. Such difference is much stronger at the ORF than at the promoter region. In addition, a protein-histone modification interaction pathway can only be inferred from the functional protein binding targets.

Conclusions

Overall, the results suggest that histone modification information can be used to distinguish the functional protein-DNA binding from the non-functional, and that the regulation of various proteins is controlled by the modification of different histone lysines such as the protein-specific histone modification levels.

Ablation of CBP in forebrain principal neurons causes modest memory and transcriptional defects and a dramatic reduction of histone acetylation but does not affect cell viability

Rubinstein-Taybi syndrome (RSTS) is an inheritable disease associated with mutations in the gene encoding the CREB (cAMP response element-binding protein)-binding protein (CBP) and characterized by growth impairment, learning disabilities, and distinctive facial and skeletal features. Studies in mouse models for RSTS first suggested a direct role for CBP and histone acetylation in cognition and memory. Here, we took advantage of the genetic tools for generating mice in which the CBP gene is specifically deleted in postmitotic principal neurons of the forebrain to investigate the consequences of the loss of CBP in the adult brain. In contrast to the conventional CBP knock-out mice, which exhibit very early embryonic lethality, postnatal forebrain-restricted CBP mutants were viable and displayed no overt abnormalities. We identified the dimer of histones H2A and H2B as the preferred substrate of the histone acetyltransferase domain of CBP. Surprisingly, the loss of CBP and subsequent histone hypoacetylation had a very modest impact in the expression of a number of immediate early genes and did not affect neuronal viability. In addition, the behavioral characterization of these mice dissociated embryonic and postnatal deficits caused by impaired CBP function, narrowed down the anatomical substrate of specific behavioral defects, and confirmed the special sensitivity of object recognition memory to CBP deficiency. Overall, our study provides novel insights into RSTS etiology and clarifies some of the standing questions concerning the role of CBP and histone acetylation in activity-driven gene expression, memory formation, and neurodegeneration.

Balancing histone methylation activities in psychiatric disorders

Alterations in histone lysine methylation and other epigenetic regulators of gene expression contribute to changes in brain transcriptomes in mood and psychosis spectrum disorders, including depression and schizophrenia. Genetic association studies and animal models implicate multiple lysine methyltransferases and demethylases in the neurobiology of emotion and cognition. Here, we review the role of histone lysine methylation and transcriptional regulation in normal and diseased neurodevelopment and discuss various methyltransferases and demethylases as potential therapeutic targets in the treatment of neuropsychiatric disease.

Requirement for Dot1l in murine postnatal hematopoiesis and leukemogenesis by MLL translocation

Disruptor of telomeric silencing 1-like (Dot1l) is a histone 3 lysine 79 methyltransferase. Studies of constitutive Dot1l knockout mice show that Dot1l is essential for embryonic development and prenatal hematopoiesis. DOT1L also interacts with translocation partners of Mixed Lineage Leukemia (MLL) gene, which is commonly translocated in human leukemia. However, the requirement of Dot1l in postnatal hematopoiesis and leukemogenesis of MLL translocation proteins has not been conclusively shown. With a conditional Dot1l knockout mouse model, we examined the consequences of Dot1l loss in postnatal hematopoiesis and MLL translocation leukemia. Deletion of Dot1l led to pancytopenia and failure of hematopoietic homeostasis, and Dot1l-deficient cells minimally reconstituted recipient bone marrow in competitive transplantation experiments. In addition, MLL-AF9 cells required Dot1l for oncogenic transformation, whereas cells with other leukemic oncogenes, such as Hoxa9/Meis1 and E2A-HLF, did not. These findings illustrate a crucial role of Dot1l in normal hematopoiesis and leukemogenesis of specific oncogenes.

A statistical framework for modeling gene expression using chromatin features and application to modENCODE datasets

We develop a statistical framework to study the relationship between chromatin features and gene expression. This can be used to predict gene expression of protein coding genes, as well as microRNAs. We demonstrate the prediction in a variety of contexts, focusing particularly on the modENCODE worm datasets. Moreover, our framework reveals the positional contribution around genes (upstream or downstream) of distinct chromatin features to the overall prediction of expression levels.

Predicting gene expression in T cell differentiation from histone modifications and transcription factor binding affinities by linear mixture models

BACKGROUND

The differentiation process from stem cells to fully differentiated cell types is controlled by the interplay of chromatin modifications and transcription factor activity. Histone modifications or transcription factors frequently act in a multi-functional manner, with a given DNA motif or histone modification conveying both transcriptional repression and activation depending on its location in the promoter and other regulatory signals surrounding it.

RESULTS

To account for the possible multi functionality of regulatory signals, we model the observed gene expression patterns by a mixture of linear regression models. We apply the approach to identify the underlying histone modifications and transcription factors guiding gene expression of differentiated CD4+ T cells. The method improves the gene expression prediction in relation to the use of a single linear model, as often used by previous approaches. Moreover, it recovered the known role of the modifications H3K4me3 and H3K27me3 in activating cell specific genes and of some transcription factors related to CD4+ T differentiation.

Gene expression differences among primates are associated with changes in a histone epigenetic modification

Changes in gene regulation are thought to play an important role in speciation and adaptation, especially in primates. However, we still know relatively little about the mechanisms underlying regulatory evolution. In particular, the extent to which epigenetic modifications underlie gene expression differences between primates is not yet known. Our study focuses on an epigenetic histone modification, H3K4me3, which is thought to promote transcription. To investigate the contribution of H3K4me3 to regulatory differences between species, we collected gene expression data and identified H3K4me3-associated genomic regions in lymphoblastoid cell lines (LCLs) from humans, chimpanzees, and rhesus macaques, using three cell lines from each species. We found strong evidence for conservation of H3K4me3 localization in primates. Moreover, regardless of species, H3K4me3 is consistently enriched near annotated transcription start sites (TSS), and highly expressed genes are more likely than lowly expressed genes to have the histone modification near their TSS. Interestingly, we observed an enrichment of interspecies differences in H3K4me3 at the TSS of genes that are differentially expressed between species. We estimate that as much as 7% of gene expression differences between the LCLs of humans, chimpanzees, and rhesus macaques may be explained, at least in part, by changes in the status of H3K4me3 histone modifications. Our results suggest a modest, yet important role for epigenetic changes in gene expression differences between primates.

Selective targeting of histone methylation

Histones are post-translationally modified by multiple histone-modifying enzymes, which in turn influences gene expression. Much of the work in the field to date has focused on genetic, biochemical and structural characterization of these enzymes. The most recent genome-wide methods provide insights into specific recruitment of histone-modifying enzymes in vivo and, therefore, onto mechanisms of establishing a differential expression pattern. Here we focus on the recruitment mechanisms of the enzymes involved in the placement of two contrasting histone marks, histone H3 lysine 4 (H3K4) methylation and histone H3 lysine 27 (H3K27) methylation. We describe distribution of their binding sites and show that recruitment of different histone-modifying proteins can be coordinated, opposed, or alternating. Specifically, genomic sites of the H3K4 histone demethylase KDM5A become accessible to its homolog KDM5B in cells with a lowered KDM5A level. The currently available data on recruitment of H3K4/H3K27 modifying enzymes suggests that the formed protein complexes are targeted in a sequential and temporal manner, but that additional, still unknown, interactions contribute to targeting specificity.

Epigenetic modifications and human disease

Epigenetics is one of the most rapidly expanding fields in biology. The recent characterization of a human DNA methylome at single nucleotide resolution, the discovery of the CpG island shores, the finding of new histone variants and modifications, and the unveiling of genome-wide nucleosome positioning maps highlight the accelerating speed of discovery over the past two years. Increasing interest in epigenetics has been accompanied by technological breakthroughs that now make it possible to undertake large-scale epigenomic studies. These allow the mapping of epigenetic marks, such as DNA methylation, histone modifications and nucleosome positioning, which are critical for regulating gene and noncoding RNA expression. In turn, we are learning how aberrant placement of these epigenetic marks and mutations in the epigenetic machinery is involved in disease. Thus, a comprehensive understanding of epigenetic mechanisms, their interactions and alterations in health and disease, has become a priority in biomedical research.

Combinatorial chromatin modification patterns in the human genome revealed by subspace clustering

Chromatin modifications, such as post-translational modification of histone proteins and incorporation of histone variants, play an important role in regulating gene expression. Joint analyses of multiple histone modification maps are starting to reveal combinatorial patterns of modifications that are associated with functional DNA elements, providing support to the ‘histone code’ hypothesis. However, due to the lack of analytical methods, only a small number of chromatin modification patterns have been discovered so far. Here, we introduce a scalable subspace clustering algorithm, coherent and shifted bicluster identification (CoSBI), to exhaustively identify the set of combinatorial modification patterns across a given epigenome. Performance comparisons demonstrate that CoSBI can generate biclusters with higher intra-cluster coherency and biological relevance. We apply our algorithm to a compendium of 39 genome-wide chromatin modification maps in human CD4⁺ T cells. We identify 843 combinatorial patterns that recur at >0.1% of the genome. A total of 19 chromatin modifications are observed in the combinatorial patterns, 10 of which occur in more than half of the patterns. We also identify combinatorial modification signatures for eight classes of functional DNA elements. Application of CoSBI to epigenome maps of different cells and developmental stages will aid in understanding how chromatin structure helps regulate gene expression.

Decoding the histone H4 lysine 20 methylation mark

The molecular biology of histone H4 lysine 20 (H4K20) methylation, like many other post-translational modifications of histones, has been the subject of intensive interest in recent years. While there is an emerging consensus linking H4K20me1, H4K20me2, and H4K20me3 to transcription, repair, and constitutive heterochromatin, respectively, the specific details of these associations and the biological mechanisms by which the methylated histones are introduced and function are now the subject of active investigation. Although a large number of methylases capable of methylating H4K20 have been identified and characterized; there is no known demethylase of H4K20, though the search is ongoing. Additionally, many recent studies have been directed at understanding the role of methylated H4K20 and other histone modifications associated with different biological processes in the context of a combinatorial histone code. It seems likely that continued study of the methylation of H4K20 will yield extremely valuable insights concerning the regulation of histone modifications before and during cell division and the impact of these modifications on subsequent gene expression.

The complexity of gene expression dynamics revealed by permutation entropy

BACKGROUND

High complexity is considered a hallmark of living systems. Here we investigate the complexity of temporal gene expression patterns using the concept of Permutation Entropy (PE) first introduced in dynamical systems theory. The analysis of gene expression data has so far focused primarily on the identification of differentially expressed genes, or on the elucidation of pathway and regulatory relationships. We aim to study gene expression time series data from the viewpoint of complexity.

RESULTS

Applying the PE complexity metric to abiotic stress response time series data in Arabidopsis thaliana, genes involved in stress response and signaling were found to be associated with the highest complexity not only under stress, but surprisingly, also under reference, non-stress conditions. Genes with house-keeping functions exhibited lower PE complexity. Compared to reference conditions, the PE of temporal gene expression patterns generally increased upon stress exposure. High-complexity genes were found to have longer upstream intergenic regions and more cis-regulatory motifs in their promoter regions indicative of a more complex regulatory apparatus needed to orchestrate their expression, and to be associated with higher correlation network connectivity degree. Arabidopsis genes also present in other plant species were observed to exhibit decreased PE complexity compared to Arabidopsis specific genes.

CONCLUSIONS

We show that Permutation Entropy is a simple yet robust and powerful approach to identify temporal gene expression profiles of varying complexity that is equally applicable to other types of molecular profile data.

Chromatin condensation via the condensin II complex is required for peripheral T-cell quiescence

Naive T cells encountering their cognate antigen become activated and acquire the ability to proliferate in response to cytokines. Stat5 is an essential component in this response. We demonstrate that Stat5 cannot access DNA in naive T cells and acquires this ability only after T-cell receptor (TCR) engagement. The transition is not associated with changes in DNA methylation or global histone modification but rather chromatin decondensation. Condensation occurs during thymocyte development and proper condensation is dependent on kleisin-β of the condensin II complex. Our findings suggest that this unique chromatin condensation, which can affect interpretations of chromatin accessibility assays, is required for proper T-cell development and maintenance of the quiescent state. This mechanism ensures that cytokine driven proliferation can only occur in the context of TCR stimulation.

Mapping dynamic histone acetylation patterns to gene expression in nanog-depleted murine embryonic stem cells

Embryonic stem cells (ESC) have the potential to self-renew indefinitely and to differentiate into any of the three germ layers. The molecular mechanisms for self-renewal, maintenance of pluripotency and lineage specification are poorly understood, but recent results point to a key role for epigenetic mechanisms. In this study, we focus on quantifying the impact of histone 3 acetylation (H3K9,14ac) on gene expression in murine embryonic stem cells. We analyze genome-wide histone acetylation patterns and gene expression profiles measured over the first five days of cell differentiation triggered by silencing Nanog, a key transcription factor in ESC regulation. We explore the temporal and spatial dynamics of histone acetylation data and its correlation with gene expression using supervised and unsupervised statistical models. On a genome-wide scale, changes in acetylation are significantly correlated to changes in mRNA expression and, surprisingly, this coherence increases over time. We quantify the predictive power of histone acetylation for gene expression changes in a balanced cross-validation procedure. In an in-depth study we focus on genes central to the regulatory network of Mouse ESC, including those identified in a recent genome-wide RNAi screen and in the PluriNet, a computationally derived stem cell signature. We find that compared to the rest of the genome, ESC-specific genes show significantly more acetylation signal and a much stronger decrease in acetylation over time, which is often not reflected in a concordant expression change. These results shed light on the complexity of the relationship between histone acetylation and gene expression and are a step forward to dissect the multilayer regulatory mechanisms that determine stem cell fate.

Stem cell differentiation and the maintenance of self-renewal are intrinsically complex processes that require coordinated regulation on many different cellular levels. Here we focus on the relationship between two important layers and follow it over the first five days of differentiation. The first layer – measured by acetylation of one of the histone proteins – describes which parts of the DNA are tightly wrapped up and which lie open. The second layer describes the activity of genes measured by their mRNA expression. Using a wide array of statistical approaches we show that changes in histone acetylation are very predictive for gene expression and that the concordance between the two levels increases over time. Concentrating on genes central to the regulatory networks in embryonic stem cells we find that key genes show very high acetylation signal in the beginning that decreases quickly over time, indicating that they lie in initially open regions that are rapidly closing down. These results are a step forward to a better understanding of the complexities of the relationship between histone acetylation and gene expression, which will help to dissect the multilayer regulatory mechanisms that determine stem cell fate.

Chromatin remodeling in the aging genome of Drosophila

Chromatin structure affects the accessibility of DNA to transcription, repair, and replication. Changes in chromatin structure occur during development, but less is known about changes during aging. We examined the state of chromatin structure and its effect on gene expression during aging in Drosophila at the whole genome and cellular level using whole-genome tiling microarrays of activation and repressive chromatin marks, whole-genome transcriptional microarrays and single-cell immunohistochemistry. We found dramatic reorganization of chromosomal regions with age. Mapping of H3K9me3 and HP1 signals to fly chromosomes reveals in young flies the expected high enrichment in the pericentric regions, the 4th chromosome, and islands of facultative heterochromatin dispersed throughout the genome. With age, there is a striking reduction in this enrichment resulting in a nearly equivalent level of H3K9me3 and HP1 in the pericentric regions, the 4th chromosome, facultative heterochromatin, and euchromatin. These extensive changes in repressive chromatin marks are associated with alterations in age-related gene expression. Large-scale changes in repressive marks with age are further substantiated by single-cell immunohistochemistry that shows changes in nuclear distribution of H3K9me3 and HP1 marks with age. Such epigenetic changes are expected to directly or indirectly impinge upon important cellular functions such as gene expression, DNA repair, and DNA replication. The combination of genome-wide approaches such as whole-genome chromatin immunoprecipitation and transcriptional studies in conjunction with single-cell immunohistochemistry as shown here provide a first step toward defining how changes in chromatin may contribute to the process of aging in metazoans.

Every methyl counts--epigenetic calculus

Histone modifications play an important role in the formation of an epigenetic memory system that maintains cellular identity. Their complex patterns have been suggested to constitute a histone code, which encodes for specific forms of chromatin. According to the histone code hypothesis these specific patterns are passed on from one cell generation to the next. This enables cells to keep a specific gene expression pattern even in absence of the specific transcription factors that initiated the expression of lineage determining genes. The methylation of specific lysine residues within the histone tails plays a particularly important role in defining the histone modification pattern as mutations of the enzymes that catalyze the formation or the removal of methyl groups have severe effects on cellular physiology. Lysines can get mono-, di- or trimethylated, but the molecular function of the different modification states is still not fully understood. In the following review we will highlight recent data that try to tackle this question and discuss their potential impact for our understanding of the role of histone methylation in epigenetic inheritance.

Hippocampal gene profiling: toward a systems biology of the hippocampus

Transcriptomics and proteomics approaches give a unique perspective for understanding brain and hippocampal functions but also pose unique challenges because of the singular complexity of the nervous system. The proliferation of genome-wide expression studies during the last decade has provided important insight into the molecular underpinnings of brain anatomy, neural plasticity, and neurological diseases. Microarray technology has dominated transcriptomics research, but this situation is rapidly changing with the recent technological advances in high-throughput sequencing. The full potential of transcriptomics in the neurosciences will be achieved as a result of its integration with other "-omics" disciplines as well as the development of novel analytical bioinformatics and systems biology tools for meta-analysis. Here, we review some of the most relevant advances in the gene profiling of the hippocampus, its relationship with proteomics approaches, and the promising perspectives for the future.

Ab initio identification of transcription start sites in the Rhesus macaque genome by histone modification and RNA-Seq

Rhesus macaque is a widely used primate model organism. Its genome annotations are however still largely comparative computational predictions derived mainly from human genes, which precludes studies on the macaque-specific genes, gene isoforms or their regulations. Here we took advantage of histone H3 lysine 4 trimethylation (H3K4me3)’s ability to mark transcription start sites (TSSs) and the recently developed ChIP-Seq and RNA-Seq technology to survey the transcript structures. We generated 14 013 757 sequence tags by H3K4me3 ChIP-Seq and obtained 17 322 358 paired end reads for mRNA, and 10 698 419 short reads for sRNA from the macaque brain. By integrating these data with genomic sequence features and extending and improving a state-of-the-art TSS prediction algorithm, we ab initio predicted and verified 17 933 of previously electronically annotated TSSs at 500-bp resolution. We also predicted approximately 10 000 novel TSSs. These provide an important rich resource for close examination of the species-specific transcript structures and transcription regulations in the Rhesus macaque genome. Our approach exemplifies a relatively inexpensive way to generate a reasonably reliable TSS map for a large genome. It may serve as a guiding example for similar genome annotation efforts targeted at other model organisms.

Variable histone modifications at the A(vy) metastable epiallele

The ability of environmental factors to shape health and disease involves epigenetic mechanisms that mediate gene-environment interactions. Metastable epiallele genes are variably expressed in genetically identical individuals due to epigenetic modifications established during early development. DNA methylation within metastable epialleles is stochastic due to probabilistic reprogramming of epigenetic marks during embryogenesis. Maternal nutrition and environment have been shown to affect metastable epiallele methylation patterns and subsequent adult phenotype. Little is known, however, about the role of histone modifications in influencing metastable epiallele expression and phenotypic variation. Utilizing chromatin immunoprecipitation followed by qPCR, we observe variable histone patterns in the 5' long terminal repeat (LTR) of the murine viable yellow agouti (A(vy)) metastable epiallele. This region contains 6 CpG sites, which are variably methylated in isogenic A(vy)/a offspring. Yellow mice, which are hypomethylated at the Avy LTR and exhibit constitutive ectopic expression of agouti (a), also display enrichment of H3 and H4 di-acetylation (p = 0.08 and 0.09, respectively). Pseudoagouti mice, in which A(vy) hypermethylation is thought to silence ectopic expression, exhibit enrichment of H4K20 tri-methylation (p = 0.01). No differences are observed for H3K4 tri-methylation (p = 0.7), a modification often enriched in the promoter of active genes. These results show for the first time the presence of variable histone modifications at a metastable epiallele, indicating that DNA methylation acts in concert with histone modifications to affect inter-individual variation of metastable epiallele expression. Therefore, the potential for environmental factors to influence histone modifications, in addition to DNA methylation, should be addressed in environmental epigenomic studies.

A unique H3K4me2 profile marks tissue-specific gene regulation

Characterization of the epigenetic landscape fundamentally contributes toward deciphering the regulatory mechanisms that govern gene expression. However, despite an increasing flow of newly generated data, no clear pattern of chromatin modifications has so far been linked to specific modes of transcriptional regulation. Here, we used high-throughput genomic data from CD4(+) T lymphocytes to provide a comprehensive analysis of histone H3 lysine 4 dimethylation (H3K4me2) enrichment in genomic regions surrounding transcriptional start sites (TSSs). We discovered that a subgroup of genes linked to T cell functions displayed high levels of H3K4me2 within their gene body, in sharp contrast to the TSS-centered profile typical of housekeeping genes. Analysis of additional chromatin modifications and DNase I hypersensitive sites (DHSS) revealed a combinatorial chromatin signature characteristic of this subgroup. We propose that this epigenetic feature reflects the activity of an as yet unrecognized, intragenic cis-regulatory platform dedicated to refining tissue-specificity in gene expression.

Development of homogeneous luminescence assays for histone demethylase catalysis and binding

Covalent modifications to histones play important roles in chromatin dynamics and the regulation of gene expression. The JumonjiC (JmjC)-containing histone demethylases (HDMs) catalyze the demethylation of methylated lysine residues on histone tails. Here we report the development of homogeneous luminescence-based assay methods for measuring the catalytic activity and the binding affinities of peptides to HDMs. The assays use amplified luminescent proximity homogeneous assay (ALPHA) technology, are sensitive and robust, and can be used for small molecule inhibitor screening of HDMs. We have profiled known inhibitors of JMJD2E and demonstrate a correlation between the inhibitor potencies determined by the ALPHA and other types of assays. Although this study focuses on the JMJD2E isoform, the catalytic turnover and binding assays described here can be used in studies on other HDMs. The assays should be useful for the development of small molecule inhibitors selective for HDM isoforms.

Kinetics of re-establishing H3K79 methylation marks in global human chromatin

We employ a stable isotope strategy wherein both histones and their methylations are labeled in synchronized human cells. This allows us to differentiate between old and new methylations on pre-existing versus newly synthesized histones. The strategy is implemented on K79 methylation in an isoform-specific manner for histones H3.1, H3.2, and H3.3. Although levels of H3.3K79 monomethylation are higher than that of H3.2/H3.1, the rate of establishing the K79 methylation is the same for all three isoforms. Surprisingly, we find that pre-existing "old" histones continue to be K79-monomethylated and -dimethylated at a rate equal to the newly synthesized histones. These observations imply that some degree of positional "scrambling" of K79 methylation occurs through the cell cycle.

Application of machine learning methods to histone methylation ChIP-Seq data reveals H4R3me2 globally represses gene expression

Background

In the last decade, biochemical studies have revealed that epigenetic modifications including histone modifications, histone variants and DNA methylation form a complex network that regulate the state of chromatin and processes that depend on it including transcription and DNA replication. Currently, a large number of these epigenetic modifications are being mapped in a variety of cell lines at different stages of development using high throughput sequencing by members of the ENCODE consortium, the NIH Roadmap Epigenomics Program and the Human Epigenome Project. An extremely promising and underexplored area of research is the application of machine learning methods, which are designed to construct predictive network models, to these large-scale epigenomic data sets.

Results

Using a ChIP-Seq data set of 20 histone lysine and arginine methylations and histone variant H2A.Z in human CD4⁺ T-cells, we built predictive models of gene expression as a function of histone modification/variant levels using Multilinear (ML) Regression and Multivariate Adaptive Regression Splines (MARS). Along with extensive crosstalk among the 20 histone methylations, we found H4R3me2 was the most and second most globally repressive histone methylation among the 20 studied in the ML and MARS models, respectively. In support of our finding, a number of experimental studies show that PRMT5-catalyzed symmetric dimethylation of H4R3 is associated with repression of gene expression. This includes a recent study, which demonstrated that H4R3me2 is required for DNMT3A-mediated DNA methylation--a known global repressor of gene expression.

Conclusion

In stark contrast to univariate analysis of the relationship between H4R3me2 and gene expression levels, our study showed that the regulatory role of some modifications like H4R3me2 is masked by confounding variables, but can be elucidated by multivariate/systems-level approaches.

Introduction into the analysis of high-throughput-sequencing based epigenome data

Sequencing-based approaches now allow high-resolution, genome-scale investigation of cellular epigenetic landscapes. For example, mapping of open chromatin regions, post-translational histone modifications and DNA methylation across a whole genome is now feasible, and new non-coding regulatory RNAs can be sensitively identified via RNA sequencing. The resulting large-scale data sets promise to contribute towards a more precise and complete understanding of gene regulation and to yield insights into the interplay between genomes and the environment. In this article, I review some of the conceptual issues and currently available software tools for the analysis of sequencing-based whole-genome epigenetics data.

Coordinated chromatin control: structural and functional linkage of DNA and histone methylation

One of the most fundamental questions in the control of gene expression in mammals is how epigenetic methylation patterns of DNA and histones are established, erased, and recognized. This central process in controlling metazoan gene expression includes coordinated covalent modifications of DNA and its associated histones. This review focuses on recent developments in characterizing the functional links between the methylation status of the DNA and of two particularly important histone marks. Mammalian DNA methylation is intricately connected to the presence of unmodified lysine 4 and methylated lysine 9 residues in histone H3. An interconnected network of methyltransferases, demethylases, and accessory proteins is responsible for changing or maintaining the modification status of specific regions of chromatin. The structural and functional interactions among members of this network are critical to processes that include imprinting and differentiation, dysregulation of which is associated with disorders ranging from inflammation to cancer.

Regulation by polycomb and trithorax group proteins in Arabidopsis

Polycomb group (PcG) and trithorax group (trxG) proteins are key regulators of homeotic genes and have crucial roles in cell proliferation, growth and development. PcG and trxG proteins form higher order protein complexes that contain SET domain proteins, with a histone methyltransferase (HMTase) activity, responsible for the different types of lysine methylation at the N-terminal tails of the core histone proteins. In recent years, genetic studies along with biochemical and cell biological analyses in Arabidopsis have enabled researchers to begin to understand how PcG and trxG proteins are recruited to chromatin and how they regulate their target genes and to elucidate their functions. This review focuses on the advances in our understanding of the biological roles of PcG and trxG proteins, their molecular mechanisms of action and further examines the role of histone marks in PcG and trxG regulation in Arabidopsis.