Linking long-term dietary patterns with gut microbial enterotypes

Gary D Wu, Jun Chen, Christian Hoffmann, Kyle Bittinger, Ying-Yu Chen, Sue A Keilbaugh, Meenakshi Bewtra, Dan Knights, William A Walters, Rob Knight, Rohini Sinha, Erin Gilroy, Kernika Gupta, Robert Baldassano, Lisa Nessel, Hongzhe Li, Frederic D Bushman, James D Lewis, Gary D Wu, Jun Chen, Christian Hoffmann, Kyle Bittinger, Ying-Yu Chen, Sue A Keilbaugh, Meenakshi Bewtra, Dan Knights, William A Walters, Rob Knight, Rohini Sinha, Erin Gilroy, Kernika Gupta, Robert Baldassano, Lisa Nessel, Hongzhe Li, Frederic D Bushman, James D Lewis

Abstract

Diet strongly affects human health, partly by modulating gut microbiome composition. We used diet inventories and 16S rDNA sequencing to characterize fecal samples from 98 individuals. Fecal communities clustered into enterotypes distinguished primarily by levels of Bacteroides and Prevotella. Enterotypes were strongly associated with long-term diets, particularly protein and animal fat (Bacteroides) versus carbohydrates (Prevotella). A controlled-feeding study of 10 subjects showed that microbiome composition changed detectably within 24 hours of initiating a high-fat/low-fiber or low-fat/high-fiber diet, but that enterotype identity remained stable during the 10-day study. Thus, alternative enterotype states are associated with long-term diet.

Figures

Fig. 1
Fig. 1
Correlation of diet and gut microbial taxa identified in the cross-sectional COMBO analysis. Columns correspond to bacterial taxa quantified using 16S rDNA tags; rows correspond to nutrients measured by dietary questionnaire. Red and blue denote positive and negative association, respectively. The intensity of the colors represents the degree of association between the taxa abundances and nutrients as measured by the Spearman's correlations. Bacterial phyla are summarized by the color code on the bottom; lower-level taxonomic assignments specified are in fig. S1. The dots indicate the associations that are significant at an FDR of 25%. The FFQ data were used for this comparison (both FFQ and Recall dietary data are shown together in fig. S1). Columns and rows are clustered by Euclidean distance, with rows separated by the predominant nutrient category.
Fig. 2
Fig. 2
Clustering of gut microbial taxa into entero-types is associated with long-term diet. (A) Clustering in the COMBO cross-sectional study using Jensen-Shannon distance. The left panel shows that the data are most naturally separated into two clusters by the PAM method. The x axis shows cluster number; the y axis shows silhouette width, a measure of cluster separation (12). The right panel shows the clustering on the first two principal components. (B) Proportions of bacterial taxa characteristic of each enterotype. Boxes represent the interquartile range (IQR) and the line inside represents the median. Whiskers denote the lowest and highest values within 1.5 × IQR. (C) The association of dietary components with each enterotype. The strength and direction of each association, as measured by the means of the standardized nutrient measurements, is shown by the color key at the lower right. Enterotype is shown at the right. Red indicates greater amounts, blue lesser amounts of each nutrient in each enterotype (complete lists of nutrients are in table S2). Columns were clustered by Euclidean distance.
Fig. 3
Fig. 3
Changes in bacterial communities during controlled feeding. Ten subjects were randomized to high-fat/low-fiber or low-fat/high-fiber diets, and microbiome composition was monitored longitudinally for 10 days by sequencing 16S rDNA gene tags (CAFE study). (A) Cluster diagram–based principal coordinates analysis using unweighted UniFrac. Colors indicate samples from each individual. (B) Day 1 samples are outliers compared to all other days, indicating change in the gut microbiome within 24 hours of initiating controlled feeding. In this analysis, weighted UniFrac distances between samples are compared within subjects in two groups. The first collection of distances compares the day 1 samples to days 2 to 10; the second group compares samples from all days to all others excluding day 1, indicating rapid change (P = 0.0003, 10,000 permutations). Error bars indicate 1 SD of the distances.

References

    1. Hooper LV, Gordon JI. Science. 2001;292:1115.
    1. Emminger A, Kahmann E, Savage DS. Cancer Lett. 1977;2:273.
    1. Gill SR, et al. Science. 2006;312:1355.
    1. Arumugam M, et al. Nature. 2011;473:174.
    1. Wu GD, et al. BMC Microbiol. 2010;10:206.
    1. Margulies M, et al. Nature. 2005;437:376.
    1. See supporting material on Science Online.

    1. Lozupone C, Hamady M, Knight R. BMC Bioinformatics. 2006;7:371.
    1. McArdle BH, Anderson MJ. Ecology. 2001;82:290.
    1. Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Nature. 2006;444:1022.
    1. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. Nucleic Acids Res. 2004;32(database issue):D277.
    1. Rousseeuw PJ. J Comput Appl Math. 1987;20:53.
    1. De Filippo C, et al. Proc Natl Acad Sci USA. 2010;107:14691.

Source: PubMed

3
Abonneren