Thank you. I started testing with STAR because it is very fast. I read that STAR could be used for DNA aassys coz it can handle gaps. Plus, using option like --alignIntronMax 1 --alignMatesGapMax <maxinsertsize-2*readlength> splicing could be prohibited. Anywa, seems I need to do everything again with bwa alignment and macs2. Btw, the reads are not paired end. Any further comments and suggestions are appreciated.
ATACseq with STAR and macs2
I am new to ATACseq data. Earlier, I have processed mainly RNAseq and few Chip seq data, thus i used same tools for atacseq as well.
I used STAR (default parameters plus --outSAMstrandField intronMotif) to generate bam and wig files, then I used macs2 (-B -q 0.01) to call peak using those bam files.
Is this a correct way to process? Also, can you recommend a tool to get the differential peaks from this ATACseq data?
• 9,865 views
•
link
updated
by
simon.vanheeringen
•
written
by
Ice P •
•
written
by
Ice P •
2 answers
You typically use the --nomodel option, as the shifting model of MACS does not really make sense for open chromatin data. I personally (and currently, edit in 10/2019) use Genrich which has an ATAC-seq mode that. For macs you can also write each read as BED (to get the cutting events which is the 5' ends of each read) and then feed this into the peak caller.
Have you tried the pipeline from the Kundaje lab?
https://github.com/kundajelab/atac_dnase_pipelines
Pretty much worked out of the box for me. For differential peaks you could do the following:
- merge all peaks into one peak set
- count reads for all experiments in this merged peak set (using bedtools multicov for instance)
- use the count table in DESeq/EdgeR
• 0 views
•
link
written
by
simon.vanheeringen
Log in to answer this question.
More posts like this
-
I do not understand how to get total peaks from MACS2 shown to me. What mistake(s) am I making here?
written by Ronin •
-
MACS2 not working in paired end mode
written by Filago •
-
Chip-seq peak calling best way usage of parameters
written by luckysardar171 •
-
Integration of ATACseq, RNAseq data with ChIPseq data
written by kuchmiyanna •
-
STAR+Cufflinks for RNA-Seq analysis
written by wangdp123
-
[MACS2] False negative peak calling problem
written by dude.dude.law •
-
Different mapping rate for tophat+mm9 vs star+mm10
written by Ice P •
-
error in macs2 usage
written by Sara
-
macs2 callpeaks genome
written by sophialovechan •
-
How can I increase the number of peaks using macs2 peak caller?
written by GK1610
We use Bowtie for mapping, MACS for peak calling, IDR filtering of peaks and DESeq2 for differential peaks. Still work in progress but our current pipeline is on github.