Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data

Overview

Journal Genomics Proteomics Bioinformatics

Specialty Biology

Date 2024 Sep 10

PMID 39255248

Authors

Siyuan Cheng

Benpeng Miao

Tiandao Li

Guoyan Zhao

Bo Zhang

Affiliations

Soon will be listed here.

Abstract

Efficient and reliable profiling methods are essential to study epigenetics. Tn5, one of the first identified prokaryotic transposases with high DNA-binding and tagmentation efficiency, is widely adopted in different genomic and epigenomic protocols for high-throughputly exploring the genome and epigenome. Based on Tn5, the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and the Cleavage Under Targets and Tagmentation (CUT&Tag) were developed to measure chromatin accessibility and detect DNA-protein interactions. These methodologies can be applied to large amounts of biological samples with low-input levels, such as rare tissues, embryos, and sorted single cells. However, fast and proper processing of these epigenomic data has become a bottleneck because massive data production continues to increase quickly. Furthermore, inappropriate data analysis can generate biased or misleading conclusions. Therefore, it is essential to evaluate the performance of Tn5-based ATAC-seq and CUT&Tag data processing bioinformatics tools, many of which were developed mostly for analyzing chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here, we conducted a comprehensive benchmarking analysis to evaluate the performance of eight popular software for processing ATAC-seq and CUT&Tag data. We compared the sensitivity, specificity, and peak width distribution for both narrow-type and broad-type peak calling. We also tested the influence of the availability of control IgG input in CUT&Tag data analysis. Finally, we evaluated the differential analysis strategies commonly used for analyzing the CUT&Tag data. Our study provided comprehensive guidance for selecting bioinformatics tools and recommended analysis strategies, which were implemented into Docker/Singularity images for streamlined data analysis.

References

Kornberg R, Lorch Y . Primary Role of the Nucleosome. Mol Cell. 2020; 79(3):371-375. DOI: 10.1016/j.molcel.2020.07.020. View

Ou J, Liu H, Yu J, Kelliher M, Castilla L, Lawson N . ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data. BMC Genomics. 2018; 19(1):169. PMC: 5831847. DOI: 10.1186/s12864-018-4559-3. View

Criscuolo A, Brisse S . AlienTrimmer removes adapter oligonucleotides with high sensitivity in short-insert paired-end reads. Commentary on Turner (2014) Assessment of insert sizes and adapter content in FASTQ data from NexteraXT libraries. Front Genet. 2014; 5:130. PMC: 4026695. DOI: 10.3389/fgene.2014.00130. View

Boyle A, Davis S, Shulha H, Meltzer P, Margulies E, Weng Z . High-resolution mapping and characterization of open chromatin across the genome. Cell. 2008; 132(2):311-22. PMC: 2669738. DOI: 10.1016/j.cell.2007.12.014. View

Li B, Carey M, Workman J . The role of chromatin during transcription. Cell. 2007; 128(4):707-19. DOI: 10.1016/j.cell.2007.01.015. View

Weinberger L, Voichek Y, Tirosh I, Hornung G, Amit I, Barkai N . Expression noise and acetylation profiles distinguish HDAC functions. Mol Cell. 2012; 47(2):193-202. PMC: 3408861. DOI: 10.1016/j.molcel.2012.05.008. View

Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25(14):1754-60. PMC: 2705234. DOI: 10.1093/bioinformatics/btp324. View

Carrozza M, Li B, Florens L, Suganuma T, Swanson S, Lee K . Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005; 123(4):581-92. DOI: 10.1016/j.cell.2005.10.023. View

Stunnenberg H, Hirst M . The International Human Epigenome Consortium: A Blueprint for Scientific Collaboration and Discovery. Cell. 2016; 167(5):1145-1149. DOI: 10.1016/j.cell.2016.11.007. View

10.

Tarbell E, Liu T . HMMRATAC: a Hidden Markov ModeleR for ATAC-seq. Nucleic Acids Res. 2019; 47(16):e91. PMC: 6895260. DOI: 10.1093/nar/gkz533. View

11.

Skene P, Henikoff S . An efficient targeted nuclease strategy for high-resolution mapping of DNA binding sites. Elife. 2017; 6. PMC: 5310842. DOI: 10.7554/eLife.21856. View

12.

Boyle A, Guinney J, Crawford G, Furey T . F-Seq: a feature density estimator for high-throughput sequence tags. Bioinformatics. 2008; 24(21):2537-8. PMC: 2732284. DOI: 10.1093/bioinformatics/btn480. View

13.

Ross-Innes C, Stark R, Teschendorff A, Holmes K, Ali H, Dunning M . Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012; 481(7381):389-93. PMC: 3272464. DOI: 10.1038/nature10730. View

14.

Meuleman W, Muratov A, Rynes E, Halow J, Lee K, Bates D . Index and biological spectrum of human DNase I hypersensitive sites. Nature. 2020; 584(7820):244-251. PMC: 7422677. DOI: 10.1038/s41586-020-2559-3. View

15.

Gilmour D, Lis J . In vivo interactions of RNA polymerase II with genes of Drosophila melanogaster. Mol Cell Biol. 1985; 5(8):2009-18. PMC: 366919. DOI: 10.1128/mcb.5.8.2009-2018.1985. View

16.

Yu F, Sankaran V, Yuan G . CUT&RUNTools 2.0: a pipeline for single-cell and bulk-level CUT&RUN and CUT&Tag data analysis. Bioinformatics. 2021; 38(1):252-254. PMC: 8696090. DOI: 10.1093/bioinformatics/btab507. View

17.

Meers M, Tenenbaum D, Henikoff S . Peak calling by Sparse Enrichment Analysis for CUT&RUN chromatin profiling. Epigenetics Chromatin. 2019; 12(1):42. PMC: 6624997. DOI: 10.1186/s13072-019-0287-4. View

18.

Anastasiadi D, Esteve-Codina A, Piferrer F . Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species. Epigenetics Chromatin. 2018; 11(1):37. PMC: 6025724. DOI: 10.1186/s13072-018-0205-1. View

19.

Albert I, Mavrich T, Tomsho L, Qi J, Zanton S, Schuster S . Translational and rotational settings of H2A.Z nucleosomes across the Saccharomyces cerevisiae genome. Nature. 2007; 446(7135):572-6. DOI: 10.1038/nature05632. View

20.

Sethi A, Gu M, Gumusgoz E, Chan L, Yan K, Rozowsky J . Supervised enhancer prediction with epigenetic pattern recognition and targeted validation. Nat Methods. 2020; 17(8):807-814. PMC: 8073243. DOI: 10.1038/s41592-020-0907-8. View