I have got pair-end metagenomic data (2 x 100bp) generated by Illumina Hiseq 2000 of several samples. Human sequences have been removed and I aligned the rest reads to bacteria reference genomes. I want to calculate organism content (namely, relative abundance) of each sample so that I can use this information in further comparative metagenomics analysis.
The question is how to calculate organism relative abundance from the rest reads. Note that different samples may have different sequencing levels, normalization may be needed. I need organism relative abundance that could correct for different sequencing levels in these samples. Is there a software which can handle this?
p.s.:Qin's paper entitled 'A human gut microbial gene catalogue established by metagenomic sequencing' talks about how relative abundance was calculted, however, there is no much detail provided. It's quite confusing.
As I am building my own analysis pipline, online resources might not fit for me.
2 answers
Check out Kraken:
http://genomebiology.com/2014/15/3/R46
https://github.com/DerrickWood/kraken
Recently published out of the Salzberg lab. It's easy to run and install, if a bit resource intensive. I would be skeptical of anything mapping below species, but to the genus level it's reasonably accurate. Not to mention blazing fast.
Edit:
Kraken won't deliver immediate relative abundance, but an identification for each reference, which will likely recaputulate the relative abundance taken together.
MetaPhilAn is another tool that is specifically for relative abundance and is worth a look.
•
written
by
Brett Kennedy
You can also try the metagenomics functions from Real Time Genomics. The RTG species tool takes SAM/BAM mapped to a reference species database and calculates abundances either by abundance of dna or by organism (taking genome lengths into account). If your reference includes taxonomic information (RTG provides pre-built ones that include bacteria and viruses) this is also taken into account, and reporting is both via tsv and krona visualization. (RTG Core is available free for non-commercial academic use).
Log in to answer this question.
-
Quality control of shotgun metagenomics data
written by yesquokkan •
-
Rescaling after relative abundance filtering of gut microbiome data
written by yesquokkan •
-
Calculation of RPM and deduplication
written by Luca Arbore •
-
How to determine the relative abundance
written by serene.s •
-
Should I calculate beta diversity with or without rarefaction?
written by dpc
-
comparative gene quantification, problem with COGs and hypothetical proteins.
written by Hansen_869 •
-
Profiling taxonomic distribution of genes in human gut microbiome
written by dazhudou1122
-
Genomic Coverage - Samtool's undocumented "depth" verses the poorly documented pileup.
written by John
-
Relative abundance (evenness) estimation/prediction in metagenomic datasets
written by Tim
-
How Pca Was Done To Relate Health Status And Abundance Of Bacteria?
written by Tianyang Li •
You can try MEGAN or MG-RAST.
Thanks for your reply.
I see, maybe you can use QIIME or Kraken.
Thanks, I will have a try.
Both Kraken and Metaphlan seem good. I will try them out.
Metaphlan seems only cover Bacteria and Archaea, which is a limiting factor for its use in Viral and Fungal metagenome.
Kraken is likely to be a resource-intensive software, but it's quite fast.
Hi,
I am also doing same kind of analysis. What do you mean by "different samples may have different sequencing levels"?
Thanks,
D