Daily Consumption of Orange Juice from Citrus sinensis L. Osbeck cv

Jan 15, 2019 - Elisa Brasili*†‡ , Neuza Mariko Aymoto Hassimotto†‡ , Federica Del Chierico§ , Federico Marini∥ , Andrea Quagliariello§ , F...
1 downloads 0 Views 529KB Size
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