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.     

Are nucleosome positions in vivo primarily determined by histone–DNA sequence preferences?

Large-scale and genome-wide studies have concluded that ∼80% of the yeast (Saccharomyces cerevisiae) genome is occupied by positioned nucleosomes. In vivo this nucleosome organization can result from a variety of mechanisms, including the intrinsic DNA sequence preferences for wrapping the DNA around the histone core. Recently, a genome-wide study was reported using massively parallel sequencing to directly compare in vivo and in vitro nucleosome positions. It was concluded that intrinsic DNA sequence preferences indeed have a dominant role in determining the in vivo nucleosome organization of the genome, consistent with a genomic code for nucleosome positioning. Some other studies disagree with this view. Using the large amount of data now available from several sources, we have attempted to clarify a fundamental question concerning the packaging of genomic DNA: to what extent are nucleosome positions in vivo determined by histone-DNA sequence preferences? We have analyzed data obtained from different laboratories in the same way, and have directly compared these data. We also identify possible problems with some of the experimental designs used and with the data analysis. Our findings suggest that DNA sequence preferences have only small effects on the positioning of individual nucleosomes throughout the genome in vivo.     

Sequence-dependent Nucleosome Positioning

Eukaryotic DNA is organized into a macromolecular structure called chromatin. The basic repeating unit of chromatin is the nucleosome, which consists of two copies of each of the four core histones and DNA. The nucleosomal organization and the positions of nucleosomes have profound effects on all DNA-dependent processes. Understanding the factors that influence nucleosome positioning is therefore of general interest. Among the many determinants of nucleosome positioning, the DNA sequence has been proposed to have a major role. Here, we analyzed more than 860,000 nucleosomal DNA sequences to identify sequence features that guide the formation of nucleosomes in vivo. We found that both a periodic enrichment of AT base pairs and an out-of-phase oscillating enrichment of GC base pairs as well as the overall preference for GC base pairs are determinants of nucleosome positioning. The preference for GC pairs can be related to a lower energetic cost required for deformation of the DNA to wrap around the histones. In line with this idea, we found that only incorporation of both signal components into a sequence model for nucleosome formation results in maximal predictive performance on a genome-wide scale. In this manner, one achieves greater predictive power than published approaches. Our results confirm the hypothesis that the DNA sequence has a major role in nucleosome positioning in vivo.     

Nucleosomal locations of dominant DNA sequence motifs for histone-DNA interactions and nucleosome positioning

DNA sequence is an important determinant of the positioning, stability, and activity of nucleosomes, yet the molecular basis of these effects remains elusive. A "consensus DNA sequence" for nucleosome positioning has not been reported and, while certain DNA sequence preferences or motifs for nucleosome positioning have been discovered, how they function is not known. Here, we report that an unexpected observation concerning the reassembly of nucleosomes during salt gradient dialysis has allowed a breakthrough in our efforts to identify the nucleosomal locations of the DNA sequence motifs that dominate histone-DNA interactions and nucleosome positioning. We conclude that a previous selection experiment for high-affinity, nucleosome-forming DNA sequences exerted selective pressure chiefly on the central stretch of the nucleosomal DNA. This observation implies that algorithms for aligning the selected DNA sequences should seek to optimize the alignment over much less than the full 147 bp of nucleosomal DNA. A new alignment calculation implemented these ideas and successfully aligned 19 of the 41 sequences in a non-redundant database of selected high-affinity, nucleosome-positioning sequences. The resulting alignment reveals strong conservation of several stretches within a central 71 bp of the nucleosomal DNA. The alignment further reveals an inherent palindromic symmetry in the selected DNAs; it makes testable predictions of nucleosome positioning on the aligned sequences and for the creation of new positioning sequences, both of which are upheld experimentally; and it suggests new signals that may be important in translational nucleosome positioning.     

Computational analysis suggests a highly bendable, fragile structure for nucleosomal DNA

Eukaryotic chromosomal DNA coils around histones to form nucleosomes. Although histone affinity for DNA depends on DNA sequence patterns, how nucleosome positioning is determined by them remains unknown. Here, we show relationships between nucleosome positioning and two structural characteristics of DNA conferred by DNA sequence. Analysis of bendability and hydroxyl radical cleavage intensity of nucleosomal DNA sequences indicated that nucleosomal DNA is bendable and fragile and that nucleosome positional stability was correlated with characteristics of DNA. This result explains how histone positioning is partially determined by nucleosomal DNA structure, illuminating the optimization of chromosomal DNA packaging that controls cellular dynamics.     

Prediction of nucleosome occupancy in Saccharomyces cerevisiae using position-correlation scoring function

Knowledge of the detailed organization of nucleosomes across genomes and the mechanisms of nucleosome positioning is critical for the understanding of gene regulation and expression. In the present work, the bias of 4-mer frequency in nucleosome and linker sequences of the S. cerevisiae genome was analyzed statistically. A novel position-correlation scoring function algorithm based on the bias of 4-mer frequency in linker sequences was presented to distinguish nucleosome vs linker sequences. Five-fold cross-validation demonstrated that the algorithm achieved a good performance with mean area under the receiver operator characteristics curve of 0.981. Next, the algorithm was used to predict nucleosome occupancy throughout the S. cerevisiae genome and relatively high correlation coefficients with experiment maps of nucleosome positioning were obtained. Besides, the distinct nucleosome depleted regions in the vicinity of regulatory sites were confirmed. The results suggest that intrinsic DNA sequence preferences in linker regions have a significant impact on the nucleosome occupancy.     

The core histone tail domains contribute to sequence-dependent nucleosome positioning

The precise positioning of nucleosomes plays a critical role in the regulation of gene expression by modulating the DNA binding activity of trans-acting factors. However, molecular determinants responsible for positioning are not well understood. We examined whether the removal of the core histone tail domains from nucleosomes reconstituted with specific DNA fragments led to alteration of translational positions. Remarkably, we find that removal of tail domains from a nucleosome assembled on a DNA fragment containing a Xenopus borealis somatic-type 5S RNA gene results in repositioning of nucleosomes along the DNA, including two related major translational positions that move about 20 bp further upstream with respect to the 5S gene. In a nucleosome reconstituted with a DNA fragment containing the promoter of a Drosophila alcohol dehydrogenase gene, several translational positions shifted by about 10 bp along the DNA upon tail removal. However, the positions of nucleosomes assembled with a DNA fragment known to have one of the highest binding affinities for core histone proteins in the mouse genome were not altered by removal of core histone tail domains. Our data support the notion that the basic tail domains bind to nucleosomal DNA and influence the selection of the translational position of nucleosomes and that once tails are removed movement between translational positions occurs in a facile manner on some sequences. However, the effect of the N-terminal tails on the positioning and movement of a nucleosome appears to be dependent on the DNA sequence such that the contribution of the tails can be masked by very high affinity DNA sequences. Our results suggest a mechanism whereby sequence-dependent nucleosome positioning can be specifically altered by regulated changes in histone tail-DNA interactions in chromatin.     

Both DNA and histone fold sequences contribute to archaeal nucleosome stability.

The roles and interdependence of DNA sequence and archaeal histone fold structure in determining archaeal nucleosome stability and positioning have been determined and quantitated. The presence of four tandem copies of TTTAAAGCCG in the polylinker region of pLITMUS28 resulted in a DNA molecule with increased affinity (DeltaDeltaG of approximately 700 cal mol(-1)) for the archaeal histone HMfB relative to the polylinker sequence, and the dominant, quantitative contribution of the helical repeats of the dinucleotide TA to this increased affinity has been established. The rotational and translational positioning of archaeal nucleosomes assembled on the (TTTAAAGCCG)(4) sequence and on DNA molecules selectively incorporated into archaeal nucleosomes by HMfB have been determined. Alternating A/T- and G/C-rich regions were located where the minor and major grooves, respectively, sequentially faced the archaeal nucleosome core, and identical positioning results were obtained using HMfA, a closely related archaeal histone also from Methanothermus fervidus. However, HMfA did not have similarly high affinities for the HMfB-selected DNA molecules, and domain-swap experiments have shown that this difference in affinity is determined by residue differences in the C-terminal region of alpha-helix 3 of the histone fold, a region that is not expected to directly interact with DNA. Rather this region is thought to participate in forming the histone dimer:dimer interface at the center of an archaeal nucleosome histone tetramer core. If differences in this interface do result in archaeal histone cores with different sequence preferences, then the assembly of alternative archaeal nucleosome tetramer cores could provide an unanticipated and novel structural mechanism to regulate gene expression.     

Nucleosome positioning on yeast plasmids is determined only by the internal signal of DNA sequence occupied by the nucleosome

The possible role of border factors in determining the nucleosome positioning on a DNA sequence was investigated. To this end a family of recombinant plasmids based on Gal10Cyc1 promoter and neomycin phosphotransferase gene NPTII were created. A DNA sequence adjoining the GalCyc promoter was varied in these plasmids. Three nearly equally represented nucleosome positions on the GalCyc promoter were found. In the basal plasmid an FRT sequence adjoins the GalCyc promoter at the right. It contains an internal signal of multiple positioning. Its replacement with different DNA sequences does not affect nucleosome positioning on the GalCyc promoter. The nucleosome positioning on the GalCyc promoter does not depend on nucleosome positioning (or its absence) on adjoining sequences. The same is true for nucleosome positioning on FRT sequence. It was found also that nucleosomes' positioning on the NPTII gene and their mutual disposition, namely the spacing between neighboring nucleosomes (linker length) are determined by the location of positioning signals only. Generally the nucleosome positioning in our experimental model is determined solely by internal DNA sequence occupied by nucleosome. On the other hand, the action of this internal positioning signal does not extend to neighboring DNA sequences.     

Histone acetylation status and DNA sequence modulate ATP-dependent nucleosome repositioning

A cell-free system derived from Drosophila embryos was used to investigate positioning of nucleosomes on specific DNA sequences. This system can be used to reconstitute differently acetylated nucleosome arrays possessing ATP-dependent dynamic properties that are not observed with chromatin assembled from pure components. Nucleosome positioning on different DNA sequences was studied by restriction endonuclease assay. The sequence of DNA and the acetylation status of histones had profound effects on the distribution of nucleosomes, suggesting their cooperative effect on nucleosome repositioning.     

Structure-based Analysis of DNA Sequence Patterns Guiding Nucleosome Positioning in vitro

Abstract

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken β-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the ‘601’ nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp—an accuracy exceeding that of earlier predictions.

Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning ‘rules’ they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL ±3.5 and ±5.5, which is similar to the α-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the ‘601’ clone) the same YR/YYRR motifs occur predominantly at the sites SHL ±1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures.

Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.     

Sequence Signatures of Nucleosome Positioning in Caenorhabditis elegans

Our recent investigation in the protist Trichomonas vaginalis suggested a DNA sequence periodicity with a unit length of 120.9 nt, which represents a sequence signature for nucleosome positioning. We now extended our observation in higher eukaryotes and identified a similar periodicity of 175 nt in length in Caenorhabditis elegans. In the process of defining the sequence compositional characteristics, we found that the 10.5-nt periodicity, the sequence signature of DNA double helix, may not be sufficient for cross-nucleosome positioning but provides essential guiding rails to facilitate positioning. We further dissected nucleosome-protected sequences and identified a strong positive purine (AG) gradient from the 5′-end to the 3′-end, and also learnt that the nucleosome-enriched regions are GC-rich as compared to the nucleosome-free sequences as purine content is positively correlated with GC content. Sequence characterization allowed us to develop a hidden Markov model (HMM) algorithm for decoding nucleosome positioning computationally, and based on a set of training data from the fifth chromosome of C. elegans, our algorithm predicted 60%-70% of the well-positioned nucleosomes, which is 15%-20% higher than random positioning. We concluded that nucleosomes are not randomly positioned on DNA sequences and yet bind to different genome regions with variable stability, well-positioned nucleosomes leave sequence signatures on DNA, and statistical positioning of nucleosomes across genome can be decoded computationally based on these sequence signatures.