Improved alignment of nucleosome DNA sequences using a mixture model

DNA sequences that are present in nucleosomes have a preferential approximately 10 bp periodicity of certain dinucleotide signals, but the overall sequence similarity of the nucleosomal DNA is weak, and traditional multiple sequence alignment tools fail to yield meaningful alignments. We develop a mixture model that characterizes the known dinucleotide periodicity probabilistically to improve the alignment of nucleosomal DNAs. We assume that a periodic dinucleotide signal of any type emits according to a probability distribution around a series of 'hot spots' that are equally spaced along nucleosomal DNA with 10 bp period, but with a 1 bp phase shift across the middle of the nucleosome. We model the three statistically most significant dinucleotide signals, AA/TT, GC and TA, simultaneously, while allowing phase shifts between the signals. The alignment is obtained by maximizing the likelihood of both Watson and Crick strands simultaneously. The resulting alignment of 177 chicken nucleosomal DNA sequences revealed that all 10 distinct dinucleotides are periodic, however, with only two distinct phases and varying intensity. By Fourier analysis, we show that our new alignment has enhanced periodicity and sequence identity compared with center alignment. The significance of the nucleosomal DNA sequence alignment is evaluated by comparing it with that obtained using the same model on non-nucleosomal sequences.     

Correlation between dinucleotide periodicities and nucleosome positioning on mouse satellite DNA

Previous studies demonstrated 16 well-defined nucleosome locations (A-P) on a tandemly repeated prototype 234 base pair (bp) mouse satellite repeat unit. We have aligned the A-P fragments to search for DNA sequence elements that might contribute to nucleosome placement at these positions. Our results demonstrate a strikingly regular, uninterrupted, periodic pattern for the AA dinucleotide occurrences along the entire length of the aligned fragments. The periodicity of the AA occurrences is about 9.7 bp. The pattern exhibits a local minimum at position 74, near the nucleosome dyad axis of symmetry. Other dinucleotides--including AC: GT, CA: TG, and CC: GG--are also placed periodically, but their patterns of occurrence are less regular and less frequent than AA. The calculated spacings between consecutive preferred nucleosome locations on mouse satellite DNA are nearly identical, corresponding to multiples of 9.7 bp. The correlation between the periodicity of dinucleotide occurrences and the average spacing of nucleosome positions suggests that the preferred nucleosome locations recur at intervals that may correspond to the DNA helical repeat in the mouse satellite nucleosomes, and that the histone octamers sample (or slip along) the duplex in steps of 9.7 bp during nucleosome formation on mouse satellite DNA.     

A translational signature for nucleosome positioning in vivo

In vivo nucleosomes often occupy well-defined preferred positions on genomic DNA. An important question is to what extent these preferred positions are directly encoded by the DNA sequence itself. We derive here from in vivo positions, accurately mapped by partial micrococcal nuclease digestion, a translational positioning signal that identifies the approximate midpoint of DNA bound by a histone octamer. This midpoint is, on average, highly A/T rich (∼73%) and, in particular, the dinucleotide TpA occurs preferentially at this and other outward-facing minor grooves. We conclude that in this set of sequences the sequence code for DNA bending and nucleosome positioning differs from the other described sets and we suggest that the enrichment of AT-containing dinucleotides at the centre is required for local untwisting. We show that this signature is preferentially associated with nucleosomes flanking promoter regions and suggest that it contributes to the establishment of gene-specific nucleosome arrays.     

Universal full-length nucleosome mapping sequence probe

For the computational sequence-directed mapping of the nucleosomes, the knowledge of the nucleosome positioning motifs – 10–11 base long sequences – and respective matrices of bendability, is not sufficient, since there is no justified way to fuse these motifs in one continuous nucleosome DNA sequence. Discovery of the strong nucleosome (SN) DNA sequences, with visible sequence periodicity allows derivation of the full-length nucleosome DNA bendability pattern as matrix or consensus sequence. The SN sequences of three species (A. thaliana, C. elegans, and H. sapiens) are aligned (512 sequences for each species), and long (115 dinucleotides) matrices of bendability derived for the species. The matrices have strong common property – alternation of runs of purine–purine (RR) and pyrimidine–pyrimidine (YY) dinucleotides, with average period 10.4 bases. On this basis the universal [R,Y] consensus of the nucleosome DNA sequence is derived, with exactly defined positions of respective penta- and hexamers RRRRR, RRRRRR, YYYYY, and YYYYYY.     

DNase I footprinting of the nucleosome in whole nuclei

DNase I was used to footprint the 147 bp DNA fragment of the nucleosome in whole chicken erythrocyte nuclei. It was found that the higher-order structure imposes an additional protection on nucleosomes at sites close to the entry and exit points of the linker DNA, around the dyad axis (site S 0). The observed protection is extended up to 20 bp on either side of S 0. It is partial (∼50%) and most probably reflects a full protection of different regions in alternatively oriented nucleosomes. These are the same regions which interact with linker histones. The results strongly support the findings by simulation of DNase I digests of unlabelled oligonucleosome fragments in the 30 nm fibre that in all nucleosomes sites S −5 to S −3 and S +3 to S +5 ara on the outside of the fibre exposed to DNase I.     

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.     

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.     

Prediction of nucleosome positions in the yeast genome based on matched mirror position filtering

Nucleosome positioning can affect the accessibility of the underlying DNA to the nuclear environment and as such plays an essential role in the regulation of cellular processes. Specific patterns have been found in the underlying DNA sequences of the nucleosome, and one of the most important patterns includes dinucleotides distributed every 10 to 11 base pairs. Based on this property, we propose to match each dinucleotide in the sequence against its mirror occurrences for 10 to 11 base pairs on both left-hand and right­hand sides. A large number of matches in a local region will then signify the existence of a nucleosome. In this paper, we propose the matched mirror position filters for efficient matching of periodic dinucleotide patterns and computationally predict the nucleosome positions. Experimental results on the Saccharomyces cerevisiae (yeast) genome show that the proposed algorithm can predict nucleosome positions effectively. More than 50% of our predicted nucleosomes are within 35 base pairs of those detected by biological experiments.     

ATP dependent histone phosphorylation and nucleosome assembly in a human cell free extract.

Physiologically spaced nucleosome formation in HeLa cell extracts is ATP dependent. ATP hydrolysis is required for chromatin assembly on both linear and covalently closed circular DNA. The link between the phosphorylation state of histones and nucleosome formation has been examined and we demonstrate that in the absence of histone phosphorylation no stable and regularly spaced nucleosomes are formed. Phosphorylated H3 stabilizes the nucleosome core; while phosphorylation of histone H2a is necessary to increase the linker length between nucleosomes from 0 to approximately 45 bp. Histone H1 alone, whether phosphorylated or unphosphorylated, does not increase the nucleosome repeat length in the absence of core histone phosphorylation. Phosphorylations of H1 and H3 correlate with condensation of chromatin. Maximum ATP hydrolysis which is necessary to increase the periodicity of nucleosomes from approximately 150 to approximately 185 bp, not only inhibits H1 and H3 phosphorylation but facilitates their dephosphorylation.     

Sequence structure of human nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.     

The histone core exerts a dominant constraint on the structure of DNA in a nucleosome.

We have examined the structures of unique sequence, A/T-rich DNAs that are predicted to be relatively rigid [oligo(dA).oligo(dT)], flexible [oligo[d(A-T)]], and curved, using the hydroxyl radical as a cleavage reagent. A 50-base-pair segment containing each of these distinct DNA sequences was placed adjacent to the T7 RNA polymerase promoter, a sequence that will strongly position nucleosomes. The final length of the DNA fragments was 142 bp, enough DNA to assemble a single nucleosome. Cleavage of DNA in solution, while bound to a calcium phosphate crystal, and after incorporation into a nucleosome is examined. We find that the distinct A/T-rich DNAs have very different structural features in solution and helical periodicities when bound to a calcium phosphate. In contrast, the organization of the different DNA sequences when associated with a histone octamer is very similar. We conclude that the histone core exerts a dominant constraint on the structure of DNA in a nucleosome and that inclusion of these various unique sequences has only a very small effect on overall nucleosome stability and structure.     

Dinucleosome DNA of human K562 cells: experimental and computational characterizations

Dinucleosome formation is the first step in the organization of the higher order chromatin structure. With the ultimate aim of elucidating the dinucleosome structure, we constructed a library of human dinucleosome DNA. The library consists of PCR-amplifiable DNA fragments obtained by treatment of nuclei of erythroid K562 cells with micrococcal nuclease followed by extraction of DNA and adaptor ligation to the blunt-ended DNA fragments. The library was then cloned using a plasmid vector and the sequences of the clones were determined. The dominating clones containing the Alu elements were removed. A total of 1002 clones, which comprised a dinucleosome database, contained 84 and 918 clones from the clones before and after removing Alu elements, respectively. Approximately 70% of the clones were between 300 and 400 bp in size and they were distributed to various locations of all chromosomes except the Y chromosome. The clones containing A(2)N(8)A(2)N(8)A(2) or T(2)N(8)T(2)N(8)T(2) sequences were classified into three types, Type I (N shape), Type II (V shape) and Type III (M shape) according to DNA curvature plots. The locations of experimentally determined curved DNA segments matched well with the calculated ones though the clones of Types I and III showed additional curved DNA segments as revealed by the curvature plots. The distributions of complementary dinucleotides in the nucleosome DNA, at the ends of the dinucleosome DNA clones, allowed us to predict the positions of the nucleosome dyad axis, and estimate the size of the nucleosome core DNA, 125nt. The distributions of AA and TT dinucleotides, as well as other RR and YY dinucleotides, showed a periodicity with an average period of 10.4 bases, close to the values observed before. Mapping of nucleosome positions in the dinucleosome database based on the observed periodicity revealed that the nucleosomes were separated by a linker of 7.5+ approximately 10 x n nt. This indicates that the nucleosome-nucleosome orientations are, typically, halfway between parallel and antiparallel. Also an important finding is that the distributions of AA/TT and other RR/YY dinucleotides, apparently, reflect both DNA curvature and DNA bendability, cooperatively contributing to the nucleosome formation.     

Human nucleosomes: special role of CG dinucleotides and Alu-nucleosomes

Background

The periodical occurrence of dinucleotides with a period of 10.4 bases now is undeniably a hallmark of nucleosome positioning. Whereas many eukaryotic genomes contain visible and even strong signals for periodic distribution of dinucleotides, the human genome is rather featureless in this respect. The exact sequence features in the human genome that govern the nucleosome positioning remain largely unknown.

Results

When analyzing the human genome sequence with the positional autocorrelation method, we found that only the dinucleotide CG shows the 10.4 base periodicity, which is indicative of the presence of nucleosomes. There is a high occurrence of CG dinucleotides that are either 31 (10.4 × 3) or 62 (10.4 × 6) base pairs apart from one another - a sequence bias known to be characteristic of Alu-sequences. In a similar analysis with repetitive sequences removed, peaks of repeating CG motifs can be seen at positions 10, 21 and 31, the nearest integers of multiples of 10.4.

Conclusions

Although the CG dinucleotides are dominant, other elements of the standard nucleosome positioning pattern are present in the human genome as well. The positional autocorrelation analysis of the human genome demonstrates that the CG dinucleotide is, indeed, one visible element of the human nucleosome positioning pattern, which appears both in Alu sequences and in sequences without repeats. The dominant role that CG dinucleotides play in organizing human chromatin is to indicate the involvement of human nucleosomes in tuning the regulation of gene expression and chromatin structure, which is very likely due to cytosine-methylation/-demethylation in CG dinucleotides contained in the human nucleosomes. This is further confirmed by the positions of CG-periodical nucleosomes on Alu sequences. Alu repeats appear as monomers, dimers and trimers, harboring two to six nucleosomes in a run. Considering the exceptional role CG dinucleotides play in the nucleosome positioning, we hypothesize that Alu-nucleosomes, especially, those that form tightly positioned runs, could serve as "anchors" in organizing the chromatin in human cells.     

A signal encoded in vertebrate DNA that influences nucleosome positioning and alignment.

Evidence is provided that the nucleotide triplet con-sensus non-T(A/T)G (abbreviated to VWG) influences nucleosome positioning and nucleosome alignment into regular arrays. This triplet consensus has been recently found to exhibit a fairly strong 10 bp periodicity in human DNA, implicating it in anisotropic DNA bendability. It is demonstrated that the experimentally determined preferences for nucleosome positioning in native SV40 chromatin can, to a large extent, be pre-dicted simply by counting the occurrences of the period-10 VWG consensus. Nucleosomes tend to form in regions of the SV40 genome that contain high counts of period-10 VWG and/or avoid regions with low counts. In contrast, periodic occurrences of the dinucleotides AA/TT, implicated in the rotational positioning of DNA in nucleosomes, did not correlate with the preferred nucleosome locations in SV40 chromatin. Periodic occurrences of AA did correlate with preferred nucleosome locations in a region of SV40 DNA where VWG occurrences are low. Regular oscillations in period-10 VWG counts with a dinucleosome period were found in vertebrate DNA regions that aligned nucleosomes into regular arrays in vitro in the presence of linker histone. Escherichia coli and plasmid DNA, which fail to align nucleosomes in vitro, lacked these regular VWG oscillations.     

Re-cracking the nucleosome positioning code

Nucleosomes, the fundamental repeating subunits of all eukaryotic chromatin, are responsible for packaging DNA into chromosomes inside the cell nucleus and controlling gene expression. While it has been well established that nucleosomes exhibit higher affinity for select DNA sequences, until recently it was unclear whether such preferences exerted a significant, genome-wide effect on nucleosome positioning in vivo. This question was seemingly and recently resolved in the affirmative: a wide-ranging series of experimental and computational analyses provided extensive evidence that the instructions for wrapping DNA around nucleosomes are contained in the DNA itself. This subsequently labeled second genetic code was based on data-driven, structural, and biophysical considerations. It was subjected to an extensive suite of validation procedures, with one conclusion being that intrinsic, genome-encoded, nucleosome organization explains approximately 50% of in vivo nucleosome positioning. Here, we revisit both the nature of the underlying sequence preferences, and the performance of the proposed code. A series of new analyses, employing spectral envelope (Fourier transform) methods for assessing key sequence periodicities, classification techniques for evaluating predictive performance, and discriminatory motif finding methods for devising alternate models, are applied. The findings from the respective analyses indicate that signature dinucleotide periodicities are absent from the bulk of the high affinity nucleosome-bound sequences, and that the predictive performance of the code is modest. We conclude that further exploration of the role of sequence-based preferences in genome-wide nucleosome positioning is warranted. This work offers a methodologic counterpart to a recent, high resolution determination of nucleosome positioning that also questions the accuracy of the proposed code and, further, provides illustrations of techniques useful in assessing sequence periodicity and predictive performance.     

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.     

Yeast nucleosome DNA pattern: deconvolution from genome sequences of S. cerevisiae

Positional correlation analysis for the complete genome of Saccharomyces cerevisiae is performed with the aim to reveal possible chromatin-related sequence features. A strong periodicity with the period 10.4 bases is detected in the distance histograms for the dinucleotides AA and TT, with the characteristic decay distance of approximately 50 base pairs. The oscillations are observed as well in the distributions of other dinucleotides. However, the respective amplitudes are small, consistent with secondary effects, due to dominant periodicity of AA and TT. The observations are in accord with earlier data on the chromatin sequence periodicities and nucleosome DNA sequence patterns. The autocorrelations of AA and TT dinucleotides in yeast include also a counter-phase component. A tentative DNA sequence pattern for the yeast nucleosomes is suggested and verified by comparison of its autocorrelation plots with the respective natural autocorrelations. The nucleosome mapping guided by the pattern is in accord with experimental data on the linker length distribution in yeast.     

The in vitro reconstitution of nucleosome and its binding patterns with HMG1／2 and HMG14／17 proteins

Using atomic force microscopy (AFM), the dynamic process of the in vitro nucleosome reconstitution followed by slow dilution from high salt to low salt was visualized. Data showed that the histone octamers were dissociatedfrom DNA at 1M NaC1. When the salt concentration was slowly reduced to 650 mM and 300 mM, the core histones bound to the naked DNA gradually. Once the salt concentration was reduced to 50 mM the classic “beads-on-a-string“ structure was clearly visualized. Furthermore, using the technique of the in vitro reconstitution of nucleosome,the mono- and di- nucleosomes were assembled in vitro with both HS2core (-10681 to -10970 bp) and NCR2 (-372to -194 bp) DNA sequences in the 5‘flanking sequence of human b-globin gene. Data revealed that HMG 1/2 and HMG 14/17 proteins binding to both DNA sequences are changeable following the assembly and disassembly of nucleosomes. We suggest that the changeable binding patterns of HMG 14/17 and HMG1/2 proteins with these regulatory elements may be critical in the process of nucleosome assembly, recruitment of chromatin-modifying activities, and the regulation of human b-globin gene expression.