Getting gene ID's from camera results
2
0
Entering edit mode
AB ▴ 100
@ab-8975
Last seen 2.9 years ago
United States

Hi,

I have used camera (limma package) to do a gene set enrichment analysis with a gene set from MSigDB. I was wondering if there is any way I could use the results from camera to get the genes belonging to that list. 

                                                                                NGenes   Correlation       Direction     PValue        FDR
HALLMARK_E2F_TARGETS                                     201      -0.0047718432        Up    0.01676313  0.5229967
HALLMARK_PI3K_AKT_MTOR_SIGNALING              96     -0.0074777132      Down  0.03990746  0.5229967
HALLMARK_IL6_JAK_STAT3_SIGNALING                 68     -0.0109702506      Down  0.04448435  0.5229967
HALLMARK_INTERFERON_ALPHA_RESPONSE      62     -0.0131317542      Down  0.04721300  0.5229967

For instance, HALLMARK_E2F_TARGETS has 201 genes. How can I find out which genes from my input gene identifiers belong to this list ?

Thanks

limma camera • 1.2k views
ADD COMMENT
2
Entering edit mode
Aaron Lun ★ 27k
@alun
Last seen 1 minute ago
The city by the bay

Well, you don't really need the results from CAMERA, because your gene sets are already defined. Assuming your list of index vectors is indices and your matrix y has the gene IDs as the row names , you could do:

rownames(y)[indices$HALLMARK_E2F_TARGETS]

If that doesn't work, you'll have to give some more details about the inputs you're using to run camera.

ADD COMMENT
0
Entering edit mode

That did work. Thank you very much. But is there anyway to do it without knowing the indices ?

I have a shiny interface that loads up the precomputed camera results from an Rdata and I need the interface to display the gene list when each row is clicked on. I have the gene set and the camera results but not the indices.

 

ADD REPLY
0
Entering edit mode

You had to know the indices to do the camera analysis in the first place. If you want to re-access to them in the future, can you not simply store them in the RData file together with the camera output?

Obviously you can't list the genes in a set without knowing the indices because the indices define which genes are in the set.

ADD REPLY
0
Entering edit mode
phil.chapman ▴ 150
@philchapman-8324
Last seen 5.9 years ago
United Kingdom

Since these are MSigDB gene set collections you could just use the .gmt files directly or one of a number of bioconductor packages that can process these for you.  I use the piano package quite a lot which has a loadGSC function, so having downloaded, for example, the hallmarks .gmt file (h.all.v5.0.symbols.gmt)  from http://software.broadinstitute.org/gsea/msigdb you can just do:

library(piano)
my_gsc <- loadGSC('h.all.v5.0.symbols.gmt')
my_gsc
my_gsc$gsc[1:5]
my_gsc$gsc[['HALLMARK_E2F_TARGETS']]

 

 

 

ADD COMMENT

Login before adding your answer.

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