SEQanswers

Go Back   SEQanswers > Applications Forums > Epigenetics

Similar Threads
Thread Thread Starter Forum Replies Last Post
mapping on to custom reference sequence PratikC Bioinformatics 13 04-03-2012 12:31 PM
bisulphite converted barcoded fragment library for SOLID ageliki Sample Prep / Library Generation 2 12-23-2011 11:33 AM
adjust solid barcoded adaptors for preparing library from bisulphite converted DNA ageliki SOLiD 1 12-22-2011 04:18 PM
Bisulphite sequencing - Ion Torrent arnaud.kr Ion Torrent 1 11-11-2011 03:06 AM
Samtools (sequence alignment/mapping) wisosonic Bioinformatics 2 03-28-2011 04:31 AM

Reply
 
Thread Tools
Old 10-28-2011, 08:03 AM   #21
kalyankpy
PostDoc
 
Location: Turku, Finland

Join Date: Mar 2010
Posts: 15
Default

Hi Yuanxin,

I am quite interested in BSMAP. I am working on RRBS data from stem cells. I am trying to install BSMAP on my computer (Ubuntu 11.10, 64 bit). But it stops with an error. I posted the error messge at http://code.google.com/p/bsmap/issues/detail?id=7 . Can you fix this!

Kalyan
kalyankpy is offline   Reply With Quote
Old 10-28-2011, 08:19 AM   #22
yxi
Junior Member
 
Location: texas

Join Date: Apr 2009
Posts: 6
Default

Hi Kalyan,

The error is in compiling samtools, bsmap used samtools API to read SAM files.

After google search I found someone experienced similar problems.

http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=606004

Hope this could help

Thanks,

Yuanxin

Quote:
Originally Posted by kalyankpy View Post
Hi Yuanxin,

I am quite interested in BSMAP. I am working on RRBS data from stem cells. I am trying to install BSMAP on my computer (Ubuntu 11.10, 64 bit). But it stops with an error. I posted the error messge at http://code.google.com/p/bsmap/issues/detail?id=7 . Can you fix this!

Kalyan
yxi is offline   Reply With Quote
Old 10-31-2011, 10:20 PM   #23
kalyankpy
PostDoc
 
Location: Turku, Finland

Join Date: Mar 2010
Posts: 15
Default BSMAP installation on ubuntu

Hi Yuanxin,

I went through the link suggested by you! I noticed that there was a bug in the samtools in earlier versions which was rectified after 0.1.13 version (I noticed that the samtools you have provided with bsmap is 0.1.7).
I tried to overcome this bug by trying to download and install samtools-0.1.17 version from ubuntu repositories. I could install it successfully from repository. I also tried by downloading the sourcecode and compiled it with 'make' command. It workd perfectly. Assuming that problem is in the samtools, I replaced the samtools folder in bsmap with the newer version of samtools source code (0.1.18) and then tried compiling (with 'make').

This option took me few step ahead. Here I am attaching the output of the 'make' command. For me the problem doesnt seem to be in samtools rather it seems to be in some code that links BSMAP and SAMTOOLS (I am not a good programmer, I am learning! Its only my assumption. As a programmer you can analyze better!).
Attached Files
File Type: txt bsmap_ubuntu_compilation_error2.txt (18.8 KB, 1 views)
kalyankpy is offline   Reply With Quote
Old 11-01-2011, 05:34 AM   #24
yxi
Junior Member
 
Location: texas

Join Date: Apr 2009
Posts: 6
Default

Good to know you solved it. We experienced a runtime error using the later versions of SAMTOOLS on Power7 CPU based system. So BSMAP included an early version SAMTOOLS.

Quote:
Originally Posted by kalyankpy View Post
Hi Yuanxin,

I went through the link suggested by you! I noticed that there was a bug in the samtools in earlier versions which was rectified after 0.1.13 version (I noticed that the samtools you have provided with bsmap is 0.1.7).
I tried to overcome this bug by trying to download and install samtools-0.1.17 version from ubuntu repositories. I could install it successfully from repository. I also tried by downloading the sourcecode and compiled it with 'make' command. It workd perfectly. Assuming that problem is in the samtools, I replaced the samtools folder in bsmap with the newer version of samtools source code (0.1.18) and then tried compiling (with 'make').

This option took me few step ahead. Here I am attaching the output of the 'make' command. For me the problem doesnt seem to be in samtools rather it seems to be in some code that links BSMAP and SAMTOOLS (I am not a good programmer, I am learning! Its only my assumption. As a programmer you can analyze better!).
yxi is offline   Reply With Quote
Old 11-01-2011, 05:46 AM   #25
kalyankpy
PostDoc
 
Location: Turku, Finland

Join Date: Mar 2010
Posts: 15
Default

Hi, The problem still xists. I dont know if my message is clear for you. I could successfully install SAMTOOLS as a standalone application. However, I fail to install it when it is part of bsmap. Even if I replace the older version of SAMTOOLS with the newer version inside the BSMAP folder, the problem continues.

I liked the logic of your programme. I am just waiting for a solution to this installation problem.

Is it possible for you to send me the executibles for a 64 bit ubuntu system that can work for me (so that I dont need to compile them on mine as this problem persists)
kalyankpy is offline   Reply With Quote
Old 11-07-2011, 03:44 AM   #26
kalyankpy
PostDoc
 
Location: Turku, Finland

Join Date: Mar 2010
Posts: 15
Default Fix to the ubuntu installation of BSMAP

Hi,

Yes finally the problem is solved. Here I explain the steps that help me resolve the installation steps in Ubuntu 11.04

tweak the makefile in the bsmap directory as follows

original line "$(CC) $(FLAGS) $(LIBS) $(REF_MODE) $(OLIGOLEN) $^ -o $@ $(THREAD) -lz -lbam"

Change it as "$(CC) $(FLAGS) $(LIBS) $(REF_MODE) $(OLIGOLEN) $^ -o $@ $(THREAD) -lbam -lz"

replace the older version of samtools source directory with the new version - 0.1.17

I dont understand why changing the order of "-lz -lbam" to "-lbam -lz" solves this issue.
kalyankpy is offline   Reply With Quote
Old 04-03-2012, 12:29 PM   #27
aniruddha.otago
Member
 
Location: New Zealand

Join Date: Jan 2010
Posts: 19
Default

Hi Guys,

Recently we have published a paper in Nucleic Acids research- "Comparison of alignment software for genome-wide bisulphite sequence data". you might find it useful to have a look at it.. http://www.ncbi.nlm.nih.gov/pubmed/22344695
aniruddha.otago is offline   Reply With Quote
Old 04-03-2012, 12:55 PM   #28
rskr
Senior Member
 
Location: Santa Fe, NM

Join Date: Oct 2010
Posts: 244
Default

Quote:
Originally Posted by aniruddha.otago View Post
Hi Guys,

Recently we have published a paper in Nucleic Acids research- "Comparison of alignment software for genome-wide bisulphite sequence data". you might find it useful to have a look at it.. http://www.ncbi.nlm.nih.gov/pubmed/22344695
Doesn't look like you compared the accuracy of the mapping algorithms say by generating randomized cpg data. Any good reason why this wasn't done?
rskr is offline   Reply With Quote
Old 04-03-2012, 01:03 PM   #29
yxibcm
Junior Member
 
Location: houston

Join Date: Jun 2010
Posts: 6
Default

The BSMAP version (v1.02/v1.2) used in this paper is a bit old. You may want to use the latest version (v2.5) to get an up-to-date comparison.

Quote:
Originally Posted by aniruddha.otago View Post
Hi Guys,

Recently we have published a paper in Nucleic Acids research- "Comparison of alignment software for genome-wide bisulphite sequence data". you might find it useful to have a look at it.. http://www.ncbi.nlm.nih.gov/pubmed/22344695
yxibcm is offline   Reply With Quote
Old 04-03-2012, 01:34 PM   #30
rskr
Senior Member
 
Location: Santa Fe, NM

Join Date: Oct 2010
Posts: 244
Default

Quote:
Originally Posted by yxibcm View Post
The BSMAP version (v1.02/v1.2) used in this paper is a bit old. You may want to use the latest version (v2.5) to get an up-to-date comparison.

I guess old funky data "The accuracy went up when they trimmed it", has to get published somehow.
rskr is offline   Reply With Quote
Old 04-03-2012, 08:07 PM   #31
aniruddha.otago
Member
 
Location: New Zealand

Join Date: Jan 2010
Posts: 19
Default

yeah when we were performing the analysis the latest version was 1.2 ( release, April, 2011). every program has newer version each week, this analysis was done at that time with latest versions of all the aligners at that particular time.

Cheers,

Aniruddha.
aniruddha.otago is offline   Reply With Quote
Old 04-18-2012, 10:42 PM   #32
aniruddha.otago
Member
 
Location: New Zealand

Join Date: Jan 2010
Posts: 19
Default

Hi Yuanxin,

Thanks, we have solved the problem and it came up with reasonable percentage of methylation. It was useful for the analysis. Thanks.
aniruddha.otago is offline   Reply With Quote
Old 04-19-2012, 06:28 PM   #33
aniruddha.otago
Member
 
Location: New Zealand

Join Date: Jan 2010
Posts: 19
Default

Hi rskr,

We never really thought of trying to do that. We were attempting
to see the performance with real RRBS data which we did. I'm not
sure that going to the trouble of trying to generate random CpG fragments
would have achieved very much at all. Any attempt to 'randomise' the
C/T conversions of somehow randomly chosen CpGs can't be guaranteed
to reflect real life methylation patterns. Until the rules for positional
methylation are more fully understood, attempts to generate randomised
data are likely to be influenced by artefacts of the process used.
aniruddha.otago is offline   Reply With Quote
Old 04-20-2012, 07:20 AM   #34
rskr
Senior Member
 
Location: Santa Fe, NM

Join Date: Oct 2010
Posts: 244
Default

Quote:
Originally Posted by aniruddha.otago View Post
Hi rskr,

We never really thought of trying to do that. We were attempting
to see the performance with real RRBS data which we did. I'm not
sure that going to the trouble of trying to generate random CpG fragments
would have achieved very much at all. Any attempt to 'randomise' the
C/T conversions of somehow randomly chosen CpGs can't be guaranteed
to reflect real life methylation patterns. Until the rules for positional
methylation are more fully understood, attempts to generate randomised
data are likely to be influenced by artefacts of the process used.
It may true that the distribution of CpG methylation is potentially not random. However simply calculating unique matches doesn't tell you anything about the performance of the mapper, which is what you claimed to have measured. It may very well be that the correct mapping was indeed ambiguous, and you don't know. If you can't model the CpG methylation, then you can't claim to have measured the performance of the mapper.
rskr is offline   Reply With Quote
Old 04-20-2012, 01:49 PM   #35
aniruddha.otago
Member
 
Location: New Zealand

Join Date: Jan 2010
Posts: 19
Default

Quote:
Originally Posted by rskr View Post
It may true that the distribution of CpG methylation is potentially not random. However simply calculating unique matches doesn't tell you anything about the performance of the mapper, which is what you claimed to have measured. It may very well be that the correct mapping was indeed ambiguous, and you don't know. If you can't model the CpG methylation, then you can't claim to have measured the performance of the mapper.
it is done on a reduced represented genome. so, say there is a fragment of 100 bp and 8 reads uniquely mapped to that. mind you each fragment and read has to have a CGG or TGG (MspI site) at the start to be able to map. say, these 8 reads has total 16 CpG site. so, we counted this as a no of unique aligned CpG site. we did make commnet with a figure that even for the same region different programs showed different mapping. However, randomising CpG site didnt seem appropriate to us in the context of this paper. "It may very well be that the correct mapping was indeed ambiguous, and you don't know". If thats the case then no point doing any other step, because simply it is not giving you correct information. The aligners choosen here are widely used in the field and are well established to yield reliable results and we have no reason to believe that the correct mapping is ambigious.
we wanted to give a feel to the molecular biologists that if you have a RRBS dataset and use this widely used aligners, this is what you can get and what are the aspects one should be aware of. and also, a biologists is more likely to will work with real sample to analyse, not on "randomised CpG data".
aniruddha.otago is offline   Reply With Quote
Old 04-20-2012, 08:27 PM   #36
rskr
Senior Member
 
Location: Santa Fe, NM

Join Date: Oct 2010
Posts: 244
Default

Quote:
Originally Posted by aniruddha.otago View Post
... we wanted to give a feel to the molecular biologists that if you have a RRBS dataset and use this widely used aligners, this is what you can get and what are the aspects one should be aware of. and also, a biologists is more likely to will work with real sample to analyse, not on "randomised CpG data".
That's fine except I gather that your data is already out of date in technology terms(GAIIx), and wasn't very good to begin with, since you ended up trimming to make the mapping rates go up, indicating the quality dropped severely near the ends. Is that because of the RRBS protocol or something else? And as you mentioned you really don't have a model for methylation so there isn't any reason to believe that your results would be reproducible.
rskr is offline   Reply With Quote
Reply

Tags
bisulphite, epigenetics, mapping

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 10:04 AM.


Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.