getting fasta sequence of a list of ensemble gene IDs

Hi,

I retrieve the Ensembl gene IDs of yeast ribisomal RNA from biomart-Ensembl

Ensembl Gene ID
RDN25-1
RDN18-2
RDN5-4
15S_rRNA
RDN37-2
RDN5-6
RDN5-3
RDN58-1
RDN18-1
RDN37-1

I need to download the fasta file of these IDs but in Ensembl I could not such a option. Do you have any suggestion please?

Thank you

gene ensembl sequence biomart

• 7,494 views • link

updated 4.1 years ago by Ram

• written 10.5 years ago by zizigolu

3 answers

You can use BioMart. Follow the sequence below.

Ensembl --> BioMart --> Choose Database --> "Ensembl Genes" --> Select Yeast Genome from list --> Filters (in the left pane) --> Gene --> Input external referenced ID list --> Paste ID's in --> Attributes (left pane) --> Sequence --> Select as needed --> Results Button at top of the page --> Export to "File" as "FASTA".

• 0 views • link

written 10.5 years ago by GenoMax

thank you so much.

sorry I searched for Arabidopsis rRNA in plant biomart, I found the IDs but the sequence was unavalaible then I tried the IDs in NCBI nucleotide but there is also nothing there.

• 0 views • link

updated 6.6 years ago by Ram

• written 10.5 years ago by zizigolu

Are your ID's from Yeast or Arabidopsis? The ones in your original post are yeast.

• 0 views • link

written 10.5 years ago by GenoMax

For Arabidopsis rRNA use the following path:

Ensembl Plant --> BioMart --> Choose Database --> Plant Mart --> Select Arabidopsis Genome from list --> Filters (in the left pane) --> Gene --> Gene type (4th option) --> Select "rRNA" --> Attributes (left pane) --> Sequence --> Select features as needed (Unspliced gene may be one option) --> Results Button at top of the page --> Export to "File" as "FASTA".

• 0 views • link

written 10.5 years ago by GenoMax

thank you soooo much genomax2,

my IDs were from yeas that you resolved my problem. then I asked you about Arabidopsis that your tip worked well as already for Arabidopsis too

• 0 views • link

updated 6.6 years ago by Ram

• written 10.5 years ago by zizigolu

I'm sorry, I'm very confused. The list are bread yeast IDs. You can easily query and get the sequence, for example. What does Arabidopsis have to do with it?

• 0 views • link

updated 6.6 years ago by Ram

• written 10.5 years ago by Emily

This is one line of code with gget seq:

pip install gget, then simply:

# Command-line
gget seq -o yeast.fasta RDN25-1 RDN18-2 RDN5-4

# Python
import gget
gget.seq(["RDN25-1", "RDN18-2", "RDN5-4"], save=True)

• 0 views • link

written 4.2 years ago by Laura Luebbert •

I could not find a decent way of automating the process through eutils (which I'm sure there are plenty of people who can show you how), but if you don't mind doing these few manually then you can find the sequence quite easily through both the

ncbi gene database:

http://www.ncbi.nlm.nih.gov/gene/9164935

or SGD's database:

http://www.yeastgenome.org/locus/S000006484/sequence

• 0 views • link

written 10.5 years ago by Jenez