SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cufflinks generated >400,000 transcripts?? IrisZhu Bioinformatics 13 05-15-2014 08:15 PM
help! samtools gave me more than 500,000 snps slowsmile Bioinformatics 1 12-15-2011 08:24 AM
PubMed: Cost-effective multiplexing before capture allows screening of 25 000 clinica Newsbot! Literature Watch 0 03-19-2011 10:41 AM
20,000 posts! ECO Site Announcements 0 07-05-2010 06:12 PM
help in building sequences in bowtie lilithdog General 1 12-17-2009 03:17 AM

Reply
 
Thread Tools
Old 12-10-2011, 01:03 PM   #1
brachysclereid
Member
 
Location: California

Join Date: Feb 2011
Posts: 32
Default building a cladogram with 100,000 sequences...

Are there any tools capable of performing multi sequence alignment and tree construction that can handle up to 100,000 sequences?
brachysclereid is offline   Reply With Quote
Old 12-12-2011, 05:13 AM   #2
mbblack
Senior Member
 
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245
Default

Quote:
Originally Posted by brachysclereid View Post
Are there any tools capable of performing multi sequence alignment and tree construction that can handle up to 100,000 sequences?
That many sequences is still pretty much in the intractable range, to my knowledge. There are tools that reportedly do alright with up to 10,000 sequences (although most of the ones I know of - SATCHMO, MUMMALS, etc) are designed for protein sequences, although MAFFT will handle protein or DNA.

In order to do what you are asking, you really should cluster the sequences first, then align from smaller sets.

MUSCLE's manual recommends using the UCLUST tool in USEARCH (http://www.drive5.com/usearch/) to reduce the number of alignments needed. And if you contact MUSCLE's author (http://www.drive5.com/muscle/about.htm) he mentions that he is working on something to leverage MUSCLE and USEARCH to deal with huge numbers of sequence alignments (without overly sacrificing accuracy).
mbblack is offline   Reply With Quote
Old 12-12-2011, 11:51 AM   #3
koadman
Member
 
Location: Sydney, Australia

Join Date: May 2010
Posts: 65
Default

Are your sequences all from the same locus, e.g. amplicon sequenced? If so, the usual approach to build a phylogeny involves first taking a representative subset, say 500 sequences, and building a profile HMM using HMMer3 or building a SCFG using Infernal if the sequences are RNA with known secondary structure. With your HMM or SCFG in hand you can then quickly construct a multiple alignment of all 100k of them. see hmmalign. You can then feed the alignment to Morgan Price's FastTree or one of the other programs that can construct a phylogeny in subquadratic time.

That approach works well if you care only about substitutions in the evolutionary history (this is what most people want since they usually hold the richest evolutionary signal). If you care about indel histories you will have to do a full multiple sequence alignment with something like MUSCLE, MAFFT or another tool. I am not aware of any indel history inference tools that can operate on such large datasets. Maybe someone else knows of such methods?
koadman is offline   Reply With Quote
Old 12-13-2011, 12:14 AM   #4
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

Can you take a subset by looking at clusters ? Some good software for this purpose is RAMMCAP.
colindaven is offline   Reply With Quote
Old 12-13-2011, 08:26 AM   #5
ssing
Member
 
Location: usa

Join Date: Jan 2009
Posts: 20
Default

RAxML should be able to handle that many tips, but it cannot align... GARLI is also efficient at making large trees from alignments, but I don't know what its upper limit is.
ssing is offline   Reply With Quote
Old 12-14-2011, 06:07 AM   #6
brachysclereid
Member
 
Location: California

Join Date: Feb 2011
Posts: 32
Default building a cladogram with 100,000 sequences

Thanks for all the advice. It will take a while to sort through all the possibilities. Perhaps clustering first will be a good approach.

I'll make a post if I can get something to work. In the mean time other suggestions would be greatly appreciated!
brachysclereid 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 04:24 PM.


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