Quantitative test of the barrier nucleosome model for statistical positioning of nucleosomes up- and downstream of transcription start sites

The positions of nucleosomes in eukaryotic genomes determine which parts of the DNA sequence are readily accessible for regulatory proteins and which are not. Genome-wide maps of nucleosome positions have revealed a salient pattern around transcription start sites, involving a nucleosome-free region (NFR) flanked by a pronounced periodic pattern in the average nucleosome density. While the periodic pattern clearly reflects well-positioned nucleosomes, the positioning mechanism is less clear. A recent experimental study by Mavrich et al. argued that the pattern observed in S. cerevisiae is qualitatively consistent with a `barrier nucleosome model', in which the oscillatory pattern is created by the statistical positioning mechanism of Kornberg and Stryer. On the other hand, there is clear evidence for intrinsic sequence preferences of nucleosomes, and it is unclear to what extent these sequence preferences affect the observed pattern. To test the barrier nucleosome model, we quantitatively analyze yeast nucleosome positioning data both up- and downstream from NFRs. Our analysis is based on the Tonks model of statistical physics which quantifies the interplay between the excluded-volume interaction of nucleosomes and their positional entropy. We find that although the typical patterns on the two sides of the NFR are different, they are both quantitatively described by the same physical model, with the same parameters, but different boundary conditions. The inferred boundary conditions suggest that the first nucleosome downstream from the NFR (the +1 nucleosome) is typically directly positioned while the first nucleosome upstream is statistically positioned via a nucleosome-repelling DNA region. These boundary conditions, which can be locally encoded into the genome sequence, significantly shape the statistical distribution of nucleosomes over a range of up to ~1000 bp to each side.Comment: includes supporting materia

Monte Carlo Simulations indicate that Chromati: Nanostructure is accessible by Light Microscopy

A long controversy exists about the structure of chromatin. Theoretically, this structure could be resolved by scattering experiments if one determines the scattering function - or equivalently the pair distribution function - of the nucleosomes. Unfortunately, scattering experiments with live cells are very difficult and limited to only a couple of nucleosomes

The Impact of Local Genome Sequence on Defining Heterochromatin Domains

Characterizing how genomic sequence interacts with trans-acting regulatory factors to implement a program of gene expression in eukaryotic organisms is critical to understanding genome function. One means by which patterns of gene expression are achieved is through the differential packaging of DNA into distinct types of chromatin. While chromatin state exerts a major influence on gene expression, the extent to which cis-acting DNA sequences contribute to the specification of chromatin state remains incompletely understood. To address this, we have used a fission yeast sequence element (L5), known to be sufficient to nucleate heterochromatin, to establish de novo heterochromatin domains in the Schizosaccharomyces pombe genome. The resulting heterochromatin domains were queried for the presence of H3K9 di-methylation and Swi6p, both hallmarks of heterochromatin, and for levels of gene expression. We describe a major effect of genomic sequences in determining the size and extent of such de novo heterochromatin domains. Heterochromatin spreading is antagonized by the presence of genes, in a manner that can occur independent of strength of transcription. Increasing the dosage of Swi6p results in increased heterochromatin proximal to the L5 element, but does not result in an expansion of the heterochromatin domain, suggesting that in this context genomic effects are dominant over trans effects. Finally, we show that the ratio of Swi6p to H3K9 di-methylation is sequence-dependent and correlates with the extent of gene repression. Taken together, these data demonstrate that the sequence content of a genomic region plays a significant role in shaping its response to encroaching heterochromatin and suggest a role of DNA sequence in specifying chromatin state

Transcriptional interaction-assisted identification of dynamic nucleosome positioning

<p>Abstract</p> <p>Background</p> <p>Nucleosomes regulate DNA accessibility and therefore play a central role in transcription control. Computational methods have been developed to predict static nucleosome positions from DNA sequences, but nucleosomes are dynamic in vivo.</p> <p>Results</p> <p>Motivated by our observation that transcriptional interaction is discriminative information for nucleosome occupancy, we developed a novel computational approach to identify dynamic nucleosome positions at promoters by combining transcriptional interaction and genomic sequence information. Our approach successfully identified experimentally determined nucleosome positioning dynamics available in three cellular conditions, and significantly improved the prediction accuracy which is based on sequence information alone. We then applied our approach to various cellular conditions and established a comprehensive landscape of dynamic nucleosome positioning in yeast.</p> <p>Conclusion</p> <p>Analysis of this landscape revealed that the majority of nucleosome positions are maintained during most conditions. However, nucleosome occupancy at most promoters fluctuates with the corresponding gene expression level and is reduced specifically at the phase of peak expression. Further investigation into properties of nucleosome occupancy identified two gene groups associated with distinct modes of nucleosome modulation. Our results suggest that both the intrinsic sequence and regulatory proteins modulate nucleosomes in an altered manner.</p

Silent chromatin at the middle and ends: lessons from yeasts

Eukaryotic centromeres and telomeres are specialized chromosomal regions that share one common characteristic: their underlying DNA sequences are assembled into heritably repressed chromatin. Silent chromatin in budding and fission yeast is composed of fundamentally divergent proteins tat assemble very different chromatin structures. However, the ultimate behaviour of silent chromatin and the pathways that assemble it seem strikingly similar among Saccharomyces cerevisiae (S. cerevisiae), Schizosaccharomyces pombe (S. pombe) and other eukaryotes. Thus, studies in both yeasts have been instrumental in dissecting the mechanisms that establish and maintain silent chromatin in eukaryotes, contributing substantially to our understanding of epigenetic processes. In this review, we discuss current models for the generation of heterochromatic domains at centromeres and telomeres in the two yeast species

The C-Terminus of Histone H2B Is Involved in Chromatin Compaction Specifically at Telomeres, Independently of Its Monoubiquitylation at Lysine 123

Telomeric heterochromatin assembly in budding yeast propagates through the association of Silent Information Regulator (SIR) proteins with nucleosomes, and the nucleosome array has been assumed to fold into a compacted structure. It is believed that the level of compaction and gene repression within heterochromatic regions can be modulated by histone modifications, such as acetylation of H3 lysine 56 and H4 lysine 16, and monoubiquitylation of H2B lysine 123. However, it remains unclear as to whether or not gene silencing is a direct consequence of the compaction of chromatin. Here, by investigating the role of the carboxy-terminus of histone H2B in heterochromatin formation, we identify that the disorderly compaction of chromatin induced by a mutation at H2B T122 specifically hinders telomeric heterochromatin formation. H2B T122 is positioned within the highly conserved AVTKY motif of the αC helix of H2B. Heterochromatin containing the T122E substitution in H2B remains inaccessible to ectopic dam methylase with dramatically increased mobility in sucrose gradients, indicating a compacted chromatin structure. Genetic studies indicate that this unique phenotype is independent of H2B K123 ubiquitylation and Sir4. In addition, using ChIP analysis, we demonstrate that telomere structure in the mutant is further disrupted by a defect in Sir2/Sir3 binding and the resulting invasion of euchromatic histone marks. Thus, we have revealed that the compaction of chromatin per se is not sufficient for heterochromatin formation. Instead, these results suggest that an appropriately arrayed chromatin mediated by H2B C-terminus is required for SIR binding and the subsequent formation of telomeric chromatin in yeast, thereby identifying an intrinsic property of the nucleosome that is required for the establishment of telomeric heterochromatin. This requirement is also likely to exist in higher eukaryotes, as the AVTKY motif of H2B is evolutionarily conserved

Differential analysis for high density tiling microarray data

<p>Abstract</p> <p>Background</p> <p>High density oligonucleotide tiling arrays are an effective and powerful platform for conducting unbiased genome-wide studies. The <it>ab initio </it>probe selection method employed in tiling arrays is unbiased, and thus ensures consistent sampling across coding and non-coding regions of the genome. These arrays are being increasingly used to study the associated processes of transcription, transcription factor binding, chromatin structure and their association. Studies of differential expression and/or regulation provide critical insight into the mechanics of transcription and regulation that occurs during the developmental program of a cell. The time-course experiment, which comprises an <it>in-vivo </it>system and the proposed analyses, is used to determine if annotated and un-annotated portions of genome manifest coordinated differential response to the induced developmental program.</p> <p>Results</p> <p>We have proposed a novel approach, based on a piece-wise function – to analyze genome-wide differential response. This enables segmentation of the response based on protein-coding and non-coding regions; for genes the methodology also partitions differential response with a 5' versus 3' versus intra-genic bias.</p> <p>Conclusion</p> <p>The algorithm built upon the framework of Significance Analysis of Microarrays, uses a generalized logic to define regions/patterns of coordinated differential change. By not adhering to the gene-centric paradigm, discordant differential expression patterns between exons and introns have been identified at a FDR of less than 12 percent. A co-localization of differential binding between RNA Polymerase II and tetra-acetylated histone has been quantified at a p-value < 0.003; it is most significant at the 5' end of genes, at a p-value < 10^-13. The prototype R code has been made available as supplementary material [see Additional file <supplr sid="S1">1</supplr>].</p> <suppl id="S1"> <title> <p>Additional file 1</p> </title> <text> <p>gsam_prototypercode.zip. File archive comprising of prototype R code for gSAM implementation including readme and examples.</p> </text> <file name="1471-2105-8-359-S1.zip"> <p>Click here for file</p> </file> </suppl

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

Histone H4 Lysine 12 Acetylation Regulates Telomeric Heterochromatin Plasticity in Saccharomyces cerevisiae

Recent studies have established that the highly condensed and transcriptionally silent heterochromatic domains in budding yeast are virtually dynamic structures. The underlying mechanisms for heterochromatin dynamics, however, remain obscure. In this study, we show that histones are dynamically acetylated on H4K12 at telomeric heterochromatin, and this acetylation regulates several of the dynamic telomere properties. Using a de novo heterochromatin formation assay, we surprisingly found that acetylated H4K12 survived the formation of telomeric heterochromatin. Consistently, the histone acetyltransferase complex NuA4 bound to silenced telomeric regions and acetylated H4K12. H4K12 acetylation prevented the over-accumulation of Sir proteins at telomeric heterochromatin and elimination of this acetylation caused defects in multiple telomere-related processes, including transcription, telomere replication, and recombination. Together, these data shed light on a potential histone acetylation mark within telomeric heterochromatin that contributes to telomere plasticity