Hostname: page-component-cc8bf7c57-n7pht Total loading time: 0 Render date: 2024-12-10T17:35:07.629Z Has data issue: false hasContentIssue false

Association of dietary fibre intake and gut microbiota in adults

Published online by Cambridge University Press:  24 October 2018

Daniel Lin
Affiliation:
Perlmutter Cancer Center, New York University Langone Medical Center, New York, NY 10016, USA
Brandilyn A. Peters
Affiliation:
Department of Population Health, New York University School of Medicine, New York, NY 10016, USA
Charles Friedlander
Affiliation:
Concorde Medical Group, New York, NY 10016, USA Department of Medicine, New York University School of Medicine, New York, NY 10016, USA
Hal J. Freiman
Affiliation:
Department of Medicine, New York University School of Medicine, New York, NY 10016, USA
James J. Goedert
Affiliation:
Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, MD 20850, USA
Rashmi Sinha
Affiliation:
Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, MD 20850, USA
George Miller
Affiliation:
Perlmutter Cancer Center, New York University Langone Medical Center, New York, NY 10016, USA Department of Surgery, New York University School of Medicine, New York, NY 10016, USA
Mitchell A. Bernstein
Affiliation:
Department of Surgery, New York University School of Medicine, New York, NY 10016, USA
Richard B. Hayes
Affiliation:
Perlmutter Cancer Center, New York University Langone Medical Center, New York, NY 10016, USA Department of Population Health, New York University School of Medicine, New York, NY 10016, USA
Jiyoung Ahn*
Affiliation:
Perlmutter Cancer Center, New York University Langone Medical Center, New York, NY 10016, USA Department of Population Health, New York University School of Medicine, New York, NY 10016, USA
*
*Corresponding author: J. Ahn, email Jiyoung.Ahn@nyumc.org
Rights & Permissions [Opens in a new window]

Abstract

Increasing evidence indicates that gut microbiota may influence colorectal cancer risk. Diet, particularly fibre intake, may modify gut microbiota composition, which may affect cancer risk. We investigated the relationship between dietary fibre intake and gut microbiota in adults. Using 16S rRNA gene sequencing, we assessed gut microbiota in faecal samples from 151 adults in two independent study populations: National Cancer Institute (NCI), n 75, and New York University (NYU), n 76. We calculated energy-adjusted fibre intake based on FFQ. For each study population with adjustment for age, sex, race, BMI and smoking, we evaluated the relationship between fibre intake and gut microbiota community composition and taxon abundance. Total fibre intake was significantly associated with overall microbial community composition in NYU (P=0·008) but not in NCI (P=0·81). In a meta-analysis of both study populations, higher fibre intake tended to be associated with genera of class Clostridia, including higher abundance of SMB53 (fold change (FC)=1·04, P=0·04), Lachnospira (FC=1·03, P=0·05) and Faecalibacterium (FC=1·03, P=0·06), and lower abundance of Actinomyces (FC=0·95, P=0·002), Odoribacter (FC=0·95, P=0·03) and Oscillospira (FC=0·96, P=0·06). A species-level meta-analysis showed that higher fibre intake was marginally associated with greater abundance of Faecalibacterium prausnitzii (FC=1·03, P=0·07) and lower abundance of Eubacterium dolichum (FC=0·96, P=0·04) and Bacteroides uniformis (FC=0·97, P=0·05). Thus, dietary fibre intake may impact gut microbiota composition, particularly class Clostridia, and may favour putatively beneficial bacteria such as F. prausnitzii. These findings warrant further understanding of diet–microbiota relationships for future development of colorectal cancer prevention strategies.

Type
Full Papers
Copyright
© The Authors 2018 

The human gastrointestinal tract hosts an estimated 100 trillion bacteria, which play a role in key physiological activities, including gastrointestinal immune system maturation and fermentation of nutrients into beneficial metabolites( 1 ). Disruption of this symbiotic relationship between human host and gut microbiota has been implicated in the development of the metabolic syndrome, as well as intestinal pathology, including inflammatory bowel disease and colorectal cancer( Reference Hold, Smith and Grange 2 Reference Han, Wang and Liu 4 ). A growing number of epidemiological studies have provided increasing evidence that perturbation of microbial community composition in the gut exists in colorectal cancer, with alterations of microbial taxa abundance in cancer cases compared with healthy controls( Reference Sobhani, Tap and Roudot-Thoraval 5 Reference Weir, Manter and Sheflin 8 ).

Colorectal cancer risk has been attributed to dietary habits, with Western-style diets – low in fibre and high in red meat and fat – associated with higher risk for colorectal cancer( Reference Bouvard, Loomis and Guyton 9 , Reference Aune, Chan and Lau 10 ). Fibre intake, in particular, has remained an appealing modifiable dietary factor, given its postulated beneficial biological effects. Fibre speeds colonic transit and may decrease exposure of colonic epithelial cells to ingested carcinogens. In addition, fibre undergoes fermentation by the microbiota to yield SCFA end products, such as butyrate, which is not only essential for colon energy metabolism and epithelial proliferation, but in mouse models also exhibits tumour-suppressive activity through histone deacetylase inhibition( Reference Donohoe, Holley and Collins 11 ).

Consequently, there has been growing interest to understand the impact of dietary fibre on gut microbiota composition, which may ultimately reduce cancer risk. A diet low in fibre – the ‘fibre gap’ – has been associated with a depletion in human gut microbiota diversity, and potentially beneficial metabolites( Reference Han, Wang and Liu 4 ). Furthermore, although short-term dietary intervention trials have demonstrated that different amounts of fibre intake can significantly alter microbiota composition in a span of a few weeks( Reference O’Keefe, Li and Lahti 12 Reference Liu, Zhao and Guo 14 ), there remain fewer studies evaluating the effect of long-term dietary habits of fibre intake on the gut microbiota in humans( Reference Wu, Chen and Hoffmann 15 ).

We investigated the association between long-term dietary fibre intake habits and gut microbiota composition in faecal samples of generally healthy adults from two independent study populations: controls from a hospital-based case–control study of colorectal cancer and gut microbiome( Reference Ahn, Sinha and Pei 6 ) and polyp-free adults from a cross-sectional colonoscopy study( Reference Peters, Dominianni and Shapiro 16 ). We sought to examine the relationship between higher dietary fibre intake and overall gut microbiota composition, as well as specific taxa abundance.

Methods

Study population

We assessed faecal samples of adults from two independent study populations: control subjects from a National Cancer Institute (NCI) hospital-based case–control study, hereafter referred to as NCI( Reference Ahn, Sinha and Pei 6 ), and polyp-free adults from a cross-sectional colonoscopy study at New York University (NYU), called the NYU Human Microbiome and Colorectal Tumor study, hereafter referred to as NYU( Reference Peters, Dominianni and Shapiro 16 ).

NCI enrolled participants from three Washington, DC area hospitals from 1985 to 1989( Reference Schiffman, Van Tassell and Robinson 17 , Reference Schiffman, Andrews and Van Tassell 18 ). We included seventy-five control subjects who were awaiting elective surgery for non-oncological, non-gastrointestinal conditions, and reported no antibiotic intake during the year before recruitment. Participants provided 2-d faecal samples that were freeze-dried, and samples with at least 100 mg of lyophilised faecal material available were included for analyses.

NYU enrolled participants from Kips Bay Endoscopy Center in New York City from 2012 to 2014. We included seventy-six polyp-free participants from this study. We excluded subjects with missing colonoscopy reports, history of inflammatory bowel disease, prior surgical anastomosis, prior history of colorectal cancer, history of familial adenomatous polyposis, and with most recent colonoscopy report >3 years before stool sample collection.

In both studies, participants provided written informed consent, reported no long-term antibiotic treatment and completed diet and demographic questionnaires. We excluded subjects with <1000 microbial sequence reads, missing or extreme energetic intake (≤2092 or >16 736 kJ/d) and with a history of other cancers, for a final sample size of 151 (n 75 in NCI, n 76 in NYU). The NYU study was approved by the institutional review board (IRB) of NYU School of Medicine, and the NCI study was approved by the IRB of NYU School of Medicine and the NCI.

Demographic information and dietary fibre assessment

Information on age, sex, height, weight, race and cigarette smoking status was collected by means of a questionnaire at stool collection. BMI was calculated by dividing weight in kg by squared height in m, and was then categorised as underweight or normal weight (<25 kg/m2), overweight (25≤BMI<30 kg/m2) or obese (BMI≥30 kg/m2), based on WHO definition( 19 ). Cigarette smoking status was defined as never, current or former smoker.

Usual dietary intake was calculated from self-administered FFQ, which queried intake frequency and portion size of food types. Nutrient values per portion were multiplied by daily frequency of intake and summed across all relevant food items. Nutrient data were standardised by total energy intake( Reference Willett, Howe and Kushi 20 ). Study-specific quartiles of total fibre intake were used (NCI: <11·3, 11·3–14·3, 14·4–16·6, ≥16·7 g/d; NYU: <19·8, 19·8–23·7, 23·8–31·6, ≥31·7 g/d). In addition, study-specific quartiles of fibre from different dietary sources were assessed: fibre from fruits and vegetables (NCI: <6·5, 6·5–8·7, 8·8–11·0, ≥11·1 g/d; NYU: <11·1, 11·1–14·9, 15·0–21·2, ≥21·3 g/d), fibre from grains (NCI: <2·7, 2·7–3·2, 3·3–4·8, ≥4·9 g/d; NYU: <9·1, 9·1–12·6, 12·7–16·6, ≥16·7 g/d) and fibre from beans (NCI: <0·9, 0·9–1·3, 1·4–2·1, ≥2·2 g/d; NYU: <1·8, 1·8–2·7, 2·8–4·4, ≥4·5 g/d).

Faecal samples

In NCI, faecal samples were collected by participants at home over a 2-d period, before hospitalisation and treatment, and stored in a plastic container in a styrofoam chest containing dry ice. Faecal samples were shipped to a USDA laboratory, lyophilised and stored at a minimum of −40°C in sealed, air-tight containers. Sample aliquots were shipped to NYU for microbiome assay. In NYU, faecal samples were collected by participants onto two sections of Beckman Coulter Hemoccult II SENSA® cards (Beckman Coulter) at home. Samples were shipped to NYU and stored immediately at −80°C.

Microbiota assay

In both NCI and NYU, DNA was extracted from faecal samples using the Mo Bio PowerSoil DNA Isolation Kit (Mo Bio) with bead-beating, as previously reported( Reference Peters, Dominianni and Shapiro 16 ). In NCI, 16S ribosomal RNA (rRNA) gene amplicons covering variable regions V3–V4 were generated using the 347F-5′GG AGGCAGCAGTRRGGAAT′-3′ and 803R 5′-CTACCRGGGTATCTAATCC-3′ primer pair( Reference Ahn, Sinha and Pei 6 , Reference Nossa, Oberdorf and Yang 21 ). Amplicons were sequenced with the 454 Roche FLX Titanium pyrosequencing system, according to the manufacturer’s protocol. In NYU, 16S rRNA gene amplicons covering the V4 region were generated using the F515-5′GTGCCAGCMGCCGCGGTAA′-3′ and R806-5′GGACTACHVGGGTWTCTAAT-3′ primer pair( Reference Caporaso, Lauber and Walters 22 ). Amplicons were sequenced with the Illumina MiSeq platform.

Sequence data processing

Because two different sequencing platforms were used, we processed the sequence data separately. Sequences were demultiplexed, and poor-quality sequences excluded, using the default parameters of Quantitative Insights into Microbial Ecology (QIIME) script split_libraries.py (for NCI) or split_libraries_fastq.py (for NYU)( Reference Caporaso, Kuczynski and Stombaugh 23 ). Filtered sequence reads were clustered into de novo operational taxonomic units (OTU) at 97 % identity, and representative sequence reads for each OTU were assigned taxonomy based on fully sequenced microbial genomes (Integrated Microbial Genomes (IMG)/Greengenes(GG))( Reference Caporaso, Kuczynski and Stombaugh 23 ). Chimeric sequences were removed with ChimeraSlayer( Reference Haas, Gevers and Earl 24 ). Blinded quality control specimens in all sequencing batches showed good reproducibility: high intraclass correlation coefficients for the Shannon diversity index and abundances of bacterial phyla and genera have been previously reported( Reference Ahn, Sinha and Pei 6 , Reference Peters, Dominianni and Shapiro 16 , Reference Love, Huber and Anders 25 ).

Statistical analysis

α-Diversity

We evaluated the association between quartiles of fibre intake and within-subject microbial diversity (α-diversity) indices of Shannon diversity and evenness( Reference Shannon 26 ). In both studies, these indices were calculated in 100 iterations of rarefied OTU tables of 1000 sequence reads per sample. We modelled the Shannon index and evenness as outcomes in linear regression, adjusting for age, sex, race, categorical BMI and cigarette smoking status.

β-Diversity

We assessed the relationship of overall gut microbiota composition and quartiles of dietary fibre intake using weighted (quantitative) and unweighted (qualitative) phylogenetic UniFrac distance matrices( Reference Lozupone, Hamady and Kelley 27 ). Permutational multivariate ANOVA (PERMANOVA) of both weighted and unweighted UniFrac distances was used to evaluate whether fibre intake is associated with overall microbial community composition, after adjusting first for age, sex, race, categorical BMI and cigarette smoking status (adonis function, ‘vegan’ package in R)( Reference McArdle and Anderson 28 ). Principal coordinate analysis plots were generated using the first two principal coordinates, and labelled according to quartile of fibre intake.

Differential abundance testing

We assessed the relationship between higher quartiles of total fibre intake and specific taxa abundance using negative binomial generalised linear models, in the ‘DESeq2’ package in R( Reference Love, Huber and Anders 25 ). Models were adjusted for age, sex, race, categorical BMI and cigarette smoking status. Nominal P values and false discovery rate (FDR) adjusted q-values were calculated( Reference Benjamini and Hochberg 29 ). DESeq2 default outlier replacement, independent filtering of low-count taxa and filtering of count outliers were turned off. Taxa models with maximum Cook’s distance >10 were removed before P value adjustment for the FDR( Reference Benjamini, Drai and Elmer 30 ).

To identify similar taxa associations in both NCI and NYU, we then performed a taxonomy-based meta-analysis to evaluate for genera and species with higher or lower abundance by fibre intake concomitantly in each of the two study populations. In addition, we performed sub-analyses to examine associations between taxon abundance and higher intake of fibre from specific sources, such as fruit and vegetables, grains and beans. We calculated nominal meta-analysis P values based on Z-score methods( Reference Evangelou and Ioannidis 31 ).

All trends were tested using the median values of each quartile of fibre intake. All analyses were performed using R, version 3.3.2.

Results

Subject characteristics

A total of 151 adult subjects were included for analysis: n 75 in NCI and n 76 in NYU. Subject characteristics are reported in Table 1. Among these participants, 73·3 % in NCI and 51·3 % in NYU were male, and 82·7 % in NCI and 85·5 % in NYU were white. Median age in both study groups generally increased with higher fibre intake, although the trends were not statistically significant.

Table 1 Participant demographic characteristics (Quartiles (Q) and medians)

* NCI, National Cancer Institute; NYU, New York University.

* NCI quartiles of fibre intake: Q1: <11·3 g; Q2: 11·3–14·3 g; Q3: 14·4–16·6 g; Q4: ≥16·7 g.

P values are calculated from Kruskal–Wallis test (continuous variables) or χ 2 test (categorical variables).

‡NYU quartiles of fibre intake: Q1: <19·8 g; Q2: 19·8–23·7 g; Q3: 23·8–31·6 g; Q4: ≥31·7 g.

Global diversity

PERMANOVA analyses of between-sample UniFrac distances demonstrated that fibre intake was significantly associated with overall microbial community composition in NYU (weighted UniFrac P=0·008; unweighted UniFrac P=0·01) but not NCI (weighted UniFrac P=0·81; unweighted UniFrac P=0·75) after adjusting for covariates of age, sex, race, BMI and cigarette smoking (Fig. 1(c)). In contrast, in both NCI and NYU, total fibre intake was not significantly associated with microbial community diversity as measured by the Shannon diversity index or evenness (Fig. 1(a)).

Fig. 1 α and β diversity in relation to quartiles of total fibre intake (a) Shannon diversity index and (b) evenness index by quartiles of total fibre intake are shown in n 151 healthy adult subjects from two independent study populations (National Cancer Institute (NCI)=75, New York University (NYU)=76). These indices were calculated in 100 iterations of rarefied operational taxonomic units tables of 1000 sequence reads per sample. Fibre 1, fibre 2, fibre 3 and fibre 4 represent increasing quartiles of fibre intake. Shannon index and evenness were modelled as outcomes in linear regression, adjusting for age, sex, race, categorical BMI and cigarette smoking status. P values of fibre variable in regression analysis are reported in the figure. (c) Principal coordinate analysis (PCoA) plots, based on weighted UniFrac phylogenetic distances, showed a difference between lowest and highest fibre intake in NYU. This relationship was not observed in NCI. PCoA plots were generated using the first two principal coordinates (PC): PC1 and PC2. P values reported in the figure are based on permutational multivariate ANOVA of weighted UniFrac distances evaluating the association between fibre intake and overall microbial community composition, after adjusting for age, sex, race, categorical BMI and cigarette smoking status.

Taxon abundance

We found fourteen genera with the same direction of association with fibre intake in NCI and NYU, out of twenty-nine total genera that overlapped (i.e. were detected) in both studies. Higher total fibre intake was associated with lower abundance of genera Actinomyces (fold change (FC)=0·95, P=0·002) of class Actinobacteria and Odoribacter (FC=0·95, P=0·03) of class Bacteroidia, and was marginally associated with lower abundance of Oscillospira (FC=0·96, P=0·06) of class Clostridia. Higher total fibre intake also tended to be associated with higher abundance of selected genera of class Clostridia: SMB53 (FC=1·04, P=0·04), Lachnospira (FC=1·03, P=0·05) and Faecalibacterium (FC=1·03, P=0·06) (Fig. 2 and Table 2).

Fig. 2 Forest plot of fold change of select genera (G) and species (S) in both the New York University and National Cancer Institute study populations with a significant or marginally significant association with higher total fibre intake, based on meta-anlysis of the two study populations. Meta-analysis nominal P values were calculated based on Z-score methods. C, class.

Table 2 Meta-analysis of the association between total fibre intake and genera and species in National Cancer Institute (NCI) and New York University (NYU) study populations* (Mean values and 95 % confidence intervals)

FC, fold change; Q, false discovery rate-adjusted q-value; NA, non-applicable.

* Taxa shown were associated with fibre intake in the same direction in both NCI and NYU. Relationship between higher quartiles of total fibre intake and differential taxon abundance was evaluated using negative binomial generalised linear models in the DESeq2 package in R. Models were adjusted for age, sex, race, categorical BMI and cigarette smoking status. Meta-analysis was conducted using Z-score methods.

At the species level, we found eight species with the same direction of association with fibre intake in the two study populations, out of seventeen total species observed overlapping in both studies. A meta-analysis at this taxonomic level showed a marginal association between higher total fibre intake and higher abundance of Faecalibacterium prausnitzii (FC=1·03, P=0·07), and lower abundance of Eubacterium dolichum (FC=0·96, P=0·04) and Bacteroides uniformis (FC=0·97, P=0·05) (Fig. 2 and Table 2).

In a sub-analysis to identify associations between taxon abundance and fibre intake from specific dietary sources, we found that higher fibre intake from fruit and vegetables was associated with lower abundance of genera Actinomyces (FC=0·97, P=0·007), Odoribacter (FC=0·96, P=0·04) and Oscillospira (FC=0·99, P=0·06). Higher abundance of genus Faecalibacterium was most significant specifically with higher fibre intake from beans (FC=1·11, P=0·01). At the species level, higher fibre intake from beans was associated with higher abundance of F. prausnitzii (FC=1·11, P=0·01) and lower abundance of B. uniformis (FC=0·87, P=0·08). A heat map representing colour-coded FC of select genera and species, by total fibre and fibre from specific sources (fruits/vegetables, grains and beans), is presented in Fig. 3. Among fibre from different food sources, there was a minimal association between fibre from grains and these specific genera and species.

Fig. 3 Heat map representing colour-coded fold changes of select genera (G) and species (S), by total fibre and fibre from specific sources (fruits/vegetables, grains, and beans). * P<0·05 and ** P<0·001 in the meta-analysis of the National Cancer Institute and New York University study populations.

Although we observed similar associations with taxon abundance and fibre intake in both the NCI and NYU studies, we also noted some inconsistent associations (online Supplementary Table S1). For example, at the phylum level, higher total fibre intake was associated with higher abundance of Proteobacteria in NCI (P=0·03) but not in NYU (P=0·53). Higher total fibre intake was also marginally associated with higher abundance of phylum Bacteroidetes in NYU (P=0·07), but not in NCI (P=0·47).

Full results of the meta-analyses examining the association between taxon abundance at each level and dietary fibre intake (total fibre and fibre from different food sources) are reported in online Supplementary Tables S1–S4.

Discussion

In the present study, we examined the association of usual dietary fibre intakes with gut microbiota composition in adults from two independent study populations. In a taxonomy-based meta-analysis, we found that higher total fibre intake is associated with specific taxon abundances, including higher abundance of select genera of Clostridia class. Some of these fibre and taxon abundance associations were consistent with specific fibre food sources, such as fruit and vegetables, grains and beans.

The usual dietary intake of participants measured in our study provides additional insight into the potential effect of longer-term dietary patterns on gut microbiota composition, compared with controlled dietary intervention studies( Reference O’Keefe, Li and Lahti 12 , Reference Tap, Furet and Bensaada 13 ). Although Wu et al.( Reference Wu, Chen and Hoffmann 15 ) demonstrated in a controlled feeding study that certain dietary modifications altered microbiome composition in humans, they did not significantly alter the overall types of microbial communities characterising individuals; the overall structure may be better correlated with long-term diet. Sonnenburg et al.( Reference Sonnenburg, Smits and Tikhonov 32 ) also reported that chronic low-fibre diets in mice may drive changes in microbial diversity. Thus, our study associations derived from long-term dietary habits are more likely to represent microbial composition in an uncontrolled, real-world setting.

We found that higher fibre intake was associated with higher abundance of select genera of Clostridia class. This finding is notable given the particular role of Clostridium spp. in colonocyte metabolism through production of SCFA, including acetate, propionate and butyrate, via fermentation( Reference Liu, Wang and He 33 ). Although butyrate is the least abundant SCFA produced, it serves as the preferred energy source for colonocytes and has garnered attention for its involvement in cell signalling, maintenance of intestinal barrier function and modulation of immune cell chemotaxis and adhesion( Reference Liu, Wang and He 33 Reference Meijer, de Vos and Priebe 37 ). Butyrate exerts an anti-inflammatory effect through inhibition of inflammatory cytokines, particularly NF-κB( Reference Aguilar, Leonel and Teixeira 38 ). Furthermore, mouse models have shown that butyrate inhibits histone deacetylases and consequently inhibits gene expression and causes tumour suppression( Reference Donohoe, Holley and Collins 11 , Reference Donohoe, Collins and Wali 39 ), carrying implications for colorectal cancer treatment.

Within the Clostridia class, we noted a marginal association between higher fibre intake and higher abundance of species F. prausnitzii. The fibre and F. prausnitzii relationship has also been observed in a cross-sectional cohort of middle-age and older adults in Spain, in which greater adherence to a Mediterranean diet, rich in higher fibre content foods such as fruits, vegetables and whole grains, correlated with higher levels of F. prausnitzii, as well as Clostridium cluster XVIa( Reference Gutierrez-Diaz, Fernandez-Navarro and Salazar 40 ). Faecalibacterium prausnitzii, one of the most abundant species found in the gut and a key producer of butyrate, has been associated with anti-inflammatory activity( Reference Arumugam, Raes and Pelletier 41 Reference Sokol, Pigneur and Watterlot 43 ). Sokol et al. ( Reference Sokol, Pigneur and Watterlot 43 ) suggested that it exerts its anti-inflammatory effects on colitis mouse models, partly owing to associated metabolites blocking production of the pro-inflammatory molecules, NF-κB and IL-8. Moreover, a reduction in F. prausnitzii was associated with a higher risk of recurrence of Crohn’s disease in patients post resection. In addition to inflammatory bowel disease, Lopez-Siles et al.( Reference Lopez-Siles, Martinez-Medina and Suris-Valls 44 ) reported that lower levels of F. prausnitzii were also found in patients with colorectal cancer compared with healthy controls. Thus, these results suggest a potential therapeutic and preventive role of F. prausnitzii in countering microbial dysbiosis in human disease. Our finding of higher abundance of F. prausnitzii particularly with higher fibre intake from beans lends support to future investigation of specific diet modifications that could affect disease states.

We also observed that higher fibre intake was associated with lower abundance of genera Odoribacter, Actinomyces and Oscillospira, and higher abundance of genus Lachnospira. In mouse models, Zackular et al.( Reference Zackular, Baxter and Iverson 45 ) reported significant microbial shifts found within stool samples of mice with colon tumours, specifically with enrichment of taxa affiliated with members of Odoribacter. Along similar lines, Thomas et al.( Reference Thomas, Jesus and Lopes 46 ) examined human tissue samples collected during colonoscopy from rectal cancer cases and non-cancer controls, and noted higher abundance of Odoribacter in the cancer cases. Kasai et al.( Reference Kasai, Sugimoto and Moritani 47 ) noted higher proportions of several genera, including Actinomyces, in faecal samples from human subjects with colorectal carcinoma. Furthermore, higher abundance of Actinomyces has been reported in both colorectal adenoma and carcinoma cases( Reference Peters, Dominianni and Shapiro 16 , Reference Kasai, Sugimoto and Moritani 47 ). Whether or not enrichment or depletion of these specific microbes reflects cause or consequence of colon tumour development remains uncertain. Nonetheless, the opposite relationship observed between abundance of these microbes and higher fibre intake in our tumour-free adult population compared with cancer cases in prior studies is notable, and merits further examination of diet, microbiota and disease associations.

Similarly, we previously reported an association between lower abundance of Clostridia and both colorectal carcinoma and adenoma( Reference Ahn, Sinha and Pei 6 , Reference Peters, Dominianni and Shapiro 16 ). In a case–control study of colorectal cancer cases and controls, stool samples from colorectal cancer cases were characterised by depletion of phylum Firmicutes, predominantly of class Clostridia, relative to controls( Reference Ahn, Sinha and Pei 6 ). In addition, we noted a depleted abundance of members of Clostridia class in stool samples of colorectal adenoma cases, compared with controls( Reference Peters, Dominianni and Shapiro 16 ). These findings suggest a potential relationship between this shift in Clostridia composition and both pre-cancerous and cancerous events. Thus, our finding that higher dietary fibre intake is associated with a higher abundance of select genera of Clostridia class encourages future investigation of specific dietary strategies and their potential protective role against colorectal cancer risk.

Our results are consistent with our previous report of dietary fibre intake and gut microbiota relationships in the NCI study( Reference Dominianni, Sinha and Goedert 48 ). Our prior analysis of these data showed higher relative abundance of Clostridia class with higher fibre from fruit and vegetable intake, as well as higher relative abundance of phylum Actinobacteria with higher fibre from bean intake. In the present meta-analysis, we confirmed a higher abundance of select genera of Clostridia class with fibre from fruit and vegetables, and a higher abundance of phylum Actinobacteria with fibre from bean intake (online Supplementary Tables S2 and S4). Differences in reported findings between these two studies may be attributed to different statistical methods. Whereas our prior study categorised taxon abundance as a binary variable, in the present study we analysed taxon abundance as a continuous variable in the DESeq2 model, which may increase statistical power.

We note the use of different sequencing platforms for the NCI and NYU studies. As described previously, amplicons in NCI were sequenced using the 454 Roche FLX Titanium pyrosequencing system, whereas the Illumina MiSeq platform was used in NYU. Certain differences in read lengths and error rates have been reported among these two sequencing platforms( Reference Liu, Li and Li 49 ). Nonetheless, comparison of these two platforms in gut microbiome samples by Allali et al.( Reference Allali, Arnold and Roach 50 ) has demonstrated that, although there may be some differences in abundance of specific taxa, similar microbiome compositional profiles were still obtained, and similar biological conclusions could therefore be drawn across both platforms. In addition, we maintained a similar system of sequence data processing using the QIIME pipeline across both of our studies, to ensure consistency in analysis.

Several potential limitations also need to be considered. Measurement error is an inherent limitation in self-reported dietary assessment. Given that the NCI and NYU studies were independent investigations conducted during different time periods, dietary questionnaires utilised differed between studies, which probably contributed to noted differences in quantity of fibre intake in the two cohorts. For this reason, we used study-specific quartile trends of higher fibre intake to evaluate relationships with global microbial diversity and specific taxon abundance. Subjects in both study cohorts were mostly white, and thus our findings may not be generalisable to more racially diverse populations. The cross-sectional design of both studies limits assessment of the temporality of the diet–microbiota relationship. Furthermore, we lacked data in our study cohorts regarding faecal metabolites. Nonetheless, strengths of this study include the relatively large sample size, excellent quality control of microbiome assays in both studies and adjustment for potential confounders. Furthermore, our study hypothesis was tested in two independent study populations, and our meta-analysis findings identified similar taxon abundance associations across both populations.

In summary, we found that fibre intake was associated with gut microbiota composition in generally healthy adults. Fibre from specific food sources may also yield differential effects on gut microbiota. In our analysis, greater fibre intake, particularly from beans, tended to be associated with greater abundance of putatively beneficial bacteria, such as F. prausnitzii, which has demonstrated anti-inflammatory activity. Fibre intake has generally been acknowledged as beneficial in relation to human diseases, such as cancer, and recommendations for sufficient intake have been promoted by the United States Department of Agriculture( Reference Fungwe, Bente and Hiza 51 ). However, there remains a substantial gap in understanding the complex interplay between specific dietary fibre sources, gut microbiota and disease state. Consequently, further investigation is warranted to study microbiota-targeted dietary interventions in both a therapeutic and preventive disease setting. Given the mounting evidence that microbial dysbiosis may affect human health and contribute to the development of colorectal cancer, it is imperative to better elucidate the association of diet–microbiota inter-relationships and their potential impact on colorectal cancer risk. A deeper understanding of diet and the gut microbiome may lay the groundwork needed for future study of diet–microbiota-based colorectal cancer prevention strategies.

Acknowledgements

Research reported in this publication was supported in part by the US National Cancer Institute under award numbers U01CA182370, R01CA164964, R03 159414 and P30CA016087. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health, the US Department of Veterans Affairs or the United States Government. The National Cancer Institute had no role in the design, analysis or writing of this article.

D. L., B. A. P., R. B. H. and J. A. designed the research, conducted the study, analysed the data, performed statistical analyses and wrote the first draft of the manuscript. C. F., H. J. F., M. A. B., R. S., J. J. G., R. B. H. and J. A. provided reagents or essential materials. J. A. had primary responsibility for final content. All authors critically revised and approved the final manuscript.

The authors declare that there are no conflicts of interest.

Supplementary material

For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114518002465

References

1. Human Microbiome Project Consortium (2012) Structure, function and diversity of the healthy human microbiome. Nature 486, 207214.Google Scholar
2. Hold, GL, Smith, M, Grange, C, et al. (2014) Role of the gut microbiota in inflammatory bowel disease pathogenesis: what have we learnt in the past 10 years? World J Gastroenterol 20, 11921210.Google Scholar
3. Marchesi, JR, Dutilh, BE, Hall, N, et al. (2011) Towards the human colorectal cancer microbiome. PLoS ONE 6, e20447.Google Scholar
4. Han, M, Wang, C, Liu, P, et al. (2017) Dietary fiber gap and host gut microbiota. Protein Pept Lett 24, 388396.Google Scholar
5. Sobhani, I, Tap, J, Roudot-Thoraval, F, et al. (2011) Microbial dysbiosis in colorectal cancer (CRC) patients. PLoS ONE 6, e16393.Google Scholar
6. Ahn, J, Sinha, R, Pei, Z, et al. (2013) Human gut microbiome and risk for colorectal cancer. J Natl Cancer Inst 105, 19071911.Google Scholar
7. Chen, W, Liu, F, Ling, Z, et al. (2012) Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLOS ONE 7, e39743.Google Scholar
8. Weir, TL, Manter, DK, Sheflin, AM, et al. (2013) Stool microbiome and metabolome differences between colorectal cancer patients and healthy adults. PLOS ONE 8, e70803.Google Scholar
9. Bouvard, V, Loomis, D, Guyton, KZ, et al. (2015) Carcinogenicity of consumption of red and processed meat. Lancet Oncol 16, 15991600.Google Scholar
10. Aune, D, Chan, DS, Lau, R, et al. (2011) Dietary fibre, whole grains, and risk of colorectal cancer: systematic review and dose-response meta-analysis of prospective studies. BMJ 343, d6617.Google Scholar
11. Donohoe, DR, Holley, D, Collins, LB, et al. (2014) A gnotobiotic mouse model demonstrates that dietary fiber protects against colorectal tumorigenesis in a microbiota- and butyrate-dependent manner. Cancer Discov 4, 13871397.Google Scholar
12. O’Keefe, SJ, Li, JV, Lahti, L, et al. (2015) Fat, fibre and cancer risk in African Americans and rural Africans. Nat Commun 6, 6342.Google Scholar
13. Tap, J, Furet, JP, Bensaada, M, et al. (2015) Gut microbiota richness promotes its stability upon increased dietary fibre intake in healthy adults. Environ Microbiol 17, 49544964.Google Scholar
14. Liu, P, Zhao, J, Guo, P, et al. (2017) Dietary corn bran fermented by Bacillus subtilis MA139 decreased gut cellulolytic bacteria and microbiota diversity in finishing pigs. Front Cell Infect Microbiol 7, 526.Google Scholar
15. Wu, GD, Chen, J, Hoffmann, C, et al. (2011) Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105108.Google Scholar
16. Peters, BA, Dominianni, C, Shapiro, JA, et al. (2016) The gut microbiota in conventional and serrated precursors of colorectal cancer. Microbiome 4, 69.Google Scholar
17. Schiffman, MH, Van Tassell, RL, Robinson, A, et al. (1989) Case–control study of colorectal cancer and fecapentaene excretion. Cancer Res 49, 13221326.Google Scholar
18. Schiffman, MH, Andrews, AW, Van Tassell, RL, et al. (1989) Case–control study of colorectal cancer and fecal mutagenicity. Cancer Res 49, 34203424.Google Scholar
19. World Health Organization (2006) BMI Classification. Global Database on Body Mass Index. Geneva: WHO.Google Scholar
20. Willett, WC, Howe, GR & Kushi, LH (1997) Adjustment for total energy intake in epidemiologic studies. Am J Clin Nutr 65, Suppl. 4, 1220S1228S (discussion 9S–31S).Google Scholar
21. Nossa, CW, Oberdorf, WE, Yang, L, et al. (2010) Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J Gastroenterol 16, 41354144.Google Scholar
22. Caporaso, JG, Lauber, CL, Walters, WA, et al. (2011) Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A 108, Suppl. 1, 45164522.Google Scholar
23. Caporaso, JG, Kuczynski, J, Stombaugh, J, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7, 335336.Google Scholar
24. Haas, BJ, Gevers, D, Earl, AM, et al. (2011) Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res 21, 494504.Google Scholar
25. Love, MI, Huber, W & Anders, S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15, 550.Google Scholar
26. Shannon, CE (1997) The mathematical theory of communication. 1963. MD Comput 14, 306317.Google Scholar
27. Lozupone, CA, Hamady, M, Kelley, ST, et al. (2007) Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol 73, 15761585.Google Scholar
28. McArdle, BH & Anderson, MJ (2001) Fitting multivariate models to community data: a comment on distance‐based redundancy analysis. Ecology 82, 290297.Google Scholar
29. Benjamini, Y & Hochberg, Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Methodol 57, 289300.Google Scholar
30. Benjamini, Y, Drai, D, Elmer, G, et al. (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125, 279284.Google Scholar
31. Evangelou, E & Ioannidis, JP (2013) Meta-analysis methods for genome-wide association studies and beyond. Nat Rev Genet 14, 379389.Google Scholar
32. Sonnenburg, ED, Smits, SA, Tikhonov, M, et al. (2016) Diet-induced extinctions in the gut microbiota compound over generations. Nature 529, 212215.Google Scholar
33. Liu, H, Wang, J, He, T, et al. (2018) Butyrate: a double-edged sword for health? Adv Nutr 9, 2129.Google Scholar
34. Pryde, SE, Duncan, SH, Hold, GL, et al. (2002) The microbiology of butyrate formation in the human colon. FEMS Microbiol Lett 217, 133139.Google Scholar
35. Donohoe, DR, Garge, N, Zhang, X, et al. (2011) The microbiome and butyrate regulate energy metabolism and autophagy in the mammalian colon. Cell Metab 13, 517526.Google Scholar
36. Chen, J, Li, Y, Tian, Y, et al. (2015) Interaction between microbes and host intestinal health: modulation by dietary nutrients and gut-brain-endocrine-immune axis. Curr Protein Pept Sci 16, 592603.Google Scholar
37. Meijer, K, de Vos, P & Priebe, MG (2010) Butyrate and other short-chain fatty acids as modulators of immunity: what relevance for health? Curr Opin Clin Nutr Metab Care 13, 715721.Google Scholar
38. Aguilar, EC, Leonel, AJ, Teixeira, LG, et al. (2014) Butyrate impairs atherogenesis by reducing plaque inflammation and vulnerability and decreasing NFkappaB activation. Nutr Metab Cardiovasc Dis 24, 606613.Google Scholar
39. Donohoe, DR, Collins, LB, Wali, A, et al. (2012) The Warburg effect dictates the mechanism of butyrate-mediated histone acetylation and cell proliferation. Mol Cell 48, 612626.Google Scholar
40. Gutierrez-Diaz, I, Fernandez-Navarro, T, Salazar, N, et al. (2017) Adherence to a Mediterranean diet influences the fecal metabolic profile of microbial-derived phenolics in a Spanish cohort of middle-age and older people. J Agric Food Chem 65, 586595.Google Scholar
41. Arumugam, M, Raes, J, Pelletier, E, et al. (2011) Enterotypes of the human gut microbiome. Nature 473, 174180.Google Scholar
42. Louis, P & Flint, HJ (2009) Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine. FEMS Microbiol Lett 294, 18.Google Scholar
43. Sokol, H, Pigneur, B, Watterlot, L, et al. (2008) Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A 105, 1673116736.Google Scholar
44. Lopez-Siles, M, Martinez-Medina, M, Suris-Valls, R, et al. (2016) Changes in the abundance of Faecalibacterium prausnitzii phylogroups I and II in the intestinal mucosa of inflammatory bowel disease and patients with colorectal cancer. Inflamm Bowel Dis 22, 2841.Google Scholar
45. Zackular, JP, Baxter, NT, Iverson, KD, et al. (2013) The gut microbiome modulates colon tumorigenesis. mBio 4, e00692e00713.Google Scholar
46. Thomas, AM, Jesus, EC, Lopes, A, et al. (2016) Tissue-associated bacterial alterations in rectal carcinoma patients revealed by 16S rRNA community profiling. Front Cell Infect Microbiol 6, 179.Google Scholar
47. Kasai, C, Sugimoto, K, Moritani, I, et al. (2016) Comparison of human gut microbiota in control subjects and patients with colorectal carcinoma in adenoma: terminal restriction fragment length polymorphism and next-generation sequencing analyses. Oncol Rep 35, 325333.Google Scholar
48. Dominianni, C, Sinha, R, Goedert, JJ, et al. (2015) Sex, body mass index, and dietary fiber intake influence the human gut microbiome. PLOS ONE 10, e0124599.Google Scholar
49. Liu, L, Li, Y, Li, S, et al. (2012) Comparison of next-generation sequencing systems. J Biomed Biotechnol 2012, 251364.Google Scholar
50. Allali, I, Arnold, JW, Roach, J, et al. (2017) A comparison of sequencing platforms and bioinformatics pipelines for compositional analysis of the gut microbiome. BMC Microbiol 17, 194.Google Scholar
51. Fungwe, TV, Bente, L & Hiza, H (2007) The Food Supply and Dietary Fiber: Its Availability and Effect on Health: Nutrition Insight 36. Alexandria, VA: USDA Center for Nutrition Policy and Promotion.Google Scholar
Figure 0

Table 1 Participant demographic characteristics (Quartiles (Q) and medians)

Figure 1

Fig. 1 α and β diversity in relation to quartiles of total fibre intake (a) Shannon diversity index and (b) evenness index by quartiles of total fibre intake are shown in n 151 healthy adult subjects from two independent study populations (National Cancer Institute (NCI)=75, New York University (NYU)=76). These indices were calculated in 100 iterations of rarefied operational taxonomic units tables of 1000 sequence reads per sample. Fibre 1, fibre 2, fibre 3 and fibre 4 represent increasing quartiles of fibre intake. Shannon index and evenness were modelled as outcomes in linear regression, adjusting for age, sex, race, categorical BMI and cigarette smoking status. P values of fibre variable in regression analysis are reported in the figure. (c) Principal coordinate analysis (PCoA) plots, based on weighted UniFrac phylogenetic distances, showed a difference between lowest and highest fibre intake in NYU. This relationship was not observed in NCI. PCoA plots were generated using the first two principal coordinates (PC): PC1 and PC2. P values reported in the figure are based on permutational multivariate ANOVA of weighted UniFrac distances evaluating the association between fibre intake and overall microbial community composition, after adjusting for age, sex, race, categorical BMI and cigarette smoking status.

Figure 2

Fig. 2 Forest plot of fold change of select genera (G) and species (S) in both the New York University and National Cancer Institute study populations with a significant or marginally significant association with higher total fibre intake, based on meta-anlysis of the two study populations. Meta-analysis nominal P values were calculated based on Z-score methods. C, class.

Figure 3

Table 2 Meta-analysis of the association between total fibre intake and genera and species in National Cancer Institute (NCI) and New York University (NYU) study populations* (Mean values and 95 % confidence intervals)

Figure 4

Fig. 3 Heat map representing colour-coded fold changes of select genera (G) and species (S), by total fibre and fibre from specific sources (fruits/vegetables, grains, and beans). * P<0·05 and ** P<0·001 in the meta-analysis of the National Cancer Institute and New York University study populations.

Supplementary material: File

Lin et al. supplementary material

Table S1

Download Lin et al. supplementary material(File)
File 127 KB
Supplementary material: File

Lin et al. supplementary material

Table S2

Download Lin et al. supplementary material(File)
File 127.7 KB
Supplementary material: File

Lin et al. supplementary material

Table S3

Download Lin et al. supplementary material(File)
File 153.5 KB
Supplementary material: File

Lin et al. supplementary material

Table S4

Download Lin et al. supplementary material(File)
File 133.1 KB