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

Physical properties of naked DNA influence nucleosome positioning and correlate with transcription start and termination sites in yeast

Abstract Background In eukaryotic organisms, DNA is packaged into chromatin structure, where most of DNA is wrapped into nucleosomes. DNA compaction and nucleosome positioning have clear functional implications, since they modulate the accessibility of genomic regions to regulatory proteins. Despite the intensive research effort focused in this area, the rules defining nucleosome positioning and the location of DNA regulatory regions still remain elusive. Results Naked (histone-free) and nucleosomal DNA from yeast were digested by microccocal nuclease (MNase) and sequenced genome-wide. MNase cutting preferences were determined for both naked and nucleosomal DNAs. Integration of their sequencing profiles with DNA conformational descriptors derived from atomistic molecular dynamic simulations enabled us to extract the physical properties of DNA on a genomic scale and to correlate them with chromatin structure and gene regulation. The local structure of DNA around regulatory regions was found to be unusually flexible and to display a unique pattern of nucleosome positioning. Ab initio physical descriptors derived from molecular dynamics were used to develop a computational method that accurately predicts nucleosome enriched and depleted regions. Conclusions Our experimental and computational analyses jointly demonstrate a clear correlation between sequence-dependent physical properties of naked DNA and regulatory signals in the chromatin structure. These results demonstrate that nucleosome positioning around TSS (Transcription Start Site) and TTS (Transcription Termination Site) (at least in yeast) is strongly dependent on DNA physical properties, which can define a basal regulatory mechanism of gene expression

Positional variations among heterogeneous nucleosome maps give dynamical information on chromatin

Although nucleosome remodeling is essential to transcriptional regulation in eukaryotes, little is known about its genome-wide behavior. Since a number of nucleosome positioning maps in vivo have been recently determined, we examined if their comparisons might be used for obtaining a genome-wide profile of nucleosome remodeling. Using seven yeast maps, the local variability of nucleosomes, measured by the entropy, was significantly higher in a set of reported unstable nucleosomes. The binding sites of four transcription factors, known as the remodeling factors, were distinctively high both in entropy and linker ratio, whereas those of Yhp1, their potential inhibitor, showed the lowest values in both of them. Taken together, our map shows the general information of nucleosome dynamics reasonably well. The “nucleosome dynamics” map provides the new significant correlation with the degree of expression variety instead of their intensity. Furthermore, the associations with gene function and histone modification were also discussed here

The Effect of Micrococcal Nuclease Digestion on Nucleosome Positioning Data

Eukaryotic genomes are packed into chromatin, whose basic repeating unit is the nucleosome. Nucleosome positioning is a widely researched area. A common experimental procedure to determine nucleosome positions involves the use of micrococcal nuclease (MNase). Here, we show that the cutting preference of MNase in combination with size selection generates a sequence-dependent bias in the resulting fragments. This strongly affects nucleosome positioning data and especially sequence-dependent models for nucleosome positioning. As a consequence we see a need to re-evaluate whether the DNA sequence is a major determinant of nucleosome positioning in vivo. More generally, our results show that data generated after MNase digestion of chromatin requires a matched control experiment in order to determine nucleosome positions

ChIPseqR: analysis of ChIP-seq experiments

<p>Abstract</p> <p>Background</p> <p>The use of high-throughput sequencing in combination with chromatin immunoprecipitation (ChIP-seq) has enabled the study of genome-wide protein binding at high resolution. While the amount of data generated from such experiments is steadily increasing, the methods available for their analysis remain limited. Although several algorithms for the analysis of ChIP-seq data have been published they focus almost exclusively on transcription factor studies and are usually not well suited for the analysis of other types of experiments.</p> <p>Results</p> <p>Here we present ChIPseqR, an algorithm for the analysis of nucleosome positioning and histone modification ChIP-seq experiments. The performance of this novel method is studied on short read sequencing data of <it>Arabidopsis thaliana </it>mononucleosomes as well as on simulated data.</p> <p>Conclusions</p> <p>ChIPseqR is shown to improve sensitivity and spatial resolution over existing methods while maintaining high specificity. Further analysis of predicted nucleosomes reveals characteristic patterns in nucleosome sequences and placement.</p

Epigenetic Transcriptional Regulation of the Growth Arrest-Specific gene 1 (Gas1) in Hepatic Cell Proliferation at Mononucleosomal Resolution

BACKGROUND: Gas1 (growth arrest-specific 1) gene is known to inhibit cell proliferation in a variety of models, but its possible implication in regulating quiescence in adult tissues has not been examined to date. The knowledge of how Gas1 is regulated in quiescence may contribute to understand the deregulation occurring in neoplastic diseases. METHODOLOGY/PRINCIPAL FINDINGS: Gas1 expression has been studied in quiescent murine liver and during the naturally synchronized cell proliferation after partial hepatectomy. Chromatin immunoprecipitation at nucleosomal resolution (Nuc-ChIP) has been used to carry out the study preserving the in vivo conditions. Transcription has been assessed at real time by quantifying the presence of RNA polymerase II in coding regions (RNApol-ChIP). It has been found that Gas1 is expressed not only in quiescent liver but also at the cell cycle G(1)/S transition. The latter expression peak had not been previously reported. Two nucleosomes, flanking a nucleosome-free region, are positioned close to the transcription start site. Both nucleosomes slide in going from the active to the inactive state and vice versa. Nuc-ChIP analysis of the acquisition of histone epigenetic marks show distinctive features in both active states: H3K9ac and H3K4me2 are characteristic of transcription in G(0) and H4R3me2 in G(1)/S transition. Sequential-ChIP analysis revealed that the "repressing" mark H3K9me2 colocalize with several "activating" marks at nucleosome N-1 when Gas1 is actively transcribed suggesting a greater plasticity of epigenetic marks than proposed until now. The recruitment of chromatin-remodeling or modifying complexes also displayed distinct characteristics in quiescence and the G(1)/S transition. CONCLUSIONS/SIGNIFICANCE: The finding that Gas1 is transcribed at the G(1)/S transition suggests that the gene may exert a novel function during cell proliferation. Transcription of this gene is modulated by specific "activating" and "repressing" epigenetic marks, and by chromatin remodeling and histone modifying complexes recruitment, at specific nucleosomes in Gas1 promoter

Contribution of histone sequence preferences to nucleosome organization: proposed definitions and methodology

We propose definitions and procedures for comparing nucleosome maps and discuss current agreement and disagreement on the effect of histone sequence preferences on nucleosome organization in vivo

Distinct Mechanisms for Induction and Tolerance Regulate the Immediate Early Genes Encoding Interleukin 1β and Tumor Necrosis Factor α

Interleukin-1β and Tumor Necrosis Factor α play related, but distinct, roles in immunity and disease. Our study revealed major mechanistic distinctions in the Toll-like receptor (TLR) signaling-dependent induction for the rapidly expressed genes (IL1B and TNF) coding for these two cytokines. Prior to induction, TNF exhibited pre-bound TATA Binding Protein (TBP) and paused RNA Polymerase II (Pol II), hallmarks of poised immediate-early (IE) genes. In contrast, unstimulated IL1B displayed very low levels of both TBP and paused Pol II, requiring the lineage-specific Spi-1/PU.1 (Spi1) transcription factor as an anchor for induction-dependent interaction with two TLR-activated transcription factors, C/EBPβ and NF-κB. Activation and DNA binding of these two pre-expressed factors resulted in de novo recruitment of TBP and Pol II to IL1B in concert with a permissive state for elongation mediated by the recruitment of elongation factor P-TEFb. This Spi1-dependent mechanism for IL1B transcription, which is unique for a rapidly-induced/poised IE gene, was more dependent upon P-TEFb than was the case for the TNF gene. Furthermore, the dependence on phosphoinositide 3-kinase for P-TEFb recruitment to IL1B paralleled a greater sensitivity to the metabolic state of the cell and a lower sensitivity to the phenomenon of endotoxin tolerance than was evident for TNF. Such differences in induction mechanisms argue against the prevailing paradigm that all IE genes possess paused Pol II and may further delineate the specific roles played by each of these rapidly expressed immune modulators. © 2013 Adamik et al

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

Oligonucleotide Sequence Motifs as Nucleosome Positioning Signals

To gain a better understanding of the sequence patterns that characterize positioned nucleosomes, we first performed an analysis of the periodicities of the 256 tetranucleotides in a yeast genome-wide library of nucleosomal DNA sequences that was prepared by in vitro reconstitution. The approach entailed the identification and analysis of 24 unique tetranucleotides that were defined by 8 consensus sequences. These consensus sequences were shown to be responsible for most if not all of the tetranucleotide and dinucleotide periodicities displayed by the entire library, demonstrating that the periodicities of dinucleotides that characterize the yeast genome are, in actuality, due primarily to the 8 consensus sequences. A novel combination of experimental and bioinformatic approaches was then used to show that these tetranucleotides are important for preferred formation of nucleosomes at specific sites along DNA in vitro. These results were then compared to tetranucleotide patterns in genome-wide in vivo libraries from yeast and C. elegans in order to assess the contributions of DNA sequence in the control of nucleosome residency in the cell. These comparisons revealed striking similarities in the tetranucleotide occurrence profiles that are likely to be involved in nucleosome positioning in both in vitro and in vivo libraries, suggesting that DNA sequence is an important factor in the control of nucleosome placement in vivo. However, the strengths of the tetranucleotide periodicities were 3–4 fold higher in the in vitro as compared to the in vivo libraries, which implies that DNA sequence plays less of a role in dictating nucleosome positions in vivo. The results of this study have important implications for models of sequence-dependent positioning since they suggest that a defined subset of tetranucleotides is involved in preferred nucleosome occupancy and that these tetranucleotides are the major source of the dinucleotide periodicities that are characteristic of positioned nucleosomes