Question: How to annotate gene list using enrichr?
gravatar for nattzy94
3 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 • 162 views
ADD COMMENTlink modified 10 weeks ago by Kevin Blighe56k • written 3 months ago by nattzy9420

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

ADD REPLYlink written 3 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 3 months ago by nattzy9420
gravatar for Kevin Blighe
10 weeks ago by
Kevin Blighe56k
Kevin Blighe56k 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 enreichment:

  • minimum genes per enrichment term
  • p-value


ADD COMMENTlink written 10 weeks ago by Kevin Blighe56k
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: 1747 users visited in the last hour