Computational integration of genomic traits into 16S rDNA microbiota sequencing studies

•Computational pipeline to retain genomic functional traits for microbiota.•Downstream data prepared as pseudo-metagenomes for follow-up software tools.•Exemplary application to single-reads and next-generation sequencing data.

Molecular sequencing techniques help to understand microbial biodiversity with regard to species richness, assembly structure and function. In this context, available methods are barcoding, metabarcoding, genomics and metagenomics. The first two are restricted to taxonomic assignments, whilst genomics only refers to functional capabilities of a single organism. Metagenomics by contrast yields information about organismal and functional diversity of a community. However currently it is very demanding regarding labour and costs and thus not applicable to most laboratories. Here, we show in a proof-of-concept that computational approaches are able to retain functional information about microbial communities assessed through 16S rDNA (meta)barcoding by referring to reference genomes. We developed an automatic pipeline to show that such integration may infer preliminary or supplementary genomic content of a community. We applied it to two biological datasets and delineated significantly overrepresented protein families between communities. The script alongside supporting data is available at http://bioapps.biozentrum.uni-wuerzburg.de.