Powered by OpenAIRE graph
Found an issue? Give us feedback
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/ Hyper Article en Lig...arrow_drop_down
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/
addClaim

This Research product is the result of merged Research products in OpenAIRE.

You have already added 0 works in your ORCID record related to the merged Research product.

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

Authors: Benoit, Gaëtan;

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

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.

Keywords

[INFO.INFO-BI] Computer Science [cs]/Bioinformatics [q-bio.QM]

Powered by OpenAIRE graph
Found an issue? Give us feedback
Related to Research communities