Comparing differentially expressed peaks to get unique values
Entering edit mode
20 months ago
rykerklie7 ▴ 20

I have done differential expression of the peaks from ATAC-seq data using DESeq2.

I have A vs B (Let us call this X), B vs C (Let us call this Y), and C vs D (Let us call this Z). I would like to compare X - (Y and Z) to get the unique peaks that are differentially expressed only in X.

I have the data in the following format

     Peak_name  Chr Start  End Width

I tried something as simple as the function merge in R but not sure if it is the right way to do it. And I am not sure if the overlap functions in diffbind works for this. I tried

      list1=merge(Y,Z, by.x = 0, by.y =0)
      list2=merge(X,list1, by.x=0, by.y=0)
ATAC-seq diffbind ChIP-Seq Genomic Ranges • 490 views
Entering edit mode
16 months ago
Rory Stark ★ 1.6k

The overlap (dba.overlap()) and Venn diagram (dba.plotVenn()) functions in DiffBind can handle these types of queries, but in R it is probably simpler to turn your peaksets into GRanges objects and use the built-in overlapping functions.


Login before adding your answer.

Traffic: 1311 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6