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.

PARMIK: PArtial Read Matching with Inexpensive K-mers

Mibianto: ultra-efficient online microbiome analysis through k -mer based metagenomics

MEM-based pangenome indexing for k-mer queries

Efficient De Novo Assembly and Recovery of Microbial Genomes from Complex Metagenomes Using a Reduced Set of k-mers

KMCP: accurate metagenomic profiling of both prokaryotic and viral populations by pseudo-mapping

Fast Detection of Maximal Exact Matches Via Fixed Sampling of Queryk-Mers and Bloom Filtering of Indexk-Mers

GPMeta: a GPU-accelerated Method for Ultrarapid Pathogen Identification from Metagenomic Sequences

PC-mer: An Ultra-fast memory-efficient tool for metagenomics profiling and classification

PanKmer: k-mer-based and reference-free pangenome analysis

Enhancing antimicrobial resistance detection with MetaGeneMiner: Targeted gene extraction from metagenomes

Parallel Read Partitioning for Concurrent Assembly of Metagenomic Data

MetaPalette: A $k$-mer painting approach for metagenomic taxonomic profiling and quantification of novel strain variation

KPop: Accurate and scalable comparative analysis of microbial genomes by sequence embeddings

Clinical PathoScope: rapid alignment and filtration for accurate pathogen identification in clinical samples using unassembled sequencing data

Exploration and retrieval of whole-metagenome sequencing samples

KmerGO: A Tool to Identify Group-Specific Sequences With k-mers

MSPKmerCounter: A Fast and Memory Efficient Approach for K-mer Counting

Mipe: A Metagenome-Based Community Structure Explorer and Ssu Primer Evaluation Tool

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

copMEM: Finding maximal exact matches via sampling both genomes

PM-profiler: a high-resolution and fast tool for taxonomy annotation of amplicon-based microbiome