Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
sequencing artifacts passing GATK's ApplyRecalibration filters efoss Bioinformatics 4 04-04-2012 07:55 AM
any pipeline for tophat with 454 sequences? feng RNA Sequencing 0 11-29-2011 05:41 PM
454 Sequences analysis!! alessandra85 Bioinformatics 2 09-14-2010 12:33 PM
Filtering and excluding reads dawe Bioinformatics 3 09-09-2009 08:47 AM
454 expected sequences and Tag for amplicon sequencing jre06 454 Pyrosequencing 3 07-22-2008 10:44 AM

Thread Tools
Old 08-22-2011, 02:03 AM   #1
Junior Member
Location: Poland

Join Date: Jul 2011
Posts: 3
Question How to collapse sequences excluding sequencing artifacts? (454)

Hi all,
I'm begginer in the next generation sequencing and recently I got my first data. I have some issues with them because there is just one study which had similar approach so far. Briefly I'm working with phylogeny of plants which includes polyploidy events. We have few markers and in some species (e.g. tetraploids) there may be few copies of one gene (marker). To obtain all possible copies we used 454 and got around 100 sequences per species. I tried to used BAPS (bayesian clustering program used for population studies mostly) to decide whether there is one copy or two of each gene. Although it is quite labourious basically it works and I couldn't figure out anything better so far. But then there is one thing which I would like to automate a little bit and make it more objective:
Even when I know that there is one or two copies in my 100 reads there are still sequences which are different only by 1 or 2 bp, usually it's easily visible that it is artifact especially when it is present only once in the dataset. I wonder wheter there is any software that could collapse sequences for me and exlude this artifacts? Let's say that I would like to have sequences that are present with minimum 20% in the dataset. The truth that in most cases in the end I need just 2 sequences out of 100...
Most programs with collapsing option just collapse everything which is identical to one haplotype and doesn't include any onformation about presence in the data or similarities. And they leave some sequences with more than one difference which makes things even more complicated.
Is there any program that I could use for my analysis or do I have to do this by hand? If anybody has any idea I would be grateful!
Thank's a lot in advance for any help!
Yevaud is offline   Reply With Quote
Old 01-29-2013, 01:08 PM   #2
Junior Member
Location: NC

Join Date: Jun 2011
Posts: 3

Hi all -
I have basically the same question. Amplicon Illumina data that's been quality filtered (fastx), duplicate reads removed (fastx), aligned to the reference (~1.2kb) with bwa. Used ShoRAH to predict haplotypes - returned 13 haps with a frequency >1%. Most of these are frequency haplotypes and are due to a single SNP in a single read. Is there a way to collapse such sequences into their nearest relatives. Have done some searching, but no luck yet.
parasitehunter is offline   Reply With Quote
Old 01-29-2013, 06:32 PM   #3
Senior Member
Location: Sydney

Join Date: Feb 2011
Posts: 149

You probably can use the program, cd-hit-est in the cd-hit suite:
You can set a threshold identity (e.g. 90%, 95%), and it clusters all smaller lengthed sequences into a longer representative one when it falls above this threshold, much like generating a unigene set.
Kennels is offline   Reply With Quote
Old 01-30-2013, 06:59 AM   #4
Junior Member
Location: NC

Join Date: Jun 2011
Posts: 3

Kennels -
Thanks for the idea - looks promising. I've been trying to get cd-hit-est to work (both locally and on their servers), but it keeps returning an error. Probably something I'm doing wrong. Or perhaps it's because all my predicted haplotypes are the same length. However, their cd_454 clusterer seems to work with my data. Hope that's legit to use ...
parasitehunter is offline   Reply With Quote

454, collapsing, phylogeny, polyploidy

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:46 PM.

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