This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Convert Gene Name to RefSeq ID

Dear all,

I have seen this post that allows gene conversion from RefSeq IDs to gene names. What I would like is a tool (command line or web-based) that:

  1. takes as input a list of gene names (shortcut name or comprehensive gene name) AND a query species
  2. outputs the list of RefSeq IDs

In my case

  1. the species would be Stegastes partitus

  2. the gene names

    LOC103370819
prtfdc1
LOC103367872
tfec
colony stimulating factor 1 receptor (csf1r)

I first thought of using blastdbcmd but it seems that blastdbcmd does not take gene names as input.

$blastdbcmd -db nt_Spar -dbtype nucl -entry prtfdc1 #nt_Spar a subset of nt with only S. partitus sequences
Error: prtfdc1: OID not found

I have tried using Batch Entrez but it does not accept gene names as input neither.

Many thanks for your help!

rna-seq ncbi sequence blast
• 5,393 views • link
updated by Sej Modha Sej Modha • written by tlorin tlorin

I didn't know this tool! But I cannot make it work B-)

• 0 views • link
written by tlorin tlorin

1 answer

You can get accession numbers for those genes by using NCBI eUtils. Here is an example: esearch -db nuccore -query "prtfdc1 [Gene] AND Stegastes partitus [ORGN]" | efetch -format docsum | xtract -pattern Caption -element Caption This produces XM_008294096 NW_007578669 You would want the NW* numbers. In that case add a pipe to grep NW* at the end of the command above.

• 0 views • link
written by GenoMax GenoMax

Didn't work at the beginning: all the commands (esearch,efetch,xtract) need to be in the path (obviously). Works perfectly now, thanks!

• 0 views • link
written by tlorin tlorin

@genomax: how would you do with a complete gene name instead of the shortcut? For instance Stegastes partitus phosphoribosyl transferase domain containing 1.

This command does not output anything: ./esearch -db nuccore -query "Stegastes partitus phosphoribosyl transferase domain [Gene] AND Stegastes partitus [ORGN]" | ./efetch -format docsum | ./xtract -pattern Caption -element Caption

• 0 views • link
written by tlorin tlorin 
esearch -db nuccore -query "Stegastes partitus phosphoribosyl transferase domain containing 1 AND Stegastes partitus [ORGN]" | efetch -format docsum | xtract -pattern Caption -element Caption | grep NW

NW_007577984
NW_007578669

While that seems to generate a result those accessions are for the genomic entries. Guess you may not be able to make some of them work.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this