Nucleosome DNA Bendability Matrix(C. elegans)

Manipulating chromatin architecture in C. elegans

BACKGROUND

Nucleosome-mediated chromatin compaction has a direct effect on the accessibility of trans-acting activators and repressors to DNA targets and serves as a primary regulatory agent of genetic expression. Understanding the nature and dynamics of chromatin is fundamental to elucidating the mechanisms and factors that epigenetically regulate gene expression. Previous work has shown that there are three types of canonical sequences that strongly regulate nucleosome positioning and thus chromatin accessibility: putative nucleosome-positioning elements, putative nucleosome-repelling sequences, and homopolymeric runs of A/T. It is postulated that these elements can be used to remodel chromatin in C. elegans. Here we show the utility of such elements in vivo, and the extreme efficacy of a newly discovered repelling sequence, PRS-322.

RESULTS

In this work, we show that it is possible to manipulate nucleosome positioning in C. elegans solely using canonical and putative positioning sequences. We have not only tested previously described sequences such as the Widom 601, but also have tested additional nucleosome-positioning sequences: the Trifonov sequence, putative repelling sequence-322 (PRS-322), and various homopolymeric runs of A and T nucleotides.

CONCLUSIONS

Using each of these types of putative nucleosome-positioning sequences, we demonstrate their ability to alter the nucleosome profile in C. elegans as evidenced by altered nucleosome occupancy and positioning in vivo. Additionally, we show the effect that PRS-322 has on nucleosome-repelling and chromatin remodeling.

Predicting the configuration and energy of DNA in a nucleosome by coarse-grain modelling

We present a novel algorithm which uses a coarse-grained model and an energy minimisation procedure to predict the sequence-dependent DNA configuration in a nucleosome together with its energetic cost.

The affinity of DNA sequences containing R5Y5 motif and TA repeats with 10.5-bp periodicity to histone octamer in vitro

Nucleosome positioning along the genome is partially determined by the intrinsic DNA sequence preferences on histone. RRRRRYYYYY (R5Y5, R = Purine and Y = Pyrimidine) motif in nucleosome DNA, which was presented based on several theoretical models by Trifonov et al., might be a facilitating sequence pattern for nucleosome assembly. However, there is not a high conformity experimental evidence to support the concept that R5Y5 motif is a key element for the determination of nucleosome positioning. In this work, the ability of the canonical, H2A.Z- and H3.3-containing octamers to assemble nucleosome on DNA templates containing R5Y5 motif and TA repeats within 10.5-bp periodicity was investigated by using salt-dialysis method in vitro. The results showed that the10.5-bp periodical distributions of both R5Y5 motif and TA repeats along DNA templates can significantly promote canonical nucleosome assembly and may be key sequence factors for canonical nucleosome assembly. Compared with TA repeats within 10.5-bp periodicity, R5Y5 motif in DNA templates did not elevate H2A.Z- and H3.3-containing nucleosome formation efficiency in vitro. This result indicates that R5Y5 motif probably isn't a pivotal factor to regulate nucleosome assembly on histone variants. It is speculated that the regulatory mechanism of nucleosome assembly is different between canonical and variant histone. These conclusions can provide a deeper insight on the mechanism of nucleosome positioning. Communicated by Ramaswamy H. Sarma.

Influence of Rotational Nucleosome Positioning on Transcription Start Site Selection in Animal Promoters

The recruitment of RNA-Pol-II to the transcription start site (TSS) is an important step in gene regulation in all organisms. Core promoter elements (CPE) are conserved sequence motifs that guide Pol-II to the TSS by interacting with specific transcription factors (TFs). However, only a minority of animal promoters contains CPEs. It is still unknown how Pol-II selects the TSS in their absence. Here we present a comparative analysis of promoters' sequence composition and chromatin architecture in five eukaryotic model organisms, which shows the presence of common and unique DNA-encoded features used to organize chromatin. Analysis of Pol-II initiation patterns uncovers that, in the absence of certain CPEs, there is a strong correlation between the spread of initiation and the intensity of the 10 bp periodic signal in the nearest downstream nucleosome. Moreover, promoters' primary and secondary initiation sites show a characteristic 10 bp periodicity in the absence of CPEs. We also show that DNA natural variants in the region immediately downstream the TSS are able to affect both the nucleosome-DNA affinity and Pol-II initiation pattern. These findings support the notion that, in addition to CPEs mediated selection, sequence-induced nucleosome positioning could be a common and conserved mechanism of TSS selection in animals.

Nucleosomal signatures impose nucleosome positioning in coding and noncoding sequences in the genome

In the yeast genome, a large proportion of nucleosomes occupy well-defined and stable positions. While the contribution of chromatin remodelers and DNA binding proteins to maintain this organization is well established, the relevance of the DNA sequence to nucleosome positioning in the genome remains controversial. Through quantitative analysis of nucleosome positioning, we show that sequence changes distort the nucleosomal pattern at the level of individual nucleosomes in three species of Schizosaccharomyces and in Saccharomyces cerevisiae. This effect is equally detected in transcribed and nontranscribed regions, suggesting the existence of sequence elements that contribute to positioning. To identify such elements, we incorporated information from nucleosomal signatures into artificial synthetic DNA molecules and found that they generated regular nucleosomal arrays indistinguishable from those of endogenous sequences. Strikingly, this information is species-specific and can be combined with coding information through the use of synonymous codons such that genes from one species can be engineered to adopt the nucleosomal organization of another. These findings open the possibility of designing coding and noncoding DNA molecules capable of directing their own nucleosomal organization.

Review fifteen years of search for strong nucleosomes

Don Crothers, Mikael Kubista, Jon Widom, and their teams have been first to look for strong nucleosomes, in a bid to reveal the nucleosome positioning pattern(s) carried by the nucleosome DNA sequences. They were first to demonstrate that the nucleosome stability correlates with 10-11 base sequence periodicity, and that the strong nucleosomes localize preferentially in centromeres. This review describes these findings and their connection to recent discovery of the strong nucleosomes (SNs) with visibly periodic nucleosome DNA sequences.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

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.

In Bacillus subtilis LutR is part of the global complex regulatory network governing the adaptation to the transition from exponential growth to stationary phase

The lutR gene, encoding a product resembling a GntR-family transcriptional regulator, has previously been identified as a gene required for the production of the dipeptide antibiotic bacilysin in Bacillus subtilis. To understand the broader regulatory roles of LutR in B. subtilis, we studied the genome-wide effects of a lutR null mutation by combining transcriptional profiling studies using DNA microarrays, reverse transcription quantitative PCR, lacZ fusion analyses and gel mobility shift assays. We report that 65 transcriptional units corresponding to 23 mono-cistronic units and 42 operons show altered expression levels in lutR mutant cells, as compared with lutR + wild-type cells in early stationary phase. Among these, 11 single genes and 25 operons are likely to be under direct control of LutR. The products of these genes are involved in a variety of physiological processes associated with the onset of stationary phase in B. subtilis, including degradative enzyme production, antibiotic production and resistance, carbohydrate utilization and transport, nitrogen metabolism, phosphate uptake, fatty acid and phospholipid biosynthesis, protein synthesis and translocation, cell-wall metabolism, energy production, transfer of mobile genetic elements, induction of phage-related genes, sporulation, delay of sporulation and cannibalism, and biofilm formation. Furthermore, an electrophoretic mobility shift assay performed in the presence of both SinR and LutR revealed a close overlap between the LutR and SinR targets. Our data also revealed a significant overlap with the AbrB regulon. Together, these findings reveal that LutR is part of the global complex, interconnected regulatory systems governing adaptation of bacteria to the transition from exponential growth to stationary phase.

Visible periodicity of strong nucleosome DNA sequences

Fifteen years ago, Lowary and Widom assembled nucleosomes on synthetic random sequence DNA molecules, selected the strongest nucleosomes and discovered that the TA dinucleotides in these strong nucleosome sequences often appear at 10-11 bases from one another or at distances which are multiples of this period. We repeated this experiment computationally, on large ensembles of natural genomic sequences, by selecting the strongest nucleosomes--i.e. those with such distances between like-named dinucleotides, multiples of 10.4 bases, the structural and sequence period of nucleosome DNA. The analysis confirmed the periodicity of TA dinucleotides in the strong nucleosomes, and revealed as well other periodic sequence elements, notably classical AA and TT dinucleotides. The matrices of DNA bendability and their simple linear forms--nucleosome positioning motifs--are calculated from the strong nucleosome DNA sequences. The motifs are in full accord with nucleosome positioning sequences derived earlier, thus confirming that the new technique, indeed, detects strong nucleosomes. Species- and isochore-specific variations of the matrices and of the positioning motifs are demonstrated. The strong nucleosome DNA sequences manifest the highest hitherto nucleosome positioning sequence signals, showing the dinucleotide periodicities in directly observable rather than in hidden form.

Periodic distribution of a putative nucleosome positioning motif in human, nonhuman primates, and archaea: mutual information analysis

Recently, Trifonov's group proposed a 10-mer DNA motif YYYYYRRRRR as a solution of the long-standing problem of sequence-based nucleosome positioning. To test whether this generic decamer represents a biological meaningful signal, we compare the distribution of this motif in primates and Archaea, which are known to contain nucleosomes, and in Eubacteria, which do not possess nucleosomes. The distribution of the motif is analyzed by the mutual information function (MIF) with a shifted version of itself (MIF profile). We found common features in the patterns of this generic decamer on MIF profiles among primate species, and interestingly we found conspicuous but dissimilar MIF profiles for each Archaea tested. The overall MIF profiles for each chromosome in each primate species also follow a similar pattern. Trifonov's generic decamer may be a highly conserved motif for the nucleosome positioning, but we argue that this is not the only motif. The distribution of this generic decamer exhibits previously unidentified periodicities, which are associated to highly repetitive sequences in the genome. Alu repetitive elements contribute to the most fundamental structure of nucleosome positioning in higher Eukaryotes. In some regions of primate chromosomes, the distribution of the decamer shows symmetrical patterns including inverted repeats.

The roles of the monomer length and nucleotide context of plant tandem repeats in nucleosome positioning

Similar to regularly spaced nucleosomes in chromatin, long tandem DNA arrays are composed of regularly alternating monomers that have almost identical primary DNA structures. Such a similarity in the structural organization makes these arrays especially interesting for studying the role of intrinsic DNA preferences in nucleosome positioning. We have studied the nucleosome formation potential of DNA tandem repeat families with different monomer lengths (ML). In total, 165 plant tandem repeat families from the PlantSat database (http://w3lamc.umbr.cas.cz/PlantSat/) were divided into two classes based on the number of nucleosome repeats in one DNA monomer. For predicting nucleosome formation potential, we developed the Phase method, which combines the advantages of multiple bioinformatics models. The Phase method was able to distinguish interfamily differences and intrafamily monomer variation and identify the influence of nucleotide context on nucleosome formation potential. Three main types of nucleosome arrangement in DNA tandem repeat arrays--regular, partially regular (partial), and flexible--were distinguished among a great variety of Phase profiles. The regular type, in which all nucleosomes of the monomer array are positioned in a context-dependent manner, is the most representative type of the class 1 families, with ML equal to or a multiple of the nucleosome repeat length (NRL). In the partially regular type, nucleotide context influences the positioning of only a subset of nucleosomes. The influence of the nucleotide context on nucleosome positioning has the least effect in the flexible type, which contains the greatest number of families (65). The majority of these families belong to class 2 and have nonmultiple ML to NRL ratios.

Chromatin organization, epigenetics and differentiation: an evolutionary perspective

Genome packaging is a universal phenomenon from prokaryotes to higher mammals. Genomic constituents and forces have however, travelled a long evolutionary route. Both DNA and protein elements constitute the genome and also aid in its dynamicity. With the evolution of organisms, these have experienced several structural and functional changes. These evolutionary changes were made to meet the challenging scenario of evolving organisms. This review discusses in detail the evolutionary perspective and functionality gain in the phenomena of genome organization and epigenetics.

Multiple levels of meaning in DNA sequences, and one more

If we define a genetic code as a widespread DNA sequence pattern that carries a message with an impact on biology, then there are multiple genetic codes. Sequences involved in these codes overlap and, thus, both interact with and constrain each other, such as for the triplet code, the intron-splicing code, the code for amphipathic alpha helices, and the chromatin code. Nucleosomes preferentially are located at the ends of exons, thus protecting splice junctions, with the N9 positions of guanines of the GT and AG junctions oriented toward the histones. Analysis of protein-coding sequences reveals numerous traces of tandem repeats, apparently formed by triplet expansion, which in effect is a genome inflation ``code''. Our data are consistent with the hypothesis that expansion of simple tandem repetition of certain aggressive triplets has been a characteristic of life from its emergence. Such expanding triplets appear to be the major factor underlying observed codon usage biases.

Predicting nucleosome positions in yeast: using the absolute frequency

Nucleosome is the basic structure of chromatin in eukaryotic cells, and they form the chromatin fiber interconnected by sections of linker DNA. Nucleosome positioning is of great significance for gene transcription regulation. In this paper, we consider the difference of absolute frequency of nucleotides between the nucleosome forming and nucleosome inhibiting sequences. Based on the 2-mer absolute frequency of nucleotides in genome, a new model is constructed to distinguish nucleosome DNA and linker DNA. When used to predict DNA potential for forming nucleosomes in S. cerevisiae, the model achieved a high accuracy of 96.05%. Thus, the model is very useful for predicting nucleosome positioning.

Apoptotic cleavage of DNA in human lymphocyte chromatin shows high sequence specificity

Apoptotic digestion of human lymphocyte chromatin results in the appearance of large amounts of nucleosome size DNA fragments. Sequencing of these fragments and analysis of the distribution of bases around the apoptotic nucleases' cutting sites revealed a rather strong consensus sequence, not observed earlier. The consensus TAAAgTAcTTTA is characterized by complementary symmetry, resembling prokaryotic restriction sites. This consensus also possesses three TA dinucleotide steps, separated by five bases (corresponding to a half-period of the DNA double helix), suggesting strong bending of the DNA at the cut sites which is perhaps required for cutting.

Methods for identifying higher-order chromatin structure

Eukaryotic genomic DNA is combined with histones, nonhistone proteins, and RNA to form chromatin, which is extensively packaged hierarchically to fit inside a cell's nucleus. The nucleosome-comprising a histone octamer with 147 base pairs of DNA wrapped around it-is the initial level and the repeating unit of chromatin packaging, which electron microscopy first made visible to the human eye as "beads on a string" nearly four decades ago. The mechanism and nature of chromatin packaging are still under intense research. Recently, classic methods like chromatin immunoprecipitation and digestion with deoxyribonuclease and micrococcal nuclease have been combined with high-throughput sequencing to provide detailed nucleosome occupancy maps, and chromosome conformation capture and its variants have revealed that higher-order chromatin structure involves long-range loop formation between distant genomic elements. This review discusses the methods for identifying higher-order chromatin structure and the information they have provided on this important topic.

Structural features based genome-wide characterization and prediction of nucleosome organization

Background

Nucleosome distribution along chromatin dictates genomic DNA accessibility and thus profoundly influences gene expression. However, the underlying mechanism of nucleosome formation remains elusive. Here, taking a structural perspective, we systematically explored nucleosome formation potential of genomic sequences and the effect on chromatin organization and gene expression in S. cerevisiae.

Results

We analyzed twelve structural features related to flexibility, curvature and energy of DNA sequences. The results showed that some structural features such as DNA denaturation, DNA-bending stiffness, Stacking energy, Z-DNA, Propeller twist and free energy, were highly correlated with in vitro and in vivo nucleosome occupancy. Specifically, they can be classified into two classes, one positively and the other negatively correlated with nucleosome occupancy. These two kinds of structural features facilitated nucleosome binding in centromere regions and repressed nucleosome formation in the promoter regions of protein-coding genes to mediate transcriptional regulation. Based on these analyses, we integrated all twelve structural features in a model to predict more accurately nucleosome occupancy in vivo than the existing methods that mainly depend on sequence compositional features. Furthermore, we developed a novel approach, named DLaNe, that located nucleosomes by detecting peaks of structural profiles, and built a meta predictor to integrate information from different structural features. As a comparison, we also constructed a hidden Markov model (HMM) to locate nucleosomes based on the profiles of these structural features. The result showed that the meta DLaNe and HMM-based method performed better than the existing methods, demonstrating the power of these structural features in predicting nucleosome positions.

Conclusions

Our analysis revealed that DNA structures significantly contribute to nucleosome organization and influence chromatin structure and gene expression regulation. The results indicated that our proposed methods are effective in predicting nucleosome occupancy and positions and that these structural features are highly predictive of nucleosome organization.

The implementation of our DLaNe method based on structural features is available online.

Thirty years of multiple sequence codes

An overview is presented on the status of studies on multiple codes in genetic sequences. Indirectly, the existence of multiple codes is recognized in the form of several rediscoveries of Second Genetic Code that is different each time. A due credit is given to earlier seminal work related to the codes often neglected in literature. The latest developments in the field of chromatin code are discussed, as well as perspectives of single-base resolution studies of nucleosome positioning, including rotational setting of DNA on the surface of the histone octamers.

"Anticipated" nucleosome positioning pattern in prokaryotes

Linguistic (word count) analysis of prokaryotic genome sequences, by Shannon N-gram extension, reveals that the dominant hidden motifs in A+T rich genomes are T(A)(T)A and G(A)(T)C with uncertain number of repeating A and T. Since prokaryotic sequences are largely protein-coding, the motifs would correspond to amphipathic alpha-helices with alternating lysine and phenylalanine as preferential polar and non-polar residues. The motifs are also known in eukaryotes, as nucleosome positioning patterns. Their existence in prokaryotes as well may serve for binding of histone-like proteins to DNA. In this case the above patterns in prokaryotes may be considered as "anticipated" nucleosome positioning patterns which, quite likely, existed in prokaryotic genomes before the evolutionary separation between eukaryotes and prokaryotes.

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.

Nucleosome DNA sequence structure of isochores

Background

Significant differences in G+C content between different isochore types suggest that the nucleosome positioning patterns in DNA of the isochores should be different as well.

Results

Extraction of the patterns from the isochore DNA sequences by Shannon N-gram extension reveals that while the general motif YRRRRRYYYYYR is characteristic for all isochore types, the dominant positioning patterns of the isochores vary between TAAAAATTTTTA and CGGGGGCCCCCG due to the large differences in G+C composition. This is observed in human, mouse and chicken isochores, demonstrating that the variations of the positioning patterns are largely G+C dependent rather than species-specific. The species-specificity of nucleosome positioning patterns is revealed by dinucleotide periodicity analyses in isochore sequences. While human sequences are showing CG periodicity, chicken isochores display AG (CT) periodicity. Mouse isochores show very weak CG periodicity only.

Conclusions

Nucleosome positioning pattern as revealed by Shannon N-gram extension is strongly dependent on G+C content and different in different isochores. Species-specificity of the pattern is subtle. It is reflected in the choice of preferentially periodical dinucleotides.

Nucleosome positioning pattern derived from oligonucleotide compositions of genomic sequences

Availability of nucleosome positioning pattern(s) is crucial for chromatin studies. The matrix form of the pattern has been recently derived (I. Gabdank, D. Barash, E. N. Trifonov. J Biomol Struct Dyn 26, 403-412 (2009), and E. N. Trifonov. J Biomol Struct Dyn 27, 741-746 (2010)). In its simplified linear form it is described by the motif CGRAAATTTYCG. Oligonucleotide components of the motif (say, triplets GRA, RAA, AAA, etc.) would be expected to appear in eukaryotic sequences more frequently. In this work we attempted the reconstruction of the bendability patterns for 13 genomes by a novel approach-extension of highest frequency trinucleotides. The consensus of the patterns reconstructed on the basis of trinucleotide frequencies in 13 eukaryotic genomes is derived: CRAAAATTTTYG. It conforms to the earlier established sequence motif. The reconstruction, thus, attests to the universality of the nucleosome DNA bendability pattern.

Predicting nucleosome positioning in genomes: physical and bioinformatic approaches

In eukaryotic genomes, nucleosomes are responsible for packaging DNA and controlling gene expression. For this reason, an increasing interest is arising on computational methods capable of predicting the nucleosome positioning along genomes. In this review we describe and compare bioinformatic and physical approaches adopted to predict nucleosome occupancy along genomes. Computational analyses attempt at decoding the experimental nucleosome maps of genomes in terms of certain dinucleotide step periodicity observed along DNA. Such investigations show that highly significant information about the occurrence of a nucleosome along DNA is intrinsic in certain features of the sequence suggesting that DNA of eukaryotic genomes encodes nucleosome organization. Besides the bioinformatic approaches, physical models were proposed based on the sequence dependent conformational features of the DNA chain, which govern the free energy needed to transform recurrent DNA tracts along the genome into the nucleosomal shape.

An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression

BACKGROUND

A nucleosome is the fundamental repeating unit of the eukaryotic chromosome. It has been shown that the positioning of a majority of nucleosomes is primarily controlled by factors other than the intrinsic preference of the DNA sequence. One of the key questions in this context is the role, if any, that can be played by the variability of nucleosomal DNA structure.

RESULTS

In this study, we have addressed this question by analysing the variability at the dinucleotide and trinucleotide as well as longer length scales in a dataset of nucleosome X-ray crystal structures. We observe that the nucleosome structure displays remarkable local level structural versatility within the B-DNA family. The nucleosomal DNA also incorporates a large number of kinks.

CONCLUSIONS

Based on our results, we propose that the local and global level versatility of B-DNA structure may be a significant factor modulating the formation of nucleosomes in the vicinity of high-plasticity genes, and in varying the probability of binding by regulatory proteins. Hence, these factors should be incorporated in the prediction algorithms and there may not be a unique 'template' for predicting putative nucleosome sequences. In addition, the multimodal distribution of dinucleotide parameters for some steps and the presence of a large number of kinks in the nucleosomal DNA structure indicate that the linear elastic model, used by several algorithms to predict the energetic cost of nucleosome formation, may lead to incorrect results.

Evaluation of elastic rod models with long range interactions for predicting nucleosome stability

The ability of a dinucleotide-step based elastic-rod model of DNA to predict nucleosome binding free energies is investigated using four available sets of elastic parameters. We compare the predicted free energies to experimental values derived from nucleosome reconstitution experiments for 84 DNA sequences. Elastic parameters (conformation and stiffnessess) obtained from MD simulations are shown to be the most reliable predictors, as compared to those obtained from analysis of base-pair step melting temperatures, or from analysis of x-ray structures. We have also studied the effect of varying the folded conformation of nucleosomal DNA by means of our Fourier - filtering knock-out and knock-in procedure. This study confirmed the above ranking of elastic parameters, and helped to reveal problems inherent in models using only a local elastic energy function. Long-range interactions were added to the elastic-rod model in an effort to improve its predictive ability. For this purpose a Debye-Huckel energy term with a single, homogenous point charge per base-pair was introduced. This term contains only three parameters, - its weight relative to the elastic energy, the Debye screening length, and a minimum sequence distance for including pairwise interactions between charges. After optimization of these parameters, our Debye-Huckel term is attractive, and yields the same level of correlation with experiment (R=0.75) as was achieved merely by varying the nucleosomal shape in the elastic-rod model. We suggest this result indicates a linker DNA - histone attraction or, possibly, entropic effects, that lead to a stabilization of a nucleosome away from the ends of DNA segments longer than 147 bp. Such effects are not accounted for by a localized elastic energy model.

A structural perspective on the where, how, why, and what of nucleosome positioning

The DNA in eukaryotic chromatin is packed by histones into arrays of repeating units called nucleosomes. Each nucleosome contains a nucleosome core, where the DNA is wrapped around a histone octamer, and a stretch of relatively unconstrained DNA called the linker DNA. Since nucleosome cores occlude the DNA from many DNA-binding factors, their positions provide important clues for understanding chromatin packing and gene regulation. Here we review the recent advances in the genome-wide mapping of nucleosome positions, the molecular and structural determinants of nucleosome positioning, and the importance of nucleosome positioning in chromatin higher order folding and transcriptional regulation.

Painting a perspective on the landscape of nucleosome positioning

DNA sequence influences the position of nucleosomes and chromatin architecture. The extent to which underlying DNA sequence affects nucleosome positioning is currently a topic of considerable discussion and active experimentation. To contribute to the discussion, I will outline a few of the methods, data and arguments that I find compelling and believe will ultimately resolve the question of what positions nucleosomes. Basically, I will give a portrait of my current perspective on what influences the landscape of nucleosome positioning and chromatin architecture.

Prediction of nucleosome positioning in genomes: limits and perspectives of physical and bioinformatic approaches

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. For this reason, an increasing interest is arising on a wide-ranging series of experimental and computational analyses capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. Based on a statistical mechanical approach, the model allows the calculation of the sequence-dependent canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for 90 single nucleosome DNA tracts and successfully compared with those obtained with nucleosome competitive reconstitution. These results, obtained for single nucleosomes, could in principle allow the calculation of the intrinsic affinity of nucleosomes along DNA sequences virtually opening the possibility of predicting the nucleosome positioning along genomes on physical basis. The theoretical nucleosome distribution was compared and validated with that of yeast and human genome experimentally determined. The results interpret on a physical basis the experimental nucleosome positioning and are comparable with those obtained adopting models based on the identification of some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.

Single-base resolution nucleosome mapping on DNA sequences

Nucleosome DNA bendability pattern extracted from large nucleosome DNA database of C. elegans is used for construction of full length (116 dinucleotide positions) nucleosome DNA bendability matrix. The matrix can be used for sequence-directed mapping of the nucleosomes on the sequences. Several alternative positions for a given nucleosome are typically predicted, separated by multiples of nucleosome DNA period. The corresponding computer program is successfully tested on best known experimental examples of accurately positioned nucleosomes. The uncertainty of the computational mapping is +/-1 base. The procedure is placed on publicly accessible server and can be applied to any DNA sequence of interest.

The DNA sequence-dependence of nucleosome positioning in vivo and in vitro

The contribution of histone-DNA interactions to nucleosome positioning in vivo is currently a matter of debate. We argue here that certain nucleosome positions, often in promoter regions, in yeast may be, at least in part, specified by the DNA sequence. In contrast other positions may be poorly specified. Positioning thus has both statistical and DNA-determined components. We further argue that the relative affinity of the octamer for different DNA sequences can vary and therefore the interaction of histones with the DNA is a 'tunable' property.

Sequence-dependent Kink-and-Slide deformations of nucleosomal DNA facilitated by histone arginines bound in the minor groove

In addition to bending and twisting deformabilities, the lateral displacements of the DNA axis (Kink-and-Slide) play an important role in DNA wrapping around the histone core (M. Y. Tolstorukov, A. V. Colasanti, D. M. McCandlish, W. K. Olson, V. B. Zhurkin, J. Mol. Biol. 371, 725-738 (2007)). Here, we show that these Kink-and-Slide deformations are likely to be stabilized by the arginine residues of histones interacting with the minor groove of DNA. The arginines are positioned asymmetrically in the minor groove, being closer to one strand. The asymmetric arginine-DNA interactions facilitate lateral displacement of base pairs across the DNA grooves, thus leading to a stepwise accumulation of the superhelical pitch of nucleosomal DNA. To understand the sequence dependence of such Kink-and-Slide deformations, we performed all-atom calculations of DNA hexamers with the YR and RY steps in the center. We found that when the unrestrained DNA deformations are allowed, the YR steps tend to bend into the major groove, and RY steps bend into the minor groove. However, when the nucleosomal Kink-and-Slide deformation is considered, the YR steps prove to be more favorable for bending into the minor groove. Overall, the Kink-and-Slide deformation energy of DNA increases in the order TA < CA < CG < GC < AC < AT. We propose a simple stereochemical model accounting for this sequence dependence. Our results agree with experimental data indicating that the TA step most frequently occurs in the minor-groove kink positions in the most stable nucleosomes. Our computations demonstrate that the Kink-and-Slide distortion is accompanied by the BI to BII transition. This fact, together with irregularities in the two-dimensional (Roll, Slide) energy contour maps, suggest that the Kink-and-Slide deformations represent a nonharmonic behavior of the duplex. This explains the difference between the two estimates of the DNA deformation energy in nucleosome - the earlier one made using knowledge-based elastic energy functions, and the current one based on all-atom calculations. Our findings are useful for refining the score functions for the prediction of nucleosome positioning. In addition, the reverse bending behavior of the YR and RY steps revealed under the Kink-and-Slide constraint is important for understanding the molecular mechanisms of binding transcription factors (such as p53) to DNA exposed on the surface of nucleosome.

Structure-based analysis of DNA sequence patterns guiding nucleosome positioning in vitro

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 beta-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 alpha-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.

Sequence periodicity in nucleosomal DNA and intrinsic curvature

BACKGROUND

Most eukaryotic DNA contained in the nucleus is packaged by wrapping DNA around histone octamers. Histones are ubiquitous and bind most regions of chromosomal DNA. In order to achieve smooth wrapping of the DNA around the histone octamer, the DNA duplex should be able to deform and should possess intrinsic curvature. The deformability of DNA is a result of the non-parallelness of base pair stacks. The stacking interaction between base pairs is sequence dependent. The higher the stacking energy the more rigid the DNA helix, thus it is natural to expect that sequences that are involved in wrapping around the histone octamer should be unstacked and possess intrinsic curvature. Intrinsic curvature has been shown to be dictated by the periodic recurrence of certain dinucleotides. Several genome-wide studies directed towards mapping of nucleosome positions have revealed periodicity associated with certain stretches of sequences. In the current study, these sequences have been analyzed with a view to understand their sequence-dependent structures.

RESULTS

Higher order DNA structures and the distribution of molecular bend loci associated with 146 base nucleosome core DNA sequence from C. elegans and chicken have been analyzed using the theoretical model for DNA curvature. The curvature dispersion calculated by cyclically permuting the sequences revealed that the molecular bend loci were delocalized throughout the nucleosome core region and had varying degrees of intrinsic curvature.

CONCLUSIONS

The higher order structures associated with nucleosomes of C.elegans and chicken calculated from the sequences revealed heterogeneity with respect to the deviation of the DNA axis. The results points to the possibility of context dependent curvature of varying degrees to be associated with nucleosomal DNA.

Identification and characterization of putative methylation targets in the MAOA locus using bioinformatic approaches

Monoamine oxidase A (MAO A) is an enzyme that catalyzes the oxidation of neurotransmitter amines. A functional polymorphism in the human MAOA gene (high- and low-MAOA) has been associated with distinct behavioral phenotypes. To investigate directly the biological mechanism whereby this polymorphism influences brain function, we recently measured the activity of the MAO A enzyme in healthy volunteers. When found no relationship between the individual's brain MAO A level and the MAOA genotype, we postulated that there are additional regulatory mechanisms that control the MAOA expression. Given that DNA methylation is linked to the regulation of gene expression, we hypothesized that epigenetic mechanisms factor into the MAOA expression. Our underplaying assumption was that the differences in an individual's genotype play a key role in the epigenetic potential of the MAOA locus and, consequently, determine the individual's level of MAO A activity in the brain. As a first step towards experimental validation of the hypothesis, we performed a comprehensive bioinformatic analysis aiming to interrogate genomic features and attributes of the MAOA locus that might modulate its epigenetic sensitivity. Major findings of our analysis are the following: (1) the extended MAOA regulatory region contains two CpG islands (CGIs), one of which overlaps with the canonical MAOA promoter and the other is located further upstream; both CGIs exhibit sensitivity to differential methylation. (2) The uVNTR's effect on the MAOA's transcriptional activity might have epigenetic nature: this polymorphic region resides within the MAOA's CGI and itself contains CpGs, thus, the number of repeating increments effectively changes the number of methylatable cytosines in the MAOA promoter. An array of in silico analyses (the nucleosome positioning, the physical properties of the local DNA, the clustering of transcription-factor binding sites) together with experimental data on histone modifications and Pol 2 sites and data from the RefSeq mRNA library suggest that the MAOA gene might have an alternative promoter. Based on our findings, we propose a regulatory mechanism for the human MAOA according to which the MAOA expression in vivo is executed by the generation of tissue-specific transcripts initiated from the alternative promoters (both CGI-associated) where transcriptional activation of a particular promoter is under epigenetic control.

A cis-regulatory signature in ascidians and flies, independent of transcription factor binding sites

BACKGROUND

Transcription initiation is controlled by cis-regulatory modules. Although these modules are usually made of clusters of short transcription factor binding sites, a small minority of such clusters in the genome have cis-regulatory activity. This paradox is currently unsolved.

RESULTS

To identify what discriminates active from inactive clusters, we focused our attention on short topologically unconstrained clusters of two ETS and two GATA binding sites, similar to the early neural enhancer of Ciona intestinalis Otx. We first computationally identified 55 such clusters, conserved between the two Ciona genomes. In vivo assay of the activity of 19 hits identified three novel early neural enhancers, all located next to genes coexpressed with Otx. Optimization of ETS and GATA binding sites was not always sufficient to confer activity to inactive clusters. Rather, a dinucleotide sequence code associated to nucleosome depletion showed a robust correlation with enhancer potential. Identification of a large collection of Ciona regulatory regions revealed that predicted nucleosome depletion constitutes a general signature of Ciona enhancers, which is conserved between orthologous loci in the two Ciona genomes and which partitions conserved noncoding sequences into a major nucleosome-bound fraction and a minor nucleosome-free fraction with higher cis-regulatory potential. We also found this signature in a large fraction of short Drosophila cis-regulatory modules.

CONCLUSION

This study indicates that a sequence-based dinucleotide signature, previously associated with nucleosome depletion and independent of transcription factor binding sites, contributes to the definition of a local cis-regulatory potential in two metazoa, Ciona intestinalis and Drosophila melanogaster.

Analysis of primary structure of chromatin with next-generation sequencing

The recent development of next-generation sequencing technology has enabled significant progress in chromatin structure analysis. Here, we review the experimental and bioinformatic approaches to studying nucleosome positioning and histone modification profiles on a genome scale using this technology. These studies advanced our knowledge of the nucleosome positioning patterns of both epigenetically modified and bulk nucleosomes and elucidated the role of such patterns in regulation of gene expression. The identification and analysis of large sets of nucleosome-bound DNA sequences allowed better understanding of the rules that govern nucleosome positioning in organisms of various complexity. We also discuss the existing challenges and prospects of using next-generation sequencing for nucleosome positioning analysis and outline the importance of such studies for the entire chromatin structure field.