Question: Functional annotation of genes on chromosome
gravatar for User000
4.9 years ago by
User000420 wrote:


I have a piece of 4Mbp of a chromosome of a plant species. I would like to see which genes are there in this piece of chromosome and what function do they have. I did a blast of a) Ensembl genes (CDS) of 10 orthologous species against my chromosome and found annotation using BioMart table, b) blast of UniprotKB SwissProt proteins of Viridiplantae against my chromosome and used the annotation table there. Do you have any suggestions if the methods I used are good? Or do you know any other ways how to annotate my 4Mbp long sequence?Any suggestions are welcomed.

blast gene • 1.3k views
ADD COMMENTlink modified 4.9 years ago by Juke345.0k • written 4.9 years ago by User000420
gravatar for Juke34
4.9 years ago by
Juke345.0k wrote:

You could use an abinitio approach as AUGUSTUS enhanced by the protein alignements coming from the dataset of the 10 close species. You could try one of the hmm profile of plant they provide (

Once done, you translate the CDSs found in Protein using gffread. Then you can perform a best blast hit approach to name the genes, and/or launch Interproscan to retrieve functional annotation. Or you could also use a clustering approach (with the protein set of close species) to name them according to the sequences present in each cluster (but it's harder to do).

ADD COMMENTlink written 4.9 years ago by Juke345.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1661 users visited in the last hour