Chemical map-based prediction of nucleosome positioning using the Bioconductor package nuCpos

Background

Assessing the nucleosome-forming potential of specific DNA sequences is important for understanding complex chromatin organization. Methods for predicting nucleosome positioning include bioinformatics and biophysical approaches. An advantage of bioinformatics methods, which are based on in vivo nucleosome maps, is the use of natural sequences that may contain previously unknown elements involved in nucleosome positioning in vivo. The accuracy of such prediction attempts reflects the genomic coordinate resolution of the nucleosome maps applied. Nucleosome maps are constructed using micrococcal nuclease digestion followed by high-throughput sequencing (MNase-seq). However, as MNase has a strong preference for A/T-rich sequences, MNase-seq may not be appropriate for this purpose. In addition to MNase-seq-based maps, base pair-resolution chemical maps of in vivo nucleosomes from three different species (budding and fission yeasts, and mice) are currently available. However, these chemical maps have yet to be integrated into publicly available computational methods.

Results

We developed a Bioconductor package (named nuCpos) to demonstrate the superiority of chemical maps in predicting nucleosome positioning. The accuracy of chemical map-based prediction in rotational settings was higher than that of the previously developed MNase-seq-based approach. With our method, predicted nucleosome occupancy reasonably matched in vivo observations and was not affected by A/T nucleotide frequency. Effects of genetic alterations on nucleosome positioning that had been observed in living yeast cells could also be predicted. nuCpos calculates individual histone binding affinity (HBA) scores for given 147-bp sequences to examine their suitability for nucleosome formation. We also established local HBA as a new parameter to predict nucleosome formation, which was calculated for 13 overlapping nucleosomal DNA subsequences. HBA and local HBA scores for various sequences agreed well with previous in vitro and in vivo studies. Furthermore, our results suggest that nucleosomal subsegments that are disfavored in different rotational settings contribute to the defined positioning of nucleosomes.

Conclusions

Our results demonstrate that chemical map-based statistical models are beneficial for studying nucleosomal DNA features. Studies employing nuCpos software can enhance understanding of chromatin regulation and the interpretation of genetic alterations and facilitate the design of artificial sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04240-2.

Sequence-dependent nucleosome formation in trinucleotide repeats evaluated by in vivo chemical mapping

Trinucleotide repeat sequences (TRSs), consisting of 10 unique classes of repeats in DNA, are members of microsatellites and abundantly and non-randomly distributed in many eukaryotic genomes. The lengths of TRSs are mutable, and the expansions of several TRSs are implicated in hereditary neurological diseases. However, the underlying causes of the biased distribution and the dynamic properties of TRSs in the genome remain elusive. Here, we examined the effects of TRSs on nucleosome formation in vivo by histone H4-S47C site-directed chemical cleavages, using well-defined yeast minichromosomes in which each of the ten TRS classes resided in the central region of a positioned nucleosome. We showed that (AAT)₁₂ and (ACT)₁₂ act as strong nucleosome-promoting sequences, while (AGG)₁₂ and (CCG)₁₂ act as nucleosome-excluding sequences in vivo. The local histone binding affinity scores support the idea that nucleosome formation in TRSs, except for (AGG)₁₂, is mainly determined by the affinity for the histone octamers. Overall, our study presents a framework for understanding the nucleosome-forming abilities of TRSs.