How to interpret WebGestalt plots for GSEA and ORA?
0
0
Entering edit mode
1 day ago
Haad • 0

Hi,

I am wondering how to interpret WebGestalt's ORA report. I did the ORA on only the significant genes from the full list (significance criteria: abs(logFC) > 1 and adjusted p-val < 0.05).

I want to know what the terms in this screenshot mean: enter image description here

I already know that enrichment ratio in this case is 14 / 0.61962 = 22.594. My questions are:

  1. What is an analyte set and how is it different from a gene set?
  2. What is the difference between "mapped input" and "gene set"?
  3. Is the expected value just the number of genes we expect to see in each analyte set?
  4. What is "overlap"? In other words, back to my second question, what is "mapped input"?

Any insight would be greatly appreciated.

I am also wondering what this highlighted sentence means in the WebGestalt job summary for ORA: enter image description here

EDIT: How is the expected value calculated for ORA here?

webgestalt gsea ora • 226 views
ADD COMMENT
0
Entering edit mode

I think it means that of the 116 IDs you supplied, 52 are in pathway_KEGG. Check out the number of entries in "Protein digestion and absorption" - that might be 103 in size.

EDIT: hsa04974 has 105 genes right now so probably, the database version that webgestalt uses could have 103. I know (because I used to work in Zhang lab) that webgestalt operates (or at least used to operate) on manually downloaded database files so that could explain the discrepancy.

ADD REPLY
0
Entering edit mode

See my screen shot with the Venn diagram above. The 14 overlapped genes were between the gene set for protein digestion and absorption (103 genes) and the "mapped input" (52 genes, out of 116 genes I supplied). My question is why is the size of mapped input not 116?

ADD REPLY
0
Entering edit mode

Are you sure that all 116 genes were found in KEGG?

ADD REPLY
0
Entering edit mode

Ah, I see what you mean. Out of 116 genes I supplied only 52 were found in KEGG. I supplied ensembl gene ids. I wonder why the number of genes found in KEGG is so low. Does this mean the results of ORA are as accurate as they would be compared to if all 116 genes were found in KEGG?

ADD REPLY

Login before adding your answer.

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