How to see if DEGs are grouped in genome?
1
0
Entering edit mode
6.1 years ago
Parham ★ 1.6k

Hi,

In my experiment I compromise replication and transcription, then look into transcriptome for changes. 

One thing that I would like to investigate is if the DEGs are significantly grouped in close-by loci through out genome or they are just scattered out in a random manner? 

I would like to have your feedback on how to approach this question and which packages do you recommend?

I prefer to work under R environment if that's possible. The organism that I work on is S. pombe. Let me know if you need more information. 

 

Thanks in advance,

Parham

DEG rna-seq overlap • 1.0k views
ADD COMMENT
1
Entering edit mode

Well, depends on how many genes you got. If you are using cufflinks, it tells you the genomic chromosome/scaffold as well as nucleotide range the DEG is located on, so if you look at the DEGs you may visually be able to tell how close they are.

ADD REPLY
0
Entering edit mode

I was previously using cufflinks but then I switched to STAR and DESeq2 for the purpose. I will look into that. Thanks!

ADD REPLY
0
Entering edit mode

Another idea is to export your annotation in BED, and then just write a quick grep script that will pull out the positions of the DEGs.

ADD REPLY
1
Entering edit mode
6.1 years ago
Christian ★ 3.0k
I ended up using the Perl script from de Preter et al. (2008), which worked quite nicely: Recommended tools for positional gene enrichment (PGE) analysis?
ADD COMMENT
0
Entering edit mode

Yeah that looks convenient. However they don't have fission yeast genome in their data base unfortunately. Thanks for recommending though!

ADD REPLY

Login before adding your answer.

Traffic: 1818 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