Removal of rare OTUs before calculating community distances / ADONIS / NMDS
0
1
Entering edit mode
5.7 years ago
dmcleanuk ▴ 10

Hi, I have a large 16s microbiome dataset with some 600 samples spanning across 5 clinical groups which consists of around 6000 OTUs. I have performed an NMDS analysis with Vegan and I would like to see if the sample distances are associated with any of my metadata such as library size/gender etc., using PERMANOVA (adonis).

I've calculated the pairwise Bray Curtis distances with vegdist:

dist.dat = vegdist(otuTable,method="bray", na.remove = TRUE)

My first question relates to whether I should be removing rare species before calculating the sample wise distances. Most of the OTUs in my OTU table are very rare, although as I understand it, bray curtis is fairly robust to this and focuses mostly on abundant species. I have included all OTUs as not to impose an arbitrary abundance threshold. Could someone comment on if this is admissible?

Next I’ve used adonis to determine if the distance between samples is associated with any of my metadata, for example library size:

adonis(dist.dat ~ otuDescription$libSize)

Call:
adonis(formula = dist.dat ~ otuDescription$libSize) 

Permutation: free
Number of permutations: 999

Terms added sequentially (first to last)

                        Df SumsOfSqs MeanSqs F.Model      R2 Pr(>F)    
otuDescription$libSize   1     3.547  3.5472  9.2281 0.01396  0.001 ***
Residuals              652   250.622  0.3844         0.98604           
Total                  653   254.169                 1.00000           
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

For which I get a significant P value. The p value is a function of sample size and thus permanova results should be discussed in context of the R2 which in my case is very low ~ 0.013. Am I correct in interpreting that this is not an important factor?

Thanks, Dave

16s microbiome Vegan NMDS adonis • 2.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 2020 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6