SEQanswers

Go Back   SEQanswers > Applications Forums > Metagenomics



Similar Threads
Thread Thread Starter Forum Replies Last Post
aligning pre-assembled transcriptomes to genomes skeffington Bioinformatics 2 02-11-2015 03:28 AM
Identify TFs in New Assembled Organism Fernas Bioinformatics 1 01-23-2015 01:17 PM
Comparing two de-novo assembled genomes calbertovj Bioinformatics 1 02-18-2014 05:08 AM
Need to identify conserved exons among 3 genomes kmkocot Bioinformatics 1 06-14-2012 07:19 AM

Reply
 
Thread Tools
Old 12-02-2015, 10:20 AM   #1
fibar
Member
 
Location: Argentina

Join Date: Feb 2013
Posts: 19
Default What's the best way to identify the taxonomy of assembled genomes?

Hi everyone,
Briefly, I applied coverage-based binning method plus pentamer-frequency ordination on metagenomic assembled contigs, to isolate individual genomes.

Now I'm looking for clever (confident, efficient) methods to answer following questions:
1) Who do those genomes belong to? (ie. identify their taxonomy)
2) How complete are they?

For (1) I've tried AmphoraNet and found it quite good. But what about blasting everything against NCBI-genomes database? Would that work fast enough on a small server (Xeon, 16 cores, 96 GB memory)?

Question (2) could probably be answered if (1) is reached succesfully.

However, it could happen that I get, e.g., a novel Proteobacteria, but how to be sure this is really the case and it is not due to misassembly issues. It's a challenging discipline.

Thanks for your insights.

Last edited by fibar; 12-02-2015 at 10:36 AM. Reason: Increase accuracy of thread.
fibar is offline   Reply With Quote
Old 03-21-2016, 04:03 AM   #2
Markiyan
Senior Member
 
Location: Cambridge

Join Date: Sep 2010
Posts: 115
Default For the course level use blastx(p), for more refined - blastn

I may suggest trying blastx or blastp (if you've got good protein predictions) for the coarse level of closest neighbour finding.

For more refined use blastn on the DNA sequence (nt or similar db).
It is a good idea to create more specific database (esp nt one), and have them on dedicated SSD (esp. nvme or PCIe one).

On the performance side of things, give blastn 900bp or 1kb chunks, turn megablast on and upp the word size. "-n T -W 32".
For blastx - raise the expect to 1e-5 or 1e-10, and give it 450 bp chunks.

Also it is good idea to play with number of cpus per job (depending on number of cores per cpu, db/ram size).

PS: also there is a blast-like tool called diamond (used by AntiSMASH), which may be a bit faster in some cases, but the developers kept changing it's output format, so not as stable as blast, if you need to write your own parser.
Markiyan is offline   Reply With Quote
Reply

Tags
assembly annotation, assembly metagenome, completeness, taxonomy

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 03:24 AM.


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