Year: 2021
Type: preprint
Abstract: ABSTRACT Optical mapping is a method for creating high resolution restriction maps of an entire genome. Optical mapping has been largely automated, and first produces single molecule restriction maps, called Rmaps, which are assembled to generate genome wide optical maps. Since the location and orientation of each Rmap is unknown, the first problem in the analysis of this data is finding related Rmaps, i.e., pairs of Rmaps that share the same orientation and have significant overlap in their genomic location. Although heuristics for identifying related Rmaps exist, they all require quantization of the data which leads to a loss in the precision. In this paper, we propose a Gaussian mixture modelling clustering based method, which we refer to as OM clust , that finds overlapping Rmaps without quantization. Using both simulated and real datasets, we show that OM clust substantially improves the precision (from 48.3% to 73.3%) over the state-of-the art methods while also reducing CPU time and memory consumption. Further, we integrated OM clust into the error correction methods (Elmeri and c OM et ) to demonstrate the increase in the performance of these methods. When OM clust was combined with c OM et to error correct Rmap data generated from human DNA, it was able to error correct close to 3x more Rmaps, and reduced the CPU time by more than 35x. Our software is written in C++ and is publicly available under GNU General Public License at https://github.com/kingufl/OMclust
Institutions: University of Florida, Morgan State University
Language: English
Subfield: Molecular Biology
Domain: Life Sciences
Open Access status: green
Funder: National Science Foundation
OpenAlex ID: W3184642780