GRanges object from gene list
1
0
Entering edit mode
rbronste ▴ 60
@rbronste-12189
Last seen 5.0 years ago

Wondering about the most straightforward way to take a list of genes and create a GRanges object, so basically pull the intervals for lets say a list of 500 RefSeq genes? Thanks!

granges bed genomicranges genomeintervals deseq2 • 1.6k views
ADD COMMENT
1
Entering edit mode
@michael-lawrence-3846
Last seen 2.9 years ago
United States

It depends on the source. For RefSeq in particular, you could use GenomicFeatures::makeTxDbFromUCSC() and ask for the "refGene" table (I think that's what it's called). Then, extract the desired features as a GRanges with e.g. genes(). There are probably resources available from AnnotationHub, as well.

ADD COMMENT
0
Entering edit mode

That would we good for refSeq for sure, how would you approach it if just using official gene symbols, such as those in DESeq2 results()?

ADD REPLY
1
Entering edit mode

You'll want to use whatever annotations were used to generate the counts. 

ADD REPLY

Login before adding your answer.

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