NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background: Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them. Results: Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs). Conclusions: The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed

BZLF1 interacts with chromatin remodelers promoting escape from latent infections with EBV.

A hallmark of EBV infections is its latent phase, when all viral lytic genes are repressed. Repression results from a high nucleosome occupancy and epigenetic silencing by cellular factors such as the Polycomb repressive complex 2 (PRC2) and DNA methyltransferases that, respectively, introduce repressive histone marks and DNA methylation. The viral transcription factor BZLF1 acts as a molecular switch to induce transition from the latent to the lytic or productive phase of EBV’s life cycle. It is unknown how BZLF1 can bind to the epigenetically silenced viral DNA and whether it directly reactivates the viral genome through chromatin remodeling. We addressed these fundamental questions and found that BZLF1 binds to nucleosomal DNA motifs both in vivo and in vitro. BZLF1 co-precipitates with cellular chromatin remodeler ATPases, and the knock-down of one of them, INO80, impaired lytic reactivation and virus synthesis. In Assay for Transposase-Accessible Chromatin-seq experiments, non-accessible chromatin opens up locally when BZLF1 binds to its cognate sequence motifs in viral DNA. We conclude that BZLF1 reactivates the EBV genome by directly binding to silenced chromatin and recruiting cellular chromatin-remodeling enzymes, which implement a permissive state for lytic viral transcription. BZLF1 shares this mode of action with a limited number of cellular pioneer factors, which are instrumental in transcriptional activation, differentiation, and reprogramming in all eukaryotic cells