Show simple item record

dc.contributor.authorHan, Renmin
dc.contributor.authorLi, Yu
dc.contributor.authorWang, Sheng
dc.contributor.authorGao, Xin
dc.date.accessioned2018-01-15T06:10:39Z
dc.date.available2018-01-15T06:10:39Z
dc.date.issued2017-12-24
dc.identifier.citationHan R, Li Y, Wang S, Gao X (2017) An accurate and rapid continuous wavelet dynamic time warping algorithm for unbalanced global mapping in nanopore sequencing. Available: http://dx.doi.org/10.1101/238857.
dc.identifier.doi10.1101/238857
dc.identifier.urihttp://hdl.handle.net/10754/626744
dc.description.abstractLong-reads, point-of-care, and PCR-free are the promises brought by nanopore sequencing. Among various steps in nanopore data analysis, the global mapping between the raw electrical current signal sequence and the expected signal sequence from the pore model serves as the key building block to base calling, reads mapping, variant identification, and methylation detection. However, the ultra-long reads of nanopore sequencing and an order of magnitude difference in the sampling speeds of the two sequences make the classical dynamic time warping (DTW) and its variants infeasible to solve the problem. Here, we propose a novel multi-level DTW algorithm, cwDTW, based on continuous wavelet transforms with different scales of the two signal sequences. Our algorithm starts from low-resolution wavelet transforms of the two sequences, such that the transformed sequences are short and have similar sampling rates. Then the peaks and nadirs of the transformed sequences are extracted to form feature sequences with similar lengths, which can be easily mapped by the original DTW. Our algorithm then recursively projects the warping path from a lower-resolution level to a higher-resolution one by building a context-dependent boundary and enabling a constrained search for the warping path in the latter. Comprehensive experiments on two real nanopore datasets on human and on Pandoraea pnomenusa, as well as two benchmark datasets from previous studies, demonstrate the efficiency and effectiveness of the proposed algorithm. In particular, cwDTW can almost always generate warping paths that are very close to the original DTW, which are remarkably more accurate than the state-of-the-art methods including FastDTW and PrunedDTW. Meanwhile, on the real nanopore datasets, cwDTW is about 440 times faster than FastDTW and 3000 times faster than the original DTW. Our program is available at https://github.com/realbigws/cwDTW.
dc.description.sponsorshipWe thank Minh Duc Cao, Lachlan J.M. Coin, Louise Roddam, and Tania Duarte for providing the nanopore sequencing data for the Pandoraea pnomenusa sample.
dc.publisherCold Spring Harbor Laboratory
dc.relation.urlhttps://www.biorxiv.org/content/early/2017/12/23/238857.1
dc.rightsThe copyright holder for this preprint is the author/funder. It is made available under a CC-BY-NC-ND 4.0 International license.
dc.rights.urihttp://creativecommons.org/licenses/by-nc-nd/4.0/
dc.titleAn accurate and rapid continuous wavelet dynamic time warping algorithm for unbalanced global mapping in nanopore sequencing
dc.typePreprint
dc.contributor.departmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
dc.contributor.departmentComputer Science Program
dc.contributor.departmentComputational Bioscience Research Center (CBRC)
dc.eprint.versionPre-print
kaust.personHan, Renmin
kaust.personLi, Yu
kaust.personWang, Sheng
kaust.personGao, Xin
refterms.dateFOA2018-06-13T19:21:11Z


Files in this item

Thumbnail
Name:
238857.1.full.pdf
Size:
2.264Mb
Format:
PDF
Description:
Preprint

This item appears in the following Collection(s)

Show simple item record

The copyright holder for this preprint is the author/funder. It is made available under a CC-BY-NC-ND 4.0 International license.
Except where otherwise noted, this item's license is described as The copyright holder for this preprint is the author/funder. It is made available under a CC-BY-NC-ND 4.0 International license.