Fast Approximation of Frequent k-Mers and Applications to Metagenomics.
Identifieur interne : 000572 ( Main/Exploration ); précédent : 000571; suivant : 000573Fast Approximation of Frequent k-Mers and Applications to Metagenomics.
Auteurs : Leonardo Pellegrina [Italie] ; Cinzia Pizzi [Italie] ; Fabio Vandin [Italie]Source :
- Journal of computational biology : a journal of computational molecular cell biology [ 1557-8666 ] ; 2019.
Abstract
Estimating the abundances of all k-mers in a set of biological sequences is a fundamental and challenging problem with many applications in biological analysis. Although several methods have been designed for the exact or approximate solution of this problem, they all require to process the entire data set, which can be extremely expensive for high-throughput sequencing data sets. Although in some applications it is crucial to estimate all k-mers and their abundances, in other situations it may be sufficient to report only frequent k-mers, which appear with relatively high frequency in a data set. This is the case, for example, in the computation of k-mers' abundance-based distances among data sets of reads, commonly used in metagenomic analyses. In this study, we develop, analyze, and test a sampling-based approach, called Sampling Algorithm for K-mErs approxIMAtion (SAKEIMA), to approximate the frequent k-mers and their frequencies in a high-throughput sequencing data set while providing rigorous guarantees on the quality of the approximation. SAKEIMA employs an advanced sampling scheme and we show how the characterization of the Vapnik-Chervonenkis dimension, a core concept from statistical learning theory, of a properly defined set of functions leads to practical bounds on the sample size required for a rigorous approximation. Our experimental evaluation shows that SAKEIMA allows to rigorously approximate frequent k-mers by processing only a fraction of a data set and that the frequencies estimated by SAKEIMA lead to accurate estimates of k-mer-based distances between high-throughput sequencing data sets. Overall, SAKEIMA is an efficient and rigorous tool to estimate k-mers' abundances providing significant speedups in the analysis of large sequencing data sets.
DOI: 10.1089/cmb.2019.0314
PubMed: 31891535
Affiliations:
Links toward previous steps (curation, corpus...)
- to stream PubMed, to step Corpus: 000301
- to stream PubMed, to step Curation: 000301
- to stream PubMed, to step Checkpoint: 000560
- to stream Ncbi, to step Merge: 002464
- to stream Ncbi, to step Curation: 002464
- to stream Ncbi, to step Checkpoint: 002464
- to stream Main, to step Merge: 000575
- to stream Main, to step Curation: 000572
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Fast Approximation of Frequent <i>k</i>
-Mers and Applications to Metagenomics.</title>
<author><name sortKey="Pellegrina, Leonardo" sort="Pellegrina, Leonardo" uniqKey="Pellegrina L" first="Leonardo" last="Pellegrina">Leonardo Pellegrina</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Pizzi, Cinzia" sort="Pizzi, Cinzia" uniqKey="Pizzi C" first="Cinzia" last="Pizzi">Cinzia Pizzi</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Vandin, Fabio" sort="Vandin, Fabio" uniqKey="Vandin F" first="Fabio" last="Vandin">Fabio Vandin</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PubMed</idno>
<date when="2019">2019</date>
<idno type="RBID">pubmed:31891535</idno>
<idno type="pmid">31891535</idno>
<idno type="doi">10.1089/cmb.2019.0314</idno>
<idno type="wicri:Area/PubMed/Corpus">000301</idno>
<idno type="wicri:explorRef" wicri:stream="PubMed" wicri:step="Corpus" wicri:corpus="PubMed">000301</idno>
<idno type="wicri:Area/PubMed/Curation">000301</idno>
<idno type="wicri:explorRef" wicri:stream="PubMed" wicri:step="Curation">000301</idno>
<idno type="wicri:Area/PubMed/Checkpoint">000560</idno>
<idno type="wicri:explorRef" wicri:stream="Checkpoint" wicri:step="PubMed">000560</idno>
<idno type="wicri:Area/Ncbi/Merge">002464</idno>
<idno type="wicri:Area/Ncbi/Curation">002464</idno>
<idno type="wicri:Area/Ncbi/Checkpoint">002464</idno>
<idno type="wicri:Area/Main/Merge">000575</idno>
<idno type="wicri:Area/Main/Curation">000572</idno>
<idno type="wicri:Area/Main/Exploration">000572</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en">Fast Approximation of Frequent <i>k</i>
-Mers and Applications to Metagenomics.</title>
<author><name sortKey="Pellegrina, Leonardo" sort="Pellegrina, Leonardo" uniqKey="Pellegrina L" first="Leonardo" last="Pellegrina">Leonardo Pellegrina</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Pizzi, Cinzia" sort="Pizzi, Cinzia" uniqKey="Pizzi C" first="Cinzia" last="Pizzi">Cinzia Pizzi</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Vandin, Fabio" sort="Vandin, Fabio" uniqKey="Vandin F" first="Fabio" last="Vandin">Fabio Vandin</name>
<affiliation wicri:level="1"><nlm:affiliation>Department of Information Engineering, University of Padova, Padova, Italy.</nlm:affiliation>
<country xml:lang="fr">Italie</country>
<wicri:regionArea>Department of Information Engineering, University of Padova, Padova</wicri:regionArea>
<wicri:noRegion>Padova</wicri:noRegion>
</affiliation>
</author>
</analytic>
<series><title level="j">Journal of computational biology : a journal of computational molecular cell biology</title>
<idno type="eISSN">1557-8666</idno>
<imprint><date when="2019" type="published">2019</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en">Estimating the abundances of all <i>k</i>
-mers in a set of biological sequences is a fundamental and challenging problem with many applications in biological analysis. Although several methods have been designed for the exact or approximate solution of this problem, they all require to process the entire data set, which can be extremely expensive for high-throughput sequencing data sets. Although in some applications it is crucial to estimate all <i>k</i>
-mers and their abundances, in other situations it may be sufficient to report only <i>frequent k</i>
-mers, which appear with relatively high frequency in a data set. This is the case, for example, in the computation of <i>k</i>
-mers' abundance-based distances among data sets of reads, commonly used in metagenomic analyses. In this study, we develop, analyze, and test a sampling-based approach, called Sampling Algorithm for <i>K</i>
-mErs approxIMAtion (SAKEIMA), to approximate the frequent <i>k</i>
-mers and their frequencies in a high-throughput sequencing data set while providing rigorous guarantees on the quality of the approximation. SAKEIMA employs an advanced sampling scheme and we show how the characterization of the Vapnik-Chervonenkis dimension, a core concept from statistical learning theory, of a properly defined set of functions leads to practical bounds on the sample size required for a rigorous approximation. Our experimental evaluation shows that SAKEIMA allows to rigorously approximate frequent <i>k</i>
-mers by processing only a fraction of a data set and that the frequencies estimated by SAKEIMA lead to accurate estimates of <i>k</i>
-mer-based distances between high-throughput sequencing data sets. Overall, SAKEIMA is an efficient and rigorous tool to estimate <i>k</i>
-mers' abundances providing significant speedups in the analysis of large sequencing data sets.</div>
</front>
</TEI>
<affiliations><list><country><li>Italie</li>
</country>
</list>
<tree><country name="Italie"><noRegion><name sortKey="Pellegrina, Leonardo" sort="Pellegrina, Leonardo" uniqKey="Pellegrina L" first="Leonardo" last="Pellegrina">Leonardo Pellegrina</name>
</noRegion>
<name sortKey="Pizzi, Cinzia" sort="Pizzi, Cinzia" uniqKey="Pizzi C" first="Cinzia" last="Pizzi">Cinzia Pizzi</name>
<name sortKey="Vandin, Fabio" sort="Vandin, Fabio" uniqKey="Vandin F" first="Fabio" last="Vandin">Fabio Vandin</name>
</country>
</tree>
</affiliations>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Main/Exploration
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000572 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Main/Exploration/biblio.hfd -nk 000572 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= MersV1 |flux= Main |étape= Exploration |type= RBID |clé= pubmed:31891535 |texte= Fast Approximation of Frequent k-Mers and Applications to Metagenomics. }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Main/Exploration/RBID.i -Sk "pubmed:31891535" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Main/Exploration/biblio.hfd \ | NlmPubMed2Wicri -a MersV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |