Summary: Genome wide optical maps are high resolution restriction maps that give a unique numeric representation to a genome. They are produced by assembling hundreds of thousands of single molecule optical maps, which are called Rmaps. Unfortunately, there exists very few choices for assembling Rmap data. There exists only one publicly-available non proprietary method for assembly and one proprietary method that is available via an executable. Furthermore, the publicly-available method, by A. Valouev et al. [Lect. Notes Comput. Sci. 3500, 489–504 (2005; Zbl 1119.92336)], follows the overlap-layout-consensus (OLC) paradigm, and therefore, is unable to scale for relatively large genomes. The algorithm behind the proprietary method, Bionano Genomics’ Solve, is largely unknown. in this paper, we extend the definition of bi-labels in the paired de Bruijn graph to the context of optical mapping data, and present the first de Bruijn graph based method for Rmap assembly. We implement our approach, which we refer to as rmapper, and compare its performance against the assembler of Valouev et al. [loc. cit.] and Solve by Bionano Genomics on data from three genomes: E. coli, human, and climbing perch fish (Anabas testudineus). Our method was the only one able to successfully run on all three genomes. The method of Valouev et al. [loc. cit.] only successfully ran on E. coli and Bionano Solve successfully ran on E. coli and human but not on the fish genome. Moreover, on the human genome rmapper was at least 130 times faster than Bionano Solve, used five times less memory and produced the highest genome fraction with zero mis-assemblies.
For the entire collection see [Zbl 1445.68019].