Compact and evenly distributed k-mer binning for genomic sequences

Published: Oct. 12, 2020, 1:01 p.m.

Link to bioRxiv paper: http://biorxiv.org/cgi/content/short/2020.10.12.335364v1?rss=1 Authors: Nyström-Persson, J. T., Keeble-Gagnere, G., Zawad, N. Abstract: The processing of k-mers (subsequences of length k) is at the foundation of many sequence processing algorithms in bioinformatics, including k-mer counting for genome size estimation, genome assembly, and taxonomic classification for metagenomics. Minimizers - ordered m-mers where m < k - are often used to group k-mers into bins as a first step in such processing. However, minimizers are known to generate bins of very different sizes, which can pose challenges for distributed and parallel processing, as well as generally increase memory requirements. Furthermore, although various minimizer orderings have been proposed, their practical value for improving tool efficiency has not yet been fully explored. Here we present a distributed k-mer counting tool based on Apache Spark, Discount, which we use to investigate the behaviour of various minimizer orderings in practice when applied to metagenomics data. Using this tool, we then introduce a new combination of frequency counted minimizers and universal k-mer hitting sets, the universal frequency ordering, which yields both evenly distributed binning and small bin sizes. We show that this new ordering allows Discount to perform distributed k-mer counting on a large dataset in as little as 1/8 of the memory of comparable approaches, making it the most efficient out-of-core distributed k-mer counting method available. Copy rights belong to original authors. Visit the link for more info