This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve Gene Sequences Via Programming From Gene Databanks

Hi dears!

I have a local database of miRNA::gene interactions for hsa, which contains just "GeneID" (Entrez ID) and "Gene Name" for Genes. I want to extract gene sequences for my genes either with the name or with the ID.

I have two different ideas:

1)I thought about finding a complete downloadable database of human genes, containing id,name and whole sequence of the genes! which I could not find!

or

2)using bioperl (or other similar stuff) and Bio::DB::EntrezGene to extract the gene sequences but it seems that "Bio::DB::EntrezGene" is broken : http://stackoverflow.com/questions/16199037/bioperl-module-biodbentrezgene-no-longer-working

what do you suggest dudes? I am really confused!!!

gene bioperl entrez
• 4,937 views • link
updated by Biostar Biostar • written by hamed2005 • hamed2005

You can fetch sequences using Entrez ID as the filter in BioMart. Please search this site for the many usage examples.

• 138 views • link
written by Neilfws Neilfws

Tnx, yes I know, Im now working on BioMart . Do you mean using EMBL Perl API ? or something else?

• 0 views • link
written by hamed2005 • hamed2005

Please note there are many BioMart installations providing many different databases, in this case you mean the Ensembl BioMart, and most likely the "Ensembl Genes" database available from that BioMart, which contains annotations using Entrez Gene.

• 80 views • link
written by Hamish Hamish

yeah bro! It was a semi-solution for me! just a question: which field in the "Ensembl BioMart" is related to the Gene Sequence? (there are lots there!)

• 0 views • link
written by hamed2005 • hamed2005

When you click "Attributes", you should see "Sequences" as an option right at the top.

• 216 views • link
written by Neilfws Neilfws

2 answers

You can probably use links like there:

http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nucleotide&rettype=fasta&retmode=text&id=**6701965,6701969,6702094,6702105,6702160**

Append the gene IDs or GI IDs or RefSeq IDs at the end separated by commas, a file is getting downloaded with all the sequences in FASTA format.

Source: FASTA Batch Download Using Entrez

• 87 views • link
written by SRKR SRKR

You can use this online tool to download fasta/genbank sequences from NCBI by submitting the GI IDs:

NCBI Sequence Retriever

Let me know if your requirement is different from what this tool provides.

• 155 views • link
written by SRKR SRKR

tnx , but it needed gi ID which i don't have, (I have EntrezGene ID and Gene Name), besides there is no way to interact with it via Perl or etc, cause I have more than 1000 genes to be extracted, and it should be done just through programming! let me know about any other options!

P.S: by the way! your tools were cool! I visited your website and I really appreciate your works!

• 0 views • link
written by hamed2005 • hamed2005

Log in to answer this question.

More posts like this