Question: Gene Set Enrichment Test Tool That Can Compare Multiple Gene Lists
5
gravatar for confusedious
5.8 years ago by
confusedious420
Australia
confusedious420 wrote:

I received some excellent recommendations of enrichment test tools for testing single gene sets here at Biostars in the past. I am hoping someone will have something similar to offer for my latest query.

I have three gene lists (~1600 genes) that are not ranked in any particular way and have no expression data (they were obtained based on the Fst of SNPs within them - I have dealt with linkage disequilibrium to the best extent I am able, so not too worried about this). I would like to compare these three gene lists for enrichment of biological processes, in particular I am interested to see if they are enriched for immunological function.

I have become familiar with DAVID, but the issue I am finding is that I have to run the three tests separately and the enrichment scores obtained for gene functional classification seem to be only internally meaningful (that is to say comparing one gene list to another isn't a simple matter of side by side comparison in terms of enrichment scores).

Does anyone know of a test that can run all three lists side by side and statistically compare them?

pathway enrichment go • 9.1k views
ADD COMMENTlink modified 5.0 years ago by azimzadeh.jam.sadegh0 • written 5.8 years ago by confusedious420
1

there is no tool that I know of. But if you are familiar with any programming language, then it is an easy task.

ADD REPLYlink written 5.8 years ago by Bharat Iyengar260

Truly? I have a novice level grasp of R - are you aware of any way this could be done simply here? I'd be very interested to hear, as this step of my research has been a bit of a bottleneck on my workflow.

ADD REPLYlink modified 5.8 years ago • written 5.8 years ago by confusedious420

Actually, as PANTHER is calling over representation from set gene lists, that may work well for my purposes. I can compare these tests side by side by p-value - I think this may do the trick.

ADD REPLYlink written 5.8 years ago by confusedious420

Dear confusedious;

I need the same thing but for several hundred tests. Can you share your current experiences? Regards

ADD REPLYlink written 5.0 years ago by azimzadeh.jam.sadegh0

You should ask this as a separate, new question, otherwise it will be buried deep into another question where very few people will see it.

ADD REPLYlink written 5.0 years ago by Istvan Albert ♦♦ 80k
3
gravatar for David Westergaard
5.8 years ago by
Copenhagen, Denmark
David Westergaard1.4k wrote:

If you are familiar with R, you could use the Bioconductor package ReactomePA (Reactome Pathway Analysis). Alternatively, you could also just use the Reactome web interface and choose an overrepresentation analysis. Both use a hypergeometric statistic to test for statistical overrepresentation. If by "Biological Processes" you are refering to the Gene Ontology category, there is also the package topGO, which offers a few other statistical tests rather than just the Hypergeometric test. For all of these methods it should be possible to save the output into a tabular format for easy comparison. If you are familiar with ggplot2, you could also produce a dot chart where dots size is proportional to -p value, to make it even easier to compare.

ADD COMMENTlink written 5.8 years ago by David Westergaard1.4k

This is straight up and down fantastic. Thank you for the recommendation. The web interface is excellent and the output it produces is straight forward to read and also lumps together the immune system in a way that works for me. I'll give the R package a try over the weekend.

ADD REPLYlink written 5.8 years ago by confusedious420
1
gravatar for Woa
5.8 years ago by
Woa2.7k
United States
Woa2.7k wrote:

Ingenuity pathway analysis and (probably Metacore GeneGO too) has the option of loading multiple gene lists for pathway enrichment analysis. They are not free, however you can download trial versions that lasts from one/two weeks.

Though you're not using gene expression data as such, but have a look at this paper if it comes handy somehow:

http://nar.oxfordjournals.org/content/38/17/e169.abstract http://systems.crump.ucla.edu/rankrank/

ADD COMMENTlink written 5.8 years ago by Woa2.7k

I have access to Metacore in my lab - I haven't tried it as of yet so I was unaware of this property. Thank you for the input, I'll take a look.

ADD REPLYlink written 5.8 years ago by confusedious420
0
gravatar for azimzadeh.jam.sadegh
5.0 years ago by
Iran, Islamic Republic Of
azimzadeh.jam.sadegh0 wrote:

Dear confusedious;

I need the same thing but for several hundred tests. Can you share your current experiences? Regards

ADD COMMENTlink written 5.0 years ago by azimzadeh.jam.sadegh0
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: 1256 users visited in the last hour