This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to retrieve nucleotid sequence from gene ids of ncbis "gene" data base?

Hello people

I would like to retrieve all sequence from a set of gene entrys of the NCBI data base "gene".

As an example I would like to retrieve all sequence of this query:

"txid511145[Organism:noexp] "

URL: https://www.ncbi.nlm.nih.gov/gene/?term=txid511145%5BOrganism%3Anoexp%5D

The only way I found so far is to download the the full genome to which the genes refer and grep all the sequence locally according to the length and starting position. Is there a better way?

Thanks

gene id ncbi
• 5,449 views • link
updated by Joseph Hughes Joseph Hughes • written by john john

I do actually, but cant figure out how to.

• 0 views • link
written by john john

in my question is the query and also the url

• 0 views • link
written by john john

I meant esearch/eutils query.. Check the link and try to build a eUtils query.

• 0 views • link
written by Santosh Anand Santosh Anand

Sorry, my bad! What I meant was to use Entrez Direct command line tools: https://www.ncbi.nlm.nih.gov/books/NBK179288/

You can make complicated queries using that and can chain queries where the results from one query are fed to the next one. See if that helps. Also a youtube video is here

• 0 views • link
written by Santosh Anand Santosh Anand

Oh thats a nice one. I didnt found that. Until now I used the r package rentrez for the calls. But it seems this one is more powerful, maybe. But unfortunately I still do not get what I want. Here is my query

esearch -db gene -query 'txid511145[Organism:noexp]' | efetch -format fasta

But this again just returns me the entries of the gene db. Example:

ID: 945651 99. dnaC DNA biosynthesis protein [Escherichia coli str. K-12 substr. MG1655] Other Aliases: b4361, ECK4351, JW4325, dnaD Annotation: NC_000913.3 (4600238..4600975, complement)

I could use the last line to get the corresponding fasta locally. But I would like to know if the server of ncbi I would do this for me or not.

• 0 views • link
written by john john

1 answer

I think something like this using the elink function of eutilities should work:

esearch -db gene -query 'txid511145[Organism:noexp]' | elink -target nuccore | efetch -format fasta

As you know what the accession number of you genome is, you are much better starting from that. The following retrieves all coding sequences for the reference genome

esearch -db nuccore -query 'NC_000913.3' | efetch -format fasta_cds_na

This gives a total of 4319 sequences. Accession number NC_000913.2 is an older version of the accession number.

• 0 views • link
written by Joseph Hughes Joseph Hughes

Oh that looked so good. But the result is really not what I hoped it would be. The returned fasta just has 38 entries and this e coli strain should have 4516 genes. Also one of the entries is the whole genome. Not really know what this results refere to any how. As all the genes map only to two entries in the nucore db "NC_000913.3" and "NC_000913.2".

• 0 views • link
written by john john

Does your starting point have to be the taxid? The problem with starting with a taxid is that it is not very precise. It sounds like you know the two full reference genomes that you want to extract genes from so why not start from those accession numbers?

• 0 views • link
written by Joseph Hughes Joseph Hughes

That works for me. Yeah the starting point is the organism so the txid. But thats okay. So I check for the best genome and work with this further. Thanks!

• 0 views • link
written by john john

Log in to answer this question.

More posts like this