Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-Seq, ATAC-Seq and NOMe-Seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However we also find a high proportion of assay specific NDRs that are often "called" by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-Seq and DNAse I-Seq provide a high NDR calling rate for relatively low sequencing costs in comparison to NOMe-Seq, NOMe-Seq singles out as it provides a multitude of information: it allows to not only detect NDRs but also endogenous DNA methylation, genome wide segmentation into heterochromatic A/B domains and local phasing of nucleosomes outside of NDRs. In summary our comparison strongly suggest to consider assay specific differences for the experimental desgin and for generalized and comparative functional interpretations.