SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Ion Torrent



Similar Threads
Thread Thread Starter Forum Replies Last Post
Ion Torrent Mappers golharam Bioinformatics 7 11-07-2012 07:06 PM
Ion Torrent $1000 Genome!? Benchtop Ion Proton Sequencer aeonsim Ion Torrent 88 10-28-2012 05:50 AM
Ion Torrent, how much to run? brianlamx Ion Torrent 18 08-31-2011 03:12 PM
ion torrent herrroaa Introductions 5 07-25-2011 06:36 AM
Ion Torrent through the roof... james hadfield Ion Torrent 14 03-21-2011 10:34 AM

Reply
 
Thread Tools
Old 08-18-2014, 10:00 AM   #1
bioWizz
Junior Member
 
Location: Manipal

Join Date: Jul 2013
Posts: 6
Default Bismark for Ion Torrent

I am trying to use Bismark on targeted bisulphite sequencing done on PGM. Average read length is about ~102 maximum being ~190.
Bismark uses bowtie2 and running it with varying parameters gives me aa lignment score of 8 - 10 % which is very low. When I use TMAP on torrent suite to map to converted reference sequence, I get decent alignment percentage.
Any suggestions on bowtie2 parameters for Ion Torrent variable read length data? Suggestion on analysis of targeted bisulphite sequencing with PGM output will also do.
Thanks
bioWizz is offline   Reply With Quote
Old 08-19-2014, 12:09 AM   #2
frozenlyse
Senior Member
 
Location: Australia

Join Date: Sep 2008
Posts: 136
Default

I wouldnt use bismark on ion torrent data as it is going to have a lot of indels due to homopolymers, and bismark doesnt run bowtie2 in local alignment mode - I would try https://github.com/brentp/bwa-meth as it uses bwa for the alignment - however it may require paired end reads!
frozenlyse is offline   Reply With Quote
Old 08-19-2014, 12:41 AM   #3
bioWizz
Junior Member
 
Location: Manipal

Join Date: Jul 2013
Posts: 6
Default

As you mentioned I had already tried running bowtie2 with local alignment mode, but it did not work.
Just looked into BWA-meth, it does say it supports single end reads. Looks like the tools I was looking for as it employs a similar method as Bismark. Will try to use it and update.
Thank you
bioWizz is offline   Reply With Quote
Old 08-19-2014, 01:37 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,476
Default

It can be a pain to tweak bowtie2's local alignment setting to get decent results (that's one of the reasons I have Bison default to end-to-end alignment). BWA-meth's defaults are pretty reasonable and might give decent enough results, though as frozenlyse mentioned, ion torrent data is often riddled with InDels. Brent's pretty good with timely feedback, so you might just drop him an email if you don't get decent results.

BTW, you can also try sending me a small test dataset and I can have a look with Bison and see if I can send you some reasonable settings, though I expect BWA-meth will work.
dpryan is offline   Reply With Quote
Old 08-19-2014, 03:58 AM   #5
IonTom
Member
 
Location: Germany

Join Date: Apr 2014
Posts: 32
Default

I don know if the original BWA works Ion Torrent data. I am currently using BWA-Mem and it works fine. Is there a methylation mapper based on Mem ?

There is also NextGenmap which works well for Ion Torrent data and also includes
a methylation mapper.
IonTom is offline   Reply With Quote
Old 08-19-2014, 04:28 AM   #6
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,476
Default

@IonTom: See the aforementioned bwa-meth package from Brent Pederson.
dpryan is offline   Reply With Quote
Old 11-07-2014, 02:41 PM   #7
vz33
Junior Member
 
Location: London

Join Date: Aug 2012
Posts: 2
Default

@bioWizz and @IonTom, which kind of bisulphite sequencing do you did by PGM? RRBS, TAB-Seq or others? Could you share a protocol to me? Thank you.
vz33 is offline   Reply With Quote
Reply

Tags
bismark, bs-seq, ion torrent

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 09:05 AM.


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