Nucleosome positioning based on the sequence word composition

The DNA of all eukaryotic organisms is packaged into nucleosomes (a basic repeating unit of chromatin). A nucleosome consists of histone octamer wrapped by core DNA and linker histone H1 associated with linker DNA. It has profound effects on all DNA-dependent processes by affecting sequence accessibility. Understanding the factors that influence nucleosome positioning has great help to the study of genomic control mechanism. Among many determinants, the inherent DNA sequence has been suggested to have a dominant role in nucleosome positioning in vivo. Here, we used the method of minimum redundancy maximum relevance (mRMR) feature selection and the nearest neighbor algorithm (NNA) combined with the incremental feature selection (IFS) method to identify the most important sequence features that either favor or inhibit nucleosome positioning. We analyzed the words of 53,021 nucleosome DNA sequences and 50,299 linker DNA sequences of Saccharomyces cerevisiae. 32 important features were abstracted from 5,460 features, and the overall prediction accuracy through jackknife cross-validation test was 76.5%. Our results support that sequence-dependent DNA flexibility plays an important role in positioning nucleosome core particles and that genome sequence facilitates the rapid nucleosome reassembly instead of nucleosome depletion. Besides, our results suggest that there exist some additional features playing a considerable role in discriminating nucleosome forming and inhibiting sequences. These results confirmed that the underlying DNA sequence plays a major role in nucleosome positioning.     

核小体定位的转录调控功能研究进展

核小体是真核生物染色质的基本组成单位,组蛋白八聚体在DNA 双螺旋上精确位置称为核小体定位．核小体定位已被证实在基因转录调控、DNA复制与修复、调控进化等过程中扮演着重要的角色．随着染色质免疫共沉淀-芯片(ChIP-chip)与染色质免疫共沉淀-测序(ChIP-seq)等高通量技术的出现,已测定了多种模式生物全基因组核小体定位图谱,掀起了一股核小体定位及其功能的研究热潮,并取得了一定的成果．本文介绍了核小体定位的概念,总结了核小体在启动子与编码区域内定位的基本模式．在此基础上,综述了核小体定位在转录起始、转录延伸、基因表达模式多样化以及可变剪接等方面的功能研究进展．     

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.     

Nucleosome destabilization in the epigenetic regulation of gene expression

Assembly, mobilization and disassembly of nucleosomes can influence the regulation of gene expression and other processes that act on eukaryotic DNA. Distinct nucleosome-assembly pathways deposit dimeric subunits behind the replication fork or at sites of active processes that mobilize pre-existing nucleosomes. Replication-coupled nucleosome assembly appears to be the default process that maintains silent chromatin, counteracted by active processes that destabilize nucleosomes. Nucleosome stability is regulated by the combined effects of nucleosome-positioning sequences, histone chaperones, ATP-dependent nucleosome remodellers, post-translational modifications and histone variants. Recent studies suggest that histone turnover helps to maintain continuous access to sequence-specific DNA-binding proteins that regulate epigenetic inheritance, providing a dynamic alternative to histone-marking models for the propagation of active chromatin.     

Nucleosomes Shape DNA Polymorphism and Divergence

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

NucleoMap: A computational tool for identifying nucleosomes in ultra-high resolution contact maps

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.     

Active Nucleosome Displacement: A Theoretical Approach

Three-quarters of eukaryotic DNA are wrapped around protein cylinders forming so-called nucleosomes that block the access to the genetic information. Nucleosomes need therefore to be repositioned, either passively (by thermal fluctuations) or actively (by molecular motors). Here we introduce a theoretical model that allows us to study the interplay between a motor protein that moves along DNA (e.g., an RNA polymerase) and a nucleosome that it encounters on its way. We aim at describing the displacement mechanisms of the nucleosome and the motor protein on a microscopic level to understand better the intricate interplay between the active step of the motor and the nucleosome-repositioning step. Different motor types (Brownian ratchet versus power-stroke mechanism) that perform very similarly under a constant load are shown to have very different nucleosome repositioning capacities.     

Nucleosome positioning: How is it established, and why does it matter?

Packaging of eukaryotic genomes into chromatin affects every process that occurs on DNA. The positioning of nucleosomes on underlying DNA plays a key role in the regulation of these processes, as the nucleosome occludes underlying DNA sequences. Here, we review the literature on mapping nucleosome positions in various organisms, and discuss how nucleosome positions are established, what effect nucleosome positioning has on control of gene expression, and touch on the correlations between chromatin packaging, sequence evolution, and the evolution of gene expression programs.     

Telomeric nucleosomes are intrinsically mobile

Nucleosomes are no longer considered only static basic units that package eukaryotic DNA but they emerge as dynamic players in all chromosomal processes. Regulatory proteins can gain access to recognition sequences hidden by the histone octamer through the action of ATP-dependent chromatin remodeling complexes that cause nucleosome sliding. In addition, it is known that nucleosomes are able to spontaneously reposition along the DNA due to intrinsic dynamic properties, but it is not clear yet to what extent sequence-dependent dynamic properties contribute to nucleosome repositioning. Here, we study mobility of nucleosomes formed on telomeric sequences as a function of temperature and ionic strength. We find that telomeric nucleosomes are highly intrinsically mobile under physiological conditions, whereas nucleosomes formed on an average DNA sequence mostly remain in the initial position. This indicates that DNA sequence affects not only the thermodynamic stability and the positioning of nucleosomes but also their dynamic properties. Moreover, our findings suggest that the high mobility of telomeric nucleosomes may be relevant to the dynamics of telomeric chromatin.     

Nucleosome positioning and nucleosome stacking: two faces of the same coin

Nucleosomes are regularly spaced along eukaryotic genomes. In the emerging model, known as "statistical positioning", this spacing is due to steric repulsion between nucleosomes and to the presence of nucleosome excluding barriers on the genome. However, new experimental evidence recently challenged the "statistical positioning" model (Z. Zhang et al., Science, 2011, 332(6032), 977-980). We propose here that the regular spacing can be better explained by adding attractive interactions between nucleosomes. In our model those attractions are due to the fact that nucleosomes are stacked in regular chromatin fibers. In a self-reinforcing mechanism, regular nucleosome spacing promotes in turn nucleosome stacking. We first show that this model can precisely account for the nucleosome spacing observed in Saccharomyces cerevisiae. We then use a simple toy model to show that attraction between nucleosomes can fasten the formation of the chromatin fiber.     

Nucleosome Interactions and Stability in an Ordered Nucleosome Array Model System

Although it is well established that the majority of eukaryotic DNA is sequestered as nucleosomes, the higher-order structure resulting from nucleosome interactions as well as the dynamics of nucleosome stability are not as well understood. To characterize the structural and functional contribution of individual nucleosomal sites, we have developed a chromatin model system containing up to four nucleosomes, where the array composition, saturation, and length can be varied via the ordered ligation of distinct mononucleosomes. Using this system we find that the ligated tetranucleosomal arrays undergo intra-array compaction. However, this compaction is less extensive than for longer arrays and is histone H4 tail-independent, suggesting that well ordered stretches of four or fewer nucleosomes do not fully compact to the 30-nm fiber. Like longer arrays, the tetranucleosomal arrays exhibit cooperative self-association to form species composed of many copies of the array. This propensity for self-association decreases when the fraction of nucleosomes lacking H4 tails is systematically increased. However, even tetranucleosomal arrays with only two octamers possessing H4 tails recapitulate most of the inter-array self-association. Varying array length shows that systems as short as dinucleosomes demonstrate significant self-association, confirming that relatively few determinants are required for inter-array interactions and suggesting that in vivo multiple interactions of short runs of nucleosomes might contribute to complex fiber-fiber interactions. Additionally, we find that the stability of nucleosomes toward octamer loss increases with array length and saturation, suggesting that in vivo stretches of ordered, saturated nucleosomes could serve to protect these regions from histone ejection.     

CENP-A Nucleosome is a Sensitive Allosteric Scaffold for DNA and Chromatin Factors

Centromeric loci of chromosomes are defined by nucleosomes containing the histone H3 variant CENP-A, which bind their DNA termini more permissively than their canonical counterpart, a feature that is critical for the mitotic fidelity. A recent cryo-EM study demonstrated that the DNA termini of CENP-A nucleosomes, reconstituted with the Widom 601 DNA sequence, are asymmetrically flexible, meaning one terminus is more clearly resolved than the other. However, an earlier work claimed that both ends could be resolved in the presence of two stabilizing single chain variable fragment (scFv) antibodies per nucleosome, and thus are likely permanently bound to the histone octamer. This suggests that the binding of scFv antibodies to the histone octamer surface would be associated with CENP-A nucleosome conformational changes, including stable binding of the DNA termini. Here, we present computational evidence that allows to explain at atomistic level the structural rearrangements of CENP-A nucleosomes resulting from the antibody binding. The antibodies, while they only bind the octamer façades, are capable of altering the dynamics of the nucleosomal core, and indirectly also the surrounding DNA. This effect has more drastic implications for the structure and the dynamics of the CENP-A nucleosome in comparison to its canonical counterpart. Furthermore, we find evidence that the antibodies bind the left and the right octamer façades at different affinities, another manifestation of the DNA sequence. We speculate that the cells could use induction of similar allosteric effects to control centromere function.