This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve RefSeq protein accession from transcript accession

Hi,

I'd like to know if anyone knows how to retrieve RefSeq protein accession # from its mRNA transcript accession # using NCBI E-utility tool?

For example: Using NM_001382556.2 to get NP_001369485.1

Thanks very much!

ncbi refseq eutils
• 2,669 views • link
updated by GenoMax GenoMax • written by speycast • speycast

1 answer

Using Entrezdirect:

$ esearch -db nuccore -query "NM_001382556" | elink -target protein | efetch -format acc
NP_001369485.1
• 0 views • link
written by GenoMax GenoMax

genomax would you happen to know how to retrieve the reference sequence genbank file for this gene (RETL1) containing the above transcript NM_001382556.2? I want to use -format gbwithparts to get the mRNA and CDS region. There are other isoform transcripts in the genbank file, how can I just get the gbwithparts with FEATURES section containing just NM_001382556.2 and CDS regions?

Thanks again in advance...

• 0 views • link
written by speycast • speycast

You should only need:

$ efetch -db nuccore  -id "NM_001382556.2" -format gb

From reference genbank file as in chromosome/genome?

• 0 views • link
written by GenoMax GenoMax

Yes as in chromosome/genome. efetch -db nuccore -id "NM_001382556.2" -format gb gives the gb in mRNA version. I would like the genome genbank format for this gene, but since it contains two isoform transcripts in the FEATURES section of full genbank file, I only want the full genbank with transcript of interest like mRNA and CDS in the FEATURES section. (like a truncated FEATURES section) with just one transcript: NM_001382556.2 along with its CDS

• 0 views • link
written by speycast • speycast

Can you tell me which specific genbank record you are looking at?

• 0 views • link
written by GenoMax GenoMax

For example the DMD full genbank record (gbwithparts) GRCh37 assembly, here's the link: https://www.ncbi.nlm.nih.gov/nuccore/NC_000023.10?report=genbank&from=31137345&to=33357726&strand=true

In the FEATURES section, there are 30 mRNA transcript_ids (first one begins with mRNA join(1..351...) and follows that is its protein and CDS section. I'm interested in getting only mRNA transcript_id NM_004006.2 and its protein_id NP_003997.1 and its CDS region for FEATURES section and all other sections the same.

• 0 views • link
written by speycast • speycast

Best I can think of is this but that is going to give you all transcripts in that range. 

$ efetch -db nuccore -id NC_000023.10 -seq_start 31137345 -seq_stop 33357726 -format gb -style withparts
• 0 views • link
written by GenoMax GenoMax

Yup, thanks very much genomax! This is super helpful already, greatly appreciate it. I guess in this case I will just use python Bio package to parse the transcript of interest. :)

• 0 views • link
written by speycast • speycast

Log in to answer this question.

More posts like this