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

Characterization of in vivo DNA-binding events of plant transcription factors by ChIP-seq : Experimental protocol and computational analysis

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is a powerful technique for genome-wide identification of in vivo binding sites of DNA-binding proteins. The technique had been used to study many DNA-binding proteins in a broad variety of species. The basis of the ChIP-seq technique is the ability to covalently cross-link DNA and proteins that are located in very close proximity. This allows the use of an antibody against the (tagged) protein of interest to specifically enrich DNAfragments bound by this protein. ChIP-seq can be performed using antibodies against the native protein or against tagged proteins. Using a specific antibody against a tag to immunoprecipitate tagged proteins eliminates the need for a specific antibody against the native protein and allows more experimental flexibility. In this chapter we present a complete workflow for experimental procedure and bioinformatic analysis that allows wet-lab biologists to perform and analyze ChIP-seq experiments.</p