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

Gaëtan Benoit 1 Pierre Peterlongo 1 Dominique Lavenier 1 Claire Lemaitre 1
1 GenScale - Scalable, Optimized and Parallel Algorithms for Genomics
Inria Rennes – Bretagne Atlantique , IRISA-D7 - GESTION DES DONNÉES ET DE LA CONNAISSANCE
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 :
Poster
JOBIM 2015, Jul 2015, Clermont-Ferrand, France. <http://www6.inra.fr/jobim2015>. <10.1093/Bioinforma-cs/btu406>
Liste complète des métadonnées

https://hal.inria.fr/hal-01180603
Contributeur : Gaëtan Benoit <>
Soumis le : lundi 27 juillet 2015 - 15:32:57
Dernière modification le : mercredi 2 août 2017 - 10:08:24
Document(s) archivé(s) le : mercredi 28 octobre 2015 - 11:02:15

Fichier

poster5.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Citation

Gaëtan Benoit, Pierre Peterlongo, Dominique Lavenier, Claire Lemaitre. Simka: fast kmer-based method for estimating the similarity between numerous metagenomic datasets. JOBIM 2015, Jul 2015, Clermont-Ferrand, France. <http://www6.inra.fr/jobim2015>. <10.1093/Bioinforma-cs/btu406>. <hal-01180603>

Partager

Métriques

Consultations de
la notice

1058

Téléchargements du document

264