12/09/2014

# Jenna M. Lang (jennomics@gmail.com). Genome Center, University of California, Davis, CA, USA Jonathan A. Eisen (jaeisen@ucdavis.edu). Genome Center, Evolution and Ecology, Medical Microbiology and Immunology, University of California, Davis, Davis, CA, USA Angela M. Zivkovic (amzivkovic@ucdavis.edu). Department of Nutrition, University of California, Davis, CA, USA

Abstract

#Abstract

Far more attention has been paid to the microbes in our feces than the microbes in our food. Research efforts dedicated to the microbes that we eat have historically been focused on a fairly narrow range of species, namely those which cause disease and those which are thought to confer some "probiotic" health benefit. Little is known about the effects of ingested microbial communities that are present in typical American diets, and even the basic questions of which microbes, how many of them, and how much they vary from diet to diet and meal to meal, have not been answered. We characterized the microbiota of three different dietary patterns in order to estimate: the average total amount of daily microbes ingested via food and beverages, and their composition in three daily meal plans representing three different dietary patterns. The three dietary patterns analyzed were: 1) the Average American (AMERICAN): focused on convenience foods, 2) USDA recommended (USDA): emphasizing fruits and vegetables, lean meat, dairy, and whole grains, and 3) Vegan (VEGAN): excluding all animal products. Meals were prepared in a home kitchen or purchased at restaurants and blended, followed by microbial analysis including aerobic, anaerobic, yeast and mold plate counts as well as 16S rRNA PCR survey analysis. Based on plate counts, the USDA meal plan had the highest total amount of microbes at $$1.3 X 10^9$$ CFU per day, followed by the VEGAN meal plan and the AMERICAN meal plan at $$6 X 10^6$$and $$1.4 X 10^6$$ CFU per day respectively. There was no significant difference in diversity among the three dietary patterns. Individual meals clustered based on taxonomic composition independent of dietary pattern. For example, meals that were abundant in Lactic Acid Bacteria were from all three dietary patterns. Some taxonomic groups were correlated with the nutritional content of the meals. Predictive metagenome analysis using PICRUSt indicated differences in some functional KEGG categories across the three dietary patterns and for meals clustered based on whether they were raw or cooked. Further studies are needed to determine the impact of ingested microbes on the intestinal microbiota, the extent of variation across foods, meals and diets, and the extent to which dietary microbes may impact human health. The answers to these questions will reveal whether dietary microbial approaches beyond probiotics taken as supplements - i.e., ingested as foods - are important contributors to the composition, inter-individual variation, and function of our gut microbiota.

#Introduction
The human gut microbiome (the total collection of microbes found in in the human gut) mediates many key biological functions and its imbalance, termed dysbiosis, is associated with a number of inflammatory and metabolic diseases from inflammatory bowel disease to asthma to obesity and insulin resistance (Machonkin 2014) (Costello 2012). How to effectively shift the microbiome and restore balance is a key question for disease prevention and treatment. The gut microbiome is influenced by a number of factors including the nature of the initial colonization at birth (e.g., vaginal vs. C-section delivery), host genotype, age, and diet. As diet is a readily modifiable factor, it is an obvious target for interventions. Several studies have confirmed high inter-individual variability in the bacterial composition of the gut microbiome in healthy individuals (Brownawell 2012), (Costello 2009). Despite this high variability at the species level, enterotypes, or distinct clusters at the genus level, were described as core microbiomes that are independent of age, gender, nationality, or BMI (Arumugam 2011). Although the concept of enterotypes is itself controversial, diet has been shown to play a key role in determining enterotype (Wu 2011, Filippo 2010, Muegge 2011). Although the core microbiota within each person are stable over longer time scales (e.g. 5 years), community composition is highly dynamic on shorter time scales (e.g. 0–50 weeks) (Faith 2013). In fact, major shifts occur within 1 day of a significant dietary change (Wu 2011, Turnbaugh 2009). “Blooms” in specific bacterial groups were observed in response to controlled feeding of different fermentable fibers (Walker 2011). Dietary changes affect both the structure and function of the gut microbiome in animals (Hildebrandt 2009), and humans under controlled feeding conditions (Wu 2011). Rapid shifts in microbiome composition are observed in response to change from a vegetarian to an animal based diet (David 2013).

An ecological perspective helps to delineate the complexity and multi-layered nature of the relationships between the microbiota, the human host, and both the nutritive and non-nutritive compounds we ingest (Costello 2012). The concept of the human gut microbiome as a distinct ecosystem or collection of micro-ecosystems allows us to identify and characterize the components of the system, including its inputs and outputs. In this case, the inputs of the system include all of the various ingested compounds that can either serve as food substrates (e.g. complex sugars) or that can be metabolized by or that affect the metabolism of the microbiota (e.g. polyphenolic compounds, environmental chemicals, medications). Some of these inputs, such as probiotics have been studied extensively. It has been well documented that certain sugars such as galactooligosaccharides, fructooligosaccharides, and oligosaccharides found in milk act as prebiotics that support the establishment and growth of certain commensal microbial species (Brownawell 2012), (de 2008), (Roberfroid 2007), (German 2008), (Zivkovic 2011). Research has also documented the effects of antibiotics, and pathogens on the microbiota composition, its recovery or lack of recovery to baseline following resolution, and the various immunological and physiological effects of these perturbations (Manichanh 2010), (Ubeda 2012), (Bien 2013)(Dethlefsen 2010).

Yet, little is known about the effects of ingested microorganisms on gut microbiota composition or function, and even the basic questions of which microbes, how many of them, and how much they vary from diet to diet and meal to meal, have not been answered. We do know about the microbial ecology of various specialty foods where fermentation, colonization, ripening, and/or aging are part of the preparation of these foods, for example pancetta (Busconi 2014) and of course cheese (Gatti 2008)(Button 2012). The microbial ecology of the surfaces of raw plant-derived foods such as fruits and vegetables has also been characterized (Leff 2013). There is a large base of literature on food-borne pathogens (Aboutaleb 2014). Furthermore, it is known that the microbial ecology of endemic microbes found on food surfaces can affect mechanisms by which pathogens colonize these foods (Critzer 2010). A recent article showed that certain ingested microbes found in foods such as cheese and deli meats were detected in the stool of individuals who consumed them, and that furthermore they were culturable and thus survived transit through the upper intestinal tract (David 2013). However, the microbial ecology or microbial assemblages of different meals and diets, as well as the total number of live microorganisms ingested in these meals and diets are largely unknown. In fact, studies of the effects of diets and foods on the gut microbiota rely on dietary recalls and other dietary reporting instruments that were not designed to capture the potential variability in aspects of foods other than their basic macronutrient and micronutrient content. Specifically, current instruments for collecting individual dietary data do not capture the provenance of foods or their preparation, both of which would likely influence certain compositional aspects of the foods, especially the microbes on those foods.

We performed a preliminary study designed to generate hypotheses about the microbes we eat, and how they vary in terms of total abundance and relative composition in different meals and dietary patterns typical of American dietary intakes. We have selected to characterize the microbiota of 15 meals that exemplify the typical meals consumed as part of three different dietary patterns in order to determine the average total amount of daily microbes ingested via food and beverages and their composition in the average American adult consuming these typical foods/diets: 1) the Average American dietary pattern (AMERICAN) focused on convenience foods, 2) the USDA recommended dietary pattern (USDA) emphasizing fruits and vegetables, lean meat, dairy, and whole grains, and 3) the Vegan (VEGAN) dietary pattern, which excludes all animal products. We used DNA sequencing, plate counting, and informatics methods to characterize microbes in these meals and dietary patterns.

#Methods
##Meal preparation
We conducted a series of experiments consisting of food preparation followed by sample preparation and microbial analysis. Food was purchased and prepared in a standard American home kitchen by the same individual using typical kitchen cleaning practices including hand washing with non-antibacterial soap between food preparation steps, washing of dishes and cooking instruments with non-antibacterial dish washing detergent, and kitchen clean-up with a combination of anti-bacterial and non-antibacterial cleaning products. Anti-bacterial products had specific anti-bacterial molecules added to them whereas "non-antibacterial" products were simple surfactant-based formulations. The goal was to simulate a typical home kitchen rather than to artificially introduce sterile practices that would be atypical of how the average American prepares their meals at home. All meals were prepared according to specific recipes (from raw ingredient preparation such as washing and chopping, to cooking and mixing).

After food preparation, meals were plated on a clean plate, weighed on a digital scale (model 157W, Escali, Minneapolis, MN), and then transferred to a blender (model 5200, Vita-Mix Corporation, Cleveland, OH) and processed until completely blended (approximately 1-3 minutes). Prepared, ready to eat foods that were purchased outside the home were simply weighed in their original packaging and then transferred to the blender. 4 mL aliquots of the blended meal composite were extracted from the blender, transported on dry ice and then stored at -80°C until analysis. The following analyses were completed using these meal composite samples: 1) total aerobic bacterial plate counts, 2) total anaerobic bacterial plate counts, 3) yeast plate counts, 4) fungal plate counts, and 5) 16S rDNA analysis for microbial ecology.

##Diet design
Diets were designed by a nutritional biologist to deliver the average number of Calories consumed by an average American per day. The average American woman is 63 inches in height and weighs 166 pounds, and the average American man is 69 inches in height and weighs 195 pounds with an average age of 35, National Health and Nutrition Examination Survey, which translates to a total daily Calorie intake range of 2,000 to 2,600 Calories per day respectively to maintain weight, as determined using the [USDA MyPlate SuperTracker tool](https://www.supertracker.usda.gov/MyWeightManager.aspx. 2010). Therefore an intermediate daily Calorie intake of about 2,200 Calories was chosen as the target.

Meal plans were created using the NutriHand program (Nutrihand Inc., Soraya, CA). Diet nutrient composition was calculated by the NutriHand program from reference nutrient data for individual foods using the USDA National Nutrient Database for Standard Reference. Three one-day meal plans were created to be representative of three typical dietary patterns that are consumed by Americans: 1) the Average American dietary pattern (AMERICAN), which includes meat and dairy and focuses on convenience foods, 2) the USDA recommended dietary pattern (USDA), which emphasizes fresh fruits and vegetables, lean meats, whole grains and whole grain products, and dairy, and 3) the Vegan dietary pattern (VEGAN), which excludes all animal products. The AMERICAN meal plan totaled 2268 Calories, which consisted of 35% fat, 53% carbohydrates of which 16.6 g was fiber, and 12% protein. The USDA meal plan totaled 2260 Calories, consisting of 25% fat, 49% carbohydrates of which 45 g was fiber, and 27% protein. The VEGAN meal plan totaled 2264 Calories and consisted of 31% fat, 54% carbohydrates of which 52 g was fiber, and 15% protein.

##Microbial community analysis
Microbial plate counts were performed by Covance Laboratories (Covance Inc., Madison, WI). Aerobic plate counts were performed according to SPCM:7, anaerobic plate counts were performed according to APCM:5 and the yeast and mold counts were performed according to Chapter 23 of the FDA's Bacteriological Analytical Manual. Plate counts were reported as colony forming units (CFU) per gram for each meal composite. The CFU/g values were multiplied by the total number of grams in each meal to obtain the CFU per meal, and the values for meals for each day were added to obtain the CFU per day for each dietary pattern (Table 1).

The taxonomic composition of each meal microbiome was assessed via amplification and sequencing of 16S rDNA from the homogenized meals. DNA was extracted from homogenized food samples with the Power Food Microbial DNA Isolation Kit (MoBio Laboratories, Inc.) according to the manufacturer’s protocol. Microbial DNA was amplified by a two-step PCR enrichment of the 16S rRNA gene (V4 region) using primers 515F and 806R, modified by addition of Illumina adaptor and barcodes sequences. All primer sequences and a detailed PCR protocol are provided in Table 2 and in a GitHub repository (https://github.com/hollybik/protocols/blob/master/16S_rRNA_twostep_PCR.tex), respectively. Libraries were sequenced using an Illumina MiSeq system, generating 250bp paired-end amplicon reads. The amplicon data was multiplexed using dual barcode combinations for each sample. We used a custom script (available in a GitHub repository (https://github.com/gjospin/scripts/blob/master/Demul_trim_prep.pl), to assign each pair of reads to their respective samples when parsing the raw data. This script allows for 1 base pair difference per barcode. The paired reads were then aligned and a consensus was computed using FLASH (Magoč 2011) with maximum overlap of 120 and a minimum overlap of 70 (other parameters were left as default). The custom script automatically demultiplexes the data into fastq files, executes FLASH, and parses its results to reformat the sequences with appropriate naming conventions for QIIME v. 1.8.0 (Caporaso 2010) in fasta format. The resulting consensus sequences were analyzed using the QIIME pipeline.

##Statistical Analyses
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 at http://nbviewer.ipython.org/gist/jennomics/c6fe5e113525c6aa8add. To explore the differences in overall microbial community composition across the 15 meals, both the phylogenetic weighted UniFrac distances (Lozupone 2011) and the taxonomic Bray-Curtis dissimilarities (Bray 1957) were calculated using the beta_diversity_through_plots.py script. This script also produced a principal coordinates analysis (PCoA) plot in which the Bray-Curtis dissimilarities between samples were used to visualize differences among groups of samples (see Figure 1 for this type of visualization for the three Diet Types.) 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 (Kruskal 1952) 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. We also used the biplot function of the make_emperor.py script to plot the family-level OTUs in PCoA space alongside each meal. To test for significant correlation between the relative abundance of a single taxonomic group and meal metadata categories (i.e. nutrient composition, whether a meal contains fermented foods, etc) at 5 taxonomic levels (phylum-genus) Pearson correlation coefficients (Pearson 1895) were calculated and tested for statistical significance using Stata (StataCorp. 2013. Stata Statistical Software: Release 13. College Station, TX: StataCorp LP.)

##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 script. 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 (Langille 2013). 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 (Parks 2014). We used STAMP to test for and visualize significant (predicted) functional differences in microbial communities between the three dietary patterns.

#Results
##Meal Composition
The detailed meal plans with all ingredients are shown in Table 3, food preparation descriptions (all steps prior to placing into blender and blending foods as described in methods section) are shown in Table 4, and nutrient values based on the USDA nutrient database are shown in Table 5. The AMERICAN meal plan consisted of a large Starbucks Mocha Frapuccino for breakfast, a McDonald’s Big Mac, French fries, and Coca Cola for lunch, Stouffer’s lasagna for dinner, and Oreo cookies for a snack. The USDA meal plan consisted of cereal with milk and raspberries for breakfast, an apple and yogurt for a morning snack, a turkey sandwich on whole wheat bread with salad (including a hard-boiled egg, grapes, parmesan cheese, and Ceasar dressing) for lunch, carrots, cottage cheese and chocolate chips for an afternoon snack, and chicken, asparagus, peas and spinach on quinoa for dinner. The VEGAN meal plan consisted of oatmeal with banana, peanut butter, and almond milk for breakfast, a protein shake (including vegetable-based protein powder, soy milk, banana and blueberries) for a morning snack, a vegetable and tofu soup (including soba noodles, spinach, carrots, celery and onions in vegetable broth) for lunch, an apple and almonds with tea for an afternoon snack, a Portobello mushroom burger (including Portobello mushroom, avocado, tomato, lettuce, and a whole wheat bun) with steamed broccoli for dinner, and popcorn, hazelnuts and fig bars for an evening snack.

The following meals contained fermented foods that contained live active cultures according to the package and were prepared without heat treatment: USDA meal plan snack #1 (yogurt), lunch (parmesan cheese), and snack #2 (cottage cheese). The following meals contained fermented foods that were cooked as part of meal preparation: VEGAN meal plan lunch (tofu), and AMERICAN meal plan lunch and dinner (cheese). Meal ingredients were purchased at local grocery stores in Saint Helena, CA, and prepared meals were purchased in restaurants in Napa, CA.

##Plate Counts
The aerobic, anaerobic, yeast and mold plate counts are shown in Table 1. The meals ranged in total numbers of microorganisms from $$2.2 X 10^4$$ CFU to $$1.2 X 10^9$$ CFU with the aerobic and anaerobic plate counts being among the highest and the yeast and mold plate counts being among the lowest across all meals. The USDA dietary pattern had the highest total microorganisms for the day at $$1.3 X 10^9$$ CFU mostly due to the higher amounts of anaerobic bacteria in the morning snack ($$5.4 X 10^7$$ CFU) and higher amounts of aerobic and anaerobic bacteria in the afternoon snack (5.5 X 108 and 6 X 108 CFU respectively). Not surprisingly, both of these meals contained fermented products, in the first case yogurt, and in the second case cottage cheese. The AMERICAN and VEGAN dietary patterns had 3 orders of magnitude less total microorganisms than the USDA dietary patterns, with total microorganisms of $$1.4 X 10^6$$ CFU and $$6 X 10^6$$ CFU respectively. Neither the AMERICAN nor the VEGAN dietary pattern meals contained fermented foods that were not heat treated as part of meal preparation. The AMERICAN lunch and dinner contained cheese that was either cooked on a grill or baked in the oven and the VEGAN lunch contained tofu, which was cooked in the vegetable broth. The USDA lunch also had the highest amounts of yeast and mold ($$1 X 10^6$$ and $$2.6 X 10^6$$ CFU respectively) of all the meals, and this meal also had relatively high amounts of aerobic bacteria ($$3.3 X 10^7$$ CFU). In the VEGAN dietary pattern, the morning snack had the highest amounts of aerobic and anaerobic bacteria ($$2 X 10^6$$ and $$1.7 X 10^6$$ CFU respectively).

##Sequence Processing and Summary Statistics
The number of high-quality sequences per sample (i.e., meal) ranged from 168,669 to 318,956 (see Table 6). Sequences were clustered and clusters were assigned to a taxonomic group (when possible) using the pick_open_reference_otus.py script with a 97% similarity cutoff and the gg_13_8_otus reference taxonomy provided by the Greengenes Database Consortium (http://greengenes.secondgenome.com.) After OTU assignment, mitochondrial and chloroplast sequences were filtered out, sequences that were observed only once across all samples were removed, and sequences that were Unassigned at the Domain taxonomic level were removed (these Unassigned sequences were verified via a manual BLAST search to be chloroplast sequences.) After this filtration, the range of sequences per sample decreased to 771-244,597. All subsequent beta diversity analyses (comparisons across samples) were performed on samples that were rarefied to 771 sequences per sample.

##Taxonomic Composition and Diversity of the Different Dietary Patterns
In terms of taxonomic alpha diversity, there was no significant difference between dietary patterns (Figure 2) (non-parametric Kruskal-Wallis test with compare_alpha_diversity.py, p>0.6). This is the case for multiple diversity metrics, including a count of the absolute number of OTUs observed, as well as the Chao1 and Shannon-Weiner (parametric and nonparametric, respectively) diversity indices, which account for the relative abundance (evenness) of the OTUs observed. We also tested for the significant variation in frequency of individual OTUs between diet types using the Kruskal-Wallis test, as implemented in the group_significance.py script. This test is appropriate for comparing independent groups, with unequal sample sizes, that may not be normally distributed. None of the OTUs were significantly different between the three diet types. The most abundant 50 OTUs (clustered at 97% similarity) belong to 25 different bacterial families, including many that are commonly found in association with plants and animals (see Figure 3).

##Factors Driving the Differences in Microbial Community Composition and Diversity of Individual Meals

There was no effect of dietary pattern on the overall community composition within individual meals (PERMANOVA with compare_categories.py, p=0.591). There was no obvious clustering based on any potentially distinguishing feature tested, including whether the meals contained fermented foods, dairy, whether they ware raw or cooked, or the calculated nutritional content (see Table 5 for complete meal metadata.) However, different meals clustered together independent of dietary pattern. For example, meals that were relatively abundant in Prevotellaceae included the USDA dinner, VEGAN dinner, AMERICAN dinner, USDA breakfast, VEGAN snack 2, and VEGAN snack 3. Prevotellaceae includes organisms that tend to be very abundant in the guts of many animals, and have been associated with Inflammatory Bowel Disease in humans (Henderson 2013)(Wu 2013). The AMERICAN snack, AMERICAN lunch, USDA snack 2 and USDA snack 1 had a high relative abundance of Streptococcaceae (Figure 4). It is difficult to know what specific features of these meals made them similar in this regard. Possible contributing factors may be provenance of ingredients and/or individual meal components such as presence of a certain fruit or vegetable.

A large amount of variation (52%) was explained by PCo1 (i.e., the eigenvector that explains the most variation)(Figure 1), and our attempts to determine the factors driving this variation led us to look at specific taxonomic groups that may be important. We did this in two ways. First, we produced a biplot with the make_emperor.py script, showing the prevalence of bacterial families in the PCoA space defined by the weighted unifrac distance between the 15 meals (Figure 5). A cluster of 4 meals, including USDA snack 1, AMERICAN snack, AMERICAN lunch, and USDA lunch, was comprised of samples that were dominated by Lactic Acid Bacteria. These are members of the order Lactobacillales, which are commonly found in association with both food products, especially in fermented milk products (Terzic-Vidojevic 2014) and human mucosal surfaces (Rizzello 2011). The Vegan Snack #1 was unique in that it was dominated (70.4%) by Xanthamonadaceae, a family containing many plant pathogens. A second cluster of 7 meals including VEGAN dinner, VEGAN breakfast, AMERICAN breakfast, AMERICAN dinner, USDA dinner, VEGAN snack 3, and USDA breakfast, was comprised of samples containing a large percentage (average = 27%) of Thermus, a clade with many heat and dessication-resistant organisms. Second, we calculated correlations between the relative abundance of a single taxon and the PCo1 value for each meal using a simple regression (see Table 7). The bacterial family most tightly correlated with PCo1 was Streptococcaceae (r = 0.852).

We also asked whether the relative abundance of any particular taxonomic group was correlated with the nutritional content of the meals via pairwise Pearson's correlations. We limited this analysis to organisms that were present in all 15 meals. Due to the exploratory nature of this study, there were no specific hypotheses tested with these correlations, and therefore no corrections for multiple hypothesis testing. Some taxa frequently abundant in human microbiome studies were found to be significantly correlated with particular nutrients (Table 7). For example, members of the genus Blautia are frequently observed in human fecal samples, and in our study, the relative abundance of this genus was found to be positively correlated with the sugar content of the meals p<0.05 (Figure 6). We emphasize here that due to the large numbers of OTUs present in this study, corrected p values were always non-significant. However, the goal of this small-scale study is to inform the development of future hypotheses, not test current ones. Nevertheless, this result suggests that there could be interesting relationships between the nutritional content of the foods that we eat, the microbes that associate with those foods, and our gut microbiome, not just because we are "feeding" our gut microbes, but because we are eating them as well (but, see Caveats section below.)

##Metagenome Prediction with PICRUSt
Because of the vast, historical effort to make the 16S rRNA gene sequence available for hundreds of thousands of organisms, we are typically able to characterize well the taxonomic diversity of most microbial communities. One might assume that each organism present in a community has some functional role to play, and the most straightforward way to predict what that role each organism might play is to use metagenomic sequencing to interrogate the genomes of all members of the community. Unfortunately, in many cases and with current sequencing technology, the amount of microbial DNA relative to host or other environmental DNA is small enough to make metagenomic sequencing infeasible. This is the case here, where the plant and animal DNA present in the food we eat is typically much more abundant than the microbial DNA. Some exceptions may exist with respect to fermented foods, but we are equally interested the microbiota associated with a wide variety of food types.

In a case like this for which metagenomic sequencing is infeasible, another approach suggests itself. There is evidence that a correlation exists between the evolutionary relatedness of two organisms and the similarity of their genomic content (Martiny 2012). This allows us to leverage the information obtained by sequencing the genome of one organism to predict the functional potential of another, even if the other genome is represented only by a 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 infer the content of the genome of the organism represented by that 16S rRNA sequence.

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 (Proctor 2011). 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. Here, we have shown the most significant KEGG functional category, for “Other N-glycan degradation" (KO 00511, p = 8.21e-3), which was highest in the VEGAN dietary pattern (Figure 7). 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 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 (Figure 8). As expected, because organisms that can form spores are more likely to survive the cooking process, Sporulation-associated genes are more abundant in cooked versus raw foods. All KEGG (Level 3) pathways that vary significantly between dietary patters are presented in Table 8.) These findings suggest that there are functional differences in bacterial populations associated with different foods and meals, and that these may be related not only to bacterial substrate preferences, but also techniques used in meal preparation.

#Summary

In this study we estimated the total numbers and kinds of microorganisms consumed in a day by an average American adult. We analyzed meals representing three typical dietary patterns, including the Average American, USDA recommended, and Vegan diet, and found that Americans likely consume in the range of $$10^6$$$$10^9$$ CFU microbes per day. The USDA meal plan included two meals with non-heat treated fermented foods, which were likely responsible for the 3-fold higher total microbes in this meal plan compared to the AMERICAN and VEGAN diets. Food preparation techniques such as heating or acid treatment can kill bacteria, however, these processes affect different bacteria to different degrees. For example, spore-forming bacteria can survive heat treatment (Stringer 2011) (also see Figure 8). Once inside the gastrointestinal tract, the low pH of the stomach, as well as bile salts also kill some bacteria, but not those that are acid and/or bile salt resistant. It is unknown what proportion of the microbes we eat make it through the hostile environment of the gastrointestinal tract. However, a recent study showed that food microbes consumed as part of fermented foods such as cheese did appear in the stool and were culturable (David 2013).

We also used PICRUSt to predict the functional potential of the microbiota associated with each meal in this study. Of course, this is not a perfect substitute for metagenomic sequencing or experimental studies, but it does allow one to develop some initial hypotheses related to the function of a microbial community. For example, between diet types, the most significant difference in KEGG functional categories was for “other N-glycan degradation.” This function was over-represented in the Vegan diet, which is perhaps not surprising given that cellulose is a glycan, and the Vegan diet is significantly higher in cellulose than the others. This suggests that when one consumes a diet that is high in cellulose, one also consumes a population of microbes that is well equipped to digest cellulose.

##Caveats
It is important to point out some caveats with regard to this study. First, the scale of this study was limited. Our objective with this preliminary study was to explore the possiblity that the microbes in our food may be meaningful contributors to the ecosystem of our gut microbiota. The current paradigm is that the harsh conditions of the human gastrointestinal tract (i.e. high acidity in the stomach, presence of bile acids and digestive enzymes in the small intestine) preclude most microbes present in and on food from playing significant roles in the gastrointestinal microbial ecosystem because they do not survive intact. Yet recent reports show that there are microbial blooms within 24 hours of large shifts in the diet (Wu 2011)(Walker 2010), attributed to changes in the available fermentable substrate, which in turn promoted the growth of specific bacteria already present in the gut. The possiblity remains that microbes present on the food itself contributed at least in part to these observed transient blooms or shifts. Very little is known about how different food matrices may promote the survival of food bacteria despite their lack of acid and bile acid resistance in vitro. Transient shifts in gut microbiome composition caused by dietary factors may be unimportant when the gut microbiome rebounds to its initial state. However, it is possible that under certain conditions, such perturbations could lead to functional and long-lasting changes. The goal of this study was not to explore the explicit effects of ingesting food microbes on changes in the gut microbiota but simply to provide preliminary data on the composition and numbers of bacteria in typical American dietary patterns.

Second, the study did not aim to be exhaustive in its exploration of diets. We did not aim to produce statistically significant differences in microbial composition and quantity in replicates across multiple days of a particular dietary pattern. Instead, the aim was to generate hypotheses about dietary microbes and their variation across meals and diets that can now be followed up with more rigorous studies. Because the purpose of this study was to generate hypotheses rather than test specific hypotheses, multiple testing corrections were not applied in our statistical analyses.

Third, the microbial counts reported here are rough estimates of the total amount of microbes consumed in a day by an average American eating meals described here. It is important to note that inherent to the plate count techniques used, not every microbe will grow under these culture conditions, and the plate counts are only estimates.

##Future Directions
It is possible that part of the high variation in gut microbiota composition observed among individuals is due to the specific and complex differences in diet beyond the nutrient composition that can be estimated from dietary records and recalls. This study begs the question: do the microbes we eat as part of our normal daily diets contribute to the composition and function of our gut microbiota? There are many questions that remain to be answered. Under what circumstances do microbes consumed as part of meals remain in the gastrointestinal tract transiently versus persistently following a meal? Do the microbes we eat affect the function of the resident gut microbiota, even if they do not affect its composition, as has been suggested by some yogurt feeding studies (McNulty 2011)? How do different cooking and preparation methods affect the microbial composition of meals and the survival characteristics of individual microbes through the gastrointestinal tract? How do specific factors such as length of transport or provenance of individual ingredients (e.g. imported vs. domestic), packaging materials, and handling of ingredients in homes alter the microbial composition of foods? The findings of this study suggest that the microbes we eat as part of normal diets vary in absolute abundance, community composition, and functional potential. This variation depends on the specific ingredients in the meals, whether and how the foods are prepared and processed, and other potential factors, not explored here, including the provenance of ingredients. The significance of this variation on the gut microbiota composition and function, and its impact on human health remains to be elucidated. In addition, much as certain gut microbes can transform and modify dietary constituents and nutrients such as polyphenolic compounds and vitamins in the gut (Tuohy 2012), it is possible that food microbes similarly modify nutritive molecules. Future studies need to explore these questions in rigorous study designs aimed at addressing key questions about the composition and content of food microbes and how these vary across diets and meals, and their impacts on the short term and long term composition and function of the gut microbiota.

All data underlying the results presented here can be found on Figshare (http://figshare.com/account/projects/2225).

Principle Coordinates Analysis (PCoA) based on Bray-Curtis dissimilarities of microbial communities found in the 15 meals, color coded according to the dietary patterns they represent. Axes are scaled to the amount of variation explained.

While some individual meals had higher alpha diversity (defined either by the number of OTUs observed or by the Chao1 and Shannon diversity measures) than others, there was no significant difference in diversity between the different dietary patterns (AMERICAN, USDA, and VEGAN).

Cumulative relative abundances, across the three diet types, of the top 50 OTUs observed in this study.
Heatmap showing relative abundance of bacterial families of individual meals. Similarities between meals are not necessarily part of the same dietary pattern. Hierarchical clustering is based on Ward clustering of the Pearson correlation coefficients, with sample by sample normalization performed using the median.

Bacterial families (light blue spheres) are displayed in a PCoA biplot based on weighted unifrac distances between the 15 meals. The size of the spheres representing taxa is correlated with the relative abundance of the labeled organism. In the interest of readability, only the bacterial families discussed in the text are labeled. Axes are scaled to amount of variation explained.

Scatterplot of the relative abundance of Blautia versus grams of sugar in each sample (i.e., meal). Pearson's r= 0.56

The abundance of the "Other glycan degradation" (KO00155) pathway found in the 15 meals when grouped by dietary pattern. Due to the exploratory nature of this data set corrections for multiple testing were not applied.

Based on the PICRUSt-predicted metagenomes, the cooked meals contained significantly more genes associated with sporulation than raw meals (ANOVA p = 0.039).

#Table 1: Microbial Plate Counts
| | | | | | | |
|------------------------------------------------------------------------------|-----------|---------------------|-----------------------|-------------|------------|----------------------|
| Dietary Pattern | Meal | Aerobic Plate Count | Anaerobic Plate Count | Yeast Count | Mold Count | Total Microorganisms |
| Average American | Breakfast | 2.15E+05 | 2.26E+05 | 5.66E+02 | 5.66E+03 | 4.48E+05 |
| | Lunch | 2.23E+05 | 1.31E+04 | 1.31E+03 | 1.31E+03 | 2.38E+05 |
| | Snack | 1.87E+04 | 2.34E+03 | 2.34E+02 | 2.34E+02 | 2.15E+04 |
| | Dinner | 1.47E+05 | 5.35E+05 | 7.75E+02 | 7.75E+02 | 6.84E+05 |
| | Total | 6.04E+05 | 7.77E+05 | 2.88E+03 | 7.98E+03 | 1.39E+06 |
| USDA Recommended | Breakfast | 1.14E+04 | 5.72E+02 | 4.29E+04 | 1.49E+06 | 1.54E+06 |
| | Snack #1 | 2.11E+05 | 5.42E+07 | 5.42E+02 | 1.19E+05 | 5.45E+07 |
| | Lunch | 3.25E+07 | 2.26E+06 | 1.06E+06 | 2.55E+06 | 3.84E+07 |
| | Snack #2 | 5.54E+08 | 6.09E+08 | 3.32E+04 | 1.39E+04 | 1.16E+09 |
| | Dinner | 3.49E+05 | 5.81E+04 | 9.69E+02 | 9.69E+03 | 4.17E+05 |
| | Total | 5.87E+08 | 6.66E+08 | 1.14E+06 | 4.18E+06 | 1.26E+09 |
| Vegan | Breakfast | 3.38E+04 | 1.99E+04 | 3.98E+02 | 9.95E+03 | 6.41E+04 |
| | Snack #1 | 1.97E+06 | 1.67E+06 | 5.41E+05 | 4.92E+05 | 4.67E+06 |
| | Lunch | 1.22E+05 | 2.34E+04 | 9.35E+02 | 9.35E+02 | 1.47E+05 |
| | Snack #2 | 9.81E+04 | 4.67E+03 | 3.50E+04 | 2.80E+05 | 4.18E+05 |
| | Dinner | 4.07E+05 | 1.45E+05 | 4.53E+02 | 9.05E+03 | 5.62E+05 |
| | Snack #3 | 1.43E+05 | 2.94E+03 | 8.40E+01 | 8.40E+01 | 1.46E+05 |
| | Total | 2.77E+06 | 1.87E+06 | 5.78E+05 | 7.92E+05 | 6.01E+06 |

#Table 2: PCR primer constructs used for Illumina MiSeq library preparation
| Name | Illumina adapter sequence | barcode | pad | linker | primer | Complete Oligo |
|---------------|-------------------------------|----------|-------------|--------|----------------------|-----------------------------------------------------------------------|
| 16S_515_bc1 | AATGATACGGCGACCACCGAGATCTACAC | AACCAGTC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACAACCAGTCTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc2 | AATGATACGGCGACCACCGAGATCTACAC | AACGCTAA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACAACGCTAATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc3 | AATGATACGGCGACCACCGAGATCTACAC | AAGACTAC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACAAGACTACTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc4 | AATGATACGGCGACCACCGAGATCTACAC | AATCGATA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACAATCGATATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc5 | AATGATACGGCGACCACCGAGATCTACAC | ACCAATTG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACACCAATTGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc6 | AATGATACGGCGACCACCGAGATCTACAC | ACTGAAGT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACACTGAAGTTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc7 | AATGATACGGCGACCACCGAGATCTACAC | ATTGCCGC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACATTGCCGCTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc8 | AATGATACGGCGACCACCGAGATCTACAC | CAACCTTA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCAACCTTATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc9 | AATGATACGGCGACCACCGAGATCTACAC | CCTAATAA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCCTAATAATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc10 | AATGATACGGCGACCACCGAGATCTACAC | CCTCTGAT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCCTCTGATTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc11 | AATGATACGGCGACCACCGAGATCTACAC | CGGTCGAG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCGGTCGAGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc12 | AATGATACGGCGACCACCGAGATCTACAC | CTAATGGC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCTAATGGCTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc13 | AATGATACGGCGACCACCGAGATCTACAC | CTCATGCG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACCTCATGCGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc14 | AATGATACGGCGACCACCGAGATCTACAC | GAACGGAG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGAACGGAGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc15 | AATGATACGGCGACCACCGAGATCTACAC | GCCTACGC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGCCTACGCTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc16 | AATGATACGGCGACCACCGAGATCTACAC | GCGTTACC | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGCGTTACCTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc17 | AATGATACGGCGACCACCGAGATCTACAC | GGAGGCTG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGGAGGCTGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc18 | AATGATACGGCGACCACCGAGATCTACAC | GGATGCCA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGGATGCCATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc19 | AATGATACGGCGACCACCGAGATCTACAC | GGATTAGG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGGATTAGGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc20 | AATGATACGGCGACCACCGAGATCTACAC | GTTGGCCG | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACGTTGGCCGTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc21 | AATGATACGGCGACCACCGAGATCTACAC | TATTAACT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACTATTAACTTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc22 | AATGATACGGCGACCACCGAGATCTACAC | TGACTGCT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACTGACTGCTTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc23 | AATGATACGGCGACCACCGAGATCTACAC | TGGCGATT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACTGGCGATTTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc24 | AATGATACGGCGACCACCGAGATCTACAC | TTCAGCGA | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACTTCAGCGATATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_515_bc25 | AATGATACGGCGACCACCGAGATCTACAC | TTGGCTAT | TATGGTAATTG | TG | TGCCAGCMGCCGCGGTAA | AATGATACGGCGACCACCGAGATCTACACTTGGCTATTATGGTAATTGTGTGCCAGCMGCCGCGGTAA |
| 16S_806R_bc1 | AATGATACGGCGACCACCGAGATCTACAC | AACCAGTC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACAACCAGTCAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc2 | AATGATACGGCGACCACCGAGATCTACAC | AACGCTAA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACAACGCTAAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc3 | AATGATACGGCGACCACCGAGATCTACAC | AAGACTAC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACAAGACTACAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc4 | AATGATACGGCGACCACCGAGATCTACAC | AATCGATA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACAATCGATAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc5 | AATGATACGGCGACCACCGAGATCTACAC | ACCAATTG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACACCAATTGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc6 | AATGATACGGCGACCACCGAGATCTACAC | ACTGAAGT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACACTGAAGTAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc7 | AATGATACGGCGACCACCGAGATCTACAC | ATTGCCGC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACATTGCCGCAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc8 | AATGATACGGCGACCACCGAGATCTACAC | CAACCTTA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCAACCTTAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc9 | AATGATACGGCGACCACCGAGATCTACAC | CCTAATAA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCCTAATAAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc10 | AATGATACGGCGACCACCGAGATCTACAC | CCTCTGAT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCCTCTGATAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc11 | AATGATACGGCGACCACCGAGATCTACAC | CGGTCGAG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCGGTCGAGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc12 | AATGATACGGCGACCACCGAGATCTACAC | CTAATGGC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCTAATGGCAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc13 | AATGATACGGCGACCACCGAGATCTACAC | CTCATGCG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACCTCATGCGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc14 | AATGATACGGCGACCACCGAGATCTACAC | GAACGGAG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGAACGGAGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc15 | AATGATACGGCGACCACCGAGATCTACAC | GCCTACGC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGCCTACGCAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc16 | AATGATACGGCGACCACCGAGATCTACAC | GCGTTACC | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGCGTTACCAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc17 | AATGATACGGCGACCACCGAGATCTACAC | GGAGGCTG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGGAGGCTGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc18 | AATGATACGGCGACCACCGAGATCTACAC | GGATGCCA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGGATGCCAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc19 | AATGATACGGCGACCACCGAGATCTACAC | GGATTAGG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGGATTAGGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc20 | AATGATACGGCGACCACCGAGATCTACAC | GTTGGCCG | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACGTTGGCCGAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc21 | AATGATACGGCGACCACCGAGATCTACAC | TATTAACT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACTATTAACTAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc22 | AATGATACGGCGACCACCGAGATCTACAC | TGACTGCT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACTGACTGCTAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc23 | AATGATACGGCGACCACCGAGATCTACAC | TGGCGATT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACTGGCGATTAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc24 | AATGATACGGCGACCACCGAGATCTACAC | TTCAGCGA | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACTTCAGCGAAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |
| 16S_806R_bc25 | AATGATACGGCGACCACCGAGATCTACAC | TTGGCTAT | AGTCAGTCAG | CC | GGACTACHVGGGTWTCTAAT | AATGATACGGCGACCACCGAGATCTACACTTGGCTATAGTCAGTCAGCCGGACTACHVGGGTWTCTAAT |

#Table 3: Food ingredients for each meal within each dietary pattern
| Average American | | | USDA | | | Vegan | | | |
|------------------ |---------- |---------------------------- |------- |---------- |----------------------------------- |------- |----------- |-------------------------------- |------- |
| | Amount | Item | grams | Amount | Item | grams | Amount | Item | grams |
| Breakfast | 20 oz | Starbucks Mocha Frappucino | 566 | 2 cups | Kashi GoLean cereal | 88 | 0.5 each | large banana | 60 |
| | | | | 1 cup | 1% milk | 232 | 1 cup | Almond Breeze almond milk | 250 |
| | | | | 0.5 cup | raspberries | 58 | 2 tsp | maple syrup | 14 |
| | | | | | | | 1 tbsp | peanut butter | 28 |
| | | | | | | | 0.5 cup | rolled oats | 46 |
| Snack #1 | | | | 1 each | small Fuji apple | 163 | 0.5 each | large banana | 62 |
| | | | | 6 oz | Yoplait strawberry yogurt | 162 | 1 cup | soy milk | 250 |
| | | | | | | | 1 scoop | Spirutein protein powder | 38 |
| | | | | | | | 1 cup | blueberries (Chile) | 142 |
| Lunch | 1 each | McDonald's Big Mac | 215 | 2 tbsp | Cesar dressing | 26 | 6 oz | firm tofu | 171 |
| | 1 large | McDonald's French Fries | 154 | 20 each | green seedless grapes | 125 | 2 oz | soba noodles | 57 |
| | 12 fl oz | McDonald's Coke | 380 | 1 each | Oroweat whole wheat burger bun | 78 | 1 cup | spinach | 28 |
| | | | | 3 cups | green leaf lettuce | 72 | 1 each | medium carrot | 71 |
| | | | | 1 each | large hard boiled egg | 52 | 2 cups | Pacific Foods vegetable broth | 480 |
| | | | | 3 tbsp | parmesan cheese, shredded | 18 | 1 stalk | medium celery | 56 |
| | | | | 2 slices | roasted turkey breast | 46 | 0.25 cup | chopped yellow onion | 65 |
| | | | | | | | 1 tsp | extra virgin olive oil | 5 |
| | | | | | | | 0.25 tsp | toasted sesame oil | 2 |
| Snack #2 | 3 each | Oreo cookies | 34 | 10 each | large baby carrots | 82 | 1 bag | green tea | 0 |
| | | | | 1 cup | 2% cottage cheese | 238 | 1 cup | water | 290 |
| | | | | 2 tbsp | semi-sweet chocolate chips | 17 | 1 each | medium Fuji apple | 153 |
| | | | | | | | 20 each | almonds | 24 |
| Dinner | 2 slices | Stouffer's Lasagna | 515 | 1 tbsp | extra virgin olive oil | 12 | 0.25 each | avocado | 38 |
| | | | | 1 cup | quinoa | 171 | 1 each | portabella mushroom | 159 |
| | | | | 0.33 cup | diced yellow onion | 35 | 1 tbsp | balsamic vinegar | 14 |
| | | | | 4 each | medium asparagus spears | 65 | 1 tbsp | Vegenaise | 14 |
| | | | | 0.5 cup | frozen green peas | 72 | 1 slice | tomato | 57 |
| | | | | 6 oz | boneless skinless chicken poached | 165 | 1 leaf | red lettuce | 14 |
| | | | | 0.5 cup | spinach | 13 | 1 cup | chopped broccoli | 80 |
| | | | | 1 tsp | lemon juice | 0.5 | 1 tsp | lemon juice | 0.5 |
| | | | | 1.5 cup | water | 435 | 1 clove | garlic | 0.5 |
| | | | | | | | 1 tbsp | chopped basil | 0.5 |
| | | | | | | | 1 bun | Oroweat whole wheat burger bun | 75 |
| Snack #3 | | | | | | | 2 cups | pop corn | 12 |
| | | | | | | | 17 each | hazelnuts | 16 |
| | | | | | | | 3 each | Newman's Own fig bars | 56 |

#Table 4: Food preparation descriptions
| | | | |
|-----------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| | Average American | USDA | Vegan |
| Breakfast | Used as purchased from Starbucks. | Cereal poured directly from box into bowl. Milk poured into measuring cup, then into cereal bowl. Raspberries washed first in colander under running water then transferred on top of cereal. | Almond milk brought to a boil, then oats added and cooked for 5 minutes on low heat. Peanut butter and maple syrup measured out then stirred into cooked oats. Banana peeled and sliced into slices on top of cooked oats. |
| Snack #1 | | Apple washed and sliced, core discarded. Yogurt used as purchased. | Soy milk measured into measuring cup, protein powder measured into scoop, banana peeled and cut in half, blueberries rinsed in colander under running water. |
| Lunch | Used as purchased from McDonald's. | Sliced roasted turkey breast deli meat take out of packaging and placed into burger bun. Lettuce rinsed in colander under running water and dried on paper towel then cut into strips and tossed with premade Cesar dressing. Egg boiled in water for 8 minutes then peeled and sliced in half and placed in top of dressed salad. Parmesan cheese shredded and added on top of salad. Grapes rinsed in colander under running water, then sliced in half and placed on top of salad. | Carrot rinsed under running water, peeled, and sliced. Celery washed under running water and sliced. Onion outer layer peeled and diced. Sliced carrot, celery and onion sauteed in olive oil for 5 minutes, then vegetable broth measured out in measuring cup and added to vegetables, brought to a boil. Tofu taken out of packaging, excess water discarded, cut into cubes, added to broth. Spinach taken out of prepackaged, prerinsed bag and added to broth. Noodles and sesame oil added to broth. Soup cooked for 8 minutes on low heat. |
| Snack #2 | Cookies taken out of packaging. | Baby carrots taken out of packaging and used. Cottage cheese measured out in measuring cup. Chocolate chips measured out and used. | Water boiled and poured into cup with tea bag, steeped for 5 minutes. Apple rinsed under running water, sliced, and core discarded. Almonds taken out of packaging. |
| Dinner | Lasagna prepared according to manufacturer instructions (taken out of freezer and baked at 400F for 1 hour and 45 minutes, cooled, then sliced. | Chicken breast taken out of plastic packaging, and placed into pot with boiling water, boiled for 3 minutes, removed from heat, covered, let stand for 18 minutes, then sliced. Quinoa rinsed in colander under running water, added to water in pan and brought to a boil, simmered covered for 20 minutes. Oil heated in large skillet over medium heat, onion peeled and diced, asparagus spears rinsed and sliced, both added to oil and cooked for 5 minutes. Peas added from frozen packaging and cooked for 1 minute. Spinach rinsed in colander under running water and added to skillet, cooked for 3 minutes. Quinoa, vegetables, and chicken combined with lemon juice. | Mushroom destemmed and peeled, soaked in vinegar, then grilled in grill pan for 5 minutes on each side. Garlic peeled and grated into Vegenaise, lemon juice added. Basil rinsed under running water, chopped and added to Vegenaise mixture. Tomato rinsed under running water, then sliced. Lettuce leaf rinsed under running water. Broccoli rinsed under running water, then steamed in colander for 3 minutes, chopped. Burger assembled: Vegenaise mixture spread onto bottom of bun, topped with mushroom, lettuce leaf, tomato slice and top of bun. |
| Snack #3 | | | Popcorn (no salt, no oil) prepared in microwave bag as directed (placed in microwave for 4 minutes). Hazelnuts taken out of packaging. Fig bars taken out of packaging. |

#Table 5: Nutrients by Meal for Each Dietary Pattern
| | | | | | | | | | | | | | | | | |
|------------------------------------------------------|---------------|-------------|---------|------------------|-----------|------------|---------|---------|---------|----------------|----------------|---------|------------------|----------------------|----------|----------|
| Dietary Pattern/Meal | Energy (kcal) | Protein (g) | Fat (g) | Carbohydrate (g) | Fiber (g) | Sugars (g) | Ca (mg) | Fe (mg) | Na (mg) | Vitamin C (mg) | Vitamin A (IU) | SFA (g) | Cholesterol (mg) | Carotene, beta (mcg) | MUFA (g) | PUFA (g) |
| AA Breakfast | 367 | 8.52 | 4.98 | 73.33 | 0 | 60 | 250 | 1.2 | 300 | 0 | 0 | 2.64 | 17 | 0 | | |
| AA Lunch | 1174 | 31.76 | 57.58 | 138.6 | 10 | 43.53 | 280 | 5.7 | 1399 | 12.1 | 412 | 11.52 | 79 | 0 | 19.63 | 7.86 |
| AA Dinner | 568 | 26.62 | 20.81 | 68.2 | 5.6 | 11.4 | 380 | 2.88 | 2102 | 7.2 | 600 | 11.4 | 56 | | | |
| AA Snack | 160 | 1 | 7 | 25 | 1 | 14 | 20 | 1.8 | 160 | 0 | 0 | 2 | 0 | | | |
| USDA Breakfast | 414 | 34.96 | 4.85 | 79.52 | 24 | 27.41 | 466 | 4.09 | 278 | 16.1 | 499 | 1.97 | 12 | 61 | 0.72 | 0.33 |
| USDA Snack #1 | 256 | 8.91 | 3.07 | 52.29 | 3.6 | 20.45 | 268 | 0.3 | 100 | 8 | 836 | 2.06 | 10 | 47 | 1.16 | 0.55 |
| USDA Lunch | 656 | 37.26 | 28.97 | 68.61 | 6 | 19.06 | 298 | 6.06 | 1516 | 110 | 5179 | 7.29 | 223 | 56 | 7.71 | 11.43 |
| USDA Snack #2 | 352 | 29.2 | 11.74 | 34.13 | 5.8 | 27.43 | 254 | 2.22 | 863 | 3.9 | 20852 | 5.99 | 23 | 9600 | 1.01 | 0.27 |
| USDA Dinner | 581 | 48.44 | 19.22 | 56.41 | 5.6 | 6.97 | 58 | 3 | 446 | 25.6 | 3271 | 1.97 | 90 | 1953 | 9.88 | 1.61 |
| VEG Breakfast | 367 | 10 | 12.78 | 58.19 | 6.7 | 28.31 | 311 | 0.9 | 312 | 6 | 544 | 2.25 | 85 | 18 | 3.9 | 2.32 |
| VEG Snack #1 | 373 | 23.76 | 4.95 | 62.8 | 7.7 | 40.46 | 373 | 6.65 | 266 | 80 | 5129 | 0.62 | 0 | 69 | 1.06 | 2.6 |
| VEG Lunch | 468 | 25.49 | 13.89 | 64.31 | 7 | 7.95 | 348 | 6.94 | 1618 | 16.2 | 13184 | 1.79 | 0 | 6850 | 5.9 | 5.82 |
| VEG Snack #2 | 233 | 5.82 | 12.59 | 30.39 | 7.3 | 19.86 | 82 | 1.11 | 9 | 8.4 | 102 | 1.48 | 0 | 50 | 7.99 | 3.48 |
| VEG Dinner | 444 | 15.4 | 20.19 | 55.62 | 16.3 | 13.74 | 176 | 3.5 | 499 | 102.3 | 1779 | 2.64 | 0 | 1018 | 10.74 | 5.56 |
| VEG Snack #3 | 378 | 7.41 | 18.69 | 50.46 | 6.8 | 23.42 | 59 | 3.02 | 169 | 1.6 | 52 | 1.69 | 0 | 17 | 12.48 | 3.51 |

#Table 6: Sequence Summary Statistics By Meal
| | | | |
|---------------------------------------------|----------------------------|-----------------------------|-----------------------------------------|
| Meal | # Sequences Pre-filtration | # Sequences Post-filtration | # OTUs (open reference, 97% similarity) |
| AMERICANbreakfast | 267254 | 226903 | 1838 |
| AMERICANdinner | 298442 | 11666 | 660 |
| AMERICANlunch | 299035 | 96898 | 622 |
| AMERICANsnack | 311311 | 279136 | 969 |
| USDAbreakfast | 318956 | 5002 | 502 |
| USDAdinner | 277213 | 6149 | 476 |
| USDAlunch | 270166 | 16456 | 607 |
| USDAsnack1 | 299998 | 226403 | 334 |
| USDAsnack2 | 238057 | 104114 | 333 |
| VEGANbreakfast | 274360 | 7310 | 399 |
| VEGANdinner | 303246 | 3576 | 417 |
| VEGANlunch | 291459 | 13874 | 480 |
| VEGANsnack1 | 244886 | 62446 | 644 |
| VEGANsnack2 | 288319 | 974 | 1053 |
| VEGANsnack3 | 168669 | 54483 | 229 |

#Table 7: Pairwise Pearson correlations between individual taxonomic groups and meal nutrient contents
| OTU Taxonomy String | Energy | Protein | Total_lipid | Carbohydrate | Fiber | Sugars | Calcium | Iron | Sodium | Vitamin_C | Vitamin_A | Fatty_acids | Cholesterol | Carotene_beta |
|-----------------------------------------------------------------------------------------------------------------|---------|---------|-------------|--------------|---------|----------|----------|---------|---------|-----------|-----------|-------------|-------------|---------------|
| k__Bacteria; p__Proteobacteria | -0.0555 | 0.1469 | -0.1615 | 0.0001 | 0.21 | 0.1548 | 0.306 | 0.5187* | -0.061 | 0.576* | -0.0015 | -0.2234 | 0.026 | -0.2217 |
| k__Bacteria; p__Proteobacteria; c__Betaproteobacteria | 0.1243 | 0.2057 | 0.2017 | -0.0193 | 0.0075 | -0.1199 | -0.0134 | 0.1572 | 0.2141 | 0.4799 | -0.1432 | 0.1704 | 0.731** | -0.2859 |
| k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria | -0.0518 | 0.1531 | -0.1693 | 0.0122 | 0.2062 | 0.1612 | 0.3298 | 0.5275* | -0.0574 | 0.5611* | 0.0223 | -0.2227 | 0 | -0.1961 |
| k__Bacteria; p__Proteobacteria; c__Betaproteobacteria; o__Burkholderiales | 0.1244 | 0.2036 | 0.2057 | -0.0206 | 0.0177 | -0.1158 | -0.0207 | 0.176 | 0.2094 | 0.5029 | -0.1219 | 0.1636 | 0.7205** | -0.2679 |
| k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Aeromonadales | 0.1881 | 0.3457 | 0.2304 | 0.0074 | -0.0007 | -0.246 | 0.0846 | 0.2563 | 0.4289 | 0.5015 | -0.0596 | 0.3302 | 0.7054** | -0.2024 |
| k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Pseudomonadales | 0.2104 | 0.3438 | 0.2531 | 0.0501 | 0.1325 | -0.1955 | 0.1255 | 0.3722 | 0.3475 | 0.6911** | 0.0597 | 0.2116 | 0.7658** | -0.1142 |
| k__Bacteria; p__Bacteroidetes; c__Bacteroidia; o__Bacteroidales; f__[Odoribacteraceae] | 0.0254 | 0.2495 | 0.0801 | -0.1591 | 0.4575 | -0.4456 | -0.5145* | -0.0743 | -0.1711 | 0.2335 | -0.2165 | -0.2105 | -0.0428 | -0.0989 |
| k__Bacteria; p__Firmicutes; c__Bacilli; o__Bacillales; f__Staphylococcaceae | -0.1049 | 0.0316 | -0.0643 | -0.1813 | 0.2204 | -0.5232* | -0.1309 | 0.1455 | 0.1666 | -0.0834 | 0.005 | -0.2458 | -0.1784 | 0.142 |
| k__Bacteria; p__Proteobacteria; c__Betaproteobacteria; o__Burkholderiales; f__Oxalobacteraceae | 0.1536 | 0.2519 | 0.2072 | 0.0201 | -0.023 | -0.0741 | 0.0677 | 0.2238 | 0.2378 | 0.5455* | -0.0402 | 0.1842 | 0.7828** | -0.2121 |
| k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Pseudomonadales; f__Pseudomonadaceae | 0.2383 | 0.3882 | 0.2629 | 0.0766 | 0.1229 | -0.1943 | 0.1457 | 0.4155 | 0.3834 | 0.7115** | 0.0968 | 0.2367 | 0.764** | -0.0822 |
| k__Bacteria; p__Firmicutes; c__Bacilli; o__Bacillales; f__Staphylococcaceae; g__Staphylococcus | -0.1031 | 0.0343 | -0.0633 | -0.1798 | 0.222 | -0.5254* | -0.1292 | 0.1482 | 0.1686 | -0.0819 | 0.0068 | -0.2455 | -0.1771 | 0.1437 |
| k__Bacteria; p__Firmicutes; c__Bacilli; o__Lactobacillales; f__Streptococcaceae; g__Lactococcus | 0.2088 | 0.2213 | 0.2812 | 0.0331 | 0.1385 | 0.142 | -0.0006 | -0.0765 | 0.136 | -0.2465 | 0.603* | 0.3762 | -0.0024 | 0.6236* |
| k__Bacteria; p__Firmicutes; c__Clostridia; o__Clostridiales; f__Lachnospiraceae; g__Blautia | -0.0198 | 0.1057 | 0.028 | -0.1283 | 0.3711 | -0.5586* | -0.2554 | 0.2255 | 0.1582 | 0.0478 | 0.0257 | -0.2386 | -0.2277 | 0.1791 |
| k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Pseudomonadales; f__Pseudomonadaceae; g__Pseudomonas | 0.2379 | 0.3875 | 0.2629 | 0.0764 | 0.1225 | -0.1935 | 0.1452 | 0.4149 | 0.3819 | 0.7113** | 0.0964 | 0.236 | 0.7645** | -0.0828 |

#Table 8: Nutrients and Meal Categories that differed significantly in KEGG (level 3) pathways as predicted by PICRUSt
| Nutrient/Descriptor | KEGG functional category (level 3) | Type of test used | p-value |
|---------------------|------------------------------------------------------------|-------------------|-------------|
| Calcium_bin | Peptidases | ANOVA | 0.039 |
| Carotene_beta_bin | Prenytransferases | Welch's T-test | 0.041 |
| Carotene_beta_bin | Vibrio cholera pathogenic cycle | 0.041 | 0.042 |
| Dairy | Other glycan degradation | Welch's T-test | 0.041 |
| Protein | Phenypropanoid biosynthesis | ANOVA | 0.021 |
| Vitamin_C_bin | Calcium signaling pathway | ANOVA | 0.023 |
| Vitamin_C_bin | Transporters | ANOVA | 0.049 |
| Vitamin_A | Cytoskeleton proteins | ANOVA | 0.017451421 |
| Vitamin_A | Peptidases | ANOVA | 0.019130989 |
| Vitamin_A | Flavonoid biosynthesis | ANOVA | 0.023972731 |
| Vitamin_A | Germination | ANOVA | 0.029007085 |
| Vitamin_A | Chaperones and folding catalysts | ANOVA | 0.031168211 |
| Total/_lipid/_bin | Secondary bile acid biosynthesis | ANOVA | 0.029181829 |
| Total/_lipid/_bin | Biosynthesis of siderophore group nonribosomal peptides | ANOVA | 0.043392615 |
| Sodium | Peptidases | ANOVA | 0.014721635 |
| Sodium | Benzoate degradation | ANOVA | 0.019865137 |
| Sodium | Limonene and pinene degradation | ANOVA | 0.031164989 |
| Sodium | Butanoate metabolism | ANOVA | 0.032537666 |
| Sodium | Nucleotide excision repair | ANOVA | 0.033766923 |
| Sodium | Phenylalanine, tyrosine and tryptophan biosynthesis | ANOVA | 0.036768304 |
| Sodium | Peroxisome | ANOVA | 0.037176174 |
| Sodium | Ethylbenzene degradation | ANOVA | 0.039923012 |
| Sodium | Naphthalene degradation | ANOVA | 0.040749468 |
| Sodium | Restriction enzyme | ANOVA | 0.043486893 |
| Sodium | Tyrosine metabolism | ANOVA | 0.047200716 |
| Iron_bin | Carbon fixation in photosynthetic organisms | ANOVA | 0.021668245 |
| Iron_bin | Protein kinases | ANOVA | 0.022263474 |
| Iron_bin | Translation proteins | ANOVA | 0.024669694 |
| Iron_bin | Pyruvate metabolism | ANOVA | 0.032434551 |
| Iron_bin | Thiamine metabolism | ANOVA | 0.038685974 |
| Iron_bin | D-Glutamine and D-glutamate metabolism | ANOVA | 0.041514757 |
| Iron_bin | One carbon pool by folate | ANOVA | 0.042928705 |
| Fiber_bin | Other glycan degradation | ANOVA | 0.014660951 |
| Fiber_bin | N-Glycan biosynthesis | ANOVA | 0.042139092 |
| Fiber_bin | Proteasome | ANOVA | 0.047590993 |
| Fiber_bin | Prostate cancer | ANOVA | 0.048468144 |
| Fiber_bin | Antigen processing and presentation | ANOVA | 0.048618504 |
| Fiber_bin | Progesterone-mediated oocyte maturation | ANOVA | 0.048618504 |
| Fiber_bin | Other transporters | ANOVA | 0.049130295 |
| Fermented | Transcription factors | Welch's T-test | 0.009468352 |
| Fermented | Phosphonate and phosphinate metabolism | Welch's T-test | 0.011059743 |
| Fermented | Cytoskeleton proteins | Welch's T-test | 0.011669915 |
| Fermented | Amoebiasis | Welch's T-test | 0.012857954 |
| Fermented | Oxidative phosphorylation | Welch's T-test | 0.013841012 |
| Fermented | Transporters | Welch's T-test | 0.015647729 |
| Fermented | Protein processing in endoplasmic reticulum | Welch's T-test | 0.017210951 |
| Fermented | Riboflavin metabolism | Welch's T-test | 0.017938599 |
| Fermented | Steroid hormone biosynthesis | Welch's T-test | 0.020976432 |
| Fermented | PPAR signaling pathway | Welch's T-test | 0.021710104 |
| Fermented | Peroxisome | Welch's T-test | 0.02186366 |
| Fermented | Citrate cycle (TCA cycle) | Welch's T-test | 0.021867149 |
| Fermented | Toluene degradation | Welch's T-test | 0.022788506 |
| Fermented | Carbon fixation pathways in prokaryotes | Welch's T-test | 0.024411154 |
| Fermented | alpha-Linolenic acid metabolism | Welch's T-test | 0.026939926 |
| Fermented | Methane metabolism | Welch's T-test | 0.030552388 |
| Fermented | Synthesis and degradation of ketone bodies | Welch's T-test | 0.031611873 |
| Fermented | RNA degradation | Welch's T-test | 0.032553878 |
| Fermented | Dioxin degradation | Welch's T-test | 0.036756638 |
| Fermented | Adipocytokine signaling pathway | Welch's T-test | 0.04018659 |
| Fermented | Benzoate degradation | Welch's T-test | 0.040925365 |
| Fermented | Chlorocyclohexane and chlorobenzene degradation | Welch's T-test | 0.04187257 |
| Fermented | Nicotinate and nicotinamide metabolism | Welch's T-test | 0.049217654 |
| Fatty_acids_bin | Other glycan degradation | ANOVA | 0.002774493 |
| Fatty_acids_bin | N-Glycan biosynthesis | ANOVA | 0.009783852 |
| Fatty_acids_bin | Proteasome | ANOVA | 0.012010853 |
| Fatty_acids_bin | Prostate cancer | ANOVA | 0.014268966 |
| Fatty_acids_bin | Antigen processing and presentation | ANOVA | 0.014432811 |
| Fatty_acids_bin | Progesterone-mediated oocyte maturation | ANOVA | 0.014432811 |
| Fatty_acids_bin | NOD-like receptor signaling pathway | ANOVA | 0.019075184 |
| Fatty_acids_bin | Chloroalkane and chloroalkene degradation | ANOVA | 0.02167892 |
| Fatty_acids_bin | Taurine and hypotaurine metabolism | ANOVA | 0.022063389 |
| Fatty_acids_bin | Sphingolipid metabolism | ANOVA | 0.028794158 |
| Fatty_acids_bin | Other transporters | ANOVA | 0.033398699 |
| Fatty_acids_bin | Primary bile acid biosynthesis | ANOVA | 0.033705405 |
| Fatty_acids_bin | Stilbenoid, diarylheptanoid and gingerol biosynthesis | ANOVA | 0.034295974 |
| Fatty_acids_bin | Glycolysis / Gluconeogenesis | ANOVA | 0.041601128 |
| Fatty_acids_bin | Amyotrophic lateral sclerosis (ALS) | ANOVA | 0.045157824 |
| Fatty_acids_bin | Purine metabolism | ANOVA | 0.048100832 |
| Energy_bin | Nucleotide excision repair | ANOVA | 0.002532996 |
| Energy_bin | Chromosome | ANOVA | 0.003249937 |
| Energy_bin | Ubiquinone and other terpenoid-quinone biosynthesis | ANOVA | 0.00500867 |
| Energy_bin | Mismatch repair | ANOVA | 0.006531889 |
| Energy_bin | Photosynthesis proteins | ANOVA | 0.00816459 |
| Energy_bin | Photosynthesis | ANOVA | 0.009638124 |
| Energy_bin | Restriction enzyme | ANOVA | 0.009805668 |
| Energy_bin | Carbohydrate metabolism | ANOVA | 0.009966615 |
| Energy_bin | Limonene and pinene degradation | ANOVA | 0.010243464 |
| Energy_bin | DNA replication proteins | ANOVA | 0.010451223 |
| Energy_bin | Lipoic acid metabolism | ANOVA | 0.010843214 |
| Energy_bin | Phenylalanine, tyrosine and tryptophan biosynthesis | ANOVA | 0.011059994 |
| Energy_bin | Peptidases | ANOVA | 0.01407109 |
| Energy_bin | DNA repair and recombination proteins | ANOVA | 0.015827001 |
| Energy_bin | Type II diabetes mellitus | ANOVA | 0.016276453 |
| Energy_bin | alpha-Linolenic acid metabolism | ANOVA | 0.018912274 |
| Energy_bin | Butanoate metabolism | ANOVA | 0.020566276 |
| Energy_bin | Homologous recombination | ANOVA | 0.025440256 |
| Energy_bin | Flavone and flavonol biosynthesis | ANOVA | 0.025662772 |
| Energy_bin | Protein export | ANOVA | 0.025889019 |
| Energy_bin | DNA replication | ANOVA | 0.026934029 |
| Energy_bin | Primary immunodeficiency | ANOVA | 0.027524778 |
| Energy_bin | Glycosphingolipid biosynthesis - lacto and neolacto series | ANOVA | 0.027782212 |
| Energy_bin | Indole alkaloid biosynthesis | ANOVA | 0.028725346 |
| Energy_bin | Amoebiasis | ANOVA | 0.028866801 |
| Energy_bin | Benzoate degradation | ANOVA | 0.029913338 |
| Energy_bin | D-Alanine metabolism | ANOVA | 0.029961727 |
| Energy_bin | C5-Branched dibasic acid metabolism | ANOVA | 0.030293739 |
| Energy_bin | Peptidoglycan biosynthesis | ANOVA | 0.031023499 |
| Energy_bin | Glycerolipid metabolism | ANOVA | 0.032020201 |
| Energy_bin | Bisphenol degradation | ANOVA | 0.032129248 |
| Energy_bin | Betalain biosynthesis | ANOVA | 0.032513208 |
| Energy_bin | Biosynthesis of siderophore group nonribosomal peptides | ANOVA | 0.033939326 |
| Energy_bin | Melanogenesis | ANOVA | 0.035235651 |
| Energy_bin | Amino acid metabolism | ANOVA | 0.037736626 |
| Energy_bin | Ribosome Biogenesis | ANOVA | 0.039545196 |
| Energy_bin | Peroxisome | ANOVA | 0.040442991 |
| Energy_bin | Steroid hormone biosynthesis | ANOVA | 0.041078867 |
| Energy_bin | Amino sugar and nucleotide sugar metabolism | ANOVA | 0.042902641 |
| Energy_bin | Phosphotransferase system (PTS) | ANOVA | 0.043431781 |
| Energy_bin | Arginine and proline metabolism | ANOVA | 0.043807314 |
| Energy_bin | Glycine, serine and threonine metabolism | ANOVA | 0.044386803 |
| Energy_bin | Riboflavin metabolism | ANOVA | 0.044908782 |
| Energy_bin | Metabolism of cofactors and vitamins | ANOVA | 0.04499208 |
| Energy_bin | Systemic lupus erythematosus | ANOVA | 0.045791434 |
| Energy_bin | Biosynthesis of type II polyketide products | ANOVA | 0.047100884 |
| DietType | Other glycan degradation | ANOVA | 0.007486572 |
| DietType | N-Glycan biosynthesis | ANOVA | 0.03160732 |
| DietType | Proteasome | ANOVA | 0.034750295 |
| DietType | Prostate cancer | ANOVA | 0.03926696 |
| DietType | Antigen processing and presentation | ANOVA | 0.039372465 |
| DietType | Progesterone-mediated oocyte maturation | ANOVA | 0.039372465 |
| Cooked | mRNA surveillance pathway | Welch's T-test | 0.006445676 |
| Cooked | Cell cycle | Welch's T-test | 0.009317081 |
| Cooked | Hepatitis C | Welch's T-test | 0.009317081 |
| Cooked | Measles | Welch's T-test | 0.009317081 |
| Cooked | mTOR signaling pathway | Welch's T-test | 0.009317081 |
| Cooked | Phagosome | Welch's T-test | 0.009317081 |
| Cooked | Transcription machinery | Welch's T-test | 0.021347218 |
| Cooked | Various types of N-glycan biosynthesis | Welch's T-test | 0.023199289 |
| Cooked | Sporulation | Welch's T-test | 0.024693602 |
| Cooked | Vibrio cholerae infection | Welch's T-test | 0.025037566 |
| Cooked | Cytoskeleton proteins | Welch's T-test | 0.04573813 |
| Cooked | Cytochrome P450 | Welch's T-test | 0.049095188 |
| Cholesterol_bin | Transcription machinery | ANOVA | 0.007389033 |
| Cholesterol_bin | Plant-pathogen interaction | ANOVA | 0.010583548 |
| Cholesterol_bin | Folate biosynthesis | ANOVA | 0.012498951 |
| Cholesterol_bin | Tetracycline biosynthesis | ANOVA | 0.034405715 |
| Cholesterol_bin | Other ion-coupled transporters | ANOVA | 0.036332221 |
| Cholesterol_bin | Proteasome | ANOVA | 0.039954688 |
| Cholesterol_bin | Valine, leucine and isoleucine biosynthesis | ANOVA | 0.040216227 |
| Cholesterol_bin | General function prediction only | ANOVA | 0.042978367 |
| Cholesterol_bin | Prostate cancer | ANOVA | 0.044525651 |
| Cholesterol_bin | Antigen processing and presentation | ANOVA | 0.044753096 |
| Cholesterol_bin | Progesterone-mediated oocyte maturation | ANOVA | 0.044753096 |
| Cholesterol_bin | NOD-like receptor signaling pathway | ANOVA | 0.047390633 |

### References

1. Nejma Aboutaleb, Edward J. Kuijper, Jaap T. van Dissel. Emerging infectious colitis. Current Opinion in Gastroenterology 30, 106-115 Ovid Technologies (Wolters Kluwer Health), 2014. Link

2. Manimozhiyan Arumugam, Jeroen Raes, Eric Pelletier, Denis Le Paslier, Takuji Yamada, Daniel R. Mende, Gabriel R. Fernandes, Julien Tap, Thomas Bruls, Jean-Michel Batto, Marcelo Bertalan, Natalia Borruel, Francesc Casellas, Leyden Fernandez, Laurent Gautier, Torben Hansen, Masahira Hattori, Tetsuya Hayashi, Michiel Kleerebezem, Ken Kurokawa, Marion Leclerc, Florence Levenez, Chaysavanh Manichanh, H. BjNielsen, Trine Nielsen, Nicolas Pons, Julie Poulain, Junjie Qin, Thomas Sicheritz-Ponten, Sebastian Tims, David Torrents, Edgardo Ugarte, Erwin G. Zoetendal, Jun Wang, Francisco Guarner, Oluf Pedersen, Willem M. de Vos, SBrunak, Joel Doré, María Antolín, François Artiguenave, Hervé M. Blottiere, Mathieu Almeida, Christian Brechot, Carlos Cara, Christian Chervaux, Antonella Cultrone, Christine Delorme, Gérard Denariaz, Rozenn Dervyn, Konrad U. Foerstner, Carsten Friss, Maarten van de Guchte, Eric Guedon, Florence Haimet, Wolfgang Huber, Johan van Hylckama-Vlieg, Alexandre Jamet, Catherine Juste, Ghalia Kaci, Jan Knol, Omar Lakhdari, Severine Layec, Karine Le Roux, Emmanuelle Maguin, Alexandre Mérieux, Raquel Melo Minardi, Christine M’rini, Jean Muller, Raish Oozeer, Julian Parkhill, Pierre Renault, Maria Rescigno, Nicolas Sanchez, Shinichi Sunagawa, Antonio Torrejon, Keith Turner, Gaetana Vandemeulebrouck, Encarna Varela, Yohanan Winogradsky, Georg Zeller, Jean Weissenbach, S. Dusko Ehrlich, Peer Bork. Enterotypes of the human gut microbiome. Nature 473, 174-180 Nature Publishing Group, 2011. Link

3. J Bien, V Palagani, P Bozko. The intestinal microbiota dysbiosis and Clostridium difficile infection: is there a relationship with inflammatory bowel disease?. 6, 53-68 (2013).

4. J. Roger Bray, J. T. Curtis. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecological Monographs 27, 325 JSTOR, 1957. Link

5. AM Brownawell, W Caers, GR Gibson, CW Kendall, KD Lewis, Y Ringel, JL Slavin. Prebiotics and the health benefits of fiber: current regulatory status, future research, and goals.. 142, 962-74 (2012).

6. Matteo Busconi, Carla Zacconi, Gianluigi Scolari. Bacterial ecology of PDO Coppa and Pancetta Piacentina at the end of ripening and after MAP storage of sliced product. International Journal of Food Microbiology 172, 13-20 Elsevier BV, 2014. Link

7. Julie E. Button, Rachel J. Dutton. Cheese microbes. Current Biology 22, R587–R589 Elsevier BV, 2012. Link

8. JG Caporaso, J Kuczynski, J Stombaugh, K Bittinger, FD Bushman, EK Costello, N Fierer, AG Peña, JK Goodrich, JI Gordon, GA Huttley, ST Kelley, D Knights, JE Koenig, RE Ley, CA Lozupone, D McDonald, BD Muegge, M Pirrung, J Reeder, JR Sevinsky, PJ Turnbaugh, WA Walters, J Widmann, T Yatsunenko, J Zaneveld, R Knight. QIIME allows analysis of high-throughput community sequencing data.. 7, 335-6 (2010).

9. EK Costello, K Stagaman, L Dethlefsen, BJ Bohannan, DA Relman. The application of ecological theory toward an understanding of the human microbiome.. 336, 1255-62 (2012).

10. E. K. Costello, C. L. Lauber, M. Hamady, N. Fierer, J. I. Gordon, R. Knight. Bacterial Community Variation in Human Body Habitats Across Space and Time. Science 326, 1694-1697 American Association for the Advancement of Science (AAAS), 2009. Link

11. Faith J Critzer, Michael P Doyle. Microbial ecology of foodborne pathogens associated with produce. Current Opinion in Biotechnology 21, 125-130 Elsevier BV, 2010. Link

12. Lawrence A. David, Corinne F. Maurice, Rachel N. Carmody, David B. Gootenberg, Julie E. Button, Benjamin E. Wolfe, Alisha V. Ling, A. Sloan Devlin, Yug Varma, Michael A. Fischbach, Sudha B. Biddinger, Rachel J. Dutton, Peter J. Turnbaugh. Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559-563 Nature Publishing Group, 2013. Link

13. L. Dethlefsen, D. A. Relman. Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proceedings of the National Academy of Sciences 108, 4554–4561 Proceedings of the National Academy of Sciences, 2010. Link

14. J. J. Faith, J. L. Guruge, M. Charbonneau, S. Subramanian, H. Seedorf, A. L. Goodman, J. C. Clemente, R. Knight, A. C. Heath, R. L. Leibel, M. Rosenbaum, J. I. Gordon. The Long-Term Stability of the Human Gut Microbiota. Science 341, 1237439-1237439 American Association for the Advancement of Science (AAAS), 2013. Link

15. C. De Filippo, D. Cavalieri, M. Di Paola, M. Ramazzotti, J. B. Poullet, S. Massart, S. Collini, G. Pieraccini, P. Lionetti. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proceedings of the National Academy of Sciences 107, 14691-14696 Proceedings of the National Academy of Sciences, 2010. Link

16. M. Gatti, J. De Dea Lindner, A. De Lorentiis, B. Bottari, M. Santarelli, V. Bernini, E. Neviani. Dynamics of Whole and Lysed Bacterial Cells during Parmigiano-Reggiano Cheese Production and Ripening. Applied and Environmental Microbiology 74, 6161-6167 American Society for Microbiology, 2008. Link

17. JB German, SL Freeman, CB Lebrilla, DA Mills. Human milk oligosaccharides: evolution, structures and bioselectivity as substrates for intestinal bacteria.. 62, 205-18; discussion 218-22 (2008).

18. G Henderson, F Cox, S Kittelmann, VH Miri, M Zethof, SJ Noel, GC Waghorn, PH Janssen. Effect of DNA extraction methods and sampling techniques on the apparent structure of cow and sheep rumen microbial communities.. 8, e74787 (2013).

19. Marie A. Hildebrandt, Christian Hoffmann, Scott A. Sherrill$$\textendash$$Mix, Sue A. Keilbaugh, Micah Hamady, Ying$$\textendash$$Yu Chen, Rob Knight, Rexford S. Ahima, Frederic Bushman, Gary D. Wu. High-Fat Diet Determines the Composition of the Murine Gut Microbiome Independently of Obesity. Gastroenterology 137, 1716-1724.e2 Elsevier BV, 2009. Link

20. William H. Kruskal, W. Allen Wallis. Use of Ranks in One-Criterion Variance Analysis. Journal of the American Statistical Association 47, 583-621 Informa UK Limited, 1952. Link

21. MG Langille, J Zaneveld, JG Caporaso, D McDonald, D Knights, JA Reyes, JC Clemente, DE Burkepile, Thurber RL Vega, R Knight, RG Beiko, C Huttenhower. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences.. 31, 814-21 (2013).

22. JW Leff, N Fierer. Bacterial communities associated with the surfaces of fresh fruits and vegetables.. 8, e59310 (2013).

23. C Lozupone, ME Lladser, D Knights, J Stombaugh, R Knight. UniFrac: an effective distance metric for microbial community comparison.. 5, 169-72 (2011).

24. TE Machonkin, MD Boshart, JA Schofield, MM Rodriguez, K Grubel, D Rokhsana, WW Brennessel, PL Holland. Structural and Spectroscopic Characterization of Iron(II), Cobalt(II), and Nickel(II) ortho-Dihalophenolate Complexes: Insights into Metal-Halogen Secondary Bonding.. (2014).

25. T Magoč, SL Salzberg. FLASH: fast length adjustment of short reads to improve genome assemblies.. 27, 2957-63 (2011).

26. C Manichanh, J Reeder, P Gibert, E Varela, M Llopis, M Antolin, R Guigo, R Knight, F Guarner. Reshaping the gut microbiome with bacterial transplantation and antibiotic intake.. 20, 1411-9 (2010).

27. Adam C Martiny, Kathleen Treseder, Gordon Pusch. Phylogenetic conservatism of functional traits in microorganisms. The ISME Journal 7, 830–838 Nature Publishing Group, 2012. Link

28. N. P. McNulty, T. Yatsunenko, A. Hsiao, J. J. Faith, B. D. Muegge, A. L. Goodman, B. Henrissat, R. Oozeer, S. Cools-Portier, G. Gobert, C. Chervaux, D. Knights, C. A. Lozupone, R. Knight, A. E. Duncan, J. R. Bain, M. J. Muehlbauer, C. B. Newgard, A. C. Heath, J. I. Gordon. The Impact of a Consortium of Fermented Milk Strains on the Gut Microbiome of Gnotobiotic Mice and Monozygotic Twins. Science Translational Medicine 3, 106ra106-106ra106 American Association for the Advancement of Science (AAAS), 2011. Link

29. B. D. Muegge, J. Kuczynski, D. Knights, J. C. Clemente, A. Gonzalez, L. Fontana, B. Henrissat, R. Knight, J. I. Gordon. Diet Drives Convergence in Gut Microbiome Functions Across Mammalian Phylogeny and Within Humans. Science 332, 970-974 American Association for the Advancement of Science (AAAS), 2011. Link

30. DH Parks, GW Tyson, P Hugenholtz, RG Beiko. STAMP: Statistical analysis of taxonomic and functional profiles.. (2014).

31. Karl Pearson. Note on Regression and Inheritance in the Case of Two Parents. Proceedings of the Royal Society of London (1854-1905) 58, 240-242 The Royal Society, 1895. Link

32. LM Proctor. The Human Microbiome Project in 2011 and beyond.. 10, 287-91 (2011).

33. Valeria Rizzello, Irene Bonaccorsi, Maria Luisa Dongarrà, Lisbeth Nielsen Fink, Guido Ferlazzo. Role of Natural Killer and Dendritic Cell Crosstalk in Immunomodulation by Commensal Bacteria Probiotics. Journal of Biomedicine and Biotechnology 2011, 1-10 Hindawi Publishing Corporation, 2011. Link

34. M Roberfroid. Prebiotics: the concept revisited.. 137, 830S-7S (2007).

35. S.C. Stringer, M.D. Webb, M.W. Peck. Lag time variability in individual spores of Clostridium botulinum. Food Microbiology 28, 228-235 Elsevier BV, 2011. Link

36. Amarela Terzic-Vidojevic, Sanja Mihajlovic, Gordana Uzelac, Katarina Veljovic, Maja Tolinacki, Milica Nikolic, Ljubisa Topisirovic, Milan Kojic. Characterization of lactic acid bacteria isolated from artisanal Travnik young cheeses sweet creams and sweet kajmaks over four seasons. Food Microbiology 39, 27-38 Elsevier BV, 2014. Link

37. Kieran M. Tuohy, Lorenza Conterno, Mattia Gasperotti, Roberto Viola. Up-regulating the Human Intestinal Microbiome Using Whole Plant Foods Polyphenols, and/or Fiber. Journal of Agricultural and Food Chemistry 60, 8776–8782 American Chemical Society (ACS), 2012. Link

38. P. J. Turnbaugh, V. K. Ridaura, J. J. Faith, F. E. Rey, R. Knight, J. I. Gordon. The Effect of Diet on the Human Gut Microbiome: A Metagenomic Analysis in Humanized Gnotobiotic Mice. Science Translational Medicine 1, 6ra14-6ra14 American Association for the Advancement of Science (AAAS), 2009. Link

39. C Ubeda, EG Pamer. Antibiotics, microbiota, and immune defense.. 33, 459-66 (2012).

40. Alan W Walker, Jennifer Ince, Sylvia H Duncan, Lucy M Webster, Grietje Holtrop, Xiaolei Ze, David Brown, Mark D Stares, Paul Scott, Aurore Bergerat, Petra Louis, Freda McIntosh, Alexandra M Johnstone, Gerald E Lobley, Julian Parkhill, Harry J Flint. Dominant and diet-responsive groups of bacteria within the human colonic microbiota. The ISME Journal 5, 220-230 Nature Publishing Group, 2011. Link

41. Alan W Walker, Jennifer Ince, Sylvia H Duncan, Lucy M Webster, Grietje Holtrop, Xiaolei Ze, David Brown, Mark D Stares, Paul Scott, Aurore Bergerat, Petra Louis, Freda McIntosh, Alexandra M Johnstone, Gerald E Lobley, Julian Parkhill, Harry J Flint. Dominant and diet-responsive groups of bacteria within the human colonic microbiota. The ISME Journal 5, 220–230 Nature Publishing Group, 2010. Link

42. G. D. Wu, J. Chen, C. Hoffmann, K. Bittinger, Y.-Y. Chen, S. A. Keilbaugh, M. Bewtra, D. Knights, W. A. Walters, R. Knight, R. Sinha, E. Gilroy, K. Gupta, R. Baldassano, L. Nessel, H. Li, F. D. Bushman, J. D. Lewis. Linking Long-Term Dietary Patterns with Gut Microbial Enterotypes. Science 334, 105–108 American Association for the Advancement of Science (AAAS), 2011. Link

43. GD Wu, FD Bushmanc, JD Lewis. Diet, the human gut microbiota, and IBD.. 24, 117-20 (2013).

44. AM Zivkovic, D Barile. Bovine milk as a source of functional oligosaccharides for improving human health.. 2, 284-9 (2011).

45. Vrese M de, J Schrezenmeir. Probiotics, prebiotics, and synbiotics.. 111, 1-66 (2008).