This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using ensembl biomart retrieving orthologous codin sequence

I have list of transcript IDs for the cow coding sequences and I want to use these IDs to retrieve orthologous sequences from other ensembl species. One such ID is ENSBTAT00000000080.

With the ensembl BioMart the procedure seemed intuitive to me....
1. Ensembl 75 as database
2. Dataset: Bos taurus
3. Filters: Gene (checked the Id list box and placed list of Ids); Multi-species comparison (checked the Homolog filter box and selected Orthologous genes (e.g Dog Orthologous genes)
4. Attributes: ticked the Sequence button at the top.
5. Click to get results.

....until the results seems not the correct ones! I am getting 100% identical sequences for whatever organism i want orthologs. For example, ENSBTAT00000000080 returns the same sequences in dog as in duck. Something is not right. What step I am missing?

gene
• 4,122 views • link
updated by Vitis Vitis • written by ishengomae ishengomae

2 answers

The sequence is coming from the database you're querying, so it'll always be the same. You need to get a list of orthologous gene/transcript IDs and then query their sequence in the database for that species. Perhaps you can do that in biomart online, but if not it's doable in the biomaRt package in R by just programmatically doing those steps.

• 0 views • link
written by Devon Ryan Devon Ryan

Take a look at ensembl compara perl API, with which you could query orthologous genes (predicted by ensembl using MCL and gene trees) in batches.

• 0 views • link
updated by Ram Ram • written by Vitis Vitis

Log in to answer this question.

More posts like this