When I run a GSEA analysis on two conditions from the same RNaseq (negative control PBS injection VS positive control CpG injection) from the same dataset/same gene list, I get results that look something like this:

Example GSEA Image

Notice in my example that many gene sets are significantly enriched in PBS VS CpG, but none are very significantly enriched in the inverse, CpG VS PBS.

I have a basic question:

If we're comparing two items, say A and B, shouldn't a up-regulated gene set that's enriched in A have a corresponding down-regulated gene set that's enriched in B?

I'm consistently confused about how to interpret the fact that one of my two conditions has many more enriched gene sets than the other. What am I missing/misunderstanding?



If you compare A vs B the genes's fold change will have the opposite sign to B vs A. So will be the gene set up or down-regulated depending on the comparison to take.

The gene set test analyze if a given group of elements is sorted in a certain way in a list (I am talking about the GSE like the one you performed or the one in the Broad Institute). Usually it is used to say that gene set X is up-regulated in the comparison A vs B, that is, gene set X (in general) is more expressed in A than in B. The enrichment is measured in a enrichment score (or a normalized enrichment score, the higher, the clearer trend/distribution group X has. To asses X is really differentially expressed a p-value is also calculated.

Despite that the number of the up-regulated gene sets sum up all the gene sets tested as pointed by juod, the number of gene sets with a a p-value below 5% and 1% is not the same in both comparisons which indicates that you are not really using the same samples in both gene sets.

Check how you performed your comparison. Are the genes in each comparison with the same value and opposite sign?
Check also the software used to run the GSEA analysis (if I recongize the type of report is from a package in Bioconductor [I don't recognize which one now]) as you (or me) might have misunderstood something or there is an error if the comparisons are done properly. For instance it could be that the report of p-values below the threshold would be from only the up-regulated gene sets.

  • $\begingroup$ Thanks for the answer! I'm still having some trouble, most likely because I'm so new to this. This came from RNAseq data from the same experiment, but there are 3 "PBS-Treated" samples and 3 "CpG-Treated" samples. I'd expect lots of upregulation in the CpG treated samples. In the GSE results, one of the top PBS groups is 'GO_NEGATIVE_REGULATION_OF_B_CELL_ACTIVATION'. This means that particular gene set is enriched in PBS and not in CpG, right? Wouldn't I expect to see B cell activation in the CpG sample, then? $\endgroup$ Sep 12 '17 at 20:29
  • $\begingroup$ It might make it easier to explain the biological answer I'm looking for. I'd want to be able to say, "When we injected with CpG (as opposed to PBS control), the following gene sets were upregulated:". Is it just that there are no such gene sets? $\endgroup$ Sep 12 '17 at 20:32
  • $\begingroup$ @JulianStanley If you use GO as gene sets you don't use the relationship and the structure of the gene ontologies, (use topGO Bioconductor packages that take that into account). I think the GSE is not correctly performed: did you check what I said? But yes GSE will answer this kind of questions when done correctly. See my edit too $\endgroup$
    – llrs
    Sep 13 '17 at 7:05
  • $\begingroup$ @Llopis I think the key here is to find out which test exactly is performed in that package. Your description sounds like Mann-Whitney or something similar, which I agree should be symmetric. However, I am more used to enrichment testing by Fisher's test (i.e. contingency table "belongs to set-does not" vs. "is upregulated-is not"). Maybe the latter approach doesn't necessarily give symmetric p-values?.. $\endgroup$
    – juod
    Sep 13 '17 at 10:45
  • $\begingroup$ @juod I agree that knowing which test OP used is important. The Mann-Whitney is different from GSEA in that you don't have up/down-regulated gene sets, but over represented (OR) gene sets. The p-values if it's a OR test won't be symmetric, if it is a GSEA they should. $\endgroup$
    – llrs
    Sep 13 '17 at 10:54

While I'm not familiar with GSEA software in particular, I believe your problem is that it only tests for upregulated gene sets. Notice:

  • 866/4408 gene sets are upregulated
  • 3542/4408 gene sets are upregulated

3542+866 = 4408. I.e., 866 sets have higher mean expression in positive condition, the rest have higher mean expression in negative condition. To avoid this confusion, I personally would reserve using the terms "up/down-regulated" only for sets where this change is at least remotely significant.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.