Question: ChIPseek online tool, errors?
0
gravatar for Samad
3.8 years ago by
Samad90
France
Samad90 wrote:

I have used online tool of analysis and annotation of ChIP-seq data "chipseek" (http://chipseek.cgu.edu.tw/ published here http://www.biomedcentral.com/1471-2164/15/539).
When I compare two lists of peaks to get VennDiag (ovelapping peaks), I get different results for the same list!

In the example bellow, from the same imported "Pif1_peaks.." list of peaks (.bed), I have a total of 1203 peaks in the 1st comparison A (338+865), however, in the second comparison the total number of peaks is 1245 (747+498) for Pif1_peaks..

In addition, in the original bed file "Pif1_..bed", I have only 1173 peaks!

Is there any explication for this?

Thank you!

                      A                                                             B



 

ADD COMMENTlink modified 3.8 years ago • written 3.8 years ago by Samad90
2
gravatar for ivivek_ngs
3.8 years ago by
ivivek_ngs4.8k
Seattle,WA, USA
ivivek_ngs4.8k wrote:

Can you please be more clear about the question? I can see the read portion for from the left in Fig A is different from that of the red in figure B. I also do not see (267+932) . Also you have to check how are the peak comparsion being done. Are they considering total end to end regions or a fraction.

ADD COMMENTlink written 3.8 years ago by ivivek_ngs4.8k

Sorry 267+932 was for another example of comparison (also different than total of initial list, and different of the two others), I updated the my post according the example in image.

The colors changes (can't control it for this online tool), but the list is the same

The comparison is done by fraction, wherever a peak ovelap another regardless of the overlapping size

Thank you for you reaction

 

 

 

ADD REPLYlink written 3.8 years ago by Samad90
1

I presume that you have bed file for both the peak files , and the tool uses interectbed and coveragebed to find the peak realtionships. However I would suggest you to do a simple one command line to check the peak relationship with intersectbed. The paper has not  shown anything about the criteria so I am assuming that the fraction they use in inetersectBed is 1.0 (100%) for overlap. If you can do that in unix you can find that easily or you can mail the developers. I am assuming they have a binning or splitting window strategy based on coverage for which you have more peaks than the usual and that is the reason you have this descreapancy. However for your testing you can use the below command if you have bedtools installed.

bedtools intersect -a file1.bed -b file2 .bed -u -f 1.0 | wc -l ## This will show the number of peaks that are overlapping with 100 % identity between regions

or

intersectBed -a file1.bed -b file2.bed -u -f 1.0 | wc -l

#Note : that you should have bedtools or bedtools2 to run the above commands.

ADD REPLYlink written 3.8 years ago by ivivek_ngs4.8k
1

Hi vchris, 

It's very useful your reply, acually I have used bedtools intersect and wanted to visualize results as Venn diagram so I can see the unique and overlapped peaks!

Your question in the previous post made me thinking & verifying in details from where they comes the additional peaks, so I realize that it's normal to have maybe more peaks in total because 1 peak in the 1st list can overlap (partially) two peaks in the 2nd, in this particular case we count 2 for the common area in the figure, this was confusing for me!

I hope it was clear my explanation

Thank you again, I appreciate

Samad

ADD REPLYlink written 3.8 years ago by Samad90
1

Yes that is obviously possible since you are overlapping region and if the overlapping fraction is not 100% percent then a a peak in A can obviously overlap with 2 or more peaks in B depending on the fraction of windows of overlap parameters applied. Good that it helps. Dont mind an upvote or an acceptance of answer if my previous queries were helpful! Good luck!

ADD REPLYlink written 3.8 years ago by ivivek_ngs4.8k
1
gravatar for Samad
3.8 years ago by
Samad90
France
Samad90 wrote:

The reply of @vchris_ngs made me thinking & verifying in details from where they comes the additional peaks, so I realize that it's normal to have maybe more peaks in total compared to the initial number, the reason is that 1 peak in the 1st list can possibly overlap (partially) two peaks in the 2nd, in this particular case we count 2 for the common area in the figure, this was confusing for me!

I hope it was clear my explanation and hope this post will be useful for others 

Thank you again vchris, I appreciate!

Samad

ADD COMMENTlink written 3.8 years ago by Samad90
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1579 users visited in the last hour