Parallel mapping with site-directed hydroxyl radicals and micrococcal nuclease reveals structural features of positioned nucleosomes in vivo

Genome-Wide Nucleosome Mapping by H3Q85C-Directed Chemical Cleavage in Saccharomyces cerevisiae

Mapping the position of nucleosomes in vivo is key to our understanding of chromatin-based processes such as gene transcription, DNA replication, and repair. Methods based on micrococcal nuclease (MNase) digestion of chromatin are widely used to map nucleosomes, but these methods suffer from some limitations. More recently, nucleosome mapping methods that take advantage of our ability to convert cysteine residues-carefully inserted in histone proteins-into nucleases provide alternatives to MNase-based assays. Here, we provide a detailed protocol for the mapping of nucleosomes via H3Q85C-directed DNA cleavage in Saccharomyces cerevisiae.

A budding yeast CAGE dataset comprising two cell types

The budding yeast Saccharomyces cerevisiae is an excellent model organism for studying chromatin regulation with high-resolution genome-wide analyses. Since newly generated genome-wide data are often compared with publicly available datasets, expanding our dataset repertoire will be beneficial for the field. Information on transcription start sites (TSSs) determined at base pair resolution is essential for elucidating mechanisms of transcription and related chromatin regulation, yet no datasets that cover two different cell types are available. Here, we present a CAGE (cap analysis of gene expression) dataset for a-cells and α-cells grown in defined and rich media. Cell type-specific genes were differentially expressed as expected, ensuring the reliability of the data. Some of the differentially expressed TSSs were medium-specific or detected due to unrecognized chromosome rearrangement. By comparing the CAGE data with a high-resolution nucleosome map, major TSSs were primarily found in +1 nucleosomes, with a peak approximately 30 bp from the promoter-proximal end of the nucleosome. The dataset is available at DDBJ/GEA.

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.