Question

GRanges object from gene list

0

Entering edit mode

rbronste ▴ 60

@rbronste-12189

Last seen 4.4 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.4k views

ADD COMMENT • link 6.3 years ago rbronste ▴ 60

score 1 · Answer 1 · 2018-01-18

1

Entering edit mode

Michael Lawrence ★ 11k

@michael-lawrence-3846

Last seen 2.4 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 • link 6.3 years ago Michael Lawrence ★ 11k

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 • link 6.3 years ago rbronste ▴ 60

1

Entering edit mode

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

ADD REPLY • link 6.3 years ago Michael Lawrence ★ 11k