Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (∼10.3 bp) at the two ends, with larger (∼11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames’ start and end sites were occupied by nucleosomes.     

The mechanics behind DNA sequence-dependent properties of the nucleosome

Chromatin organization and composition impart sophisticated regulatory features critical to eukaryotic genomic function. Although DNA sequence-dependent histone octamer binding is important for nucleosome activity, many aspects of this phenomenon have remained elusive. We studied nucleosome structure and stability with diverse DNA sequences, including Widom 601 derivatives with the highest known octamer affinities, to establish a simple model behind the mechanics of sequence dependency. This uncovers the unique but unexpected role of TA dinucleotides and a propensity for G|C-rich sequence elements to conform energetically favourably at most locations around the histone octamer, which rationalizes G|C% as the most predictive factor for nucleosome occupancy in vivo. In addition, our findings reveal dominant constraints on double helix conformation by H3-H4 relative to H2A-H2B binding and DNA sequence context-dependency underlying nucleosome structure, positioning and stability. This provides a basis for improved prediction of nucleosomal properties and the design of tailored DNA constructs for chromatin investigations.     

DNA methylation, nucleosome formation and positioning.

Recent mapping of nucleosome positioning on several long gene regions subject to DNA methylation has identified instances of nucleosome repositioning by this base modification. The evidence for an effect of CpG methylation on nucleosome formation and positioning in chromatin is reviewed here in the context of the complex sequence-structure requirements of DNA wrapping around the histone octamer and the role of this epigenetic mark in gene repression.     

Effects of DNA sequence and conformation on nucleosome formation

A simple theoretical analysis of the free energy balance controlling nucleosome formation shows that the specific effects of different DNA sequences and/or conformations observed in vitro are mainly due to their different elastic properties. A calculation of the elastic free energy required to fold DNA on histone octamers yields quantitative results rationalizing the experimental findings provided that: (i) the average helical repeat of DNA on nucleosomes is greater than 10.2 bp per turn, and (ii) poly[dG.dC] adopts an A-type conformation.     

The influence of DNA stiffness upon nucleosome formation

The rotational and translational positioning of nucleosomes on DNA is dependent to a significant extent on the physicochemical properties of the double helix. We have investigated the influence of the axial flexibility of the molecule on the affinity for the histone octamer by substituting selected DNA sequences with either inosine for guanosine or diaminopurine for adenine. These substitutions, respectively, remove or add a purine 2-amino group exposed in the minor groove and, respectively, decrease and increase the apparent persistence length. We observe that for all sequences tested inosine substitution, with one exception, increases the affinity for histone binding. Conversely diaminopurine substitution decreases the affinity. In the sole example where replacement of guanosine with inosine decreases the persistence length as well as the affinity for histones, the substitution concomitantly removes an intrinsic curvature of the DNA molecule. We show that, to a first approximation, the binding energy of DNA to histones at 1M NaCl is directly proportional to the persistence length. The data also indicate that a high local flexibility of DNA can favour strong rotational positioning.     

Evolutionary footprints of nucleosome positions in yeast

Using genome-wide maps of nucleosome positions in yeast, we have analyzed the influence of chromatin structure on the molecular evolution of genomic DNA. We have observed, on average, 10-15% lower substitution rates in linker regions than in nucleosomal DNA. This widespread local rate heterogeneity represents an evolutionary footprint of nucleosome positions and reveals that nucleosome organization is a genomic feature conserved over evolutionary timescales.     

Different mechanism for in vitro formation of nucleosome core particles

The interaction of different histone oligomers with nucleosomes has been investigated by using nondenaturing gel electrophoresis. In the presence of 0.2 M NaCl, the addition of the pairs H2A,H2B or H3,H4 or the four core histones to nucleosome core particles produces a decrease in the intensity of the core particle band and the appearance of aggregated material at the top of the gel, indicating that all these histone oligomers are able to associate with nucleosomes. Equivalent results were obtained by using oligonucleosome core particles. Additional electrophoretic results, together with second-dimension analysis of histone composition and fluorescence and solubility studies, indicate that H2A,H2B, H3,H4, and the four core histones can migrate spontaneously from the aggregated nucleosomes containing excess histones to free core DNA. In all cases the estimated yield of histone transfer is very high. Furthermore, the results obtained from electron microscopy, solubility, and supercoiling assays demonstrate the transfer of excess histones from oligonucleosomes to free circular DNA. However, the extent of solubilization obtained in this case is lower than that observed with core DNA as histone acceptor. Our results demonstrate that nucleosome core particles can be formed in 0.2 M NaCl by the following mechanisms: (1) transfer of excess core histones from oligonucleosomes of free DNA, (2) transfer to excess H2A,H2B and H3,H4 associated separately with oligonucleosomes to free DNA, (3) transfer to excess H2A,H2B initially associated with oligonucleosomes to DNA, followed by the reaction of the resulting DNA-(H2A,H2B) complex with oligonucleosomes containing excess H3,H4, and (4) a two-step transfer reaction similar to that indicated in (3), in which excess histones H3,H4 are transferred to DNA before the reaction with oligonucleosomes containing excess H2A,H2B. The possible biological implications of these spontaneous reactions are discussed in the context of the present knowledge of the nucleosome function.     

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.     

Noncoding DNA, isochores and gene expression: nucleosome formation potential

The nucleosome formation potential of introns, intergenic spacers and exons of human genes is shown here to negatively correlate with among-tissues breadth of gene expression. The nucleosome formation potential is also found to negatively correlate with the GC content of genomic sequences; the slope of regression line is steeper in exons compared with noncoding DNA (introns and intergenic spacers). The correlation with GC content is independent of sequence length; in turn, the nucleosome formation potential of introns and intergenic spacers positively (albeit weakly) correlates with sequence length independently of GC content. These findings help explain the functional significance of the isochores (regions differing in GC content) in the human genome as a result of optimization of genomic structure for epigenetic complexity and support the notion that noncoding DNA is important for orderly chromatin condensation and chromatin-mediated suppression of tissue-specific genes.     

DNA sequence plays a major role in determining nucleosome positions in yeast CUP1 chromatin.

The role of DNA sequence in determining nucleosome positions in vivo was investigated by comparing the positions adopted by nucleosomes reconstituted on a yeast plasmid in vitro using purified core histones with those in native chromatin containing the same DNA, described previously. Nucleosomes were reconstituted on a 2.5 kilobase pair DNA sequence containing the yeast TRP1ARS1 plasmid with CUP1 as an insert (TAC-DNA). Multiple, alternative, overlapping nucleosome positions were mapped on TAC-DNA. For the 58 positioned nucleosomes identified, the relative positioning strengths and the stabilities to salt and temperature were determined. These positions were, with a few exceptions, identical to those observed in native, remodeled TAC chromatin containing an activated CUP1 gene. Only some of these positions are utilized in native, unremodeled chromatin. These observations suggest that DNA sequence is likely to play a very important role in positioning nucleosomes in vivo. We suggest that events occurring in yeast CUP1 chromatin determine which positions are occupied in vivo and when they are occupied.     

DNA bending and its relation to nucleosome positioning

X-ray and solution studies have shown that the conformation of a DNA double helix depends strongly on its base sequence. Here we show that certain sequence-dependent modulations in structure appear to determine the rotational positioning of DNA about the nucleosome. Three different experiments are described. First, a piece of DNA of defined sequence (169 base-pairs long) is closed into a circle, and its structure examined by digestion with DNAase I: the helix adopts a highly preferred configuration, with short runs of (A, T) facing in and runs of (G, C) facing out. Secondly, the same sequence is reconstituted with a histone octamer: the angular orientation around the histone core remains conserved, apart from a small uniform increase in helix twist. Finally, it is shown that the average sequence content of DNA molecules isolated from chicken nucleosome cores is non-random, as in a reconstituted nucleosome: short runs of (A, T) are preferentially positioned with minor grooves facing in, while runs of (G, C) tend to have their minor grooves facing out. The periodicity of this modulation in sequence content (10.17 base-pairs) corresponds to the helix twist in a local frame of reference (a result that bears on the change in linking number upon nucleosome formation). The determinants of translational positioning have not been identified, but one possibility is that long runs of homopolymer (dA) X (dT) or (dG) X (dC) will be excluded from the central region of the supercoil on account of their resistance to curvature.     

A 145-base pair DNA sequence that positions itself precisely and asymmetrically on the nucleosome core.

A 145-bp DNA sequence, cloned from Escherichia coli, was reconstituted into nucleosome core particles by a number of methods. The behaviour of the resulting complex upon sucrose gradient sedimentation and nucleoprotein gel electrophoresis closely resembled that of control bulk nucleosome core particles. DNase I digestion of the 32P-end-labelled complex revealed the 10-bp periodicity of cleavages expected for DNA bound on a histone surface. The narrow cleavage sites observed (1 bp wide) imply that the sequence occupies a single preferred position on the nucleosome core, accurate to the level of single base pairs. By relating the digestion pattern observed to the pattern of site protection found for random sequence nucleosomes, the DNA position was found to be offset by 17 bp from that in the normal core particle. A number of experiments argue against the involvement of length or end effects and suggest that it is some feature of the DNA sequence itself that determines this precise positioning of DNA on the nucleosome.     

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.     

DNA sequence alterations affect nucleosome array formation of the chicken ovalbumin gene

The role of the large amount (more than half of the genome) of noncoding DNA in higher organisms is not well understood. DNA evolved to function in the context of chromatin, and the possibility exists that some of the noncoding DNA serves to influence chromatin structure and function. In this age of genomics and bioinformatics, genomic DNA sequences are being searched for informational content beyond the known genetic code. The discovery that period-10 non-T, A/T, G (VWG) triplets are among the most abundant motifs in human genomic DNA suggests that they may serve some function in higher organisms. In this paper, we provide direct evidence that the regular oscillation of period-10 VWG that occurs in the chicken ovalbumin gene sequence with a dinucleosome-like period facilitates nucleosome array formation. Using a linker histone-dependent in vitro chromatin assembly system that spontaneously aligns nucleosomes into a physiological array, we show that nucleosomes tend to avoid DNA regions with low period-10 VWG counts. This avoidance leads to the formation of an array with a nucleosome repeat equal to half the period value of the oscillation in period-10 VWG, as determined by Fourier analysis. Two different half-period deletions in the wild-type DNA sequence altered the nucleosome array, as predicted computationally. In contrast, a full-period deletion had an insignificant effect on the nucleosome array formed, also consistent with the prediction. An inversion mutation, with no DNA sequences deleted, again altered the nucleosome array formed, as predicted computationally. Hence, a VWG dinucleosome signal is plausible.     

DNA accessibility to minor groove ligands in core nucleosome and chromatosome

Using the circular dichroism spectra induced in the visible by the binding to the minor groove of DNA, we found that Hoechst 33258 is able to occupy its specific sites even when these are located inside the nucleosome structure. This high accessibility of the DNA in the nucleosome is not modified by the removal of the amino-terminal domains of the octamer histones and is not reduced by the presence of linker histone. Interesting and reasonable differences were found in the association constants.