This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Map a list of gene names in UniProt

Hi all, I have a long list of gene names (>7000) which is coming from more than 50 organisms. I am not worried about where it is coming from, however, I need to get the most accurate GO terms or at least UniProtKB of these genes so that I can categorize them. The ID mapping of UniProt requires the organism name and this is taking huge time and data overlap. My aim is to be able to categorize these genes based on GO terms and pathways they are associated to. I would appreciate any suggestion. Thank you.

mapping uniprot id
• 3,103 views • link
updated by jv jv • written by Space_Life • Space_Life

Thank you for replying. The genes are from UniProt. I tried its REST API but it just gets stuck with no output.

<import urllib.parse
import urllib.request

url = 'https://www.uniprot.org/uploadlists/'

params = {
'from': 'GENENAME',    
'to': 'ACC',
'format': 'tab',
'columns': 'id,entry_name,reviewed',
'query': 'tsf, pyrH, frr, ispH, btuD_1, rpoC, rpoB, rplL, rplJ, rplA, rplK',
'taxon': '391774'
}

data = urllib.parse.urlencode(params)
data = data.encode('utf-8')
req = urllib.request.Request(url, data)
with urllib.request.urlopen(req) as f:
   response = f.read()
print(response.decode('utf-8'))

I just tested this. If I remove 'taxon' from this code, the output is nothing. Thank you.

• 0 views • link
written by Space_Life • Space_Life

First, I think there is some confusion. Your gene names are not accession ids which is what I was inquiring about in my first question. What you have are generic gene symbols and the associated taxonomy but which should be enough to get you the info you want from UniProt.

I'm not quite sure why your above query did not return any results as I am not familiar with python urllib functions. But I suspect that your url is improperly formatted. For an example, the following url will produce the results that you desire for gene pyrH from taxonomy 391774

https://www.uniprot.org/uniprot/?query=gene:pyrH+and+taxonomy:391774&&columns=id,entry_name,reviewed&format=tab

You'll want to make sure that your code is producing a url that is formatted as shown above. To query for multiple genes, either loop this for each gene or see the below example:

https://www.uniprot.org/uniprot/?query=gene:(pyrH+or+rpoC+or+tsf)+and+taxonomy:391774&&columns=id,entry_name,reviewed&format=tab

• 0 views • link
written by jv jv

0 answers

No answers yet.

Log in to answer this question.

More posts like this