Peak caller for RIP-seq analysis - required or not?
3
1
Entering edit mode
3.6 years ago
priya.dalvi ▴ 10

Hello everyone! I am a biologist trying to analyse my generated RIP-seq data (+/- treatment, pulldown with histone biding proteins). I have performed genome mapping using HiSat2 and not sure which package to use for Peak calling. I tried Piranha, but the output from it doesn't seem right. I tried to find published work but there are not many people who have performed RIP-seq. Any suggestions regarding the peak caller would be really helpful. Also what would be the appropriate step after peak calling? Just peak annotation and motif analysis? Or no peak calling is needed. In this case how could I proceed further. I would highly appreciate any help in this direction. Thank you very much in advance.

RNA-Seq sequencing • 2.6k views
0
Entering edit mode

Thanks for your comment. I tried RIPSeeker, but it a bit complicated for me as a beginner to compute it,. As is it is RNA-IP I assume there are many split reads so I performed Hisat2 mapping. But, I am still unsure of MACS2 after Hisat2

0
Entering edit mode
3.6 years ago

Not really sure, but I'll be facing this problem soon enough too.

Try MACS2 for peak calling as one of the classics. There are a few bioconductor R packages on RIP-seq. They should give you an overview on further steps.

https://bioconductor.org/packages/release/bioc/html/RIPSeeker.html

Was HiSat2 necessary, i.e. did you see many split reads ? What was your read length? Also, just manually converting your BAMs to bigwig and looking at a few prominent peaks will get you a long way.

0
Entering edit mode

Hi, I used STAR instead Hisat2. Is it Correct?

0
Entering edit mode
3.1 years ago

Hi, Priya. Have solved the problem? I have RIPseq data to analysis. I tried RIPSkeer and MACS, but the results seem not good. Could you tell me what software you use and the parameters? Thank you!

1
Entering edit mode

RIPSkeer doesen't work in the latest R and Bioconductor Versions. How Can I solve this? Thanks

0
Entering edit mode
3.1 years ago

I used macs and the parameters are: macs2 callpeak -t ${sample[0]}.sorted.bam -c${control[0]}.sorted.bam -f BAM -g dm -n ${sample[0]} -B -q 0.01 --nomodel --nolambda --extsize 200 >${sample[0]}.log

But the results are weird. There are almost 15000 peaks. My data came from fly.

0
Entering edit mode

Hi xiaocong3333, have you found a solution? I too have RIP-seq data of the fly model. Looking for a peak caller for motif analysis.

Traffic: 1410 users visited in the last hour
FAQ
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.