Jenna M. Lang edited Methods.md  about 10 years ago

Commit id: a0931345bef0b9f06d5964b95b72c416950163b0

deletions | additions      

       

Unless otherwise noted, all statistical analyses were performed using python scripts implemented in QIIME v. 1.8.0, and all python scripts referenced here are QIIME scripts. The IPython notebook file used for all QIIME analyses is available HERE. To explore the differences in overall microbial community composition across the 15 meals, both the phylogenetic weighted UniFrac distances (REF) and the taxonomic Bray-Curtis dissimilarities (REF) were calculated using the beta\_diversity\_through\_plots.py script. To test for the significance of dietary pattern on the overall microbial community composition, we used a permutational multivariate ANOVA as implemented in the compare\_categories.py script. To test for significant differences in taxonomic richness across dietary patterns, we used the non-parametric Kruskal-Wallis test with the FDR (false discovery rate) correction as implemented in compare\_alpha\_diversity.py. To test for the significant variation in frequency of individual OTUs across dietary patterns, we used the Kruskal-Wallis test with the FDR correction as implemented in the group_significance.py script. To test for significant correlation between the relative abundance of a single taxonomic group and meal metadata categories at 5 taxonomic levels (phylum-genus,) Pearson correlation coefficients were calculated and tested for statistical significance using Stata  ##Metagenome Prediction with PICRUSt  A synthetic metagenome was generated based on the observed 16S rDNA sequences for each meal. To do this, the 16S rDNA sequences were clustered into a collection of OTUs sharing 99% sequence identity, using the pick\_closed\_reference\_otus.py. The resultant OTU table was normalized with respect to inferred 16S rRNA gene copy numbers using the normalize\_by\_copy\_number.py script distributed with PICRUSt v. 1.0.0. The normalized OTU table was used to predict meal microbial metagenomes with PICRUSt's predict\_metagenomes.py script. The final predicted metagenome is output as a .biom table, which is suitable for analysis with a tool such as STAMP (REF). We used STAMP to test for and visualize significant (predicted) functional differences in microbial communities between the three dietary patterns.