Visible periodicity of strong nucleosome DNA sequences

A study of strong nucleosomes in the human genome

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.

•

Naked DNA correction improved the accuracy of nucleosome map in partial digestion

•

Level of MNase digestion has effects on nucleosome organization

•

A type of strong nucleosomes (SNs) exist across different nucleosome maps

•

Nucleosome stability may be related to its location and the DNA sequence

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

Protein intrinsic disorder negatively associates with gene age in different eukaryotic lineages

The emergence of new protein-coding genes in a specific lineage or species provides raw materials for evolutionary adaptations. Until recently, the biology of new genes emerging particularly from non-genic sequences remained unexplored. Although the new genes are subjected to variable selection pressure and face rapid deletion, some of them become functional and are retained in the gene pool. To acquire functional novelties, new genes often get integrated into the pre-existing ancestral networks. However, the mechanism by which young proteins acquire novel interactions remains unanswered till date. Since structural orientation contributes hugely to the mode of proteins' physical interactions, in this regard, we put forward an interesting question - Do new genes encode proteins with stable folds? Addressing the question, we demonstrated that the intrinsic disorder inversely correlates with the evolutionary gene ages - i.e. young proteins are richer in intrinsic disorder than the ancient ones. We further noted that young proteins, which are initially poorly connected hubs, prefer to be structurally more disordered than well-connected ancient proteins. The phenomenon strikingly defies the usual trend of well-connected proteins being highly disordered in structure. We justified that structural disorder might help poorly connected young proteins to undergo promiscuous interactions, which provides the foundation for novel protein interactions. The study focuses on the evolutionary perspectives of young proteins in the light of structural adaptations.

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.

Transcription factors operate TATA switches via rotational remodeling of local columnar chromatin structure

Our earlier study on the nucleosomes containing TBP binding sites (TATA boxes) indicated that generally the same sequence, which harbors the TATA box, encodes simultaneously an alternative rotational setting of the box, so that the TATA element is either exposed (position "minor groove out") or hidden in position "minor groove in". The sequence elements (dinucleotides) residing on the inner surface of DNA in contact with histone octamers are identified by calculating YR tracks in the promoter regions of the genes - periodically reappearing YR elements, at distances of 10-11 bases from one another. Non-YR elements of the YR tracks are also verified by nucleosome mapping procedure based on alternation of runs of purines with runs of pyrimidines. The tracks observed in yeast promoter regions are found to split into two, passing further downstream either through the TATA element or 4-6 bases toward the tail of the box. The points of splitting which play the role of TATA switches are located in close vicinity or within the transcription factor binding sites. This suggests the regulatory function of transcription factors changing the YR tracks by over- or under-twisting DNA induced by their binding.

Nucleosome repeat lengths and columnar chromatin structure

Thorough quantitative study of nucleosome repeat length (NRL) distributions, conducted in 1992 by J. Widom, resulted in a striking observation that the linker lengths between the nucleosomes are quantized. Comparison of the NRL average values with the MNase cut distances predicted from the hypothetical columnar structure of chromatin (this work) shows a close correspondence between the two. This strongly suggests that the NRL distribution, actually, reflects the dominant role of columnar chromatin structure common for all eukaryotes.

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.

Sequence structure of Lowary/Widom clones forming strong nucleosomes

Lowary and Widom selected from random sequences those which form exceptionally stable nucleosomes, including clone 601, the current champion of strong nucleosome (SN) sequences. This unique sequence database (LW sequences) carries sequence elements which confer stability on the nucleosomes formed on the sequences, and, thus, may serve as source of information on the structure of "ideal" or close to ideal nucleosome DNA sequence. An important clue is also provided by crystallographic study of Vasudevan and coauthors on clone 601 nucleosomes. It demonstrated that YR·YR dinucleotide stacks (primarily TA·TA) follow one another at distances 10 or 11 bases or multiples thereof, such that they all are located on the interface between DNA and histone octamer. Combining this important information with alignment of the YR-containing 10-mers and 11-mers from LW sequences, the bendability matrices of the stable nucleosome DNA are derived. The matrices suggest that the periodically repeated TA (YR), RR, and YY dinucleotides are the main sequence features of the SNs. This consensus coincides with the one for recently discovered SNs with visibly periodic DNA sequences. Thus, the experimentally observed stable LW nucleosomes and SNs derived computationally appear to represent the same entity - exceptionally stable SNs.

Strong nucleosomes of yeasts

Yeast genome lacks visibly periodic sequences characteristic of strong nucleosomes (SNs) originally discovered in A. thaliana, C. elegans, and H. sapiens. Yet, the sequences with good match to the (RRRRRYYYYY)n consensus of the SNs do show preference to centromere regions of Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Cryptococcus neoformans - property characteristic of SNs of higher eukaryotes. Candida albicans is the first exception detected so far, where their SNs do not have any affinity to the centromeres, nor pericentromeric regions. Three of the four yeast genomes analyzed possess unique repeating centromere-specific SN sequences (C. albicans, again, is an exception). The results firmly indicate that centromeres of plants, animals, and yeasts in general have special chromatin structure, favoring SNs.

TA-periodic ("601"-like) centromeric nucleosomes of A. thaliana

The bulk of strong nucleosomes (SNs, with visibly periodic DNA sequences) is described by consensus pattern of 5 or 6 base runs of purines alternating with similar runs of pyrimidines - RR/YY SNs. Yet, the strongest known nucleosome positioning sequence, the 601 clone of Lowary and Widom, is rather periodic repetition of TA dinucleotides following one another every 10 bases. We located "601"-like TA-periodic sequences in the genome of A. thaliana. Several families of such sequences are discovered repeating almost exclusively in centromeres. Thus, while A. thaliana SNs of RR/YY type have strong affinity to pericentromeric regions, as it has been previously found, the SNs of TA periodic type concentrate rather in centromeres.

Reading sequence-directed computational nucleosome maps

Recently developed latest version of the sequence-directed single-base resolution nucleosome mapping reveals existence of strong nucleosomes and chromatin columnar structures (columns). Broad application of this simple technique for further studies of chromatin and chromosome structure requires some basic understanding as to how it works and what information it affords. The paper provides such an introduction to the method. The oscillating maps of singular nucleosomes, of short and long oligonucleosome columns, are explained, as well as maps of chromatin on satellite DNA and occurrences of counter-phase (antiparallel) nucleosome neighbors.

Strong nucleosomes of mouse genome including recovered centromeric sequences

Recently discovered strong nucleosomes (SNs) characterized by visibly periodical DNA sequences have been found to concentrate in centromeres of Arabidopsis thaliana and in transient meiotic centromeres of Caenorhabditis elegans. To find out whether such affiliation of SNs to centromeres is a more general phenomenon, we studied SNs of the Mus musculus. The publicly available genome sequences of mouse, as well as of practically all other eukaryotes do not include the centromere regions which are difficult to assemble because of a large amount of repeat sequences in the centromeres and pericentromeric regions. We recovered those missing sequences using the data from MNase-seq experiments in mouse embryonic stem cells, where the sequence of DNA inside nucleosomes, including missing regions, was determined by 100-bp paired-end sequencing. Those nucleosome sequences, which are not matching to the published genome sequence, would largely belong to the centromeres. By evaluating SN densities in centromeres and in non-centromeric regions, we conclude that mouse SNs concentrate in the centromeres of telocentric mouse chromosomes, with ~3.9 times excess compared to their density in the rest of the genome. The remaining non-centromeric SNs are harbored mainly by introns and intergenic regions, by retro-transposons, in particular. The centromeric involvement of the SNs opens new horizons for the chromosome and centromere structure studies.

Regulation of the nucleosome repeat length in vivo by the DNA sequence, protein concentrations and long-range interactions

The nucleosome repeat length (NRL) is an integral chromatin property important for its biological functions. Recent experiments revealed several conflicting trends of the NRL dependence on the concentrations of histones and other architectural chromatin proteins, both in vitro and in vivo, but a systematic theoretical description of NRL as a function of DNA sequence and epigenetic determinants is currently lacking. To address this problem, we have performed an integrative biophysical and bioinformatics analysis in species ranging from yeast to frog to mouse where NRL was studied as a function of various parameters. We show that in simple eukaryotes such as yeast, a lower limit for the NRL value exists, determined by internucleosome interactions and remodeler action. For higher eukaryotes, also the upper limit exists since NRL is an increasing but saturating function of the linker histone concentration. Counterintuitively, smaller H1 variants or non-histone architectural proteins can initiate larger effects on the NRL due to entropic reasons. Furthermore, we demonstrate that different regimes of the NRL dependence on histone concentrations exist depending on whether DNA sequence-specific effects dominate over boundary effects or vice versa. We consider several classes of genomic regions with apparently different regimes of the NRL variation. As one extreme, our analysis reveals that the period of oscillations of the nucleosome density around bound RNA polymerase coincides with the period of oscillations of positioning sites of the corresponding DNA sequence. At another extreme, we show that although mouse major satellite repeats intrinsically encode well-defined nucleosome preferences, they have no unique nucleosome arrangement and can undergo a switch between two distinct types of nucleosome positioning.

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.

Strong nucleosomes reside in meiotic centromeres of C. elegans

Recently discovered strong nucleosomes (SNs) are characterized by strongly periodical DNA sequence, with visible rather than hidden sequence periodicity. In a quest for possible functions of the SNs, it has been found that the SNs concentrate within centromere regions of A. thaliana chromosomes . They, however, have been detected in Caenorhabditis elegans as well, although the holocentric chromosomes of this species do not have centromeres. Scrutinizing the SNs of C. elegans and their distributions along the DNA sequences of the chromosomes, we have discovered that the SNs are located mainly at the ends of the chromosomes of C. elegans. This suggests that, perhaps, the ends of the chromosomes fulfill some function(s) of centromeres in this species, as also indicated by the cytogenetic studies on meiotic chromosomes in spermatocytes of C. elegans, where the end-to-end association is observed. The centromeric involvement of the SNs, also found in A. thaliana, opens new horizons for the chromosome and centromere structure studies.

Strong nucleosomes of A. thaliana concentrate in centromere regions

Earlier identified strongest nucleosome DNA sequences of A. thaliana, those with visible 10-11 base sequence periodicity, are mapped along chromosomes. Resulting positional distributions reveal distinct maxima, one per chromosome, located in the centromere regions. Sequence-directed nucleosome mapping demonstrates that the strong nucleosomes (SNs) make tight arrays, several 'parallel' nucleosomes each, suggesting a columnar chromatin structure. The SNs represent a new class of centromeric nucleosomes, presumably, participating in synapsis of chromatids and securing the centromere architecture.