Tool stats & trends
Looking to identify usage trends or leading experts?
|Alternative names||Model-based Analysis for ChIP-Seq, MACS2, macs2|
|Interface||Command line interface|
|Restrictions to use||None|
|Input data||A tag file, a treatment file.|
|Output data||A file with information about called peaks, peak locations, peak summits locations, negative peaks, a script to produce a PDF image, a file that can be viewed through the UCSC genome browser, a diagnosis report and an optional file for the subpeaks option.|
|Programming languages||C, Python|
|License||BSD 3-clause “New” or “Revised” License|
Numpy, GCC, Cython
No version available
- Issues: https://github.com/taoliu/MACS/
- person_outline Wei Li
- person_outline Tao Liu
Publication for Model-based Analysis for ChIP-Seq
Distinct epigenetic landscapes underlie the pathobiology of pancreatic cancer subtypes
[…] alyzed using the HiChIP pipeline. Briefly, paired-end reads were mapped by BWA and pairs with one or both ends uniquely mapped were retained. H3K4me3, H3K4me1, and H3K27ac peaks were called using the MACS2 software package at false discovery rate (FDR) ≤ 1%. SICER was used to identify enriched domains for H3K27me3 and H3K9me3. For data visualization, BEDTools in combination with in-house scripts w […]
Methionine metabolism influences genomic architecture and gene expression through H3K4me3 peak width
[…] th were computed on this combinational peak set using the extended reads coverage files (that is, number of fragments extended from ChIP-seq reads mapped to each base pair on the genome) generated by MACS2. Height was computed by searching for position with highest read coverage. Area was computed by integrating the reads coverage file over the peak region. Width was quantified by first normalizin […]
Comprehensive epigenetic landscape of rheumatoid arthritis fibroblast like synoviocytes
[…] OA samples, each with six histone modifications (H3K4me1, H3K4me3, H3K9me3, H3K27ac, H3K27me3, and H3K36me3) and input data, were mapped to hg19 genome using bwa (0.7.7-r441). Peaks were called using MACS2 (184.108.40.20650420.1) using a p-value cutoff of 0.1 for broad marks (9me3, 27me3, 36me3) and 0.01 for narrow marks (4me1, 4me3, 27ac). To assess the library complexity and the enrichment of the dat […]
PCGF5 is required for neural differentiation of embryonic stem cells
[…] hen, the aligned reads were converted to bam format using samtools and duplicates were removed by Picard (http://broadinstitute.github.io/picard). Finally, the histone-enriched regions were called by MACS2 with default parameters. For PCGF5 ChIP-seq, enriched peaks were called by MACS2 with p-value <10−4 as cutoff, then peaks with q-value less than 0.01 were chosen for further analysis. Peaks were […]
TGF β induces miR 100 and miR 125b but blocks let 7a through LIN28B controlling PDAC progression
[…] ocus (see Fig. ) with the aim of disrupting part or the entire miRNA locus. For deletion of MIR100HG promoter region (MIR100HG∆P), pairs of sgRNAs were chosen to remove the SMAD2/3 peaks predicted by MACS2 from our ChIP-seq experiments. Finally, to generate LIN28B KO clones a single sgRNA targeting the genomic region downstream of the AUG translation start site codon was used. Oligonucleotides con […]
NOTCH mediated non cell autonomous regulation of chromatin structure during senescence
[…] d reads mapping to the ‘blacklisted’ regions identified by ENCODE were further removed. Average fragment size was determined using the ChIPQC Bioconductor package, and peak calling was performed with MACS2 (version 2.1.0), using fragment size as an extension size (--extsize) parameter. High-confidence peak sets for each condition were identified separately using only those peak regions that were p […]
Looking to check out a full list of citations?
I think it's important to understand that it is designed for point-source data, and for that it works perfectly (that's why I gave maximum rating, despite the constraints described in the following section, which I don't see as a fault - it's really not intended for the tasks below.) Numerous tests show (even a paper published by the MACS authors) that when it comes to diffuse enrichments, most notably, ChIP-seq with histone marks, then MACS fails to identify the peaks correctly. Usually it manages to detect some parts of it, some local summits within the dispersed signal, if it detects something at all. Although in the literature it is used for various datasets, including histone marks, ATAC-seq, MeDIP-seq etc., I wouldn't recommend it for detecting other enrichments than point-source peaks. There are much superior tools out there for diffuse signals (like SICER). I have to mention though that maybe H3K4me3 is an exception, it is a very "TF-like" histone mark, with sharp, clear, high enrichments, and for that mark only MACS might be acceptable. But I wouldn't trust it for any other HMs (and I've been working with ChIP-seq data of all histone marks you can imagine for ~10 years). I would also take the results of MeDIP-seq. ATAC-seq and similar data with a grain of salt if it was processed with MACS, because while it can certainly detect some enrichment, not all ofthem can be considered point-source, many methylated regions, open chromatin etc. are way broader than that.
With that said, I think it's a great tool for what it was designed.
Be aware that the MACS website is still saying that the newest version is version 1.4.2, which is not true, it has long been superceded by the major version 2, which you can find on GitHub (at the time of writing this review the latest stable version is 2.1.0).
Easy to install and use.