A workflow for simplified analysis of ATAC-cap-seq data in R.
A workflow for simplified analysis of ATAC-cap-seq data in R.
复制标题
DOI:
10.1093/gigascience/giy080
复制
发表时间:
2018-07-01
期刊:
影响因子:
9.2
通讯作者:
MacLean D
中科院分区:
文献类型:
--
作者:
Shrestha RK;Ding P;Jones JDG;MacLean D
Assay for Transposase-Accessible Chromatin (ATAC)-cap-seq is a high-throughput sequencing method that combines ATAC-seq with targeted nucleic acid enrichment of precipitated DNA fragments. There are increased analytical difficulties arising from working with a set of regions of interest that may be small in number and biologically dependent. Common statistical pipelines for RNA sequencing might be assumed to apply but can give misleading results on ATAC-cap-seq data. A tool is needed to allow a nonspecialist user to quickly and easily summarize data and apply sensible and effective normalization and analysis. We developed atacR to allow a user to easily analyze their ATAC enrichment experiment. It provides comprehensive summary functions and diagnostic plots for studying enriched tag abundance. Application of between-sample normalization is made straightforward. Functions for normalizing based on user-defined control regions, whole library size, and regions selected from the least variable regions in a dataset are provided. Three methods for detecting differential abundance of tags from enriched methods are provided, including bootstrap t, Bayes factor, and a wrapped version of the standard exact test in the edgeR package. We compared the precision, recall, and F-score of each detection method on resampled datasets at varying replicate, significance threshold, and genes changed and found that the Bayes factor method had the greatest overall detection power, though edgeR was slightly stronger in simulations with lower numbers of genes changed. Our package allows a nonspecialist user to easily and effectively apply methods appropriate to the analysis of ATAC-cap-seq in a reproducible manner. The package is implemented in pure R and is fully interoperable with common workflows in Bioconductor.
登录
查看更多内容
影响因子:
48
作者:
Huber W;Carey VJ;Gentleman R;Anders S;Carlson M;Carvalho BS;Bravo HC;Davis S;Gatto L;Girke T;Gottardo R;Hahne F;Hansen KD;Irizarry RA;Lawrence M;Love MI;MacDonald J;Obenchain V;Oleś AK;Pagès H;Reyes A;Shannon P;Smyth GK;Tenenbaum D;Waldron L;Morgan M
通讯作者:
Morgan M
DOI:
10.1111/tpj.12307
发表时间:
2013-11
期刊:
The Plant journal : for cell and molecular biology
影响因子:
--
作者:
Jupe F;Witek K;Verweij W;Sliwka J;Pritchard L;Etherington GJ;Maclean D;Cock PJ;Leggett RM;Bryan GJ;Cardle L;Hein I;Jones JD
通讯作者:
Jones JD
影响因子:
14.9
作者:
Lun AT;Smyth GK
通讯作者:
Smyth GK
影响因子:
12.3
作者:
Litzenburger UM;Buenrostro JD;Wu B;Shen Y;Sheffield NC;Kathiria A;Greenleaf WJ;Chang HY
通讯作者:
Chang HY
影响因子:
64.8
作者:
通讯作者:
--