Bracken: Estimating species abundance in metagenomics data

Jennifer Lu, Florian P. Breitwieser, Peter Thielen, Steven L. Salzberg

Research output: Contribution to journalArticle

Abstract

Metagenomic experiments attempt to characterize microbial communities using high-throughput DNA sequencing. Identification of the microorganisms in a sample provides information about the genetic profile, population structure, and role of microorganisms within an environment. Until recently, most metagenomics studies focused on high-level characterization at the level of phyla, or alternatively sequenced the 16S ribosomalRNAgene that is present in bacterial species. As the cost of sequencing has fallen, though, metagenomics experiments have increasingly used unbiased shotgun sequencing to capture all the organisms in a sample. This approach requires a method for estimating abundance directly from the raw read data. Here we describe a fast, accurate new method that computes the abundance at the species level using the reads collected in a metagenomics experiment. Bracken (Bayesian Reestimation of Abundance after Classification with KrakEN) uses the taxonomic assignments made by Kraken, a very fast read-level classifier, along with information about the genomes themselves to estimate abundance at the species level, the genus level, or above. We demonstrate that Bracken can produce accurate species- and genus-level abundance estimates even when a sample contains multiple near-identical species.

Original languageEnglish (US)
Article numbere104
JournalPeerJ Computer Science
Volume2017
Issue number1
DOIs
StatePublished - Jan 1 2017

    Fingerprint

Keywords

  • Bayesian estimation
  • Metagenomics
  • Microbiome
  • Species abundance

ASJC Scopus subject areas

  • Computer Science(all)

Cite this