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.

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

Identifying Group-Specific Sequences for Microbial Communities Using Long K-Mer Sequence Signatures.

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

These are not the k-mers you are looking for: efficient online k-mer counting using a probabilistic data structure

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

A computational toolset for rapid identification of SARS-CoV-2, other viruses and microorganisms from sequencing data

Kmerlight: fast and accurate k-mer abundance estimation

GECKO is a genetic algorithm to classify and explore high throughput sequencing data

Klumpy: A Tool to Evaluate the Integrity of Long-Read Genome Assemblies and Illusive Sequence Motifs

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

Kmer2SNP: Reference-Free Heterozygous SNP Calling Using k-mer Frequency Distributions

A high-precision genome size estimator based on the k-mer histogram correction

KAT: a K-mer analysis toolkit to quality control NGS datasets and genome assemblies

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

A Survey of K-mer Methods and Applications in Bioinformatics

Gencore: an efficient tool to generate consensus reads for error suppressing and duplicate removing of NGS data

Klumpy: A tool to evaluate the integrity of long‐read genome assemblies and illusive sequence motifs

Kerfuffle: a web tool for multi-species gene colocalization analysis

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

EGM: encapsulated gene-by-gene matching to identify gene orthologs and homologous segments in genomes

Analytical Performance Validation of Next-Generation Sequencing Based Clinical Microbiology Assays Using a K-mer Analysis Workflow