More posts like this
-
Should I calculate beta diversity with or without rarefaction?
written by dpcHi there!!! I am doing metagenomic analysis with metaphlan output data that show relative abundance s of the species. I am going to calculate the …
-
Can I use pseudocounts for differential abundance of OTUs?
written by dpcHi !!! I am working with WGS metagenome data profiled by MetaPhlAn software which gives relative abundance of taxa (no read count data of the …
-
How to convert merged_abundance-table into a biom file?
written by dpcHi there!!! I have profile around 60 WGS metagenomic samples with `MetaPhlAn 3.0`. In the next step I have merged all the output profile files …
-
Import metaphlan output to phyloseq
written by dpcHi there Can anyone please tell me how can I import MetaPhlAn biom or profile file to phyloseq for subsequent statistical analysis? (I have analysed …
-
Process to merge biom files
written by dpcHi, I am analysing shotgun metagenomic sequencing data with MetaPhlAN 3.0 software. I want to do statistical analysis from the metaphlan abundance or the biom …
-
Softwares for statistitical analysis of MetaPhlAn output
written by dpcHi community!!! I have profiled shotgun metagenome data with MetaPhlAn. Now, I want to do statistical analysis (e.g. alpha-, beta- diversity comparison, rarefaction graph, boxplots, …
-
Visualise metaphlan output with KRONA
written by dpcHi biostars!!! Is there anyone using krona with metaphlan taxonomic profile output file? I am struggling with krona. The interactive pie chart is not giving …
-
Number of OTUs obtained from mothur analysis
written by dpcHi friends!!! I have analysed gut metagenome 16S rRNA data with MOTHUR package. I am getting aroung 23,000 OTUs (cut off= 0.03) in my output …
-
Mock community in MOTHUR
written by dpcHi community!!! I am running MOTHUR with already submitted metagenomic sequences in SRA. In that case I am not getting any "MOCK" community sequences (as …
-
454 GS FLX WGS data analysis
written by dpcHi friends!!! - Can single end WGS metagenome sequence data sequenced by 454 GS FLX instrument be analysed by MetaPhlAn software? - I am seeing …
Hi,
According to one of the authors, yes you should rarefy to even sampling depth. Though there are other opinions on the field. Check the discussion here: https://groups.google.com/forum/#!topic/metaphlan-users/e2tets6NHK8
António
Thanks antonioggsousa... I've already checked the discussion. I also have read this article and this one also. they are articles from Cell and plos computational biology journals. And, I am also getting some articles where rarefaction analysis has been done. So, overall I'm confused. However, Do you know. if I can do this step after the Bowtie2 alignment step or after the profiling step (in MetaPhlAn analysis)?
Thanks, DC7
Yes, there isn't a consensus. This is one of these things that there is a lot of discussion and you probably will never find a consensus solution.
I'm not familiar with MetaPhlAn, but I would say after the profiling step. I worked (and still work) with 16S rRNA gene amplicon data analysis and data transformation or rarefaction or normalization is done after you map and annotate (or profile) your OTUs/ASVs.
The rarefaction here (or normalization or transformation) is important to reduce the difference between the sequencing depth/coverage per sample if you want to perform statistical comparisons between different samples.
I hope this answers your question,
António
Thanks a lot Antonio. Yes. During analysis with MOTHUR, I also have done that after profiling. But, in MetaPhlAn, I am really getting confused when this step to be performed. And I also have to do statistical analysis here.
DC7