The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: GRanges object from gene list
0
gravatar for rbronste
13 months ago by
rbronste60
rbronste60 wrote:

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!

ADD COMMENTlink modified 13 months ago • written 13 months ago by rbronste60
Answer: GRanges object from gene list
1
gravatar for Michael Lawrence
13 months ago by
United States
Michael Lawrence10k wrote:

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 COMMENTlink written 13 months ago by Michael Lawrence10k

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 REPLYlink written 13 months ago by rbronste60
1

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

ADD REPLYlink written 13 months ago by Michael Lawrence10k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 270 users visited in the last hour