Subscriber access provided by Nottingham Trent University
Biotechnology and Biological Transformations
Modeling and regulation of higher alcohol production through the combined effects of the C/N ratio and microbial interaction Jian Jiang, Yuancai Liu, Huanhuan Li, Qiang Yang, Qun Wu, Shenxi Chen, Jie Tang, and Yan Xu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b04545 • Publication Date (Web): 03 Sep 2019 Downloaded from pubs.acs.org on September 5, 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 33
Journal of Agricultural and Food Chemistry
1
Modeling and regulation of higher alcohol production through the combined
2
effects of the C/N ratio and microbial interaction
3 4
Jian Jiang a, Yuancai Liu b, Huanhuan Li a, Qiang Yang b, Qun Wu a,*, Shenxi Chen
5
b,
6
a Key
7
Key Laboratory of Food Science and Technology, School of Biotechnology,
8
Jiangnan University, Wuxi, Jiangsu 214122, China
9
b Hubei
Jie Tang b, and Yan Xu a,*. Laboratory of Industrial Biotechnology of Ministry of Education, State
Provincial Key Laboratory for Quality and Safety of Traditional Chinese
10
Medicine Health Food, Jing Brand Research Institute, Jing Brand Co., Ltd,
11
Daye, Hubei 435100, China
12 13
* Corresponding author: Yan Xu, Qun Wu
14
Postal address: Key Laboratory of Industrial Biotechnology of Ministry of Education,
15
State Key Laboratory of Food Science and Technology, School of Biotechnology,
16
Jiangnan University, Wuxi, Jiangsu 214122, China
17
Tel: +86 510 85918201
18
Fax: +86 510 85864112
19
E-mail:
[email protected] or
[email protected] 20
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 33
21
ABSTRACT
22
Too large of a higher alcohol content has negative effects on the liquor taste and
23
health. Revealing the key microbes and their key driving forces is essential to regulate
24
the higher alcohol content in spontaneous liquor fermentation. Herein, we used high-
25
throughput sequencing associated with a multivariate statistical algorithm to reveal
26
the contributing microbes for higher alcohol production in Chinese light-aroma-type
27
liquor and identified that Saccharomyces and Pichia were the main contributors. In
28
addition, the C/N ratio and microbial interaction were found to significantly affect the
29
production of higher alcohols. Herein, we used response surface methodology to
30
establish a predictive model for higher alcohol production with the regulating factors,
31
and the content of total higher alcohols decreased significantly from 328.80 ± 24.83
32
mg/L to 114.88 ± 5.02 mg/L with the optimized levels of the regulators. This work
33
would facilitate the control of flavor production via regulating microbial communities
34
in food fermentation.
35 36
KEYWORDS: higher alcohols, microbial interaction, C/N ratio, Chinese Baijiu,
37
predictive modeling
38
2 ACS Paragon Plus Environment
Page 3 of 33
Journal of Agricultural and Food Chemistry
39
INTRODUCTION
40
Chinese Baijiu is one of the six famous distilled liquors in the world. The quality of
41
Chinese Baijiu is critically related to the complex flavor compounds.1 Higher alcohols
42
are a group of flavor compounds that coexist in fermented alcoholic beverages such as
43
brandy,2 whiskey,3 Chinese Baijiu,4, 5 wine,6 Huangjiu,7 beer8. The typical higher
44
alcohols in Chinese Baijiu are 2-methylpropanol and 3-methylbutanol.4, 9, 10 The
45
proper higher alcohol contents contribute to the formation of the rich flavor of
46
Chinese Baijiu.4, 11 However, too larger of a higher alcohol contents in liquor will
47
cause a bitter taste and off-flavors as well as neurotoxic effects.12-14 It is vital to
48
control the content of higher alcohols in Chinese Baijiu to improve the liquor quality.
49
The higher alcohols are generally considered to be produced by Saccharomyces
50
cerevisiae during the liquor fermentation process,15 and via the amino acid catabolic
51
pathway (Ehrlich pathway)16, 17 along with the sugar metabolism synthetic pathway
52
(Harris pathway).18, 19 Currently, methods for regulating higher alcohols in fermented
53
alcoholic beverages mainly focus on genetic modification of a single strain,20-22 such
54
as knockout of ILV1, LEU1, LEU2, BAT2 and IAH1 and overexpression of BAT1 and
55
ATF1 in S. cerevisiae. However, these single strain modification methods are not
56
applicable in spontaneous food fermentation containing multiple species. As a result,
57
it would be challenging to control the higher alcohols in a multi-species food
58
fermentation process.
59 60
Chinese Baijiu is produced by a typical solid-state spontaneous fermentation process.1 The contributing microbes for the production of higher alcohols are unclear. 3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
61
In addition, the microbial community and its metabolism are strongly influenced by
62
abiotic and biotic factors such as physiochemical factors and microbial interaction.23,
63
24
64
the microbial community, as well as the yield of metabolites during
65
biotransformation.25-27 However, the key controlling factors and their effects on
66
regulating higher alcohol production in multi-species communities are still unknown.
67
Thus, identifying the key higher alcohol-producing microbes and regulating their
68
metabolisms based on the key controlling factors are vital to regulate the higher
69
alcohols in liquor fermentation.
70
Page 4 of 33
For example, the C/N ratio in the substrate is an important factor which can affect
In this study, we revealed the key microbes and the key factors associated with the
71
higher alcohol production. Furthermore, we established a predictive model to regulate
72
the higher alcohols by using a response surface methodology with the key factors in a
73
synthetic microbial community. This modeling strategy is beneficial for the control of
74
flavor metabolism in food fermentation to improve the food quality.
75 76 77
MATERIALS AND METHODS Samples and Analytical Reagents. Samples were collected in May 2018 from a
78
famous mechanized liquor factory in Huangshi, Hubei province, China (29.76 N,
79
115.42 E). A mixture of the starter (Qu), raw material (sorghum), and fermented
80
grains from the last fermentation round (volume ratio = 1:100:300) was put into a
81
truck-shaped tank (1.7 m × 1.0 m× 1.0 m) and then sealed for a 15-day fermentation.
82
Samples (250 g) were collected from three separate fermentation tanks in the middle 4 ACS Paragon Plus Environment
Page 5 of 33
Journal of Agricultural and Food Chemistry
83
layer (0.5 m deep) on days 0, 1, 2, 3, 5, 7, 9, 11, and 14. These samples were stored at
84
-20 °C for DNA extraction and physicochemical analysis. Additional samples (100g)
85
from the sampling position were stored at 4 °C for strain screening.
86
Ethanol (99.9%, HPLC grade) and L-menthol (99.5%) were purchased from the
87
J&K Scientific Co., Ltd. (Beijing, China). Ethyl acetate, 2-methylpropanol, and 3-
88
methylbutanol were obtained from the Sigma-Aldrich Co., Ltd. (Shanghai, China)
89
with at least 99.0% purity.
90
Yeast Strains and Culture Procedure. One S. cerevisiae strain JN-JY-1 and one
91
Pichia kudriavzevii strain JN-JY-2 were used in this study. The yeast strains were
92
isolated from Chinese light-aroma-type liquor fermentation and deposited in the
93
China General Microbiological Culture Collection Center with the accession numbers
94
CGMCC 8130 and 12418. A single colony was obtained using Wallerstein-
95
laboratory-nutrient-agar (WL) medium (30 °C). The strains were activated using yeast
96
extract-peptone-dextrose (YPD) medium at 30 °C, 200 rpm, for 24 h. The mono- and
97
coculture were carried out using sorghum extract medium. Fermentation was
98
performed at 20 °C, 200 rpm, for 72 h.
99
Physicochemical Analysis. The temperature was measured in real time at the
100
sampling location during fermentation. The moisture of samples was measured by
101
drying 10 g of the fermented grains at 105 °C for 2 h to achieve a constant weight.
102
The acidity and amino acid nitrogen content of fermented grains were analyzed by
103
acid-base titration. To analyze the glucose content in the fermented grains, 10 g
104
samples were mixed with 25 mL of distilled water, ultrasonically treated at 0 °C for 5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
105
30 min, and then centrifuged at 4 °C and 12857 × g for 10 min. The obtained
106
supernatant was filtered through a 0.22 μm syringe filter (Nylon Acrodisc, Waters
107
Co., Milford, MA). The glucose contents were determined via high-performance
108
liquid chromatography (HPLC) (Agilent 1200 HPLC, Agilent Technology, Santa
109
Clara, CA) and a refractive index detector (RI) (WGE, Germany). The detection
110
conditions were as previously reported.28
111
Determination of the Higher Alcohol Concentration. The pretreatment method
112
for determining higher alcohol concentration in the fermented grains was the same as
113
determining glucose content. To analyze the higher alcohol content in the fermented
114
broth, 10 mL samples were centrifuged at 4 °C and 12857 × g for 10 min. The
115
supernatant was then collected and stored at -20 °C prior to its analysis.
116
Page 6 of 33
Higher alcohols were determined by gas chromatography-flame ionization
117
detection (GC-FID) as reported.29 The sampling conditions were as follows:
118
supernatant (1 mL), ethyl acetate (1 mL) and L-menthol (10 μL, 10.10 g/L, internal
119
standard) were mixed with 0.4 g of NaCl. After 3 minutes of oscillation, samples were
120
centrifuged at 4 °C and 3214 × g for 2 min. The obtained supernatant was filtered
121
through a 0.22 μm syringe filter (Nylon Acrodisc, Waters Co., Milford, MA).
122
GC-FID analysis was conducted on an Agilent 7890 GC system coupled to an
123
Agilent FID detector (Agilent, Folsom, CA). The supernatant (1 μL) was analyzed on
124
a CP-Wax column (30 m × 0.25 mm inner diameter, 0.25 μm film thickness; J&W
125
Scientific, Folsom, CA). The injector temperature was set at 250 °C, and the split
126
ratio was 5:1. The oven temperature was held at 60 °C for 3 min, then raised to 6 ACS Paragon Plus Environment
Page 7 of 33
Journal of Agricultural and Food Chemistry
127
135 °C at a rate of 6 °C/min, then raised to 200 °C at a rate of 30 °C/min, and finally
128
held at 200 °C for 3 min. Nitrogen was used as the column carrier gas at a constant
129
flow rate of 9 mL/min. The identification of compounds was carried out by comparing
130
their retention indices and higher alcohols standards.
131
DNA Extraction, Amplification, and Sequencing. Total genomic DNA was
132
extracted via the E.Z.N.A.® Soil DNA Kit (Omega Bio-tek, Norcross, GA) from 7 g
133
of each sample according to the manufacturer’s instructions. For fungi, the ITS2 was
134
amplified by polymerase chain reactions with the primers of ITS3 and ITS4.30 For
135
bacteria, the V3-V4 hypervariable region of the 16S rRNA gene was amplified by
136
PCR with the universal primers of the forward 338F and the reverse 806R.31 These
137
primers contained a set of barcode sequences unique to each sample. PCR products
138
were purified with a PCR purification kit (TaKaRa, Dalian, China), and their
139
concentrations were measured by the Thermo Scientific NanoDrop 8000 UV-Vis
140
Spectrophotometer (NanoDrop Technologies, Wilmington, DE). The barcoded PCR
141
products were merged into equimolar quantities and subjected to high-throughput
142
sequencing using a MiSeq Benchtop Sequencer for 2 × 300 bp paired-end sequencing
143
(Illumina, San Diego, CA) at Beijing Allwegene Tech, Ltd. (Beijing, China).32 All of
144
the sequence data were submitted in the DNA Data Bank of Japan (DDBJ) database
145
with the accession numbers of DRA008439 (16S) and DRA008440 (ITS2).
146
Sequence Processing. All of the raw sequences generated from MiSeq were
147
processed via the QIIME pipeline (v 1.8.0).33 Quality trimming was conducted by
148
removing the sequences with quality scores < 20 or lengths < 200 bp. Sequences that 7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 33
149
did not perfectly match the PCR primer or had non-assigned tags were also
150
removed.33 Chimeras were removed using the UCHIME software.34 Next, the
151
trimmed sequences were clustered into operational taxonomic units (OTUs) at 97%
152
similarity via UCLUST under the confidence threshold of 90%.35 Chao1 richness and
153
Shannon diversity indices were calculated by QIIME as previously reported.33
154
Real-Time Quantitative PCR (qPCR). To estimate the populations of P.
155
kudriavzevii and S. cerevisiae in mono- and co-cultures, absolute quantification qPCR
156
was performed on Real-Time PCR System (Applied Biosystems, Foster City, CA).
157
The genomic DNA was used as the template to amplify P. kudriavzevii by primers
158
PK1 and PK2, and S. cerevisiae by primers SC1 and SC2.36 Details about the primers
159
were listed in Table S4. Each reaction was performed in 20.0 μL containing 10.0 μL
160
SYBR Green Supermix (SYBR Premix ExTaq II, Takara, Dalian, China), 0.4 μL of
161
each primer (20 μM), 1.0 μL of DNA template and 8.2 μL ddH2O. The amplification
162
conditions and calibration curves were constructed as followed: preheating at 98 °C
163
for 3 mins, 40 cycles of 98 °C for 30 s, 60 °C for 30 s, and an increase of 0.5 °C every
164
5 s from 72 °C to 95 °C for melting curve analysis to confirm the specificity of the
165
amplification.
166
Box-Behnken Design. To model the effects of independent variables, including the
167
C/N ratio, pH and P. kudriavzevii/S. cerevisiae ratio, a response surface methodology
168
combined with Box-Behnken design was employed. The initial concentration of S.
169
cerevisiae was 1 × 105 CFU/mL. The initial concentration of P. kudriavzevii ranged
170
from 1 × 104 CFU/mL to 1 × 106 CFU/mL to explore the impact of the P. 8 ACS Paragon Plus Environment
Page 9 of 33
Journal of Agricultural and Food Chemistry
171
kudriavzevii/S. cerevisiae ratio. The total higher alcohol production was selected as
172
the response variable. In total, 17 runs were carried out to optimize the three
173
variables. The response surface (3D) and corresponding contour (2D) plots were
174
generated from the models. The optimum values of the variables were calculated from
175
the response surface (Design Expert software, version 8.0.6.1, Minneapolis,
176
Minnesota).
177
Data Analysis. Analysis of variance was conducted using one-way analysis of
178
variance (ANOVA). To establish the relationship between microbial communities and
179
higher alcohols, all possible Pearson correlation coefficients (ρ) among the dominant
180
genera (maximum relative abundance ≥ 1%) were calculated via SPSS Statistics 22.0
181
(IBM SPSS Statistics, Chicago, IL), |ρ| > 0.6 with statistically significant (P < 0.05)
182
was considered as a robust correlation.37 Gephi (Web Atlas, Paris, France) was used
183
to create the network to sort through and visualize the correlations.38 Redundancy
184
analysis (RDA) and the Monte Carlo permutation test were conducted based on
185
metabolites and the microbial communities via the vegan package in R (http://vegan.r-
186
forge.r-project.org/). ANOVA was carried out to examine the statistical importance of
187
regression coefficients. The coefficient of regression (R2) was calculated to find the
188
goodness of fit of the model. The F–test was used to assess the significance of the
189
model terms and equation.
190 191
RESULTS AND DISCUSSION
192
Dynamics of the higher alcohol contents during the liquor fermentation process. 9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
193
Figure 1 shows the dynamics of 2-methylpropanol and 3-methylbutanol in the
194
fermented grains. These two higher alcohols exhibited similar dynamic profiles, and
195
both reached the highest level on the 3rd day. The total higher alcohols increased
196
from 51.93 ± 3.06 mg/kg to 280.34 ± 33.82 mg/kg during the liquor fermentation
197
process. The accumulation of 3-methylbutanol reached 203.27 ± 23.00 mg/kg,
198
accounting for 72.5% of the total higher alcohol content.
Page 10 of 33
199
The production rate of the two higher alcohols increased rapidly during the initial
200
two days and then decreased until the end. The highest rates of 2-methylpropanol and
201
3-methylbutanol production were up to 29.74 mg/(kgday) and 79.57 mg/(kgday). As
202
a result, the whole fermentation process can be divided into two stages, stage I (0-2 d)
203
and stage II (2-14 d), based on the production rate of higher alcohols. Higher alcohols
204
mainly accumulated in stage I.
205 206
Identification of the key higher alcohol-producing microbes during the liquor
207
fermentation process. The microbial community diversity in the fermented grains
208
was studied via 16S rRNA and ITS sequence analysis. Across all samples, 995,309
209
reads (36,863 ± 12,221 on average) and 777,911 reads (28,812 ± 8,938 on average)
210
were identified for fungi and bacteria after quality control, respectively. Meanwhile,
211
we obtained 1,473 and 4,316 OTUs for fungi and bacteria, respectively, with 97%
212
similarity. The Good’s coverage for ITS2 and 16S rRNA genes was over 98.8%,
213
which indicated adequate sequencing depth of most of the samples (Table S1). The
214
fungal and bacterial α-diversity were quantified by the Chao1 richness and Shannon 10 ACS Paragon Plus Environment
Page 11 of 33
215 216
Journal of Agricultural and Food Chemistry
diversity (Figure S1). A total of 51 bacterial genera and 17 fungal genera were identified in the
217
fermentation process. Across all of the samples, only 17 bacterial genera and 9 fungal
218
genera were dominant (maximum relative abundance ≥ 1%) (Figure 2A). For fungal
219
communities, Rhizopus, Pichia and Saccharomyces were the predominant genera
220
(average abundance > 10%) in stage I. In stage II, Pichia and Saccharomyces were
221
still the predominant genera, while Rhizopus and Wickerhamomyces were the
222
subdominant genera (1% ≤ average abundance ≤ 10%).
223
For the bacterial community, Acetobacter, Gluconobacter and Lactobacillus were
224
the predominant genera in stage I, while Bacillus, Geobacillus, Klebsiella,
225
Leuconostoc, Pediococcus and Weissella were the subdominant genera. In stage II,
226
Lactobacillus was the only predominant genus, while Acetobacter, Gluconobacter,
227
Leuconostoc, Pediococcus and Weissella became the subdominant genera.
228
To identify the key higher alcohol-contributing microbes, we calculated all of the
229
possible Pearson correlation coefficients among the dominant genera and the contents
230
of 2-methylpropanol and 3-methylbutanol. Herein, seven genera were identified that
231
had significant correlations with the contents of the two higher alcohols (P < 0.05,
232
|ρ| > 0.6), including five fungal and two bacterial genera (Figure 2B). Among them,
233
Saccharomyces and Pichia had positive correlations with 2-methylpropanol and 3-
234
methylbutanol, indicating they might be the key higher alcohol-producers.
235 236
Identification of the key driving forces for microbial community succession 11 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 33
237
during the liquor fermentation process. The dynamics of six physicochemical
238
factors in the fermented grains during the liquor fermentation process were shown in
239
Figure S2, including the temperature, moisture, acidity, glucose, amino acid nitrogen,
240
and C/N ratio.
241
We used redundancy analysis to uncover the major abiotic factors driving the
242
microbial community succession (Figure 3). The correlations between
243
physicochemical factors and microbes were calculated via variation partitioning
244
analysis. Focusing on the main accumulation stage of higher alcohols (stage I), the six
245
factors accounted for 89.47% of the variation in the dominant microbes (Table S2).
246
The Monte Carlo permutation test indicated that the C/N ratio and acidity had
247
significant influences on the dominant microbes in both stage I and stage II (Table
248
S3). The factors, C/N ratio and acidity played vital roles in the succession of the
249
microbial community during the liquor fermentation process, which might be related
250
to the higher alcohol production.
251 252
Effect of microbial interaction on the higher alcohol production. As mentioned
253
above, we discovered that Saccharomyces and Pichia had the potential to produce
254
higher alcohols, and here we evaluated the higher alcohol production capacities of
255
these two potential producers based on culture-dependent methods. S. cerevisiae and
256
P. kudriavzevii were the dominant species of these two genera in liquor fermentation.
257
Therefore, S. cerevisiae JN-JY-1 and P. kudriavzevii JN-JY-2 were both isolated from
258
Chinese light-aroma-type liquor fermentation and used in the fermentation 12 ACS Paragon Plus Environment
Page 13 of 33
259 260
Journal of Agricultural and Food Chemistry
experiments. The biomass of the two yeasts tended to be steady after 24 h, whereas the
261
production of total higher alcohols stabilized after 60 h. After 72 h of mono-culture, S.
262
cerevisiae produced a total of 370.98 ± 18.75 mg/L higher alcohols. P. kudriavzevii
263
produced 38.82 ± 1.35 mg/L, accounting for approximately 10% of the production of
264
S. cerevisiae. S. cerevisiae was considered to be related to the production of higher
265
alcohols.17 Here we confirmed that P. kudriavzevii also had the ability to produce
266
higher alcohols. This would be in line with some other studies that some non-
267
conventional yeasts also have the potential to produce higher alcohols.39, 40
268
Microbial interaction is an important factor that affects the quality of fermented
269
foods.1 Next, we explored the impact of microbial interaction on higher alcohol
270
production. We developed a S. cerevisiae-P. kudriavzevii co-culture system to study
271
their interaction with respect to higher alcohol production.
272
The biomass of P. kudriavzevii was nearly the same in mono- and co-cultures,
273
whereas the biomass of S. cerevisiae decreased significantly in the co-cultures after
274
48 h compared with that in a mono-culture. In addition, the production of higher
275
alcohols decreased by 19.8% in co-cultures compared to the sum of the mono-cultures
276
of S. cerevisiae and P. kudriavzevii. Co-culturing affected the specific rate of
277
production of higher alcohols (Fig. 4D). These results indicated that the interaction
278
between S. cerevisiae and P. kudriavzevii could regulate the microbial metabolic
279
activity of higher alcohols in the artificial simulating system, which may deepen our
280
understanding of the regulation of higher alcohols. 13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 33
281 282
Predictive modeling for the regulation of higher alcohols. To further analyze the
283
regulation of biotic and abiotic factors on the production of higher alcohols, we
284
established a predictive model for the production of higher alcohols in the co-culture
285
system with S. cerevisiae and P. kudriavzevii using the C/N ratio, pH, and the
286
inoculation ratio as independent variables. Based on the response surface methodology
287
combined with Box-Behnken design, the empirical quadratic model was selected to
288
develop the correlation between the response of the higher alcohol contents and the
289
independent variables as follows:
290
Y (mg/L)=148.74+25.63A-9.39B-39.57C-16.87AB-
291
8.02AC+14.61BC+17.83A2+2.50B2+27.36C2.
292
In the equation, Y represents the total content of higher alcohols (mg/L), A is lg Pichia
293
(CFU/mL), B is the pH, and C is the C/N ratio. The F-value of the model was 13.66,
294
which implied that the model was significant. The P-value of the lack of fit (P = 0.1297)
295
and the model (P = 0.0012) revealed that the actual corresponding higher alcohol values
296
exhibited a good fit with this model (Fig. 5A). The correlation coefficient R2 was used
297
to assess the accuracy of the regression models established, with a closer value to unity
298
indicating a more precise response value estimated by the models.41 In this case, the R2
299
was 0.946, indicating the accuracy of the regression model.
300
The effects of the tested factors on the control of higher alcohols were visualized in
301
the response surfaces (Figure 6). The effects of the pH may be limited, as the P-value
302
of the pH was over 0.05 (Table 2), whereas the factors lg Pichia and C/N ratio had 14 ACS Paragon Plus Environment
Page 15 of 33
Journal of Agricultural and Food Chemistry
303
significant effects on the higher alcohol production (P < 0.05). The total content of
304
higher alcohols decreased to a minimal level (114.88 ± 5.02 mg/L) when the C/N ratio,
305
pH and P. kudriavzevii/S. cerevisiae ratio were 271.53:1, 3.59 and 104.08:105,
306
respectively. The content of total higher alcohols decreased by 65.1%, compared with
307
the co-culture of P. kudriavzevii and S. cerevisiae before optimization. This indicated
308
that this established predictive model would be efficient to predict and control higher
309
alcohol production. Besides, the model created in this study was based on a synthetic
310
microbial community including two yeasts. In the real liquor fermentation system, the
311
microbial species and interactions are more complex, deeper research is needed in the
312
future to explain the metabolism mechanisms of higher alcohols, which can facilitate
313
the regulation the higher alcohol contents more efficient in the liquor fermentation.
314
In this work, we revealed that S. cerevisiae and P. kudriavzevii were both higher
315
alcohol-contributing microbes. Moreover, the C/N ratio and microbial interaction were
316
found to significantly affect the higher alcohol production. Hence, we established a
317
predictive model for the production of higher alcohols based on these regulators and
318
efficiently controlled the higher alcohol production. This work provides an efficient
319
strategy for controlling flavor metabolism in food fermentations and is beneficial for
320
improving food quality.
321 322
ACKNOWLEDGEMENTS
323
We thank Bowen Wang and Shilei Wang for the use of R and Design Expert.
324 325
Funding 15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 33
326
This work was supported by the National Key R&D Program of China
327
(2018YFD0400402), the National Natural Science Foundation of China (31530055),
328
the Jiangsu Province Science and Technology Project (BE2017705), the Postgraduate
329
Research & Practice Innovation Program of Jiangsu Province (SJCX19_0772), the
330
China Postdoctoral Science Foundation (2017M611702), and the National First-Class
331
Discipline Program of Light Industry Technology and Engineering (LITE2018-12),
332
and the Priority Academic Program Development of Jiangsu Higher Education
333
Institutions, the 111 Project (No. 111-2-06).
334
Notes
335
We declare no conflicts of interest.
336 337
ASSOCIATED CONTENT
338
Supporting Information
339
Chao1 indices, Shannon indices of microbial communities, dynamics of
340
physicochemical factors during the liquor fermentation process, mash sample
341
statistical information about Illumina sequencing results, the interpretation rate of
342
single environment factor, the Monte Carlo permutation test of the dominant
343
microbes, and details about the primers are showed in the Supporting Information.
344
16 ACS Paragon Plus Environment
Page 17 of 33
Journal of Agricultural and Food Chemistry
345
REFERENCES
346
1.
347
Food Sci. Technol. 2017, 63, 18-28.
348
2.
349
Giarnpaoli, P. Evolution of volatile compounds during the distillation of cognac spirit.
350
J. Agric. Food Chem. 2017, 65, 7736-7748.
351
3.
352
composition analysis and authentication of whisky. J. Sci. Food Agric. 2015, 95,
353
2159-2166.
354
4.
355
type Chinese liquor by gas chromatography-olfactometry, quantitative measurements,
356
aroma recombination, and omission studies. J. Agric. Food Chem. 2014, 62, 5796-
357
5804.
358
5.
359
potent odorants in Chinese roasted sesame-like flavor type liquor by headspace solid
360
phase microextraction-aroma extract dilution analysis, with special emphasis on
361
sulfur-containing odorants. J. Agric. Food Chem. 2017, 65, 123-131.
362
6.
363
Y.; Reynolds, A. G.; Duan, C. Characterization and differentiation of key odor-active
364
compounds of ‘Beibinghong’ icewine and dry wine by gas chromatography-
365
olfactometry and aroma reconstitution. Food Chem. 2019, 287, 186-196.
366
7.
Jin, G.; Zhu, Y.; Xu, Y. Mystery behind Chinese liquor fermentation. Trends
Awad, P.; Athès, V.; Decloux, M. E.; Ferrari, G.; Snakkers, G.; Raguenaud, P.;
Wiśniewska, P.; Dymerski, T.; Wardencki, W.; Namieśnik, J. Chemical
Gao, W.; Fan, W.; Xu, Y. Characterization of the key odorants in light aroma
Sha, S.; Chen, S.; Qian, M.; Wang, C.; Xu, Y. Characterization of the typical
Lan, Y.; Xiang, X.; Qian, X.; Wang, J.; Ling, M.; Zhu, B.; Liu, T.; Sun, L.; Shi,
Yang, Y.; Xia, Y.; Lin, X.; Wang, G.; Zhang, H.; Xiong, Z.; Yu, H.; Yu, J.; Ai, L. 17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 33
367
Improvement of flavor profiles in Chinese rice wine by creating fermenting yeast with
368
superior ethanol tolerance and fermentation activity. Food Res. Int. 2018, 108, 83-92.
369
8.
370
reaction products and yeast strain on the synthesis of key higher alcohols and esters in
371
beer fermentations. Food Chem. 2017, 232, 595-601.
372
9.
373
aroma-type liquor by gas chromatography-olfactometry, quantitative measurements,
374
aroma recombination, and omission studies. J. Agric. Food Chem. 2015, 63, 3660-
375
3668.
376
10. Wang, X.; Fan, W.; Xu, Y. Comparison on aroma compounds in Chinese soy
377
sauce and strong aroma type liquors by gas chromatography-olfactometry, chemical
378
quantitative and odor activity values analysis. Eur. Food Res. Technol. 2014, 239,
379
813-825.
380
11. Fan, W.; Qian, M. C. Characterization of aroma compounds of Chinese
381
“Wuliangye” and “Jiannanchun” liquors by aroma extract dilution analysis. J. Agric.
382
Food Chem. 2006, 54, 2695-2704.
383
12. McKarns, S. C.; Hansch, C.; Caldwell, W. S.; Morgan, W. T.; Moore, S. K.;
384
Doolittle, D. J. Correlation between hydrophobicity of short-chain aliphatic alcohols
385
and their ability to alter plasma membrane integrity. Fundam. Appl. Toxicol. 1997, 36,
386
62-70.
387
13. Strubelt, O.; Deters, M.; Pentz, R.; Siegers, C. P.; Younes, M. The toxic and
388
metabolic effects of 23 aliphatic alcohols in the isolated perfused rat liver. Toxicol.
Dack, R. E.; Black, G. W.; Koutsidis, G.; Usher, S. J. The effect of Maillard
Fan, H.; Fan, W.; Xu, Y. Characterization of key odorants in Chinese chixiang
18 ACS Paragon Plus Environment
Page 19 of 33
Journal of Agricultural and Food Chemistry
389
Sci. 1999, 49, 133-142.
390
14. Yang, D.; Luo, X.; Wang, X. Characteristics of traditional Chinese shanlan wine
391
fermentation. J. Biosci. Bioeng. 2014, 117, 203-207.
392
15. Wu, Q.; Kong, Y.; Xu, Y. Flavor profile of Chinese liquor is altered by
393
interactions of intrinsic and extrinsic microbes. Appl. Environ. Microbiol. 2016, 82,
394
422-430.
395
16. Sentheshanuganathan, S. The mechanism of the formation of higher alcohols
396
from amino acids by Saccharomyces cerevisiae. Biochem. J. 1960, 74, 568-576.
397
17. Hazelwood, L. A.; Daran, J.-M.; van Maris, A. J. A.; Pronk, J. T.; Dickinson, J.
398
R. The ehrlich pathway for fusel alcohol production: a century of research on
399
Saccharomyces cerevisiae metabolism. Appl. Environ. Microbiol. 2008, 74, 2259-
400
2266.
401
18. Giudici, P.; Romano, P.; Zambonelli, C. A biometric study of higher alcohol
402
production in Saccharomyces cerevisiae. Can. J. Microbiol. 1990, 36, 61-64.
403
19. Choi, Y. J.; Lee, J.; Jang, Y.-S.; Lee, S. Y. Metabolic engineering of
404
microorganisms for the production of higher alcohols. mBio 2014, 5, e01524-14.
405
20. Li, W.; Chen, S.; Wang, J.; Zhang, C.; Shi, Y.; Guo, X.; Chen, Y.; Xiao, D.
406
Genetic engineering to alter carbon flux for various higher alcohol productions by
407
Saccharomyces cerevisiae for Chinese Baijiu fermentation. Appl. Microbiol.
408
Biotechnol. 2018, 102, 1783-1795.
409
21. Ma, L.; Huang, S.; Du, L.; Tang, P.; Xiao, D. Reduced production of higher
410
alcohols by Saccharomyces cerevisiae in red wine fermentation by simultaneously 19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 33
411
overexpressing BAT1 and deleting BAT2. J. Agric. Food Chem. 2017, 65, 6936-6942.
412
22. Li, W.; Wang, J.; Zhang, C.; Ma, H.; Xiao, D. Regulation of Saccharomyces
413
cerevisiae genetic engineering on the production of acetate esters and higher alcohols
414
during Chinese Baijiu fermentation. J. Ind. Microbiol. Biotechnol. 2017, 44, 949-960.
415
23. Frindte, K.; Pape, R.; Werner, K.; Löffler, J.; Knief, C. Temperature and soil
416
moisture control microbial community composition in an arctic-alpine ecosystem
417
along elevational and micro-topographic gradients. ISME J. 2019.
418
24. Yu, K.; Yi, S.; Li, B.; Guo, F.; Peng, X.; Wang, Z.; Wu, Y.; Alvarez-Cohen, L.;
419
Zhang, T. An integrated meta-omics approach reveals substrates involved in
420
synergistic interactions in a bisphenol A (BPA)-degrading microbial community.
421
Microbiome 2019, 7, 16.
422
25. Carsanba, E.; Papanikolaou, S.; Erten, H. Production of oils and fats by
423
oleaginous microorganisms with an emphasis given to the potential of the
424
nonconventional yeast Yarrowia lipolytica. Crit. Rev. Biotechnol. 2018, 38, 1230-
425
1243.
426
26. Bahram, M.; Hildebrand, F.; Forslund, S. K.; Anderson, J. L.; Soudzilovskaia, N.
427
A.; Bodegom, P. M.; Bengtsson-Palme, J.; Anslan, S.; Coelho, L. P.; Harend, H.;
428
Huerta-Cepas, J.; Medema, M. H.; Maltz, M. R.; Mundra, S.; Olsson, P. A.; Pent, M.;
429
Põlme, S.; Sunagawa, S.; Ryberg, M.; Tedersoo, L.; Bork, P. Structure and function
430
of the global topsoil microbiome. Nature 2018, 560, 233–237.
431
27. Syaichurrozi, I. Biogas production from co-digestion Salvinia molesta and rice
432
straw and kinetics. Renewable Energy 2018, 115, 76-86. 20 ACS Paragon Plus Environment
Page 21 of 33
Journal of Agricultural and Food Chemistry
433
28. Wu, Q.; Chen, L.; Xu, Y. Yeast community associated with the solid state
434
fermentation of traditional Chinese Maotai-flavor liquor. Int. J. Food Microbiol.
435
2013, 166, 323-330.
436
29. López-Vázquez, C.; Bollaín, M. H.; Berstsch, K.; Orriols, I. Fast determination of
437
principal volatile compounds in distilled spirits. Food Control 2010, 21, 1436-1441.
438
30. Levin, R. A.; Myers, N. R.; Bohs, L. Phylogenetic relationships among the “spiny
439
solanums” (Solanum subgenus Leptostemonum, Solanaceae). Am. J. Bot. 2006, 93,
440
157-169.
441
31. Zhang, X.; Tian, X.; Ma, L.; Feng, B.; Liu, Q.; Yuan, L.; Fan, C.; Huang, H.;
442
Huang, H.; Yang, Q. Biodiversity of the symbiotic bacteria associated with toxic
443
marine dinoflagellate Alexandrium tamarense. J. Biosci. Med. 2015, 3, 23-28.
444
32. Ren, G.; Ren, W.; Teng, Y.; Li, Z. Evident bacterial community changes but only
445
slight degradation when polluted with pyrene in a red soil. Front. Microbiol. 2015, 6,
446
22.
447
33. Caporaso, J. G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F. D.;
448
Costello, E. K.; Fierer, N.; Peña, A. G.; Goodrich, J. K.; Gordon, J. I.; Huttley, G. A.;
449
Kelley, S. T.; Knights, D.; Koenig, J. E.; Ley, R. E.; Lozupone, C. A.; McDonald, D.;
450
Muegge, B. D.; Pirrung, M.; Reeder, J.; Sevinsky, J. R.; Tumbaugh, P. J.; Walters, W.
451
A.; Widmann, J.; Yatsunenko, T.; Zaneveld, J.; Knight, R. QIIME allows analysis of
452
high-throughput community sequencing data. Nat. Methods 2010, 7, 335-336.
453
34. Edgar, R. C.; Haas, B. J.; Clemente, J. C.; Quince, C.; Knight, R. UCHIME
454
improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 219421 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 22 of 33
455
2200.
456
35. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST.
457
Bioinformatics 2010, 26, 2460-2461.
458
36. Zott, K.; Claisse, O.; Lucas, P.; Coulon, J.; Lonvaud-Funel, A.; Masneuf-
459
Pomarede, I. Characterization of the yeast ecosystem in grape must and wine using
460
real-time PCR. Food Microbiol. 2010, 27, 559-567.
461
37. De Pasquale, I.; Di Cagno, R.; Buchin, S.; De Angelis, M.; Gobbetti, M.
462
Microbial ecology dynamics reveal a succession in the core microbiota involved in
463
the ripening of pasta filata Caciocavallo Pugliese cheese. Appl. Environ. Microbiol.
464
2014, 80, 6243-6255.
465
38. Wang, X.; Du, H.; Xu, Y. Source tracking of prokaryotic communities in
466
fermented grain of Chinese strong-flavor liquor. Int. J. Food Microbiol. 2017, 244,
467
27-35.
468
39. Gamero, A.; Quintilla, R.; Groenewald, M.; Alkema, W.; Boekhout, T.;
469
Hazelwood, L. High-throughput screening of a large collection of non-conventional
470
yeasts reveals their potential for aroma formation in food fermentation. Food
471
Microbiol. 2016, 60, 147-159.
472
40. Liu, S.; Laaksonen, O.; Kortesniemi, M.; Kalpio, M.; Yang, B. Chemical
473
composition of bilberry wine fermented with non-Saccharomyces yeasts (Torulaspora
474
delbrueckii and Schizosaccharomyces pombe) and Saccharomyces cerevisiae in pure,
475
sequential and mixed fermentations. Food Chem. 2018, 266, 262-274.
476
41. Qu, J.; Meng, X.; You, H.; Ye, X.; Du, Z. Utilization of rice husks functionalized 22 ACS Paragon Plus Environment
Page 23 of 33
Journal of Agricultural and Food Chemistry
477
with xanthates as cost-effective biosorbents for optimal Cd(II) removal from aqueous
478
solution via response surface methodology. Bioresour. Technol. 2017, 241, 1036-
479
1042.
480
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 24 of 33
481
FIGURE CAPTIONS
482
Figure 1. Dynamics of the two main higher alcohol contents during the liquor
483
fermentation process: (A) 2-methylpropanol and (B) 3-methylbutanol.
484
Figure 2. Dynamics of the microbial composition during the liquor fermentation
485
process. (A) Relative abundances of fungal and bacterial genera during Chinese light-
486
aroma-type liquor fermentation. (B) Correlation network between microbes and
487
higher alcohols calculated statistically (P < 0.05) by the Pearson correlation
488
coefficient (|ρ| > 0.6). Green modules represent fungal genera, and orange modules
489
represent bacterial genera. The thickness of the edges shows the value of the Pearson
490
correlation coefficient. Red and blue edges indicate negative and positive correlations
491
between the microbes and higher alcohols.
492
Figure 3. Redundancy analysis of microbes and the physicochemical factors during
493
the liquor fermentation process: stage I (A) and stage II (B).
494
Figure 4. Biomass of Pichia kudriavzevii (A) and Saccharomyces cerevisiae (B), and
495
the total production of higher alcohols in mono- and co-cultures in sorghum media
496
(C). Relationship between the specific rate of formation higher alcohols and the
497
growth rate of microbes (D).
498
Figure 5. (A) Predicted and actual values of the total higher alcohols. (B) Total higher
499
alcohol contents in the 17 runs carried out in simulated fermentations.
500
Figure 6. Response surface (3D) and corresponding contour (2D) plots for the total
501
higher alcohols. (A and B) Effects of the Pichia/Saccharomyces ratio and pH on the
502
total higher alcohol contents. (C and D) Effects of the Pichia/Saccharomyces ratio
503
and C/N ratio on the total higher alcohol contents. (E and F) Effects of the pH and
504
C/N ratio on the total higher alcohol contents.
505 24 ACS Paragon Plus Environment
Page 25 of 33
Journal of Agricultural and Food Chemistry
Table 1. Experimental design matrix for the total higher alcohols. Variables lg Pichia (CFU/mL) a pH C/N a The
Symbol coded A B C
Range and levels -1 0 4 5 3.5 4.5 50:1 237.5:1
1 6 5.5 425:1
initial concentration of Saccharomyces was 1 × 105 CFU/mL. Thus, the
Pichia/Saccharomyces ratio ranged from 1:10 to 10:1.
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 33
Table 2. ANOVA for the response surface quadratic model. Sum of df a Squares 25561.77 9
Mean Squared 2840.20
5255.56
1
705.14 12526.23 1138.78 257.21 853.71 1338.05 26.41 3151.35 1455.18
Lack of Fit
1052.25
Pure Error Cor Total
402.92 4 27016.95 16
Source Model A-lg Pichia B-pH C-C/N AB AC BC A2 B2 C2 Residual
a df,
13.66
P-value Prob > F 0.0012
5255.56
25.28
0.0015
1 1 1 1 1 1 1 1 7
705.14 12526.23 1138.78 257.21 853.71 1338.05 26.41 3151.35 207.88
3.39 60.26 5.48 1.24 4.11 6.44 0.13 15.16
0.1081 0.0001 0.0518 0.3027 0.0823 0.0388 0.7320 0.0059
3
350.75
3.48
0.1297
F-value
Significant
Not significant
100.73
degree of freedom.
26 ACS Paragon Plus Environment
Page 27 of 33
Journal of Agricultural and Food Chemistry
Figure 1
27 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 28 of 33
Figure 2
28 ACS Paragon Plus Environment
Page 29 of 33
Journal of Agricultural and Food Chemistry
Figure 3
29 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 30 of 33
Figure 4
30 ACS Paragon Plus Environment
Page 31 of 33
Journal of Agricultural and Food Chemistry
Figure 5
31 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 32 of 33
Figure 6
32 ACS Paragon Plus Environment
Page 33 of 33
Journal of Agricultural and Food Chemistry
TOC Graphic
33 ACS Paragon Plus Environment