Simka: fast kmer-based method for estimating the similarity between numerous metagenomic datasets

Gaëtan Benoit 1, 2
1 GenScale - Scalable, Optimized and Parallel Algorithms for Genomics
Abstract : Comparative metagenomics aims to provide high-level information based on DNA material sequenced from different environments. The purpose is mainly to estimate proximity between two or more environmental sites at the genomic level. One way to estimate similarity is to count the number of similar DNA fragments. From a computational point of view, the problem is thus to calculate the intersections between datasets of reads. Resorting to traditional methods such as all-versus-all sequence alignment is not possible on current metagenomic projects. For instance, the Tara Oceans project involves hundreds of datasets of more than 100M reads each. Maillet et al. defined the following heuristic in their method called Commet[1]. Two reads are considered similar if they share t non-overlapping kmers (words of length k). This method is currently the fastest but still does not scale on Tara Oceans samples. To tackle this issue, we introduce a new similarity function between two datasets, called Simka, based on their amount of shared kmers. To scale on large metagenomic projects, we use a new technique which is able to count the kmers of N datasets simultaneously. This method also offers new possibilities such as filtering low frequency kmers which potentially contain sequencing errors. Simka was tested and compared to Commet on 21 Tara Oceans samples. This shows that our kmer-based similarity function is very close to the read-based one of Commet. Regarding sample proximity, both methods identify the same clusters of datasets. Commet required a few weeks to compute all the intersections whereas Simka took only 4 hours.
Type de document :
Communication dans un congrès
RCAM, Oct 2015, Paris, France
Liste complète des métadonnées
Contributeur : Gaëtan Benoit <>
Soumis le : vendredi 18 décembre 2015 - 16:31:24
Dernière modification le : jeudi 17 mai 2018 - 09:59:09
Document(s) archivé(s) le : samedi 19 mars 2016 - 10:10:47


  • HAL Id : hal-01231795, version 1


Gaëtan Benoit. Simka: fast kmer-based method for estimating the similarity between numerous metagenomic datasets. RCAM, Oct 2015, Paris, France. 〈hal-01231795〉



Consultations de la notice


Téléchargements de fichiers