Jenna M. Lang edited Results PICRUSt.md  over 9 years ago

Commit id: 8e27c7c771fc65bd57c2fc5f3e8440ed118f4dd2

deletions | additions      

       

In a case like this for which metagenomic sequencing is infeasible, another approach suggests itself. There is good evidence that a correlation exists between the evolutionary relatedness of two organisms and the similarity of their genomic content \cite{Martiny_2012}. This means that we can leverage the information obtained by sequencing the genome of one organism to predict the functional potential of another, even if that other genome is represented by a single 16S rRNA sequence. The power of this approach is increased when very many, very closely-related genome sequences are available. This predictive approach has recently been implemented in the software package PICRUSt. PICRUSt uses the phylogenetic placement of a 16S rRNA sequence within a phylogeny of sequenced genomes to “reconstruct” what the genome of the organism containing that 16S rRNA sequence might look like.  With PICRUSt one can calculate a metric (NSTI) that measures how closely related the average 16S rDNA sequence in an environmental sample is to an available sequenced genome. When this number is low, PICRUSt is likely to perform well in predicting the genomes of the organisms in an environmental sample (*i.e.*, a metagenome). The average NSTI for our 15 meals was 0.038, which is on par with the NSTI for the Human Microbiome Project samples (mean NSTI = 0.03 ± 0.02 s.d.), for which a massive effort has been made to obtain reference genome sequences \cite{22018227}. This low NSTI metric suggests that PICRUSt may perform well when predicting the metaboloc potential of the microbial communities found in the meals prepared for this study. The Here, we have shown the  most significantdifference in  KEGG functional categories was category,  for “Other N-glycan degradation" (KO 00511, *p* = 8.21e-3), which was highest in the VEGAN dietary pattern (Figure 7). Also, Again, this is not a significant result  when a p-value correction is applied, but is nevertheless highlighted as a potential source of information when using a pilot study like this to inform future research questions. As a simple sanity check for the PICRUSt predictions, we compared the relative abundance of genes present in the KEGG functional category "Sporulation" between  meals that were cooked were compared to those that were raw, raw (Figure 8). As expected, because organisms that can form spores are more likely to survive  the cooked meals had cooking process, Sporulation-associated genes are  more functional categories for sporulation (Fig. 8). abundant in cooked versus raw foods.  These findings suggest that there are functional differences in bacterial populations associated with different foods and meals, and that these may be related to bacterial substrate preferences and/or techniques used in meal preparation.