Question: How to annotate gene list using enrichr?
gravatar for nattzy94
10 months ago by
nattzy9420 wrote:

I am analysing a published rnaseq read counts table that is upper quartile normalised. I have correlated a transcript of interest (transcript A) against all other transcripts and obtained a 2 column table with transcript/gene names in the first column and Pearson correlation (ranging from -1 to 1) in the second column. I want to annotate the gene list so that I can get the gene function in a 3rd column. I've passed the list of genes into enrichr but the list of gene functions that it returns is not the same length as the list of genes that I passed in. Any ideas on how else I can do this?

rna-seq annotation gene • 271 views
ADD COMMENTlink modified 9 months ago by Kevin Blighe66k • written 10 months ago by nattzy9420

not all the genes have annotated functions? Can you give an example of what enrichr has skipped?

ADD REPLYlink written 10 months ago by b.bearmi0

It's skipped about 3000 genes. I passed ~8000 genes into queryMany using the GoBiologicalProcess db. The resulting GO_Biological_Process_2018$Genes only has 5000 entries.

Should I include the correlation values when I input the data into queryMany?

ADD REPLYlink written 10 months ago by nattzy9420
gravatar for Kevin Blighe
9 months ago by
Kevin Blighe66k
Kevin Blighe66k wrote:

There may be an issue with ID matching, or it may be that it simply could not enrich those 3000 'missing' genes. No need to input the correlation values, as it won't use them and may throw a syntax error.

Another thing to check is thresholds used for enrichment:

  • minimum genes per enrichment term
  • p-value


ADD COMMENTlink modified 4 months ago • written 9 months ago by Kevin Blighe66k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1369 users visited in the last hour