Podbat: A Novel Genomic Tool Reveals Swr1-Independent H2A.Z Incorporation at Gene Coding Sequences through Epigenetic Meta-Analysis

Epigenetic regulation consists of a multitude of different modifications that determine active and inactive states of chromatin. Conditions such as cell differentiation or exposure to environmental stress require concerted changes in gene expression. To interpret epigenomics data, a spectrum of different interconnected datasets is needed, ranging from the genome sequence and positions of histones, together with their modifications and variants, to the transcriptional output of genomic regions. Here we present a tool, Podbat (Positioning database and analysis tool), that incorporates data from various sources and allows detailed dissection of the entire range of chromatin modifications simultaneously. Podbat can be used to analyze, visualize, store and share epigenomics data. Among other functions, Podbat allows data-driven determination of genome regions of differential protein occupancy or RNA expression using Hidden Markov Models. Comparisons between datasets are facilitated to enable the study of the comprehensive chromatin modification system simultaneously, irrespective of data-generating technique. Any organism with a sequenced genome can be accommodated. We exemplify the power of Podbat by reanalyzing all to-date published genome-wide data for the histone variant H2A.Z in fission yeast together with other histone marks and also phenotypic response data from several sources. This meta-analysis led to the unexpected finding of H2A.Z incorporation in the coding regions of genes encoding proteins involved in the regulation of meiosis and genotoxic stress responses. This incorporation was partly independent of the H2A.Z-incorporating remodeller Swr1. We verified an Swr1-independent role for H2A.Z following genotoxic stress in vivo. Podbat is open source software freely downloadable from www.podbat.org, distributed under the GNU LGPL license. User manuals, test data and instructions are available at the website, as well as a repository for third party–developed plug-in modules. Podbat requires Java version 1.6 or higher

The Insulator Binding Protein CTCF Positions 20 Nucleosomes around Its Binding Sites across the Human Genome

Chromatin structure plays an important role in modulating the accessibility of genomic DNA to regulatory proteins in eukaryotic cells. We performed an integrative analysis on dozens of recent datasets generated by deep-sequencing and high-density tiling arrays, and we discovered an array of well-positioned nucleosomes flanking sites occupied by the insulator binding protein CTCF across the human genome. These nucleosomes are highly enriched for the histone variant H2A.Z and 11 histone modifications. The distances between the center positions of the neighboring nucleosomes are largely invariant, and we estimate them to be 185 bp on average. Surprisingly, subsets of nucleosomes that are enriched in different histone modifications vary greatly in the lengths of DNA protected from micrococcal nuclease cleavage (106–164 bp). The nucleosomes enriched in those histone modifications previously implicated to be correlated with active transcription tend to contain less protected DNA, indicating that these modifications are correlated with greater DNA accessibility. Another striking result obtained from our analysis is that nucleosomes flanking CTCF sites are much better positioned than those downstream of transcription start sites, the only genomic feature previously known to position nucleosomes genome-wide. This nucleosome-positioning phenomenon is not observed for other transcriptional factors for which we had genome-wide binding data. We suggest that binding of CTCF provides an anchor point for positioning nucleosomes, and chromatin remodeling is an important component of CTCF function

The SWR1 Histone Replacement Complex Causes Genetic Instability and Genome-Wide Transcription Misregulation in the Absence of H2A.Z

The SWR1 complex replaces the canonical histone H2A with the variant H2A.Z (Htz1 in yeast) at specific chromatin regions. This dynamic alteration in nucleosome structure provides a molecular mechanism to regulate transcription, gene silencing, chromosome segregation and DNA repair. Here we show that genetic instability, sensitivity to drugs impairing different cellular processes and genome-wide transcriptional misregulation in htz1Δ can be partially or totally suppressed if SWR1 is not formed (swr1Δ), if it forms but cannot bind to chromatin (swc2Δ) or if it binds to chromatin but lacks histone replacement activity (swc5Δ and the ATPase-dead swr1-K727G). These results suggest that in htz1Δ the nucleosome remodelling activity of SWR1 affects chromatin integrity because of an attempt to replace H2A with Htz1 in the absence of the latter. This would impair transcription and, either directly or indirectly, other cellular processes. Specifically, we show that in htz1Δ, the SWR1 complex causes an accumulation of recombinogenic DNA damage by a mechanism dependent on phosphorylation of H2A at Ser129, a modification that occurs in response to DNA damage, suggesting that the SWR1 complex impairs the repair of spontaneous DNA damage in htz1Δ. In addition, SWR1 causes DSBs sensitivity in htz1Δ; consistently, in the absence of Htz1 the SWR1 complex bound near an endonuclease HO-induced DSB at the mating-type (MAT) locus impairs DSB-induced checkpoint activation. Our results support a stepwise mechanism for the replacement of H2A with Htz1 and demonstrate that a tight control of this mechanism is essential to regulate chromatin dynamics but also to prevent the deleterious consequences of an incomplete nucleosome remodelling

The Euchromatic and Heterochromatic Landscapes Are Shaped by Antagonizing Effects of Transcription on H2A.Z Deposition

A role for variant histone H2A.Z in gene expression is now well established but little is known about the mechanisms by which it operates. Using a combination of ChIP–chip, knockdown and expression profiling experiments, we show that upon gene induction, human H2A.Z associates with gene promoters and helps in recruiting the transcriptional machinery. Surprisingly, we also found that H2A.Z is randomly incorporated in the genome at low levels and that active transcription antagonizes this incorporation in transcribed regions. After cessation of transcription, random H2A.Z quickly reappears on genes, demonstrating that this incorporation utilizes an active mechanism. Within facultative heterochromatin, we observe a hyper accumulation of the variant histone, which might be due to the lack of transcription in these regions. These results show how chromatin structure and transcription can antagonize each other, therefore shaping chromatin and controlling gene expression

Predicting Human Nucleosome Occupancy from Primary Sequence

Nucleosomes are the fundamental repeating unit of chromatin and comprise the structural building blocks of the living eukaryotic genome. Micrococcal nuclease (MNase) has long been used to delineate nucleosomal organization. Microarray-based nucleosome mapping experiments in yeast chromatin have revealed regularly-spaced translational phasing of nucleosomes. These data have been used to train computational models of sequence-directed nuclesosome positioning, which have identified ubiquitous strong intrinsic nucleosome positioning signals. Here, we successfully apply this approach to nucleosome positioning experiments from human chromatin. The predictions made by the human-trained and yeast-trained models are strongly correlated, suggesting a shared mechanism for sequence-based determination of nucleosome occupancy. In addition, we observed striking complementarity between classifiers trained on experimental data from weakly versus heavily digested MNase samples. In the former case, the resulting model accurately identifies nucleosome-forming sequences; in the latter, the classifier excels at identifying nucleosome-free regions. Using this model we are able to identify several characteristics of nucleosome-forming and nucleosome-disfavoring sequences. First, by combining results from each classifier applied de novo across the human ENCODE regions, the classifier reveals distinct sequence composition and periodicity features of nucleosome-forming and nucleosome-disfavoring sequences. Short runs of dinucleotide repeat appear as a hallmark of nucleosome-disfavoring sequences, while nucleosome-forming sequences contain short periodic runs of GC base pairs. Second, we show that nucleosome phasing is most frequently predicted flanking nucleosome-free regions. The results suggest that the major mechanism of nucleosome positioning in vivo is boundary-event-driven and affirm the classical statistical positioning theory of nucleosome organization

Transcription Initiation Patterns Indicate Divergent Strategies for Gene Regulation at the Chromatin Level

The application of deep sequencing to map 5′ capped transcripts has confirmed the existence of at least two distinct promoter classes in metazoans: “focused” promoters with transcription start sites (TSSs) that occur in a narrowly defined genomic span and “dispersed” promoters with TSSs that are spread over a larger window. Previous studies have explored the presence of genomic features, such as CpG islands and sequence motifs, in these promoter classes, but virtually no studies have directly investigated the relationship with chromatin features. Here, we show that promoter classes are significantly differentiated by nucleosome organization and chromatin structure. Dispersed promoters display higher associations with well-positioned nucleosomes downstream of the TSS and a more clearly defined nucleosome free region upstream, while focused promoters have a less organized nucleosome structure, yet higher presence of RNA polymerase II. These differences extend to histone variants (H2A.Z) and marks (H3K4 methylation), as well as insulator binding (such as CTCF), independent of the expression levels of affected genes. Notably, differences are conserved across mammals and flies, and they provide for a clearer separation of promoter architectures than the presence and absence of CpG islands or the occurrence of stalled RNA polymerase. Computational models support the stronger contribution of chromatin features to the definition of dispersed promoters compared to focused start sites. Our results show that promoter classes defined from 5′ capped transcripts not only reflect differences in the initiation process at the core promoter but also are indicative of divergent transcriptional programs established within gene-proximal nucleosome organization

Genomic Organization of H2Av Containing Nucleosomes in Drosophila Heterochromatin

H2Av is a versatile histone variant that plays both positive and negative roles in transcription, DNA repair, and chromatin structure in Drosophila. H2Av, and its broader homolog H2A.Z, tend to be enriched toward 5′ ends of genes, and exist in both euchromatin and heterochromatin. Its organization around euchromatin genes and other features have been described in many eukaryotic model organisms. However, less is known about H2Av nucleosome organization in heterochromatin. Here we report the properties and organization of individual H2Av nucleosomes around genes and transposable elements located in Drosophila heterochromatic regions. We compare the similarity and differences with that found in euchromatic regions. Our analyses suggest that nucleosomes are intrinsically positioned on inverted repeats of DNA transposable elements such as those related to the “1360” element, but are not intrinsically positioned on retrotransposon-related elements

In Vivo Expression of MHC Class I Genes Depends on the Presence of a Downstream Barrier Element

Regulation of MHC class I gene expression is critical to achieve proper immune surveillance. In this work, we identify elements downstream of the MHC class I promoter that are necessary for appropriate in vivo regulation: a novel barrier element that protects the MHC class I gene from silencing and elements within the first two introns that contribute to tissue specific transcription. The barrier element is located in intergenic sequences 3′ to the polyA addition site. It is necessary for stable expression in vivo, but has no effect in transient transfection assays. Accordingly, in both transgenic mice and stably transfected cell lines, truncation of the barrier resulted in transcriptional gene silencing, increased nucleosomal density and decreased histone H3K9/K14 acetylation and H3K4 di-methylation across the gene. Significantly, distinct sequences within the barrier element govern anti-silencing and chromatin modifications. Thus, this novel barrier element functions to maintain transcriptionally permissive chromatin organization and prevent transcriptional silencing of the MHC class I gene, ensuring it is poised to respond to immune signaling

PRC1 and PRC2 Are Not Required for Targeting of H2A.Z to Developmental Genes in Embryonic Stem Cells

The essential histone variant H2A.Z localises to both active and silent chromatin sites. In embryonic stem cells (ESCs), H2A.Z is also reported to co-localise with polycomb repressive complex 2 (PRC2) at developmentally silenced genes. The mechanism of H2A.Z targeting is not clear, but a role for the PRC2 component Suz12 has been suggested. Given this association, we wished to determine if polycomb functionally directs H2A.Z incorporation in ESCs. We demonstrate that the PRC1 component Ring1B interacts with multiple complexes in ESCs. Moreover, we show that although the genomic distribution of H2A.Z co-localises with PRC2, Ring1B and with the presence of CpG islands, H2A.Z still blankets polycomb target loci in the absence of Suz12, Eed (PRC2) or Ring1B (PRC1). Therefore we conclude that H2A.Z accumulates at developmentally silenced genes in ESCs in a polycomb independent manner

Dissecting Epigenetic Silencing Complexity in the Mouse Lung Cancer Suppressor Gene Cadm1

Disease-oriented functional analysis of epigenetic factors and their regulatory mechanisms in aberrant silencing is a prerequisite for better diagnostics and therapy. Yet, the precise mechanisms are still unclear and complex, involving the interplay of several effectors including nucleosome positioning, DNA methylation, histone variants and histone modifications. We investigated the epigenetic silencing complexity in the tumor suppressor gene Cadm1 in mouse lung cancer progenitor cell lines, exhibiting promoter hypermethylation associated with transcriptional repression, but mostly unresponsive to demethylating drug treatments. After predicting nucleosome positions and transcription factor binding sites along the Cadm1 promoter, we carried out single-molecule mapping with DNA methyltransferase M.SssI, which revealed in silent promoters high nucleosome occupancy and occlusion of transcription factor binding sites. Furthermore, M.SssI maps of promoters varied within and among the different lung cancer cell lines. Chromatin analysis with micrococcal nuclease also indicated variations in nucleosome positioning to have implications in the binding of transcription factors near nucleosome borders. Chromatin immunoprecipitation showed that histone variants (H2A.Z and H3.3), and opposing histone modification marks (H3K4me3 and H3K27me3) all colocalized in the same nucleosome positions that is reminiscent of epigenetic plasticity in embryonic stem cells. Altogether, epigenetic silencing complexity in the promoter region of Cadm1 is not only defined by DNA hypermethylation, but high nucleosome occupancy, altered nucleosome positioning, and ‘bivalent’ histone modifications, also likely contributed in the transcriptional repression of this gene in the lung cancer cells. Our results will help define therapeutic intervention strategies using epigenetic drugs in lung cancer