how to select peaks at not transcribed genes
0
0
Entering edit mode
20 days ago
Lila M ▴ 910

Hi there, Maybe this is a silly question. I have several peaks called with macss (n = 12808) and I want to create two groups based peaks at transcribed genes and peaks at not transcribed genes. My first approach is to annotate peaks with Chipseeker , which gaves me a total of n = 12576 peaks, so I can consider those peaks at transcribed genes. In order to identify the peaks at not trasncribed genes, my idea is use bedtools intersect -v and the remaining 232 peaks would be consider as not trasncribed. Is this a good approach? or do you guys have any better idea? Thanks!

bedtools transcribed peaks not_trasncribed • 329 views
0
Entering edit mode

Which type data are we talking about? What was the experiment? There are literally no details in this question right now. It the assay measuring transcription or is this a DNA-based genomic assay such as ChIP-seq?

0
Entering edit mode

rigth, it is a ChIP seq analysis

0
Entering edit mode

What makes you think genes with a peak are transcribed while those without are not? Unless your IP was for Pol2, this is an overly simplistic assumption that I very much doubt is true. Especially since ChIPseeker just assigns to the closest TSS, so some of these peaks will likely not even impact the gene they are assigned to.

0
Entering edit mode

indeed, my IP was for Pol2

0
Entering edit mode

In that case, wouldn't you consider all peaks to be associated with transcription? Even if your 232 peaks are not associated with genes, they may still mark transcription of unannotated genes, eRNAs, etc.

0
Entering edit mode

That is my concern... this is why I asked in here :)

1
Entering edit mode

Shouldn't you assess transcription by an actual RNA-based assaay such as CAGE or RNA-seq rather than guessing based on PolII?