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 an 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.Comment: accepted at PLoS Computational Biolog

Patterns and Mechanisms of Ancestral Histone Protein Inheritance in Budding Yeast

Tracking of ancestral histone proteins over multiple generations of genome replication in yeast reveals that old histones move along genes from 3′ toward 5′ over time, and that maternal histones move up to around 400 bp during genomic replication

The Inheritance of Histone Modifications Depends upon the Location in the Chromosome in Saccharomyces cerevisiae

Histone modifications are important epigenetic features of chromatin that must be replicated faithfully. However, the molecular mechanisms required to duplicate and maintain histone modification patterns in chromatin remain to be determined. Here, we show that the introduction of histone modifications into newly deposited nucleosomes depends upon their location in the chromosome. In Saccharomyces cerevisiae, newly deposited nucleosomes consisting of newly synthesized histone H3-H4 tetramers are distributed throughout the entire chromosome. Methylation of lysine 4 on histone H3 (H3-K4), a hallmark of euchromatin, is introduced into these newly deposited nucleosomes, regardless of whether the neighboring preexisting nucleosomes harbor the K4 mutation in histone H3. Furthermore, if the heterochromatin-binding protein Sir3 is unavailable during DNA replication, histone H3-K4 methylation is introduced onto newly deposited nucleosomes in telomeric heterochromatin. Thus, a conservative distribution model most accurately explains the inheritance of histone modifications because the location of histones within euchromatin or heterochromatin determines which histone modifications are introduced

The Set2/Rpd3S Pathway Suppresses Cryptic Transcription without Regard to Gene Length or Transcription Frequency

In cells lacking the histone methyltransferase Set2, initiation of RNA polymerase II transcription occurs inappropriately within the protein-coding regions of genes, rather than being restricted to the proximal promoter. It was previously reported that this “cryptic” transcription occurs preferentially in long genes, and in genes that are infrequently transcribed. Here, we mapped the transcripts produced in an S. cerevisiae strain lacking Set2, and applied rigorous statistical methods to identify sites of cryptic transcription at high resolution. We find that suppression of cryptic transcription occurs independent of gene length or transcriptional frequency. Our conclusions differ with those reported previously because we obtained a higher-resolution dataset, we accounted for the fact that gene length and transcriptional frequency are not independent variables, and we accounted for several ascertainment biases that make cryptic transcription easier to detect in long, infrequently transcribed genes. These new results and conclusions have implications for many commonly used genomic analysis approaches, and for the evolution of high-fidelity RNA polymerase II transcriptional initiation in eukaryotes

A complete set of nascent transcription rates for yeast genes

The amount of mRNA in a cell is the result of two opposite reactions: transcription and mRNA degradation. These reactions are governed by kinetics laws, and the most regulated step for many genes is the transcription rate. The transcription rate, which is assumed to be exercised mainly at the RNA polymerase recruitment level, can be calculated using the RNA polymerase densities determined either by run-on or immunoprecipitation using specific antibodies. The yeast Saccharomyces cerevisiae is the ideal model organism to generate a complete set of nascent transcription rates that will prove useful for many gene regulation studies. By combining genomic data from both the GRO (Genomic Run-on) and the RNA pol ChIP-on-chip methods we generated a new, more accurate nascent transcription rate dataset. By comparing this dataset with the indirect ones obtained from the mRNA stabilities and mRNA amount datasets, we are able to obtain biological information about posttranscriptional regulation processes and a genomic snapshot of the location of the active transcriptional machinery. We have obtained nascent transcription rates for 4,670 yeast genes. The median RNA polymerase II density in the genes is 0.078 molecules/kb, which corresponds to an average of 0.096 molecules/gene. Most genes have transcription rates of between 2 and 30 mRNAs/hour and less than 1% of yeast genes have >1 RNA polymerase molecule/gene. Histone and ribosomal protein genes are the highest transcribed groups of genes and other than these exceptions the transcription of genes is an infrequent phenomenon in a yeast cell

Replication and Active Demethylation Represent Partially Overlapping Mechanisms for Erasure of H3K4me3 in Budding Yeast

Histone modifications affect DNA–templated processes ranging from transcription to genomic replication. In this study, we examine the cell cycle dynamics of the trimethylated form of histone H3 lysine 4 (H3K4me3), a mark of active chromatin that is viewed as “long-lived” and that is involved in memory during cell state inheritance in metazoans. We synchronized yeast using two different protocols, then followed H3K4me3 patterns as yeast passed through subsequent cell cycles. While most H3K4me3 patterns were conserved from one generation to the next, we found that methylation patterns induced by alpha factor or high temperature were erased within one cell cycle, during S phase. Early-replicating regions were erased before late-replicating regions, implicating replication in H3K4me3 loss. However, nearly complete H3K4me3 erasure occurred at the majority of loci even when replication was prevented, suggesting that most erasure results from an active process. Indeed, deletion of the demethylase Jhd2 slowed erasure at most loci. Together, these results indicate overlapping roles for passive dilution and active enzymatic demethylation in erasing ancestral histone methylation states in yeast

The Genomic Distribution and Function of Histone Variant HTZ-1 during C. elegans Embryogenesis

In all eukaryotes, histone variants are incorporated into a subset of nucleosomes to create functionally specialized regions of chromatin. One such variant, H2A.Z, replaces histone H2A and is required for development and viability in all animals tested to date. However, the function of H2A.Z in development remains unclear. Here, we use ChIP-chip, genetic mutation, RNAi, and immunofluorescence microscopy to interrogate the function of H2A.Z (HTZ-1) during embryogenesis in Caenorhabditis elegans, a key model of metazoan development. We find that HTZ-1 is expressed in every cell of the developing embryo and is essential for normal development. The sites of HTZ-1 incorporation during embryogenesis reveal a genome wrought by developmental processes. HTZ-1 is incorporated upstream of 23% of C. elegans genes. While these genes tend to be required for development and occupied by RNA polymerase II, HTZ-1 incorporation does not specify a stereotypic transcription program. The data also provide evidence for unexpectedly widespread independent regulation of genes within operons during development; in 37% of operons, HTZ-1 is incorporated upstream of internally encoded genes. Fewer sites of HTZ-1 incorporation occur on the X chromosome relative to autosomes, which our data suggest is due to a paucity of developmentally important genes on X, rather than a direct function for HTZ-1 in dosage compensation. Our experiments indicate that HTZ-1 functions in establishing or maintaining an essential chromatin state at promoters regulated dynamically during C. elegans embryogenesis

Recruitment of Histone Deacetylase 3 to the Interferon-A Gene Promoters Attenuates Interferon Expression

Induction of Type I Interferon (IFN) genes constitutes an essential step leading to innate immune responses during virus infection. Sendai virus (SeV) infection of B lymphoid Namalwa cells transiently induces the transcriptional expression of multiple IFN-A genes. Although transcriptional activation of IFN-A genes has been extensively studied, the mechanism responsible for the attenuation of their expression remains to be determined.In this study, we demonstrate that virus infection of Namalwa cells induces transient recruitment of HDAC3 (histone deacetylase 3) to IFN-A promoters. Analysis of chromatin-protein association by Chip-QPCR demonstrated that recruitment of interferon regulatory factor (IRF)3 and IRF7, as well as TBP correlated with enhanced histone H3K9 and H3K14 acetylation, whereas recruitment of HDAC3 correlated with inhibition of histone H3K9/K14 acetylation, removal of IRF7 and TATA-binding protein (TBP) from IFN-A promoters and inhibition of virus-induced IFN-A gene transcription. Additionally, HDAC3 overexpression reduced, and HDAC3 depletion by siRNA enhanced IFN-A gene expression. Furthermore, activation of IRF7 enhanced histone H3K9/K14 acetylation and IFN-A gene expression, whereas activation of both IRF7 and IRF3 led to recruitment of HDAC3 to the IFN-A gene promoters, resulting in impaired histone H3K9 acetylation and attenuation of IFN-A gene transcription.Altogether these data indicate that reversal of histone H3K9/K14 acetylation by HDAC3 is required for attenuation of IFN-A gene transcription during viral infection

Dynamic histone H3 methylation during gene induction: HYPB/Setd2 mediates all H3K36 trimethylation

Understanding the function of histone modifications across inducible genes in mammalian cells requires quantitative, comparative analysis of their fate during gene activation and identification of enzymes responsible. We produced high-resolution comparative maps of the distribution and dynamics of H3K4me3, H3K36me3, H3K79me2 and H3K9ac across c-fos and c-jun upon gene induction in murine fibroblasts. In unstimulated cells, continuous turnover of H3K9 acetylation occurs on all K4-trimethylated histone H3 tails; distribution of both modifications coincides across promoter and 5′ part of the coding region. In contrast, K36- and K79-methylated H3 tails, which are not dynamically acetylated, are restricted to the coding regions of these genes. Upon stimulation, transcription-dependent increases in H3K4 and H3K36 trimethylation are seen across coding regions, peaking at 5′ and 3′ ends, respectively. Addressing molecular mechanisms involved, we find that Huntingtin-interacting protein HYPB/Setd2 is responsible for virtually all global and transcription-dependent H3K36 trimethylation, but not H3K36-mono- or dimethylation, in these cells. These studies reveal four distinct layers of histone modification across inducible mammalian genes and show that HYPB/Setd2 is responsible for H3K36 trimethylation throughout the mouse nucleus

Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data

<p>Abstract</p> <p>Background</p> <p>Initiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models.</p> <p>Results</p> <p>Here, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).</p> <p>As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling.</p> <p>Conclusions</p> <p>In this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.</p