Periodicity of dinucleotides in nucleosomes derived from simian virus 40 chromatin

ADAPTIVE MACHINE LEARNING TECHNIQUE FOR PERIODICITY DETECTION IN BIOLOGICAL SEQUENCES

Researchers devoted considerable effort to detect the periodicity in DNA sequences, namely, the DNA segments that wrap the Histone protein. It is anticipated that periodic dinucleotide signals are indicators of certain important spots (like binding regions) within a DNA sequence; they are equally spaced along nucleosomal DNA with ~ 10 base-pair period. Positioned nucleosomes are believed to play an important role in transcriptional regulation and for the organization of chromatin in cell nuclei. In this paper, we describe and apply a dynamic periodicity detection algorithm to discover the periodicity of certain dinucleotides in DNA and Protein sequences. Our algorithm is based on suffix tree as the underlying data structure. The proposed approach is suitable to analyze different kinds of data and can serve different targets. It considers the periodicity of alternative substrings like the three dinucleotides AA/TA/TT, in addition to considering dynamic window to detect the periodicity of certain instances of substrings. We tested the applicability, effectiveness and resilience of the proposed approach to noise as compared to the other existing algorithms described in the literature.

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.

High-resolution mapping of sequence-directed nucleosome positioning on genomic DNA

We have mapped in vitro nucleosome positioning on the sheep beta-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.

Nucleosome Positioning Determinants

A previous report demonstrated that one site in a nucleosome assembled onto a synthetic positioning sequence known as Fragment 67 is hypersensitive to permanganate. The site is required for positioning activity and is located 1.5 turns from the dyad, which is a region of high DNA curvature in the nucleosome. Here, the permanganate sensitivity of the nucleosome positioning Fragment 601 was examined in order to expand the dataset of nucleosome sequences containing KMnO(4) hypersensitive sites. The hyperreactive T residue in the six sites detected as well as the one in Fragment 67 and three in the 5 S rDNA positioning sequence were contained within a TA step. Seven of the ten sequences were of the form CTAGPuG or the related sequence TTAAPu. These motifs were also found in the binding sites of several transcriptional regulatory proteins that kink DNA. In order to assess the significance of these sites, the 10 bp positioning determinant in Fragment 67 was removed and replaced by the nine sequences from the 5 S rDNA and Fragment 601. The results demonstrated that these derivative fragments promoted high nucleosome stability and positioning as compared to a control sequence that contained an AT step in place of the TA step. The importance of the TA step was further tested by making single base-pair substitutions in Fragment 67 and the results revealed that stability and positioning activity followed the order: TA>TG>TT>/=TC approximately GG approximately GA approximately AT. Sequences flanking the TA step were also shown to be critical for nucleosome stability and positioning. Nucleosome positioning was restored to near wild-type levels with (CTG)(3), which can form slipped stranded structures and with one base bulges that kink DNA. The results of this study suggest that local DNA structures are important for positioning and that single base-pair changes at these sites could have profound effects on those genomic functions that depend on ordered nucleosomes.

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.

Wavelet Analysis of DNA Bending Profiles reveals Structural Constraints on the Evolution of Genomic Sequences

Analyses of genomic DNA sequences have shown in previous works that base pairs are correlated at large distances with scale-invariant statistical properties. We show in the present study that these correlations between nucleotides (letters) result in fact from long-range correlations (LRC) between sequence-dependent DNA structural elements (words) involved in the packaging of DNA in chromatin. Using the wavelet transform technique, we perform a comparative analysis of the DNA text and of the corresponding bending profiles generated with curvature tables based on nucleosome positioning data. This exploration through the optics of the so-called `wavelet transform microscope' reveals a characteristic scale of 100-200 bp that separates two regimes of different LRC. We focus here on the existence of LRC in the small-scale regime (≲ 200 bp). Analysis of genomes in the three kingdoms reveals that this regime is specifically associated to the presence of nucleosomes. Indeed, small scale LRC are observed in eukaryotic genomes and to a less extent in archaeal genomes, in contrast with their absence in eubacterial genomes. Similarly, this regime is observed in eukaryotic but not in bacterial viral DNA genomes. There is one exception for genomes of Poxviruses, the only animal DNA viruses that do not replicate in the cell nucleus and do not present small scale LRC. Furthermore, no small scale LRC are detected in the genomes of all examined RNA viruses, with one exception in the case of retroviruses. Altogether, these results strongly suggest that small-scale LRC are a signature of the nucleosomal structure. Finally, we discuss possible interpretations of these small-scale LRC in terms of the mechanisms that govern the positioning, the stability and the dynamics of the nucleosomes along the DNA chain. This paper is maily devoted to a pedagogical presentation of the theoretical concepts and physical methods which are well suited to perform a statistical analysis of genomic sequences. We review the results obtained with the so-called wavelet-based multifractal analysis when investigating the DNA sequences of various organisms in the three kingdoms. Some of these results have been announced in B. Audit et al. [1, 2].

Long-range correlations between DNA bending sites: relation to the structure and dynamics of nucleosomes

It has been established that the precise positioning of nucleosomes on genomic DNA can be achieved, at least for a minority of them, through sequence-dependent processes. However, to what extent DNA sequences play a role in the positioning of the major part of nucleosomes is still debated. The aim of the present study is to examine to what extent long-range correlations (LRC) are related to the presence of nucleosomes. Using the wavelet transform technique, we perform a comparative analysis of the DNA text and of the corresponding bending profiles generated with curvature tables based on nucleosome positioning data. The exploration of a number of eukaryotic and bacterial genomes through the optics of the so-called "wavelet transform microscope" reveals a characteristic scale of 100-200 bp that separates two regimes of different LRC. Here, we focus on the existence of LRC in the small-scale regime (10-200 bp) which are actually observed in eukaryotic genomes, in contrast to their absence in eubacterial genomes. Analysis of viral DNA genomes shows that, like their host's genomes, eukaryotic viruses present LRC but eubacterial viruses do not. There is one exception for genomes of poxviruses (Vaccinia and Melamoplus sanguinipes) which do not replicate in the cell nucleus and do not exhibit LRC. No small-scale LRC are detected in the genomes of all examined RNA viruses, with the exception of retroviruses. These results together with the observation of LRC between particular sequence motifs known to participate in the formation of nucleosomes (e.g. AA dinucleotides) strongly suggest that the 10-200 bp LRC are a signature of the sequence-dependence of nucleosome positioning. Finally, we discuss possible interpretations of these LRC in terms of the physical mechanisms that might govern the positioning and the dynamics of the nucleosomes along the DNA chain through cooperative processes.

Sequence patterns defining the 5' boundary of human genes

In studies of both short and relatively long human genomic DNA, we found a clustering of the consensus site for the transcription factor GCF at the 5' boundary of a subset of human genes. In studies of promoter regions with known transcription initiation site, the cluster of consensus GCF site appeared near the transcription initiation site and in some sequences it extended into the transcribed region defining the leader mRNA. We also found a detectable correlation between the 5' boundary of human genes and recognition motifs for other transcription factors that bind to GC-rich sequences. But in these cases, the correlation was not as general as the correlation observed for the consensus GCF site.

Archaeal histone selection of nucleosome positioning sequences and the procaryotic origin of histone-dependent genome evolution

Archaeal histones and the eucaryal (eucaryotic) nucleosome core histones have almost identical histone folds. Here, we show that DNA molecules selectively incorporated by rHMfB (recombinant archaeal histone B from Methanothermus fervidus) into archaeal nucleosomes from a mixture of approximately 10(14) random sequence molecules contain sequence motifs shown previously to direct eucaryal nucleosome positioning. The dinucleotides GC, AA (=TT) and TA are repeated at approximately 10 bp intervals, with the GC harmonic displaced approximately 5 bp from the AA and TA harmonics [(GCN(3)AA or TA)(n)]. AT and CG were not strongly selected, indicating that TA not equalAT and GC not equalCG in terms of facilitating archaeal nucleosome assembly. The selected molecules have affinities for rHMfB ranging from approximately 9 to 18-fold higher than the level of affinity of the starting population, and direct the positioned assembly of archaeal nucleosomes. Fourier-transform analyses have revealed that AA dinucleotides are much enriched at approximately 10. 1 bp intervals, the helical repeat of DNA wrapped around a nucleosome, in the genomes of Eucarya and the histone-containing Euryarchaeota, but not in the genomes of Bacteria and Crenarchaeota, procaryotes that do not have histones. Facilitating histone packaging of genomic DNA has apparently therefore imposed constraints on genome sequence evolution, and since archaeal histones have no structure in addition to the histone fold, these constraints must result predominantly from histone fold-DNA contacts. Based on the three-domain universal phylogeny, histones and histone-dependent genome sequence evolution most likely evolved after the bacterial-archaeal divergence but before the archaeal-eucaryal divergence, and were subsequently lost in the Crenarchaeota. However, with lateral gene transfer, the first histone fold could alternatively have evolved after the archaeal-eucaryal divergence, early in either the euryarchaeal or eucaryal lineages.

DNA repeats and archaeal nucleosome positioning

Archaeal histones, homologs of the eucaryal nucleosome core histones, have been identified in the Euryarchaeota. They assemble as tetramers with dsDNA to form archaeal nucleosomes that resemble the central structure of the eucaryal nucleosome formed by the histone (H3-H4)2 tetramer. Eucaryal and archaeal nucleosomes assemble preferentially on DNA molecules that best accommodate the severe bends found within these structures, and here we discuss the relationships between archaeal and eucaryal nucleosomes, repeating DNA sequences, and nucleosome positioning.

Sequence motifs and free energies of selected natural and non-natural nucleosome positioning DNA sequences

Our laboratories recently completed SELEX experiments to isolate DNA sequences that most-strongly favor or disfavor nucleosome formation and positioning, from the entire mouse genome or from even more diverse pools of chemically synthetic random sequence DNA. Here we directly compare these selected natural and non-natural sequences. We find that the strongest natural positioning sequences have affinities for histone binding and nucleosome formation that are sixfold or more lower than those possessed by many of the selected non-natural sequences. We conclude that even the highest-affinity sequence regions of eukaryotic genomes are not evolved for the highest affinity or nucleosome positioning power. Fourier transform calculations on the selected natural sequences reveal a special significance for nucleosome positioning of a motif consisting of approximately 10 bp periodic placement of TA dinucleotide steps. Contributions to histone binding and nucleosome formation from periodic TA steps are more significant than those from other periodic steps such as AA (=TT), CC (=GG) and more important than those from the other YR steps (CA (=TG) and CG), which are reported to have greater conformational flexibility in protein-DNA complexes even than TA. We report the development of improved procedures for measuring the free energies of even stronger positioning sequences that may be isolated in the future, and show that when the favorable free energy of histone-DNA interactions becomes sufficiently large, measurements based on the widely used exchange method become unreliable.

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.

A topological approach to nucleosome structure and dynamics: the linking number paradox and other issues

The linking number paradox of DNA in chromatin (two negative crossings around the octamer, associated with a unit linking number reduction), which is 21 years old this year, has come of age. After stirring much debate in the past, the initially hypothetical explanation of the paradox by DNA overtwisting on the nucleosome surface is now presented as a hard fact in recent textbooks. The first part of this article presents a historical perspective of the problem and details the numerous attempts to measure DNA local periodicity, which in one remarkable example sowed the seeds for the discovery of DNA bending. The second part is devoted to the DNA minicircle system, which has been developed in the author's laboratory as an alternative to the local-periodicity-measurement approach. It offers a simple proposal: a unit linking number reduction associated with a single crossing. This conclusion is contrasted with the latest high-resolution crystallographic data of the nucleosome in the third part of the article, and the fourth part examines the available evidence supporting an extension of these results to nucleosomes in chromatin. The last part addresses another basic question pertaining to nucleosome dynamics, the conformational flexibility of the histone tetramer.

New DNA sequence rules for high affinity binding to histone octamer and sequence-directed nucleosome positioning

DNA sequences that position nucleosomes are of increasing interest because of their relationship to gene regulation in vivo and because of their utility in studies of nucleosome structure and function in vitro. However, at present our understanding of the rules for DNA sequence-directed nucleosome positioning is fragmentary, and existing positioning sequences have many limitations. We carried out a SELEX experiment starting with a large pool of chemically synthetic random. DNA molecules to identify those individuals having the highest affinity for histone octamer. A set of highest-affinity molecules were selected, cloned, and sequenced, their affinities (free energies) for histone octamer in nucleosome reconstitution measured, and their ability to position nucleosomes in vitro assessed by native gel electrophoresis. The selected sequences have higher affinity than previously known natural or non-natural sequences, and have a correspondingly strong nucleosome positioning ability. A variety of analyses including Fourier transform, real-space correlation, and direct counting computations were carried out to assess non-random features in the selected sequences. The results reveal sequence rules that were already identified in earlier studies of natural nucleosomal DNA, together with a large set of new rules having even stronger statistical significance. Possible physical origins of the selected molecules' high affinities are discussed. The sequences isolated in this study should prove valuable for studies of chromatin structure and function in vitro and, potentially, for studies in vivo.

Nucleosome packaging and nucleosome positioning of genomic DNA

The goals of this study were to assess the extent to which bulk genomic DNA sequences contribute to their own packaging in nucleosomes and to reveal the relationship between nucleosome packaging and positioning. Using a competitive nucleosome reconstitution assay, we found that at least 95% of bulk DNA sequences have an affinity for histone octamer in nucleosomes that is similar to that of randomly synthesized DNA; they contribute little to their own packaging at the level of individual nucleosomes. An equation was developed that relates the measured free energy to the fractional occupancy of specific nucleosome positions. Evidently, the bulk of eukaryotic genomic DNA is also not evolved or constrained for significant sequence-directed nucleosome positioning at the level of individual nucleosomes. Implications for gene regulation in vivo are discussed.