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.     

DNA Physical Properties and Nucleosome Positions Are Major Determinants of HIV-1 Integrase Selectivity

Retroviral integrases (INs) catalyse the integration of the reverse transcribed viral DNA into the host cell genome. This process is selective, and chromatin has been proposed to be a major factor regulating this step in the viral life cycle. However, the precise underlying mechanisms are still under investigation. We have developed a new in vitro integration assay using physiologically-relevant, reconstituted genomic acceptor chromatin and high-throughput determination of nucleosome positions and integration sites, in parallel. A quantitative analysis of the resulting data reveals a chromatin-dependent redistribution of the integration sites and establishes a link between integration sites and nucleosome positions. The co-activator LEDGF/p75 enhanced integration but did not modify the integration sites under these conditions. We also conducted an in cellulo genome-wide comparative study of nucleosome positions and human immunodeficiency virus type-1 (HIV-1) integration sites identified experimentally in vivo. These studies confirm a preferential integration in nucleosome-covered regions. Using a DNA mechanical energy model, we show that the physical properties of DNA probed by IN binding are important in determining IN selectivity. These novel in vitro and in vivo approaches confirm that IN has a preference for integration into a nucleosome, and suggest the existence of two levels of IN selectivity. The first depends on the physical properties of the target DNA and notably, the energy required to fit DNA into the IN catalytic pocket. The second depends on the DNA deformation associated with DNA wrapping around a nucleosome. Taken together, these results indicate that HIV-1 IN is a shape-readout DNA binding protein.     

Direct measurements of the nucleosome-forming preferences of periodic DNA motifs challenge established models

Several periodic motifs have been implicated in facilitating the bending of DNA around the histone core of the nucleosome. For example, di-nucleotides AA/TT/TA and GC at ∼10-bp periods, but offset by 5 bp, are found with higher-than-expected occurrences in aligned nucleosomal DNAs in vitro and in vivo. Additionally, regularly oscillating period-10 trinucleotide motifs non-T, A/T, G and their complements have been implicated in the formation of regular nucleosome arrays. The effects of these periodic motifs on nucleosome formation have not been systematically tested directly by competitive reconstitution assays. We show that, in general, none of these period-10 motifs, except TA, in certain sequence contexts, facilitates nucleosome formation. The influence of periodic TAs on nucleosome formation is appreciable; with some of the 200-bp DNAs out-competing bulk nucleosomal DNA by more than 400-fold. Only the nucleotides immediately flanking TA influence its nucleosome-forming ability. Period-10 TA, when flanked by a pair of permissive nucleotides, facilitates DNA bending through compression of the minor groove. The free energy change for nucleosome formation decreases linearly with the number of consecutive TAs, up to eight. We suggest how these data can be reconciled with previous findings.     

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.     

Reproducibility and Consistency of In Vitro Nucleosome Reconstitutions Demonstrated by Invitrosome Isolation and Sequencing

Nucleosomes and their positions in the eukaryotic genome play an important role in regulating gene expression by influencing accessibility to DNA. Many factors influence a nucleosome''s final position in the chromatin landscape including the underlying genomic sequence. One of the primary reasons for performing in vitro nucleosome reconstitution experiments is to identify how the underlying DNA sequence will influence a nucleosome''s position in the absence of other compounding cellular factors. However, concerns have been raised about the reproducibility of data generated from these kinds of experiments. Here we present data for in vitro nucleosome reconstitution experiments performed on linear plasmid DNA that demonstrate that, when coverage is deep enough, these reconstitution experiments are exquisitely reproducible and highly consistent. Our data also suggests that a coverage depth of 35X be maintained for maximal confidence when assaying nucleosome positions, but lower coverage levels may be generally sufficient. These coverage depth recommendations are sufficient in the experimental system and conditions used in this study, but may vary depending on the exact parameters used in other systems.     

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.     

Predicting nucleosome positions on the DNA: combining intrinsic sequence preferences and remodeler activities

Nucleosome positions on the DNA are determined by the intrinsic affinities of histone proteins to a given DNA sequence and by the ATP-dependent activities of chromatin remodeling complexes that can translocate nucleosomes with respect to the DNA. Here, we report a theoretical approach that takes into account both contributions. In the theoretical analysis two types of experiments have been considered: in vitro experiments with a single reconstituted nucleosome and in vivo genome-scale mapping of nucleosome positions. The effect of chromatin remodelers was described by iteratively redistributing the nucleosomes according to certain rules until a new steady state was reached. Three major classes of remodeler activities were identified: (i) the establishment of a regular nucleosome spacing in the vicinity of a strong positioning signal acting as a boundary, (ii) the enrichment/depletion of nucleosomes through amplification of intrinsic DNA-sequence-encoded signals and (iii) the removal of nucleosomes from high-affinity binding sites. From an analysis of data for nucleosome positions in resting and activated human CD4⁺ T cells [Schones et al., Cell 132, p. 887] it was concluded that the redistribution of a nucleosome map to a new state is greatly facilitated if the remodeler complex translocates the nucleosome with a preferred directionality.     

Long-range oscillation in a periodic DNA sequence motif may influence nucleosome array formation

We have experimentally examined the characteristics of nucleosome array formation in different regions of mouse liver chromatin, and have computationally analyzed the corresponding genomic DNA sequences. We have shown that the mouse adenosine deaminase (MADA) gene locus is packaged into an exceptionally regular nucleosome array with a shortened repeat, consistent with our computational prediction based on the DNA sequence. A survey of the mouse genome indicates that <10% of 70 kb windows possess a nucleosome-ordering signal, consisting of regular long-range oscillations in the period-10 triplet motif non-T, A/T, G (VWG), which is as strong as the signal in the MADA locus. A strong signal in the center of this locus, confirmed by in vitro chromatin assembly experiments, appears to cooperate with weaker, in-phase signals throughout the locus. In contrast, the mouse odorant receptor (MOR) locus, which lacks locus-wide signals, was representative of ~40% of the mouse genomic DNA surveyed. Within this locus, nucleosome arrays were similar to those of bulk chromatin. Genomic DNA sequences which were computationally similar to MADA or MOR resulted in MADA- or MOR-like nucleosome ladders experimentally. Overall, we provide evidence that computationally predictable information in the DNA sequence may affect nucleosome array formation in animal tissue.     

High-Resolution Mapping of Sequence-Directed Nucleosome Positioning on Genomic DNA

We have mapped in vitro nucleosome positioning on the sheep β-lactoglobulin gene using high-throughput sequencing to characterise the DNA sequences recovered from reconstituted nucleosomes. This methodology surpasses previous approaches for coverage, accuracy and resolution and, most importantly, offers a simple yet rapid and relatively inexpensive method to characterise genomic DNA sequences in terms of nucleosome positioning capacity. We demonstrate an unambiguous correspondence between in vitro and in vivo nucleosome positioning around the promoter of the gene; identify discrete, sequence-specific nucleosomal structures above the level of the canonical core particle—a feature that has implications for regulatory protein access and higher-order chromatin packing; and reveal new insights into the involvement of periodically organised dinucleotide sequence motifs of the type GG and CC and not AA and TT, as determinants of nucleosome positioning—an observation that supports the idea that the core histone octamer can exploit different patterns of sequence organisation, or structural potential, in the DNA to bring about nucleosome positioning.     

The main role of the sequence-dependent DNA elasticity in determining the free energy of nucleosome formation on telomeric DNAs

Using a competitive reconstitution assay, we measured the free energy spent in nucleosome formation of eight telomeric DNAs, differing in sequence and/or in length. The obtained values are in satisfactorily good agreement with those derived from a theoretical model that allows the calculation of the free energy of nucleosome formation on the basis of sequence-dependent DNA elasticity, using a statistical thermodynamic approach. Both theoretical and experimental evaluations show that telomeres are characterized by the highest free energies of nucleosome formation among all the DNA sequences so far studied. The free energy of nucleosome formation varies according to the different telomeric sequences and the length of the fragments. Theoretical analysis and experimental mapping by lambda exonuclease show that telomeric nucleosomes occupy multiple positions spaced every telomeric repeat. Sequence-dependent DNA elasticity appears as the main determinant of the stability of telomeric nucleosomes and their multiple translational positioning.     

Multiple Aspects of ATP-Dependent Nucleosome Translocation by RSC and Mi-2 Are Directed by the Underlying DNA Sequence

Background

Chromosome structure, DNA metabolic processes and cell type identity can all be affected by changing the positions of nucleosomes along chromosomal DNA, a reaction that is catalysed by SNF2-type ATP-driven chromatin remodelers. Recently it was suggested that in vivo, more than 50% of the nucleosome positions can be predicted simply by DNA sequence, especially within promoter regions. This seemingly contrasts with remodeler induced nucleosome mobility. The ability of remodeling enzymes to mobilise nucleosomes over short DNA distances is well documented. However, the nucleosome translocation processivity along DNA remains elusive. Furthermore, it is unknown what determines the initial direction of movement and how new nucleosome positions are adopted.

Methodology/Principal Findings

We have used AFM imaging and high resolution PAGE of mononucleosomes on 600 and 2500 bp DNA molecules to analyze ATP-dependent nucleosome repositioning by native and recombinant SNF2-type enzymes. We report that the underlying DNA sequence can control the initial direction of translocation, translocation distance, as well as the new positions adopted by nucleosomes upon enzymatic mobilization. Within a strong nucleosomal positioning sequence both recombinant Drosophila Mi-2 (CHD-type) and native RSC from yeast (SWI/SNF-type) repositioned the nucleosome at 10 bp intervals, which are intrinsic to the positioning sequence. Furthermore, RSC-catalyzed nucleosome translocation was noticeably more efficient when beyond the influence of this sequence. Interestingly, under limiting ATP conditions RSC preferred to position the nucleosome with 20 bp intervals within the positioning sequence, suggesting that native RSC preferentially translocates nucleosomes with 15 to 25 bp DNA steps.

Conclusions/Significance

Nucleosome repositioning thus appears to be influenced by both remodeler intrinsic and DNA sequence specific properties that interplay to define ATPase-catalyzed repositioning. Here we propose a successive three-step framework consisting of initiation, translocation and release steps to describe SNF2-type enzyme mediated nucleosome translocation along DNA. This conceptual framework helps resolve the apparent paradox between the high abundance of ATP-dependent remodelers per nucleus and the relative success of sequence-based predictions of nucleosome positioning in vivo.     

Crystal Structures of Nucleosome Core Particles Containing the ‘601’ Strong Positioning Sequence

Nucleosome positioning plays a key role in genomic regulation by defining histone-DNA context and by modulating access to specific sites. Moreover, the histone-DNA register influences the double-helix structure, which in turn can affect the association of small molecules and protein factors. Analysis of genomic and synthetic DNA has revealed sequence motifs that direct nucleosome positioning in vitro; thus, establishing the basis for the DNA sequence dependence of positioning would shed light on the mechanics of the double helix and its contribution to chromatin structure in vivo. However, acquisition of well-diffracting nucleosome core particle (NCP) crystals is extremely dependent on the DNA fragment used for assembly, and all previous NCP crystal structures have been based on human α-satellite sequences. Here, we describe the crystal structures of Xenopus NCPs containing one of the strongest known histone octamer binding and positioning sequences, the so-called ‘601’ DNA.Two distinct 145-bp 601 crystal forms display the same histone-DNA register, which coincides with the occurrence of DNA stretching-overtwisting in both halves of the particle around five double-helical turns from the nucleosome center, giving the DNA an ‘effective length’ of 147 bp. As we have found previously with stretching around two turns from the nucleosome center for a centromere-based sequence, the terminal stretching observed in the 601 constructs is associated with extreme kinking into the minor groove at purine-purine (pyrimidine-pyrimidine) dinucleotide steps. In other contexts, these step types display an overall nonflexible behavior, which raises the possibility that DNA stretching in the nucleosome or extreme distortions in general have unique sequence dependency characteristics. Our findings indicate that DNA stretching is an intrinsically predisposed site-specific property of the nucleosome and suggest how NCP crystal structures with diverse DNA sequences can be obtained.     

cgDNA: a software package for the prediction of sequence-dependent coarse-grain free energies of B-form DNA

cgDNA is a package for the prediction of sequence-dependent configuration-space free energies for B-form DNA at the coarse-grain level of rigid bases. For a fragment of any given length and sequence, cgDNA calculates the configuration of the associated free energy minimizer, i.e. the relative positions and orientations of each base, along with a stiffness matrix, which together govern differences in free energies. The model predicts non-local (i.e. beyond base-pair step) sequence dependence of the free energy minimizer. Configurations can be input or output in either the Curves+ definition of the usual helical DNA structural variables, or as a PDB file of coordinates of base atoms. We illustrate the cgDNA package by comparing predictions of free energy minimizers from (a) the cgDNA model, (b) time-averaged atomistic molecular dynamics (or MD) simulations, and (c) NMR or X-ray experimental observation, for (i) the Dickerson–Drew dodecamer and (ii) three oligomers containing A-tracts. The cgDNA predictions are rather close to those of the MD simulations, but many orders of magnitude faster to compute. Both the cgDNA and MD predictions are in reasonable agreement with the available experimental data. Our conclusion is that cgDNA can serve as a highly efficient tool for studying structural variations in B-form DNA over a wide range of sequences.     

Dynamics of nucleosome invasion by DNA binding proteins

Nucleosomes sterically occlude their wrapped DNA from interacting with many large protein complexes. How proteins gain access to nucleosomal DNA target sites in vivo is not known. Outer stretches of nucleosomal DNA spontaneously unwrap and rewrap with high frequency, providing rapid and efficient access to regulatory DNA target sites located there; however, rates for access to the nucleosome interior have not been measured. Here we show that for a selected high-affinity nucleosome positioning sequence, the spontaneous DNA unwrapping rate decreases dramatically with distance inside the nucleosome. The rewrapping rate also decreases, but only slightly. Our results explain the previously known strong position dependence on the equilibrium accessibility of nucleosomal DNA, which is characteristic of both selected and natural sequences. Our results point to slow nucleosome conformational fluctuations as a potential source of cell-cell variability in gene activation dynamics, and they reveal the dominant kinetic path by which multiple DNA binding proteins cooperatively invade a nucleosome.     

Learning a Weighted Sequence Model of the Nucleosome Core and Linker Yields More Accurate Predictions in Saccharomyces cerevisiae and Homo sapiens

DNA in eukaryotes is packaged into a chromatin complex, the most basic element of which is the nucleosome. The precise positioning of the nucleosome cores allows for selective access to the DNA, and the mechanisms that control this positioning are important pieces of the gene expression puzzle. We describe a large-scale nucleosome pattern that jointly characterizes the nucleosome core and the adjacent linkers and is predominantly characterized by long-range oscillations in the mono, di- and tri-nucleotide content of the DNA sequence, and we show that this pattern can be used to predict nucleosome positions in both Homo sapiens and Saccharomyces cerevisiae more accurately than previously published methods. Surprisingly, in both H. sapiens and S. cerevisiae, the most informative individual features are the mono-nucleotide patterns, although the inclusion of di- and tri-nucleotide features results in improved performance. Our approach combines a much longer pattern than has been previously used to predict nucleosome positioning from sequence—301 base pairs, centered at the position to be scored—with a novel discriminative classification approach that selectively weights the contributions from each of the input features. The resulting scores are relatively insensitive to local AT-content and can be used to accurately discriminate putative dyad positions from adjacent linker regions without requiring an additional dynamic programming step and without the attendant edge effects and assumptions about linker length modeling and overall nucleosome density. Our approach produces the best dyad-linker classification results published to date in H. sapiens, and outperforms two recently published models on a large set of S. cerevisiae nucleosome positions. Our results suggest that in both genomes, a comparable and relatively small fraction of nucleosomes are well-positioned and that these positions are predictable based on sequence alone. We believe that the bulk of the remaining nucleosomes follow a statistical positioning model.