Subscriber access provided by Iowa State University | Library
Bioactive Constituents, Metabolites, and Functions
Daily consumption of orange juice from Citrus sinensis L. Osbeck cv. Cara Cara and cv. Bahia differently affects Gut Microbiota Profiling as unveiled by an integrated meta-omics approach Elisa Brasili, Neuza Mariko Aymoto Hassimotto, Federica Del Chierico, Frederico Marini, Andrea Quagliarello, Fabio Sciubba, Alfredo Miccheli, Lorenza Putignani, and Franco Maria Lajolo J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b05408 • Publication Date (Web): 15 Jan 2019 Downloaded from http://pubs.acs.org on January 16, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 38
Journal of Agricultural and Food Chemistry
Daily consumption of orange juice from Citrus sinensis L. Osbeck cv. Cara Cara and cv. Bahia differently affects Gut Microbiota Profiling as unveiled by an integrated metaomics approach Elisa Brasiliab*, Neuza Mariko Aymoto Hassimottoab, Federica Del Chiericoc, Federico Marinid, Andrea Quagliarielloc, Fabio Sciubbad, Alfredo Micchelid, Lorenza Putignanice, Franco Lajoloab
aDepartment
of Food Science and Experimental Nutrition, School of Pharmaceutical Science,
University of São Paulo, São Paulo, Brazil bFood
Research Center (FoRC), CEPID-FAPESP (Research Innovation and Dissemination
Centers Sao Paulo Research Foundation), São Paulo, Brazil cUnit
of Human Microbiome, Children's Hospital and Research Institute Bambino Gesù,
Rome, Italy dDepartment
eUnit
of Chemistry, University of Rome "La Sapienza", Rome, Italy
of Parasitology, Children's Hospital and Research Institute Bambino Gesù, Rome, Italy
Corresponding Author *(E.B.) Phone: (+55) 11-30910128. Fax: (+55) 11-3815 4410. E-mail:
[email protected] https://orcid.org/0000-0002-0623-2705
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
1
ABSTRACT
2
We have investigated the effect of intake of two different orange juices from Citrus sinensis
3
cv. ‘Cara Cara’ and cv. ‘Bahia’ on faecal microbiota and metabolome using an integrated
4
meta-omics approach. Following a randomized cross-over design, healthy subjects daily
5
consumed 500 mL of orange juice from Cara Cara or Bahia juices or an isocaloric control
6
drink. Stools were collected at baseline (T0) and after a week (T7) of intervention.
7
Operational taxonomic units (OTUs) were pyrosequenced targeting 16S rRNA and faecal
8
metabolites were analysed by an untargeted metabolomics approach based on 1H-NMR-
9
Spectroscopy. The major shift observed in microbiota composition after orange juice intake
Page 2 of 38
10
was the increased abundance of a network of Clostridia OTUs from Mogibacteriaceae,
11
Tissierellaceae, Veillonellaceae, Odoribacteraceae and Ruminococcaceae families, whose
12
members were differently affected by Cara Cara or Bahia juice consumption. A core of six
13
metabolites such as inositol, choline, lysine, arginine, urocanic acid and formate significantly
14
increased in Cara Cara compared to Bahia group.
15 16
Keywords: ‘Cara Cara’ orange juice, ‘Bahia’ orange juice, gut microbiota, Clostridia,
17
pyrosequencing, 1H-NMR-based metabolomics
18 19 20 21 22 23 24 25 2
ACS Paragon Plus Environment
Page 3 of 38
Journal of Agricultural and Food Chemistry
26
INTRODUCTION
27
The gut microbiota is a metabolically active community of microorganisms that interact with
28
one another and with the host organism influencing many aspects of human health.1,2
29
Genetic, lifestyle and environmental factors, including diet exert the greatest influence on gut
30
microbiota by modifying not only its composition but also its functionality to a largely
31
unknown extent. Emerging evidence suggests that exists a high interindividual variability in
32
gut microbiota composition and functions, making it difficult to predict the response of
33
individual’s gut microbiota to a given dietary intervention.3,4
34
Citrus sinensis L. Osbeck orange juice is an excellent dietary source of bioactive compounds
35
including vitamin C, flavonoids such as hesperidin and narirutin, carotenoids, sugars and
36
fibres that contribute to general health status and reduce the risk of metabolic and chronic
37
diseases.5-6 There are evidences that the bioactive components present in orange juice are
38
associated with the metabolism of gut microbiota.7 Most of ingested orange juice flavanones
39
reach the colon where are subjected to the action of the microorganisms able to hydrolyse
40
flavanone glycosides. A part of the released aglycones are absorbed and transformed to
41
glucuronide and sulphate metabolites and a further portion is excreted in the faeces.8,9 In this
42
context, Pereira-Caro et al.10, demonstrated that the aglycone hesperetin could be converted
43
into 3-(3'-hydroxy-4'-methoxyphenyl)propionic acid, 3-(3',4'-dihydroxyphenyl)propionic
44
acid, and 3-(3'-hydroxyphenyl)propionic acid, while 3-(phenyl)propionic acid is the major
45
end product derived from naringenin. In a recent in vitro study, it was demonstrated that fresh
46
and pasteurized orange juices increased Lactobacillus spp., Enterococcus spp.,
47
Bifidobacterium spp., and Clostridium spp. and reduced enterobacteria.11 These variations in
48
microbiota composition were also associated to an increase in butyric, acetic, and propionic
49
acid concentrations, and a decrease in ammonium production.11
3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 38
50
The metabolic relationships between nutrition and host-microbe co-metabolism are still
51
poorly characterized. Additionally, most of the studies focused on a single bioactive
52
compound and selected bacterial populations and there are not reports examining the impact
53
of orange juice on the human gut microbiota composition and metabolism.
54
Major advances in metabolomics and metagenomics technologies have opened a new
55
scenario in the comprehension of the gut ecosystem by shedding light on its composition,
56
modulation and interplay with microorganisms.12,13
57
Recently, we have analyzed the chemical composition of two different Brazilian orange
58
juices from Citrus sinensis L. Osbeck, namely ‘Bahia’ and ‘Cara Cara’. A higher ascorbic
59
acid content was detected in Bahia as compared to Cara Cara juice.14 Total flavanones content
60
as well as hesperidin and didymin levels were higher in Cara Cara with respect to Bahia
61
juice. Cara Cara was also characterized by a significantly higher sucrose content and a higher
62
and diversified carotenoid content compared to Bahia juice with a mixture of (Z)-isomers of
63
lycopene, all-E-β-carotene, phytoene, and phytofluene isomers accounting for the highest
64
carotenoid proportion. In addition, Cara Cara showed a higher content of insoluble fibre
65
compared to Bahia juice, whereas a higher content of soluble fibre was found in Bahia
66
compared to Cara Cara juice.14
67
It should be also considered that little is known about how the gut microbiota varies across
68
different countries and that Brazil presents a unique combination of large demographic
69
distribution and nutrition transition.
70
On account of these evidence, we investigated the effect of two different orange juices from
71
Citrus sinensis L. Osbeck ‘Bahia’ and ‘Cara Cara’ on Brazilian healthy volunteer gut
72
microbiota composition and activities using an integrated meta-omics approach. OTUs were
73
identified by pyrosequencing of 16S ribosomal RNA (rRNA) and metabolic profiling was
74
characterized by an untargeted NMR-based metabolomics approach. 4
ACS Paragon Plus Environment
Page 5 of 38
Journal of Agricultural and Food Chemistry
75 76
MATERIAL AND METHODS
77
Subjects and study design
78
The study protocol was approved by the Ethics Committee of University of São Paulo (CEP
79
FCF 1097907) and was formally registered as clinical trial (NCT02685124).
80
In this study, healthy volunteers, aged 18-45 years old and concordant for normal weight
81
(BMI= 18.5-24.9 kg/m2) were included. Before enrolling in the study, participants underwent
82
a health screening examination to exclude those with a history of gastrointestinal diseases,
83
chronic illness, taking medicaments or nutritional supplements, eating disorder or following
84
special diets, smokers and pregnant, which could interfere with the results of the study. A
85
complete blood exam of all participants was performed before they started the interventional
86
study.
87
A convenience sample of 21 participants (n=11 men and n=10 women) who have completed
88
a food-frequency questionnaire was included. All participants gave a written informed
89
consent.
90
Following a randomized controlled cross-over design (Figure S1), each participant consumed
91
500 mL of orange juice from Citrus sinensis L. Osbeck cv. Cara Cara or Bahia or an
92
isocaloric control drink containing water, sucrose and vitamin C (CTRLs), daily for a period
93
of 7 days with a 1-week washout period between consecutive interventions.
94
During the trial, participants recorded their dietary intake and were asked to follow a diet low
95
in polyphenolic compounds by avoiding fruits and vegetables, especially tomatoes and other
96
citrus fruits, chocolate, high-fibre products, and beverages, such as tea and coffee. A day
97
before the beginning and the end of each intervention period, participants were given a
98
standard low polyphenol meal (roasted chicken breast with potato puree). They were
99
instructed to continue the low polyphenol diet throughout the intervention period. 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
100
To analyse the gut microbiota composition and faecal metabolome, stool samples were
101
collected from each subject in sterile plastic pots at the time 0 (T0) and after 7 days (T7) of
102
intervention. Samples were stored at −80 °C until further analyses.
Page 6 of 38
103 104
Gut microbiota profiling
105
The genomic DNA was isolated from faecal samples using the QIAamp DNA Stool Mini Kit
106
(Qiagen, Germany) according to the protocol described by Del Chierico et al.15 Briefly, the
107
amplification of the 16S rRNA V1-V3 region was performed using a 454- Junior Genome
108
Sequencer (Roche 454 Life Sciences, Branford, CT). Reads were analysed with Quantitative
109
Insights into Microbial Ecology (QIIME v.1.8.0)16. To characterized the taxonomical
110
structure of the faecal samples, the sequences were organized into OTUs by clustering at a
111
threshold of 97% similarity. OTUs were matched against Greengenes database (v. 13.8) for
112
taxonomic assignment. The representative sequences were submitted to PyNAST for
113
sequence alignment16, and to UCLUST for sequence clustering.17
114 115
Faecal metabolites analysis by 1H NMR Spectroscopy
116
Faecal water samples were obtained by modifying published procedures.18,19 Briefly, 500 mg
117
of frozen faeces were mixed with ice-cold PBS/D2O (1.0 mL, pH 7.12, deuterium content
118
99.96%, PBS 0.1 mM). After centrifugation at 10,000×g for 25 min at 4° C, supernatants
119
were filtered through a cell strainer (40 μm pore size, BD bioscience) and further centrifuged.
120
PBS-D2O (200 μL) containing 2 mM trimethylsilyl-sodium-deuterated propionate (TSP), as
121
internal standard, and sodium azide solution 0.1% (v/v) were added to the supernatant (500
122
μL). Finally, samples were transferred into 5-mm NMR tubes and analysed by high-
123
resolution NMR spectroscopy. 1H-NMR spectra of faecal water were acquired at 298 K using
124
a Bruker AVANCE 400 spectrometer (Bruker BioSpin GmbH, Rheinstetten, Germany) 6
ACS Paragon Plus Environment
Page 7 of 38
Journal of Agricultural and Food Chemistry
125
equipped with a magnet operating at 9.4 T and 400.13 MHz for 1H frequency, as described
126
elsewhere.18,20 1H spectra were acquired employing the presat pulse sequence for solvent
127
suppression and resonance assignments were determined by comparison with literature
128
data21,22 as well as with online (Human Metabolome Database) 23 and in-house databases.
129
Two-dimensional (TOCSY and HSQC) NMR experiments were carried out to confirm the
130
assignments. 1H- and 2D-NMR spectra were processed and quantified as previously
131
described by Calvani et al.18 The quantification of metabolites was obtained by comparison of
132
the integrals of specific signals with the internal standard TSP integral. Concentrations of
133
metabolites are expressed in μmol/L.
134 135
Statistical Analysis
136
Ecological diversity for each intervention group was assessed by: i) number of obtained
137
OTUs; ii) Shannon index providing the entropy information of the OTU abundances and
138
accounts for both richness and evenness; iii) Chao1 metric estimating species richness; iv)
139
phylogenetic distance (PD_whole_tree) to assess quantitative measure of phylogenetic
140
diversity; v) Good’s coverage measuring what percent of the total species is represented in a
141
sample.
142
To measure the differences in taxonomic abundance profiles from different samples (β-
143
diversity), Bray-Curtis algorithm was used. Adonis test (R vegan package) was performed to
144
test for significant differences between groups (999 permutations).
145
The α and β diversity were performed using both QIIME software and phyloseq package in
146
R.24 Data normality was checked using the Shapiro-Wilk test by Sigma plot (V13.0). Since
147
the data distribution was not normal, nonparametric test (Mann-Whitney) was applied to
148
compare the data. Significantly different OTUs were selected by false discovery rate (FDR)-
149
adjusted p values. The meta-omic data from metagenomic and metabolomic were analysed 7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 38
150
using partial least squares discriminant analysis (PLS-DA). PLS-DA was performed using in-
151
house routines running under Matlab R2015b environment (The MathWorks, Natick, MA,
152
USA).
153
In order to obtain an unbiased estimate of the prediction performances, since the available
154
number of observation wouldn't allow for splitting the data into a single representative
155
training/test set pair, a repeated double cross-validation procedure was used. In detail,
156
cancelation groups were defined so to have an outer and an inner cross-validation loop, the
157
latter for model selection and the former to evaluate the performances of the optimized model
158
on a set of samples not used during the optimization stage; in order not to rely on a single
159
definition of the cross-validation splits, the procedure was repeated for 30 times and the
160
classification results have been averaged. Moreover, as a further validation stage, in order to
161
show that the classification results were significantly higher than those which could be
162
obtained by chance, the double-cross-validated predictions were compared to their respective
163
distribution under the null hypothesis, which in turn were estimated non-parametrically by
164
means of permutation tests with 1000 randomizations.
165
The significance of variables was evaluated using the variable importance in the projection
166
(VIP) method. Metabolites with VIP > 1.5 were considered significantly different and
167
selected for further analysis.
168 169
RESULTS
170
Gut microbiota profiling
171
A total of 63,627 16S rRNA sequencing reads were obtained from the 53 faecal samples. The
172
effect of consumption of orange juice from two cultivars on faecal microbiota was assessed
173
based on α-diversity indices and β-diversity (Bray-Curtis).
8
ACS Paragon Plus Environment
Page 9 of 38
Journal of Agricultural and Food Chemistry
174
The Mann-Whitney test was used to highlight significant differences in the α-diversity term
175
for the number of OTUs, Shannon index, and Chao1 index and the entire phylogenetic
176
distance tree for all three groups (Cara Cara, Bahia and CTRLs).
177
Significantly different (p< 0.05) Shannon index (Figure S2) and phylogenetic distance tree
178
(Table S1) were observed in CTRLs group between T0 and T7. Conversely, the mean
179
differences in the α-shannon index, Chao1 and entire phylogenetic distance tree values
180
between the T0 and T7 in Cara Cara and Bahia groups were not significant.
181
The β-diversity analysis highlighted the existence of a significant difference between T0 and
182
T7 for both Bahia and Cara Cara groups (p-value =0.037 and 0.024, respectively), while no
183
difference was detected among CTRLs samples (Figure 1).
184
The most of the OTU sequences were attributed to ten main phyla: Firmicutes, Bacteroidetes,
185
Cyanobacteria, Synergistetes, TM7, Fusobacteria, Tenericutes, Verrucomicrobia,
186
Actinobacteria and Proteobacteria. The most abundant phyla in all the volunteers were
187
Firmicutes, Bacteroidetes and Tenericutes (Figure 2). At time 0, all groups preserved the
188
same proportion between Firmicutes, Bacteroidetes and Tenericutes. Cara Cara group was
189
characterized by Firmicutes (67%), Bacteroidetes (29%) and Tenericutes (2%). Similarly,
190
Bahia group showed Firmicutes (70%), Bacteroidetes (23%) and Tenericutes (4%) levels.
191
Although the CTRLs group presented the higher Tenericutes level (5%), this one was not
192
significantly different from Cara Cara and Bahia group and all the observed little variations at
193
phylum level did not attain the statistical significance.
194
When we compared T0 and T7 within each intervention group, a significantly increase in
195
Firmicutes level was observed in CTRLs group (0.6960 versus 0.7995; p=0.0379) after intake
196
of the drink containing sucrose and vitamin C (Figure 2).
197
To evaluate the gut microbiota composition at the final time of each intervention period, a
198
comparison among the T7 of three groups was carried out. The Proteobacteria significantly 9
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 38
199
increased at T7 in Cara Cara compared to CTRLs group (0.01368 versus 0.00976; p =0.034)
200
(Figure 2).
201
Within the Firmicutes, the Mogibacteriaceae and Tissierellaceae families increased after
202
intake of Cara Cara orange juice, whereas among Bacteroidetes, the Odoribacteraceae family
203
and the Odoribacter genus decreased. The Mogibacteriaceae family statistically increased
204
also after intake of Bahia orange juice. In Bahia group, the Enterococcaceae family and
205
Veillonellaceae family increased while Ruminococcaceae family and the Faecalibacterium
206
prausnitzii decreased. Ruminococcaceae family and Oscillospira genus as well as
207
Porphyromonadaceae family and Parabacteroides distasonis significantly decreased in
208
CTRLs group (Table 1).
209
To highlight other differences among the treatments, we compared the final time (T7) of all
210
the groups (Table 2). Among the Actinobacteria, the relative abundance of the
211
Coriobacteriaceae family and Adlercreutzia genus were significantly higher in Bahia
212
compared to CTRLs group. Within the Bacteroidetes, the Porphyromonadaceae family and
213
Parabacteroides genus alike the Odoribacteraceae family and Butyricimonas genus were
214
increased in Cara Cara compared to CTRLs group. Among Firmicutes, a significant increase
215
of Enterococcaceae family and Enterococcus genus was observed in Bahia compared to
216
CTRLs group. Veillonellaceae family showed a significant increase across the groups CTRLs
217
< Cara Cara < Bahia. After Bahia juice intake, a significantly higher abundance of
218
Clostridiaceae family and Clostridium genus as well as Ruminococcaceae family and
219
Anaerotruncus genus compared to CTRLs group was also observed. In addition,
220
Anaerotruncus genus appeared to decrease across the groups Bahia > Cara Cara > CTRLs.
221
On the other hand, Cara Cara group showed a higher relative abundance of
222
Christensenellaceae, Lachnospiraceae and Ruminococcaceae families compared to CTRLs
10
ACS Paragon Plus Environment
Page 11 of 38
Journal of Agricultural and Food Chemistry
223
group. Interestingly, at species level, only Eubacterium dolichum showed evidence of
224
variation between groups, being detected exclusively in Cara Cara group.
225
Regarding the Proteobacteria, a higher abundance of the Enterobacteriaceae family was
226
observed in Cara Cara compared to CTRLs group, while a significant increase of
227
Pasteurellaceae was observed in Bahia compared to CTRLs group (Table 2).
228 229
1H-NMR-based
230
Typical 1H NMR and TOCSY spectra of faecal water at baseline (T0) are illustrated in
231
Figure S3 A and B respectively. We identified and quantified a series of faecal metabolites
232
including amino acids (glycine, phenylalanine, lysine, arginine), microbial-related
233
metabolites (isovalerate, acetate, butyrate, isobutyrate, propionate), alcohols (methanol,
234
ethanol, isopropanol), phenolic-derived compounds (3-(4-Hydroxyphenyl)-propionic acid),
235
histidine metabolism (urocanate) and others (oxoglutarate, inositol, dihydroxyacetone,
236
choline, p-cresol, formate, uracil, niacin) (Table S2).
237
In order to obtain a balanced data matrix for statistical analysis, we selected faecal water
238
samples at each experimental time point from volunteers that completed the study.
239
A preliminary analysis was conducted to characterize the differences among subjects (inter-
240
individual variation) and separate it from the treatment-related (intra-individual variation). To
241
investigate the effect of each treatment on participants (intra-individual variation), a PLS-DA
242
classification model was constructed and validated using post-intervention-baseline values
243
for pairs of classes. The relevance of individual metabolites was identified by inspecting the
244
values of the variable importance in projection (VIP) index. The first PLS-DA model “Bahia
245
versus Cara Cara” was built using three LVs that accounted for more than 65% of the
246
variance originally present in the X block. The model allowed to correctly predict the effect
247
of intake of Bahia and Cara Cara juices in 100.0% of participants in the calibration phase and
metabolic profiles of faecal water
11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 38
248
92.9% of participants in double cross-validation (outer loop; 100% Bahia and 87.5% Cara
249
Cara). The classification capacity of the model is also evident by inspecting the projection of
250
participants into the space spanned by the first three LVs of the PLS-DA model (Figure 3,
251
Figure S4), which shows a clear separation between “Bahia” and “Cara Cara” groups. A total
252
of six variables contributed significantly to the discrimination model, in particular choline,
253
inositol, urocanate, lysine, formate and arginine as shown in Table 3. The second PLS-DA
254
model “Bahia versus Control” was built using the first LV that accounted 41.46% of the
255
variance originally present in the X block (Figure 4). As indicated by the cross-validation
256
procedure, the model predicted the effect of Bahia juice and control drink intake in 86.7% of
257
participants in calibration (100% Bahia and 77.8% Control) and 80.0% in double cross-
258
validation (83.3% Bahia and 77.8% Control). Relevant variables that contributed
259
significantly to the discrimination model were isovalerate, choline, p-cresol, isopropanol,
260
phenylalanine + phenylacetate, formate and α-ketoglutarate (Table 3). The PLS-DA model
261
for Cara Cara and Control groups did not show a clear discrimination among the classes,
262
indicating an unsatisfactory model with low predictive ability (data not shown).
263 264
Integrated profiling of the metagenomics and metabolomic results
265
The data referring to OTUs and faecal metabolites were composed into a single data matrix
266
and processed using multivariate analysis. To prevent the misleading results potentially due
267
to the large quantity of null data, a new data matrix was assembled considering variables and
268
subjects populated more than 70%. The final data set included 56 variables (32 OTUs and 24
269
fecal metabolites). PLS-DA was used to identify an integrated metagenomics/metabolomics
270
profile for Cara Cara, Bahia and Control group. The PLS-DA analysis showed that only the
271
Cara Cara group could be significantly distinguished from CTRLs group within OTUs
272
variables (Table S3, Figure S5). Based on PLS-DA model, Cara Cara group had 12
ACS Paragon Plus Environment
Page 13 of 38
Journal of Agricultural and Food Chemistry
273
significantly higher levels of Enterobacteriaceae, Lachnospiraceae, Lachnospira,
274
Bacteroidaceae, Bacteroides ovatus, Porphyromonadaceae, Parabacteroides,
275
Ruminococcaceae, Oscillospira, Erysipelotrichaceae compared to CTRLs group. The PLS-
276
DA models built to compare Cara Cara versus Bahia group and Bahia versus CTRL group
277
did not display significant discriminating features between groups, indicating that these PLS-
278
DA models were not able to display statistically robust covariations of metagenomics and
279
metabolomics variables.
280
In order to explore the relationships between the bacteria and faecal metabolites after orange
281
juice intake, the Spearman rank test was performed in Cara Cara and Bahia groups. The
282
results are reported in Table 4.
283
In Cara Cara group, the Mogibacteriaceae family showed consistent positive correlations with
284
isobutyrate, isovalerate, α-ketoglutarate, inositol, choline, lysine, arginine, glycine and p-
285
cresol.
286
Butyrate positively correlated with Lachnospiraceae and Enterobacteriaceae. Similarly,
287
propionate showed a positive correlation with Lachnospiraceae family. Conversely,
288
Odoribacter genus negatively correlated with butyrate, isobutyrate and propionate in Cara
289
Cara group. A series of alcohols were also inversely correlated to Odoribacter after Cara
290
Cara juice intake, including ethanol, isopropanol and methanol.
291
In Bahia group, butyrate, propionate and acetate negatively correlated with Bacteroides.
292
Interestingly, desaminotyrosine showed consistent positive correlations with Bacteroides
293
uniformis and Odoribacter and a negative correlation with Lachnospira. Faecalibacterium
294
prausnitzii showed only a strong negative correlation with fumarate. Among the bacterial
295
products, urocanate was positively correlated with Lachnospiraceae in Cara Cara group and
296
with Coriobacteriaceae in Bahia group.
297 13
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 38
298
DISCUSSION
299
The most of dietary interventions in human cohorts have investigated microbial community
300
variations on timescales of weeks25 to months26. A limited number of researches examined
301
the effect of short-term consumption of diets on microbial community structure and
302
activities.27 Further studies are needed to establish how rapidly and reproducibly the human
303
gut microbiome responds to short-term diet change.
304
Orange juice is an important source of dietary bioactive compounds including polyphenols,
305
carbohydrates and fibres responsible for various health effects.
306
The association of citrus fruit bioactive compounds to gut microbiota metabolism has been
307
previously demonstrated using in vitro human cell lines and animal models.4,10,28-29 The
308
clinical trials in which citrus fruit bioactive compounds were associated to gut microbiota
309
community and metabolism are, however, very limited.30
310
In vivo studies carried out to investigate the effect of orange juice flavanone bioavailability in
311
healthy humans highlighted a large interindividual variability in the biological response to
312
(poly)phenols, probably due to differences in the gut microbiota or to epigenetic differences
313
between the different volunteers.31,32 Flavanone bioavailability seems to be dependent on the
314
presence of specific microorganisms able to remove the rutinosides from the juice glycosides,
315
which results in aglycones that are then absorbed from the gut. 29,31,33,34 In view of the above,
316
It is essential to study the correlations between gut microbiota composition and the metabolic
317
signature generated by citrus (poly)phenols assumption in both health and illness states.
318
The impact of two distinguished cultivars of orange juice from Citrus sinensis Osbeck ‘Cara
319
Cara’ and ‘Bahia’ characterized by a different content of vitamin C, flavanones and
320
carbohydrates on the human faecal microbiota and metabolome was investigated.
321
We observed that daily consumption of Cara Cara or Bahia orange juices differently affects
322
the gut microbiota profiling. 14
ACS Paragon Plus Environment
Page 15 of 38
Journal of Agricultural and Food Chemistry
323
The major shift observed in microbiota composition after orange juice intake was the
324
increased abundance of a network of Clostridia OTUs from Mogibacteriaceae,
325
Tissierellaceae, Veillonellaceae, Lachnospiraceae and Ruminococcaceae families, whose
326
members were differently affected by Cara Cara or Bahia juice consumption.
327
There are evidences that hesperidin and naringenin flavanones are specifically exposed to
328
anaerobic degradation by Clostridia bacteria, since they are able to synthesize enzymes that
329
cleaved the C ring of flavonoids.35
330
Although the Mogibacteriaceae family was not characterized at genus and species level, we
331
observed that it was statistically increased both after Cara Cara and after Bahia juice intake.
332
Not much is known about the function of the Mogibacteriaceae family in the gut. However,
333
in a recent study about faecal microbiota of healthy Japanese, it was observed that the
334
abundance of Mogibacteriaceae was significantly higher in lean than in obese subjects and
335
was also negatively associated with bowel movement frequency.36
336
In addition to unclassified members of Tissierellaceae that significantly discriminated T0 and
337
T7, the intervention with ‘Cara Cara’ juice also resulted in increased numbers of
338
Lachnospiraceae, Lachnospira and Dorea, as well as numbers of Ruminococcaceae,
339
Anaerotruncus and Oscillospira. Lachnospiraceae and Ruminococcaceae are the most
340
abundant Firmicutes families in gut environment, accounting for roughly 50% and 30% of
341
phylotypes, respectively.37 The decrease in the relative abundance of Lachnospiraceae
342
members in the gut microbiota was associated with compromised health status of patients
343
suffering from colorectal cancer 38, ulcerative colitis39, type 140 and type 2 diabetes.41
344
Likewise, the importance of the members of Ruminococcaceae family for the gut health is
345
indicated by their reduced abundance in faeces of Crohn's disease42 and ulcerative colitis
346
patients.39
15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 38
347
Members of Lachnospiraceae and Ruminococcaceae have in common the capacity to
348
generate SCFAs, mainly butyrate from fermentation of non-digestible plant fibers.43 Butyrate
349
exerts a key role in maintaining gut homeostasis and epithelial integrity, since it is used as the
350
primary energy source by colonic epithelial cells and then is transported to the mitochondria
351
to undergo β-oxidation to acetyl-CoA that enters the TCA cycle and results in ATP
352
production.43 Interestingly, after Cara Cara juice intake consistent positive correlations were
353
observed between butyrate and Lachnospiraceae family, as well as among the most abundant
354
SCFAs present in the colon, including acetate, butyrate and propionate.
355
The consumption of Bahia juice statistically increased the relative abundance of
356
Veillonellaceae but also Enterococcaceae and Coriobacteriaceae families.
357
Veillonella spp. are an essential component of the upper-gastrointestinal microbiota44 where
358
they form a trophic chain with Streptococcus spp able to convert lactate into propionate.45 Up
359
to now, it is not clear the role of Veillonella in human health, although several studies have
360
found an increase of Veillonella in faecal samples of patients suffering from irritable bowel
361
syndrome.46,47 The relative abundance of Faecalibacterium prausnitzii was lower after Bahia
362
juice intake, although no significant change was observed comparing T7 of all groups,
363
suggesting that its level was maintained constant among the subjects independently from the
364
dietary intervention. It is well known that F. prausnitzii is able to utilize simple carbohydrates
365
as glucose and fructose, as well as complex molecules as pectin and N-acetylglucosamine to
366
synthetize butyrate, formate and lactate.48
367
F. prausnitzii has also found to exhibit anti-inflammatory properties in mice 49 and was
368
already associated with a wide range of metabolic processes in the human gastrointestinal
369
tract.50 However, in this study, we observed that the selected orange juices were not potential
370
candidate foods to stimulate F. prausnitzii growth.
16
ACS Paragon Plus Environment
Page 17 of 38
Journal of Agricultural and Food Chemistry
371
The increase of Enterococcus (Enterococcaceae) and Adlercreutzia (Coriobacteriaceae)
372
indicates that Bahia orange juice, differently from Cara Cara one, is able to modulate also
373
Actinobacteria and Proteobacteria phyla. The data on Enterococcus spp. in human health are
374
contradictory. Enterococcus spp. are widely described as opportunistic pathogens, although
375
these species can exhibit probiotic abilities.51 Adlercreutzia spp. are well known gut species
376
able to produce equol, an isoflavandiol metabolized from daidzein.52 The ability to produce
377
equol is significantly higher in Asian as compared to Caucasian gut microbiota53 probably
378
due to the adaptation of the microbiota to the higher availability of isoflavones mainly
379
derived from soy beans. However, the consumption of soy and soy-based foods is also widely
380
spread in Brazil.54
381
After the Cara Cara juice intake, a modulation of faecal members of Bacteroidetes from
382
Bacteroidaceae, Barnesiellaceae and Porphyromonadaceae families was also observed. In this
383
group, the increase of Parabacteroides genus and Bacteroides ovatus was associated to the
384
decrease of Odoribacter genus. In previous studies, it was observed that tumour-bearing mice
385
showed enrichment in OTUs associated with members of the Odoribacter genus.55 In
386
humans, Bacteroides species are able to degrade plant fibre and complex polysaccharides,
387
including pectin and xylans derived from fruit and vegetable.56 It has been demonstrated that
388
Bacteroides ovatus, B. thetaiotaomicron and B. uniformis ferment a large number of complex
389
polysaccharides. In particular, B. ovatus metabolizes plant cell components, specifically
390
hemicelluloses of which the orange juice is rich, targeting the β1,4-glucosidic linkages.57
391
The importance of Parabacteroides genus is evinced by its ability to decrease the severity of
392
intestinal inflammation in animal models of acute and chronic colitis.58 Among species of
393
Parabacteroides genus, Parabacteroides distasonis, owns the ability to reduce the intestinal
394
inflammation by inducing the anti-inflammatory cytokine IL-10 and suppressing the secretion
395
of inflammatory cytokines such as IL-17, IL-6, and IFN-γ.59 Unexpectedly, in our study a 17
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 38
396
significant increase of this bacterium was observed after intake of control drink containing
397
vitamin C and sucrose.
398
The Bacteroides species exert also an important role in protein metabolism, since several
399
species show proteolytic activity, attributed to protease enzymes linked to the cell wall.60
400
As revealed by the PLS-DA model, a core of six host and microbial-derived metabolites
401
including choline, inositol, urocanate, lysine, arginine and formate significantly distinguished
402
Cara Cara from Bahia group. It is well known that colonic bacteria utilize the amino acids
403
including lysine, arginine, glycine as well as leucine, valine, and isoleucine (BCAA) to
404
generate a complex mixture of metabolic end products including acetate, butyrate, propionate
405
(SCFAs) and branched-chain fatty acids (BCFA; valerate, isobutyrate, and isovalerate).61
406
Importantly, these bacterial metabolites have been found to influence epithelial physiology by
407
influencing the signalling pathways in epithelial cells and by modulating the mucosal
408
immune system of the host.62 Although, no significant increase in SCFAs level was observed
409
after intake of Cara Cara or Bahia juice, we detected a significant increase in isovalerate and
410
a positive correlation between Mogibacteriaceae and isovalerate, α-ketoglutarate, choline,
411
lysine, arginine and glycine in Cara Cara group suggesting a link between these bacteria and
412
their metabolic products with chemical composition of Cara Cara orange juice. As microbial
413
metabolic product, urocanate is the first intermediate of the histidine degradation pathway
414
and can be produced by the bacterial histidine-ammonia lyase enzyme, which convert
415
histidine in urocanate and ammonia. The observed positive correlation between urocanate
416
and Lachnospiraceae or Coriobacteriaceae family, respectively in Cara Cara and Bahia group,
417
suggests that the two orange juices could affect the modulation of different bacteria involved
418
in histidine degradation reactions leading to the urocanate synthesis.
419
The orange juice consumption from Citrus sinensis Osbeck ‘Cara Cara’ and ‘Bahia’ induced
420
significant changes in the gut microbiota which were mainly associated to a core of Clostridia 18
ACS Paragon Plus Environment
Page 19 of 38
Journal of Agricultural and Food Chemistry
421
OTUs potentially involved in several host physiological processes including glucose and
422
protein metabolism, immune homeostasis and body energy balance.
423
Together, these findings demonstrate that the orange juice intake could modulate the host-
424
microbe interactions heralding a promising future therapeutic approach in nutrition studies
425
targeting chronic human diseases, including inflammatory bowel disease, obesity, type 2
426
diabetes, cardiovascular disease, and cancer. In this scenario, further research using large,
427
long-term intervention studies to evaluate the effect of Cara Cara or Bahia orange juice
428
consumption on gut microbiota composition and function would be helpful in making
429
specific dietary recommendations to population. As far as we know, this is the first study of
430
the Brazilian healthy volunteer faecal microbiota and metabolome after orange juice intake.
431 432
Funding
433
Funding has been provided by the Sao Paulo Research Foundation (FAPESP) (Grant
434
2016/08796-7) and the National Counsel of Technological and Scientific Development
435
(CNPq Grant 150033/2015-0).
436 437
The authors declare no competing financial interest.
438 439
ACKNOWLEDGMENTS
440
We thank Fundecitrus company in Araraquara (São Paulo, Brazil) for providing the orange
441
juices.
442 443 444
Supporting Information description 19
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 38
445
Figure S1. Experimental design.
446
Figure S2. Alpha diversity (Shannon index).
447
Figure S3. (A) Typical 1H-NMR spectra of faecal water at baseline (T0). (B) Typical TOCSY
448
spectrum of faecal water at baseline (T0).
449
Figure S4. PLS-DA loading plot of Bahia versus Cara Cara group obtained from NMR
450
metabolomics data.
451
Figure S5. Distribution of number of misclassifications (NMC), area under the ROC curve
452
(AUROC), and discriminant Q2 (DQ2) values under their respective null hypothesis as
453
estimated by permutation tests with 1,000 randomizations (blue histograms) and corresponding
454
values obtained by the PLS-DA model on unpermuted data (red dashed lines).
455
Table S1. Ecological diversity for each intervention group assessed by (i) number of obtained
456
OTUs; ii) Shannon index; iii) Chao1 metric estimating species richness; iv) phylogenetic
457
distance (PD_whole_tree); v) Good’s coverage.
458
Table S2. Assignments list: 1H chemical shifts of fecal metabolite signals.
459
Table S3. List of discriminant OTUs between Cara Cara and CTRLs groups as revealed by
460
PLS-DA model.
461 462 463
References 1) Holmes, E., Li, J. V., Athanasiou, T., Ashrafian, H., Nicholson, J. K. Understanding
464
the role of gut microbiome – host metabolic signal disruption in health and disease.
465
Trends Microbiol. 2011, 19, 349-59.
466 467 468 469
2) Ventura, M., Toole, P. W. O., Vos, W. M. De, Sinderen, D. Van. Selected aspects of the human gut microbiota. Cell Mol Life Sci. 2018, 75, 81-82. 3) Flint, H. J., Duncan, S. H., Louis, P. The impact of nutrition on intestinal bacterial communities. Curr Opin Microbiol. 2017, 38, 59-65. 20
ACS Paragon Plus Environment
Page 21 of 38
470
Journal of Agricultural and Food Chemistry
4) Espín, J. C., González-sarrías, A., Tomás-barberán, F. A., Espín, J. C., González-
471
sarrías, A., Tomás-barberán, F. A. The gut microbiota : A key factor in the therapeutic
472
effects of ( poly ) phenols. Biochem Pharmacol. 2017, 139, 82-93.
473
5) Crowe-white, K., Parrott, J. S., Stote, K. S., Gutschall, M. Benson-davies, S.
474
Metabolic Impact of 100 % Fruit Juice Consumption on Antioxidant / oxidant Status
475
and Lipid Profiles of Adults: An Evidence-based Review. Crit Rev Food Sci
476
Nutr. 2015, 57, 152-162.
477
6) Chaves, D. F. S., Carvalho, P. C., Brasili, E., Rogero, M. M., Mariko, N., Hassimotto,
478
A., Diedrich, J. K., Moresco, J. J., Yates, J. R., Lajolo, F. M. Proteomic Analysis of
479
Peripheral Blood Mononuclear Cells after a High-Fat, High-Carbohydrate Meal With
480
Orange Juice. J Proteome Res. 2017, 16, 4086-4092.
481 482 483
7) Laparra, J. M., Sanz, Y. Interactions of gut microbiota with functional food components and nutraceuticals. Pharmacol Res. 2010, 61, 219-25. 8) Borges, G., Mullen, W., Mullan, A., Lean, M. E., Roberts, S.A., Crozier, A.
484
Bioavailability of multiple components following acute ingestion of a polyphenol-rich
485
drink. Mol. Nutr. Food Res. 2010, 54, S268–S277.
486
9) Borges, G., Lean, M.E.J., Roberts, S.A., CrozierA. Bioavailability of dietary
487
(poly)phenols: a study with ileostomists to discriminate between absorption in small
488
and large intestine. Food Funct. 2013, 4, 754–762.
489
10) Pereira-caro, G., Borges, G., Ky, I., Ribas, A., Calani, L., Rio, D. Del, Clifford, M.
490
N., Roberts, S. A., Crozier, A. In vitro colonic catabo lism of orange juice ( poly )
491
phenols. Mol Nutr Food Res. 2015, 59, 465-75.
492
11) Rocha Faria Duque, A.L., Monteiro, M., Tallarico Adorno, M.A., Kimiko Sakamoto,
493
I., Sivieri, K. An exploratory study on the influence of orange juice on gut microbiota
494
using a dynamic colonic model. Food Res. Int. 2016, 84, 160-169. 21
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
495 496 497
12) Moco, S., Martin, F.P., Rezzi, S. Metabolomics View on Gut Microbiome Modulation by Polyphenol- rich Foods. J Proteome Res. 2012, 11, 4781-90. 13) Vernocchi, P., Del Chierico, F., Putignani, L. Gut Microbiota Profiling :
498
Metabolomics Based Approach to Unravel Compounds Affecting Human Health.
499
Front Microbiol. 2016, 7, 1144.
500
Page 22 of 38
14) Brasili, E., Chaves, D.F, Xavier, A.A., Mercadante, A.Z., Hassimotto, N.M., Lajolo,
501
F.M. Effect of Pasteurization on Flavonoids and Carotenoids in Citrus sinensis (L.)
502
Osbeck cv. 'Cara Cara' and 'Bahia' Juices. J Agric Food Chem. 2017, 65, 1371-1377.
503
15) Del Chierico, F., Nobili, V.; Vernocchi, P.; Russo, A.; Stefanis, C., Gnani, D. et al.
504
Gut microbiota profiling of pediatric nonalcoholic fatty liver disease and obese
505
patients unveiled by an integrated meta-omics-based approach. Hepatology. 2017, 65,
506
451-464.
507
16) Caporaso, J.G., Bittinger, K., Bushman, F.D., DeSantis, T.Z., Andersen, G.L. Knight,
508
R. PyNAST: a flexible tool for aligning sequences to a template alignment.
509
Bioinformatics, 2010, 26, 266-7.
510 511 512
17) Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26, 2460-1. 18) Calvani, R., Brasili, E., Praticò, G., Capuani, G., Tomassini, A., Marini, F., ... &
513
Miccheli, A. Fecal and urinary NMR-based metabolomics unveil an aging signature
514
in mice. Exp gerontol. 2014, 49, 5-11.
515
19) Saric, J., Wang, Y., Li, J., Coen, M., Utzinger, J., Marchesi, J. R., et al. Species
516
Variation in the Fecal Metabolome Gives Insight into Differential Gastrointestinal
517
Function research articles. J Proteome Res. 2008, 7, 352-60.
518
20) Miccheli, A., Marini, F., Capuani, G., Tomassini, A., Delfini, M., Di Cocco, M.E., et 22
ACS Paragon Plus Environment
Page 23 of 38
Journal of Agricultural and Food Chemistry
519
al. The Influence of a Sports Drink on the Postexercise Metabolism of Elite Athletes
520
as Investigated by NMR- Based Metabolomics The Influence of a Sports Drink on the
521
Postexercise Metabolism of Elite Athletes as Investigated by NMR-Based
522
Metabolomics. J Am Coll Nutr. 2009, 28, 553-64.
523
21) Martin, F.P., Wang, Y., Yap, I.K., Sprenger, N., Lindon, J.C., Rezzi, S., et al.
524
Topographical variation in murine intestinal metabolic
525
profiles in relationto microbiome speciation and functional ecological activity. J
526
Proteome Res. 2009, 8, 3464-74.
527
22) Tian, Y., Zhang, L., Wang, Y., Tang, H. Age-
528
related topographical metabolic signatures for the ratgastrointestinal contents. J
529
Proteome Res. 2012, 11, 1397-411.
530
23) Wishart, D.S; Jewison, T; Guo, A.C.; Wilson, M.; Knox, C.; et al. HMDB 3.0- The
531
Human Metabolome Database in 2013. Nucleic Acids Res. 2013, 41, D801-7.
532
24) McMurdie, P.J., Holmes, S. phyloseq: an R package for reproducible interactive
533
analysis and graphics of microbiome census data. PLoS One. 2013, 8, e61217.
534
25) Russell, W. R., Gratz, S. W., Duncan, S. H., Holtrop, G., Ince, J., Scobbie, L., ... &
535
Duthie, G. G. High-protein, reduced-carbohydrate weight-loss diets promote
536
metabolite profiles likely to be detrimental to colonic health–. Am J Clin Nutr. 2011,
537
93, 1062-1072.
538 539 540
26) Ley, R.E., Turnbaugh, P.J., Klein, S., Gordon, J.I. Microbial ecology: human gut microbes associated with obesity. Nature. 2006; 444:1022–1023. 27) David, L. A., Maurice, C. F., Carmody, R. N., Gootenberg, D. B., Button, J. E.,
541
Wolfe, B. E., ... & Biddinger, S. B. Diet rapidly and reproducibly alters the human gut
542
microbiome. Nature, 2014, 505, 559.
23
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
543
28) Unno, T., Hisada, T., Takahashi, S. Hesperetin Modifies the Composition of Fecal
544
Microbiota and Increases Cecal Levels of Short-Chain Fatty Acids in Rats. J Agric
545
Food Chem. 2015, 63, 7952-7
546
Page 24 of 38
29) Pereira-Caro, G., Borges, G., van der Hooft, J., Clifford, M.N., Del Rio D., Lean,
547
M.E. et al. Orange juice (poly)phenols are highly bioavailable in humans. Am J Clin
548
Nutr. 2014, 100, 1378-8.
549
30) Cuervo, A., Hevia, A., López, P., Suárez, A., Sánchez, B., Margolles, A., González,
550
S. Association of polyphenols from oranges and apples with specific intestinal
551
microorganisms in systemic lupus erythematosus patients. Nutrients. 2015, 7, 1301-
552
17.
553
31) Vallejo, F., Larrosa, M., Escudero, E., Zafrilla, M.P., Cerdá, B., Boza, J., García-
554
Conesa, M.T., Espín, J.C., Tomás-Barberán, F.A. Concentration and solubility of
555
flavanones in orange beverages affect their bioavailability in humans. J Agric Food
556
Chem. 2010, 58, 6516-24.
557 558 559
32) Tomas-Barberan, F. A., Selma, M. V., & Espín, J. C. Polyphenols’ Gut Microbiota Metabolites: Bioactives or Biomarkers?. J. Agric. Food Chem. 2018, 66, 3593-3594. 33) Pereira-Caro, G., Polyviou, T., Ludwig, I.A., Nastase, A.M., Moreno-Rojas, J.M.,
560
Garcia, A.L., Malkova, D., Crozier, A. Bioavailability of orange juice (poly)phenols:
561
the impact of short-term cessation of training by male endurance athletes. Am J Clin
562
Nutr. 2017, 106, 791-800.
563
34) Pereira-Caro, G., Fernández-Quirós, B., Ludwig, I.A., Pradas, I., Crozier,
564
A.,, Moreno-Rojas, J.M. Catabolism of citrus flavanones by
565
the probiotics Bifidobacteriumlongum and Lactobacillus rhamnosus. Eur J Nutr.
566
2018, 57, 231-242.
567
35) Schoefer, L., Mohan, R., Schwiertz, A., Braune, A., Blaut, M. 24
ACS Paragon Plus Environment
Page 25 of 38
Journal of Agricultural and Food Chemistry
568
Anaerobic degradation of flavonoids by Clostridium orbiscindens. Appl Environ
569
Microbiol. 2003, 69, 5849-54.
570
36) Oki, K., Toyama, M., Banno, T., Chonan, O., Benno Y., Watanabe K. Comprehensive
571
analysis of the fecal microbiota of healthy Japanese adults reveals a new bacterial
572
lineage associated with a phenotype characterized by a high frequency of bowel
573
movements and a lean body type. BMC Microbiol. 2016, 16:284.
574
37) Jalanka-Tuovinen, J., Salonen, A., Nikkilä, J., Immonen, O., Kekkonen, R., Lahti, L.,
575
et al. Intestinal microbiota in healthy adults: temporal analysis reveals individual and
576
common core and relation to intestinal symptoms. PLoS One, 2011, 6:e23035.
577
38) Wang, T., Cai, G., Qiu, Y., Fei, N., Zhang, M., Pang, X., et al. Structural segregation
578
of gut microbiota between colorectal cancer patients and healthy volunteers. ISME J.
579
2012, 6:320-9.
580
39) Rajilić-Stojanović, M., Shanahan, F., Guarner, F., de Vos, W.M. Phylogenetic
581
analysis of dysbiosis in ulcerative colitis during remission. Inflamm Bowel Dis. 2013,
582
19, 481-8.
583
40) Murri, M., Leiva, I., Gomez-Zumaquero, J.M., Tinahones, F.J., Cardona, F., Soriguer,
584
F. et al. (2013). Gut microbiota in children with type 1 diabetes differs from that in
585
healthy children: a case-control study. BMC Med. 2013, 11:46.
586
41) Qin, J.; Li, Y.; Cai, Z.; Li, S.; Zhu, J.; Zhang, F.; Liang, S.; Zhang, W.; Hansen, T.;
587
Sanchez, G.; et al. A metagenome-wide association study of gut microbiota in type 2
588
diabetes. Nature, 2012, 490, 55-60.
589
42) Kang, S., Denman, S.E., Morrison, M., Yu, Z., Dore, J., Leclerc, M. et al. Dysbiosis
590
of fecal microbiota in Crohn's disease patients as revealed by a custom phylogenetic
591
microarray. Inflamm Bowel Dis. 2010, 16, 2034-42.
592
43) Vital, M., Howe, C., Tiedje, M. 25
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 38
593
Revealing the Bacterial Butyrate Synthesis Pathways by Analyzing (Meta)genomic D
594
ata mBio. 2014, 5: e00889-14.
595
44) Van den Bogert, B., de Vos, W.M., Zoetendal, E.G., Kleerebezem, M. Microarray
596
analysis and barcoded pyrosequencing provide consistent microbial profiles
597
depending on the source of human intestinal samples. Appl Environ Microbiol. 2011,
598
77, 2071-80.
599
45) Zoetendal, E.G., Raes, J., van den Bogert, B., Arumugam
600
,M., Booijink, C.C., Troost, F.J., et al. The human small intestinal microbiota is
601
driven by rapid uptake and conversion of simple carbohydrates. ISME J. 2012, 6,
602
1415-26.
603
46) Tana, C., Umesaki, Y., Imaoka, A., Handa, T., Kanazawa, M., Fukudo, S. Altered
604
profiles of intestinal microbiota and organic acids may be the origin of symptoms in
605
irritable bowel syndrome. Neurogastroenterol Motil. 2010, 22, 512-9.
606
47) Saulnier, D. M., Riehle, K., Mistretta, T.-A., Diaz, M.-A., Mandal, D., Raza, S., …
607
Versalovic, J. Gastrointstinal microbiome signatures of pediatric patients with
608
irritable bowel syndrome. Gastroenterology, 2011, 141, 1782–1791.
609
48) Lopez-Siles, M., Khan, T.M., Duncan, S.H., Harmsen, H.J., Garcia-Gil, L.J., Flint,
610
H.J. Cultured representatives of two majo phylogroups of human colonic
611
Faecalibacterium prausnitzii can utilize pectin, uronic acids, and host-derived
612
substrates for growth. Appl Environ Microbiol. 2012, 78, 420-8.
613
49) Sokol, H., Pigneur, B., Watterlot, L., Lakhdari, O., Bermúdez-Humarán,
614
L.G., Gratadoux, J.J., et al. Faecalibacterium prausnitzii is an anti-
615
inflammatory commensal bacterium identified by gut microbiota analysis of Crohn
616
disease patients. Proc Natl Acad Sci U S A. 2008, 105, 16731-6.
617
50) Lepage, P.; Häsler, R.; Spehlmann, M. E.; Rehman, A.; Zvirbliene, A.; Begun, A., et 26
ACS Paragon Plus Environment
Page 27 of 38
Journal of Agricultural and Food Chemistry
618
al. Twin Study Indicates Loss of Interaction Between Microbiota and Mucosa of
619
Patients With Ulcerative Colitis. YGAST, 2011, 141, 227–236.
620
51) Cuív, P. Ó.; Klaassens, E. S.; Smith, W. J.; Mondot, S.; Durkin, A. S.; Harkins, D.
621
M.; Foster, L.; Mccorrison, J.; Torralba, M.; Nelson, K. E. Probiotic Strain Isolated
622
from Human Feces. 2013, 1, 12–13.
623
52) Maruo, T., Sakamoto, M., Ito, C., Toda, T. & Benno, Y. (2008). Adlercreutzia
624
equolifaciens gen. nov., sp. nov., an equol-producing bacterium isolated from human
625
faeces, and emended description of the genus Eggerthella. Int J Syst Evol Microbiol.
626
2008, 58, 1221–1227.
627
53) Song, K.B., Atkinson, C., Frankenfeld, C.L., Jokela, T., Wähälä, K., Thomas, W.T., et
628
al. Prevalence of Daidzein-Metabolizing Phenotypes Differs between Caucasian and
629
Korean American Women and Girls. J Nutr. 2006, 136, 1347–1351.
630
54) Genovese, M.I.; Lajolo, F.M. Isoflavones in soy-based foods consumed in
631
Brazil:levels, distribution, and estimated intake. J Agric Food Chem. 2002, 50, 5987–
632
93.
633 634 635
55) Louis, P., Hold, G.L., Flint, H.J. The gut microbiota, bacterial metabolites and colorectal cancer. Nat Rev Microbiol. 2014, 12, 661-72. 56) Henning, S.M.; Yang, J.; Shao, P.; Lee, R.P.; Huang, J.; Ly, A.; Hsu, M.; Lu,
636
Q.Y.; Thames, G.; Heber, D. Li, Z. Health benefit of vegetable/fruit juice-based diet:
637
Role of microbiome. Sci Rep. 2017, 7:2167.
638
57) Martens, E.C., Lowe, E.C., Chiang, H., Pudlo, N.A., Wu, M., McNulty, N.P et al.
639
Recognition and degradation of plant cell wall polysaccharides by two human gut
640
symbionts. PLoS Biol. 2011, 9:e1001221.
641
58) Kverka, M., Zakostelska, Z., Klimesova, K., Sokol, D., Hudcovic, T., Hrncir, T., et al.
642
Oral administration of Parabacteroides distasonis antigens attenuates experimental 27
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 28 of 38
643
murine colitis through modulation of immunity and microbiota composition. Clin Exp
644
Immunol. 2011, 163, 250-9.
645
59) Li, J., Sung, C.Y., Lee, N., Ni, Y., Pihlajamäki, J., Panagiotou, G., et al. Probiotics
646
modulated gut microbiota suppresses hepatocellular carcinoma growth in mice. Proc
647
Natl Acad Sci U S A. 2016, 113:E1306-15.
648 649 650 651 652 653
60) Macfarlane, G.T., Cummings, J.H., Allison, C. Protein degradation by human intestinal bacteria. J Gen Microbiol. 1986, 132, 1647-56. 61) Dai, Z.L., Wu, G., Zhu, W.Y. Amino acid metabolism in intestinal bacteria: links between gut ecology and host health. Front Biosci (Landmark Ed). 2011, 16, 1768-86. 62) Schaible, U.E., Kaufmann, S.H. A nutritive view on the host-pathogen interplay. Trends Microbiol. 2005, 13, 373-80.
654 655
28
ACS Paragon Plus Environment
Page 29 of 38
Journal of Agricultural and Food Chemistry
FIGURE CAPTIONS 656
Figure 1. Metagenomics beta-diversity plots. Principal component plots of the diversity in
657
the microbial communities between T0 and T7 faecal samples in Cara Cara (A) Bahia (B)
658
and CTRLs (C) groups.
659
Figure 2. Pie charts of OTU distribution. The distribution of the OTUs was compared
660
between T0 and T7 (*) within Cara Cara (A), Bahia (B) and CTRLs (C) groups, and between
661
T7 (**) of different groups (P-value