SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Help regg. tools for comparing two sets of contigs from de-novo assembly jsreddy82 Bioinformatics 1 03-04-2014 06:13 AM
How are K-mers used to form contigs in de novo assembly? kdarbs Bioinformatics 2 09-26-2013 05:11 PM
Miseq de novo assembly : Ambigous base pairs (NNs) in the contigs ndeshpan Bioinformatics 2 07-21-2013 03:59 PM
de novo assembly with MIRA and 454 single-end reads. Too much contigs fgajardoe De novo discovery 6 04-17-2013 05:03 AM
How to allow mapped contigs to be used in de novo assembly projects edleloth 454 Pyrosequencing 1 03-08-2011 07:13 AM

Reply
 
Thread Tools
Old 02-06-2014, 10:16 AM   #1
witty
Junior Member
 
Location: RI

Join Date: Nov 2013
Posts: 5
Default Annotation for contigs from de novo assembly

Hi,

I want to annotate my assembled contigs (from de novo assembly). I used BLASTX and only got 10~20% percentage of hits(evalue=1e-5). Now all my differentially expressed contigs (genes) have no annotation. At least I want to know what these genes are, e.g, signaling, transmembrane etc.

Thanks a lot!
Victoria
witty is offline   Reply With Quote
Old 02-06-2014, 10:33 AM   #2
themerlin
Member
 
Location: Flagstaff, AZ

Join Date: Feb 2010
Posts: 51
Default

I'd give Prokka a try:

http://www.vicbioinformatics.com/software.prokka.shtml
themerlin is offline   Reply With Quote
Old 02-06-2014, 11:00 AM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,077
Default

Provided Victoria is working with a prokaryotic genome

NCBI has a eukaryotic annotation pipeline: http://www.ncbi.nlm.nih.gov/genome/a...n_euk/process/ and a prokaryotic one: https://www.ncbi.nlm.nih.gov/genome/annotation_prok/ If I recall right, you will have to make the sequence public though at some point in time if you use these.

Other eukaryotic options (have not used myself):

Pasa: http://pasa.sourceforge.net/
Maker: http://www.gmod.org/wiki/MAKER
GenoMax is offline   Reply With Quote
Old 02-06-2014, 10:45 PM   #4
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

I think Blast2GO would also be useful
dnusol is offline   Reply With Quote
Old 02-07-2014, 01:16 AM   #5
ChrisAU
Junior Member
 
Location: Tübingen

Join Date: Jan 2012
Posts: 9
Default

I've also had good experience with Blast2GO, it doesn't require installation and is quite easy to handle. Also, they updated the quite ugly colours of their pie charts
ChrisAU is offline   Reply With Quote
Old 02-07-2014, 11:30 AM   #6
witty
Junior Member
 
Location: RI

Join Date: Nov 2013
Posts: 5
Default

Hi,

Thank you for your reply. I understand that blast2go (see the below link) just used blast result so basically it won't provide more annotated contigs than BLASTX that I did, is it correct?

https://sites.google.com/a/brown.edu...d/b2g4pipe-2-5

The organism I want to annotate is the protist, Oxyrrhis Marina.

Thank you!
Victoria
witty is offline   Reply With Quote
Old 02-11-2014, 02:10 AM   #7
Krish_143
Member
 
Location: Sweden

Join Date: Jan 2012
Posts: 45
Default

RAST annotation.
http://rast.nmpdr.org/
__________________
Krishna
Krish_143 is offline   Reply With Quote
Old 02-25-2014, 12:26 AM   #8
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

Hi Victoria, I guess you could use several databases to increase your chances of annotation. What databases have you used? I don't have experience with protists but in general a good start could be to compare against GenBank and Uniprot's Swiss-Prot and TrEMBL protein databases. Have you tried a less conservative e-value? Also try to download similar species that are annotated to compare directly. This reference may help you

http://www.plosone.org/article/info%...l.pone.0014202

Dave
dnusol is offline   Reply With Quote
Old 02-25-2014, 05:44 AM   #9
Birdman
Member
 
Location: Montreal

Join Date: Jan 2014
Posts: 21
Default

You can try the Trinotate pipeline. It involves several tools (TransDecoder to get plausible ORFs, PFAM, HMMER, signalIP, tmHMM, RNAmmer) to obtain a quite complete annotation report. They give a lot of details on the website on how to use it.
Birdman is offline   Reply With Quote
Old 03-02-2014, 09:43 AM   #10
bastianwur
Member
 
Location: Germany/Netherlands

Join Date: Feb 2014
Posts: 98
Default

Run a gene prediction tool (e.g. prodigal) over it, throw the proteins in InterproScan, and check if you get anything interesting for your analysis.

Might as well be good to know how long the contigs are.
Will not be of much use to annotate stuff, which is considerable less long than 900 bp.
bastianwur is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:50 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO