Question: Thousands of significant differentially expressed genes but no significant enriched pathways by GSEA?
0
gravatar for 88adavis
24 days ago by
88adavis0
88adavis0 wrote:

Hello everyone,

I've recently performed a GSEA analysis (using fGSEA in R) on an RNA-seq data set, consisting of 4 pairwise comparisons: 1.) drug 1 vs. untreated, 2.) drug 2 vs. untreated, and 3.) drug 1 + 2 vs. untreated).

Using edgeR, I found that there were ~2000 differentially expressed genes (DEGs) for drug 1 vs. untreated, ~7000 DEGs for drug 2 vs. untreated, and >9000 genes for drug 1+2 vs. untreated. After ranking the full list of genes in each of the 3 signatures (by the following formula: -log10(adjusted p-value)/sign(logFC)), I performed GSEA using fGSEA (v. 1.10.1) in R, against the Human GO and Reactome databases (~3700 pathways). After applying a Benjamini-Hochberg p-value adjustment, I obtained the following number of ‘significant’ gene sets/pathways (adj. p < 0.05): 1.) 70 for drug 1 vs. no tx, 2.) 0 for drug 2 vs. no tx and 3.) 106 for drug 1 + 2.

I find it strange that there aren’t any ‘significant’ gene sets obtained for the ‘drug 2 vs. no tx’ signature, given the large number of DEGs found. In fact, this signature shared more than 5000 significant DEGs (with same sign(logFC) direction) with the combo drug 1 + 2 treatment.

I have 2 questions: 1.) Is it possible to have no significant pathways despite having such a large fraction of the transcriptome being differentially expressed? If so, what conditions would account for that? 2.) I ranked my genes using -log10(adjusted p-value)/sign(logFC) (instead of raw p-value, as noted here), as I only had access to logFC and B-H adjusted p-values for each of the gene signatures. Is it appropriate to rank on adj. p-value? Could this account for the few enriched pathways observed for ‘drug 2 vs. no tx’?

Thanks, Andrew

ADD COMMENTlink modified 16 days ago by alserg300 • written 24 days ago by 88adavis0
2
gravatar for colin.kern
24 days ago by
colin.kern750
United States
colin.kern750 wrote:

1.) Is it possible to have no significant pathways despite having such a large fraction of the transcriptome being differentially expressed? If so, what conditions would account for that?

Yes, absolutely. Let's say, for easy math, that your species has 1000 genes and contains a pathway consisting of 10 of those genes. You do an experiment and find 100 DEGs. If you were to randomly choose 100 genes, you would expect 1 of these on average to be part of that pathway. So if your DEGs include, say, 7 of those 10 genes, that's clearly significant. Even 3 or 4 might be good enough for significance. But if you have 500 DEGs, now a random sampling of 500 genes is going to have 5 of the 10 genes from that pathway in it. It's much harder to get a significant enrichment for that pathway now. If your DEG set contains 6 or 7 of those genes, that might not be significant.

You could try a more stringent significance cut-off for your DEGs, or look for something like a batch effect or other technical bias that could be causing so many DEGs. Do you expect the drugs you're testing to have such a huge effect on gene expression?

ADD COMMENTlink written 24 days ago by colin.kern750
1
gravatar for kristoffer.vittingseerup
24 days ago by
European Union
kristoffer.vittingseerup2.9k wrote:

Did you try a functional-class scorring gene-set analysis? Aka one that does not rely on the cutoff on you previous DE analysis. Since you are using edgeR you can do that directly with the camera() function (technically it is from limma - but edgeR loads limma so it is redily available and accepts edgeR objects).

ADD COMMENTlink written 24 days ago by kristoffer.vittingseerup2.9k
0
gravatar for alserg
16 days ago by
alserg300
alserg300 wrote:

What value of nperm did you use for fgsea? Low values, like 1000, or sometimes even 10000 will limit the minial possible nominal P-value, which after BH-correction can lead to non-significant results.

ADD COMMENTlink written 16 days ago by alserg300
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 2.3.0
Traffic: 1002 users visited in the last hour