map methylome data with GKNO mosaik program

Hi,

I found the program mosaik in GKNO tool sets (Lee W-P, Stromberg MP, Ward A, Stewart C, Garrison EP, et al. (2014) MOSAIK: A Hash-Based Algorithm for Accurate Next-Generation Sequencing Short- Read Mapping. PLoS ONE 9(3): e90581. doi:10.1371/journal.pone.0090581) is a powerful mapping method. For pair-end reads, it utilizes the knowledge of approximate fragment length in library construction to rescue the ambiguously/multiple anchored read if its mate is uniquely aligned.

I am wondering whether this progam is able to process methylome data generated by MethylC-seq, and whether the above strategy is deployed in processing such data. My datasets were sequenced on HiSeq2000 with PE100 (pair-end 100bp).

Do anyone has any experience with this problem?

Thanks,

Zhiguang Li

methylc-seq gkno mosaik

• 3,718 views • link

updated 4.7 years ago by Ram

• written 12.0 years ago by zhiguangli88 •

3 answers

My I guess is no, it can't handle BS-Seq data.

If you want to align BS-Seq reads, from directional protocol, with an aligner not designed for BS-Seq data you could follow these steps (but you would probably reinvent the wheel...):

Convert all the C in reference genome to T. Append this converted genome to the one where all the G are converted to A. I.e. each original chromosome is present twice: as C to T converted and G to A converted. Make chromosome names unique (e.g. by appending to the original name _CT or _GA). Index this reference with your favourite aligner.
Convert all the C in your reads mate 1 to T. Convert all the G in mate 2 to A
Align converted reads to the converted reference.
Parse chromosome names in output SAM/BAM to have the original chromosome names (e.g. by removing _CT or _GA)

(Hope I got it right...)

• 0 views • link

updated 4.7 years ago by Ram

• written 12.0 years ago by dariober

Dear Dariober,

Thanks very much for your explanation. But looks like it is pretty complicated. Do you know any public available program that can handle pair-end methylome data?

Thanks again.

Zhiguang Li

• 0 views • link

updated 4.7 years ago by Ram

• written 12.0 years ago by zhiguangli88 •

bismark, which uses bowtie or bowtie2, it's probably the most popular (and for good reason!). Recently I got pretty good results with bwameth.py, which is backed by bwa mem (see also this Bwa-Meth: Align And Tabulate Bs-Seq Reads). Definitively, unless you have a good reason to use a specific aligner, don't reinvent the wheel!

PS: Your comment should have been a "comment" not an answer.

• 0 views • link

written 12.0 years ago by dariober

looks right. at the end, you also need to recover the original reads without the in silico conversions and report those instead of the converted.

• 0 views • link

updated 4.7 years ago by Ram

• written 12.0 years ago by brentp

Sure... Thanks for pointing it out. (And thanks for bwameth!)

• 0 views • link

updated 4.7 years ago by Ram

• written 12.0 years ago by dariober

sure, glad someone is using it.

• 0 views • link

written 12.0 years ago by brentp

As mentioned above, Mosaik is not ideal for your task. I am not familiar with processing methylome data, but you would definitely want to employ a tool that has been developed specifically for handling this type of data. Mosaik was not developed with this in mind. We would be happy to host tools to help parcel together a pipeline, making it a more straightforward task to handle this problem, but the current gkno toolset does not currently contain tools or pipelines for methylome processing.

• 0 views • link

written 12.0 years ago by alistairnward