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.     

Human centromere protein B induces translational positioning of nucleosomes on alpha-satellite sequences

The human centromere proteins A (CENP-A) and B (CENP-B) are the fundamental centromere components of chromosomes. CENP-A is the centromere-specific histone H3 variant, and CENP-B specifically binds a 17-base pair sequence (the CENP-B box), which appears within every other alpha-satellite DNA repeat. In the present study, we demonstrated centromere-specific nucleosome formation in vitro with recombinant proteins, including histones H2A, H2B, H4, CENP-A, and the DNA-binding domain of CENP-B. The CENP-A nucleosome wraps 147 base pairs of the alpha-satellite sequence within its nucleosome core particle, like the canonical H3 nucleosome. Surprisingly, CENP-B binds to nucleosomal DNA when the CENP-B box is wrapped within the nucleosome core particle and induces translational positioning of the nucleosome without affecting its rotational setting. This CENP-B-induced translational positioning only occurs when the CENP-B box sequence is settled in the proper rotational setting with respect to the histone octamer surface. Therefore, CENP-B may be a determinant for translational positioning of the centromere-specific nucleosomes through its binding to the nucleosomal CENP-B box.     

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.     

Spontaneous access to DNA target sites in folded chromatin fibers

DNA wrapped in nucleosomes is sterically occluded from many protein complexes that must act on it; how such complexes gain access to nucleosomal DNA is not known. In vitro studies on isolated nucleosomes show that they undergo spontaneous partial unwrapping conformational transitions, which make the wrapped nucleosomal DNA transiently accessible. Thus, site exposure might provide a general mechanism allowing access of protein complexes to nucleosomal DNA. However, existing quantitative analyses of site exposure focused on single nucleosomes, while the presence of neighbor nucleosomes and concomitant chromatin folding might significantly influence site exposure. In this work, we carried out quantitative studies on the accessibility of nucleosomal DNA in homogeneous nucleosome arrays. Two striking findings emerged. Organization into chromatin fibers changes the accessibility of nucleosomal DNA only modestly, from ∼ 3-fold decreases to ∼ 8-fold increases in accessibility. This means that nucleosome arrays are intrinsically dynamic and accessible even when they are visibly condensed. In contrast, chromatin folding decreases the accessibility of linker DNA by as much as ∼ 50-fold. Thus, nucleosome positioning dramatically influences the accessibility of target sites located inside nucleosomes, while chromatin folding dramatically regulates access to target sites in linker DNA.     

Heterogeneity of chromatin subunits in vitro and location of histone H1.

Chromatin subunits ("nucleosomes") which were purified by sucrose gradient centrifugation of a staphylococcal nuclease digest of chromatin have been studied. We found that such a preparation contains nucleosomes of two discrete types which can be separated from each other by polyacrylamide gel electrophoresis. Nucleosome of the first type contains all five histones and a DNA segment of approximately 200 base pairs long, whereas nucleosome of the second type lacks histone H1 and its DNA segment is approximately 170 base pairs long, i.e., about 30 base pairs shorter than the DNA segment of the nucleosome of the first type. Purified dimer of the nucleosome also can be fractionated by gel electrophoresis into three discrete bands which correspond to dinucleosomes containing two molecules of histone H1, one and no H1. These and related findings strongly suggest that the H1 molecule is bound to a short (approximately 30 base pairs) terminal stretch of the nucleosomal DNA segment which can be removed by nuclease (possibly in the form of H1-DNA complex) without any significant disturbance of main structural features of the nucleosome.     

HMG-D and histone H1 alter the local accessibility of nucleosomal DNA

There is evidence that HMGB proteins facilitate, while linker histones inhibit chromatin remodelling, respectively. We have examined the effects of HMG-D and histone H1/H5 on accessibility of nucleosomal DNA. Using the 601.2 nucleosome positioning sequence designed by Widom and colleagues we assembled nucleosomes in vitro and probed DNA accessibility with restriction enzymes in the presence or absence of HMG-D and histone H1/H5. For HMG-D our results show increased digestion at two spatially adjacent sites, the dyad and one terminus of nucleosomal DNA. Elsewhere varying degrees of protection from digestion were observed. The C-terminal acidic tail of HMG-D is essential for this pattern of accessibility. Neither the HMG domain by itself nor in combination with the adjacent basic region is sufficient. Histone H1/H5 binding produces two sites of increased digestion on opposite faces of the nucleosome and decreased digestion at all other sites. Our results provide the first evidence of local changes in the accessibility of nucleosomal DNA upon separate interaction with two linker binding proteins.     

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.     

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.     

Nucleosomes Shape DNA Polymorphism and Divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.     

Fragmentation of chromatin by endogenous nucleases, accumulation of the subnucleosomal fragments with high electrophoretic mobility]

Digestion of chromatin by endogenous nucleases to nucleosomes (140-160 base pairs of DNA) is accompanied by the accumulation of subnucleosomal DNP particles with high electrophoretic mobility (20-40 base pairs of DNA). All histones associate with the 140-160 base pairs fragment. The production of subnucleosomal DNP particles does not correlate with the degradation of histone H1 and the appearance of nucleosomes lacking histone H1. Degradation of the protein in this fragment is accompanied by the appearance of free DNA. The data obtained are in agreement with the hypothesis on the origin of subnucleosomes from the nucleosomal locus preferentially associated with the non-histone proteins and on the autonomy of these loci and of the loci associated with histone H1 in the nucleosome.     

Using DNA mechanics to predict in vitro nucleosome positions and formation energies

In eukaryotic genomes, nucleosomes function to compact DNA and to regulate access to it both by simple physical occlusion and by providing the substrate for numerous covalent epigenetic tags. While competition with other DNA-binding factors and action of chromatin remodeling enzymes significantly affect nucleosome formation in vivo, nucleosome positions in vitro are determined by steric exclusion and sequence alone. We have developed a biophysical model, DNABEND, for the sequence dependence of DNA bending energies, and validated it against a collection of in vitro free energies of nucleosome formation and a set of in vitro nucleosome positions mapped at high resolution. We have also made a first ab initio prediction of nucleosomal DNA geometries, and checked its accuracy against the nucleosome crystal structure. We have used DNABEND to design both strong and weak histone- binding sequences, and measured the corresponding free energies of nucleosome formation. We find that DNABEND can successfully predict in vitro nucleosome positions and free energies, providing a physical explanation for the intrinsic sequence dependence of histone–DNA interactions.     

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.     

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.     

Characterization of Dnmt1 Binding and DNA Methylation on Nucleosomes and Nucleosomal Arrays

The packaging of DNA into nucleosomes and the organisation into higher order structures of chromatin limits the access of sequence specific DNA binding factors to DNA. In cells, DNA methylation is preferentially occuring in the linker region of nucleosomes, suggesting a structural impact of chromatin on DNA methylation. These observations raise the question whether DNA methyltransferases are capable to recognize the nucleosomal substrates and to modify the packaged DNA. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the maintenance DNA methyltransferase Dnmt1. Our binding studies show that Dnmt1 has a DNA length sensing activity, binding cooperatively to DNA, and requiring a minimal DNA length of 20 bp. Dnmt1 needs linker DNA to bind to nucleosomes and most efficiently recognizes nucleosomes with symmetric DNA linkers. Footprinting experiments reveal that Dnmt1 binds to both DNA linkers exiting the nucleosome core. The binding pattern correlates with the efficient methylation of DNA linkers. However, the enzyme lacks the ability to methylate nucleosomal CpG sites on mononucleosomes and nucleosomal arrays, unless chromatin remodeling enzymes create a dynamic chromatin state. In addition, our results show that Dnmt1 functionally interacts with specific chromatin remodeling enzymes to enable complete methylation of hemi-methylated DNA in chromatin.     

A statistical thermodynamic model applied to experimental AFM population and location data is able to quantify DNA-histone binding strength and internucleosomal interaction differences between acetylated and unacetylated nucleosomal arrays

Imaging of nucleosomal arrays by atomic force microscopy allows a determination of the exact statistical distributions for the numbers of nucleosomes per array and the locations of nucleosomes on the arrays. This precision makes such data an excellent reference for testing models of nucleosome occupation on multisite DNA templates. The approach presented here uses a simple statistical thermodynamic model to calculate theoretical population and positional distributions and compares them to experimental distributions previously determined for 5S rDNA nucleosomal arrays (208-12,172-12). The model considers the possible locations of nucleosomes on the template, and takes as principal parameters an average free energy of interaction between histone octamers and DNA, and an average wrapping length of DNA around the octamers. Analysis of positional statistics shows that it is possible to consider interactions between nucleosomes and positioning effects as perturbations on a random positioning noninteracting model. Analysis of the population statistics is used to determine histone-DNA association constants and to test for differences in the free energies of nucleosome formation with different types of histone octamers, namely acetylated or unacetylated, and different DNA templates, namely 172-12 or 208-12 5S rDNA multisite templates. The results show that the two template DNAs bind histones with similar affinities but histone acetylation weakens the association of histones with both templates. Analysis of locational statistics is used to determine the strength of specific nucleosome positioning tendencies by the DNA templates, and the strength of the interactions between neighboring nucleosomes. The results show only weak positioning tendencies and that unacetylated nucleosomes interact much more strongly with one another than acetylated nucleosomes; in fact acetylation appears to induce a small anticooperative occupation effect between neighboring nucleosomes.     

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.