Abstract:By decomposing genome sequences into k -mers, it is possible to estimate genome differences without alignment. Techniques such as k- mer minimisers, for example MinHash, have been developed and are often accurate approximations of distances based on full k -mer sets. These and other alignment-free methods avoid the large temporal and computational expense of alignment. However, these k -mer set comparisons are not entirely accurate within-species and can be completely inaccurate within-lineage. This is due, in part, to their inability to distinguish core polymorphism from accessory differences. Here we present a new approach, KmerAperture , which uses information on the k -mer relative genomic positions to determine the type of polymorphism causing differences in k -mer presence and absence between pairs of genomes. Single SNPs are expected to result in contiguous of k unique k -mers per genome. On the other hand, contiguous series > k may be caused by accessory differences of length S - k +1; when the start and end of the sequence are contiguous with homologous sequence. Alternatively, they may be caused by multiple SNPs within k bp from each other and KmerAperture can determine whether that is the case. To demonstrate use cases KmerAperture was benchmarked using datasets including a very low diversity simulated population with accessory content independent from the number of SNPs, a simulated population were SNPs are spatially dense, a moderately diverse real cluster of genomes ( Escherichia coli ST1193) with a large accessory genome and a low diversity real genome cluster ( Salmonella Typhimurium ST34). We show that KmerAperture can accurately distinguish both core and accessory sequence diversity without alignment, outperforming other k -mer based tools. The KmerAperture algorithm provides a substantive progression in alignment-free methodologies in bacterial comparative genomics. The utility of bacterial genomes for epidemiology and uncovering the genetic basis for phenotypic diversity and adaptation lies in our ability to compare them at scale. Perhaps the most limiting feature of genetic analysis workflows is sequence alignment. A number of k -mer based alignment-free methodologies have been developed to avoid the temporal and compute cost of alignment. However, in closely related bacteria the signal of polymorphisms, as they are represented by unique k -mers, is lost to even a small amount of unshared sequence, whether real or artefactual. Here, we show that it's possible to discern SNP diversity without alignment, including when SNPs are within k of one another, a perennial problem for k -mer analyses.

Diminishing Return for Increased Mappability with Longer Sequencing Reads: Implications of the k-mer Distributions in the Human Genome

Sequence Repetitiveness Quantification and De Novo Repeat Detection by Weighted K-Mer Coverage.

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

k-mer-based approaches to bridging pangenomics and population genetics

KmerAperture : Retaining k -mer synteny for alignment-free extraction of core and accessory differences between bacterial genomes

pathMap: a path-based mapping tool for long noisy reads with high sensitivity

Perm: Efficient Mapping of Short Sequencing Reads with Periodic Full Sensitive Spaced Seeds

Sequencing depth and coverage: key considerations in genomic analyses

Assessment of Kmer Degeneration Method for Complicated Genomes.

Personalizing pangenome graphs with k -mers

A Short Reads Alignment Algorithm Oriented To Massive Data

Estimating the total genome length of a metagenomic sample using k-mers

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Guide to k-mer approaches for genomics across the tree of life

A Survey of K-mer Methods and Applications in Bioinformatics

Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects

KMAP: Kmer Manifold Approximation and Projection for visualizing DNA sequences

Estimation of genome size using k-mer frequencies from corrected long reads

COPE: an Accurate K-Mer-based Pair-End Reads Connection Tool to Facilitate Genome Assembly

Fleximer: Accurate Quantification of RNA-Seq via Variable-Length k-mers

Revisiting K-mer Profile for Effective and Scalable Genome Representation Learning