This is a test version of Biostars. For the public version, visit https://www.biostars.org.
calculating abundance of metagenomes

Hi, We performed shotgun metagenomics for dozens of stool samples. Next, we assembled the reads into contigs and binned the contigs to reconstruct the genomes. We want to use the reconstructed metagenomes in order to calculate the abundance of the microbes they represent in all samples.

  1. Is there any recommended aligner (and parameters) for such task? (considering that we align against not-full genomes, multimapping etc.)
  2. I came across some studies that calculated relative abundance by dividing the number of reads aligned to a metagenome by the total number of reads of a sample. Any Idea why genome length was not also taken into account?
  3. Is there any reason to calculate absolute abundance in such case? (maybe for downstream analyses/tools?)

Thanks!

assembly alignment
• 6,627 views • link
updated by colindaven colindaven • written by biobiu biobiu

Thanks. Yes I've seen, we think on a way to estimate abundance without using marker genes.

• 0 views • link
written by biobiu biobiu

Why not to use marker genes methods?

• 0 views • link
written by Costadia • Costadia

We assume that we have un-characterized species of bacteria and non-bacterial genomes.

• 0 views • link
written by biobiu biobiu

2 answers

https://github.com/HRGV/phyloFlashphyloFlash which uses SSU 16S against SILVA db.

• 0 views • link
written by Naren Naren

It's not really intended to answer your question, and you've probably solved this by now, but our pipeline can do the alignment and your intended normalization (including number of reads and genome length):

https://github.com/MHH-RCUG/Wochenende

Seeing as you have a custom output, you'll need to create a reference sequence (combine contigs to "genomes" by removing headers) and index it.

• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this