eCLIP peak calling
2
1
Entering edit mode
3.2 years ago
GiusiG ▴ 10

Hi everyone,

I'm facing a problem calling peaks from a eCLIP experiment. I have 3 bam files: 2 replicates for one condition and a negative control. I don't know if this is the best option but I'm calling the peaks using MACS (from Galaxy), which is supposed to be for ChiP-seq and, for this reason, it doesn't give me the strand of the signal in the output file but I really need this information. Is there a way to tell MACS that I want the strand information to be kept? Or maybe is there a better way to call peaks for a CLIP experiment?

Once I have the bed files with the peaks I will use bedtools to extract the RNAs were I can find a peak for different factors.

Thank you in advance for your help

CLIP peak calling MACs Galaxy • 2.3k views
ADD COMMENT
2
Entering edit mode

We have had success with the CTK toolkits across various CLIP data. Installation can be a bit painful but I think it's worth the investment. https://zhanglab.c2b2.columbia.edu/index.php/CTK_Documentation

ADD REPLY
0
Entering edit mode
3.0 years ago

Hi, we developed a tool called PureCLIP, which calls crosslink sites (or high resolution binding regions) from such data. It is available on GitHub https://github.com/skrakau/PureCLIP and can be installed via Bioconda. There is a branch on GitHub to integrate replicates, which will be released within the next month.

ADD COMMENT
0
Entering edit mode
3.0 years ago
Ashley ▴ 90

Hi, I think Piranha is also a good strategy for analyzing CLIP-seq data. You can download from: http://smithlabresearch.org/software/piranha/ Thanks. Ashley

ADD COMMENT

Login before adding your answer.

Traffic: 1737 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6