Subscriber access provided by NEW YORK UNIV
Article
Metabolic Effects of the pksCT Gene on Monascus aurantiacus Li As3.4384 Using Gas ChromatographyTime-of-Flight Mass Spectrometry-Based Metabolomics Zhibing Huang, Shuyun Zhang, Yang Xu, Laisheng Li, and Yanping Li J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.5b06082 • Publication Date (Web): 29 Jan 2016 Downloaded from http://pubs.acs.org on January 31, 2016
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 free 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 accessible to all readers and 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.
Journal of Agricultural and Food Chemistry 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
Metabolic Effects of the pksCT Gene on Monascus aurantiacus Li As3.4384
1 2
Using Gas Chromatography-Time-of-Flight Mass Spectrometry-Based
3
Metabolomics
4
Zhibing Huang,*,† Shuyun Zhang,† Yang Xu,† Laisheng Li,§ and Yanping Li†
5 6 7
†
State Key Laboratory of Food Science and Technology, Sino–German Joint
8
Research Institute, Nanchang University, No. 235 Nanjing East Road, Nanchang
9
330047, China
10 11
§
Center of Analysis and Testing, Nanchang University, No. 235 Nanjing East Road,
Nanchang 330047, China
12 13
Running head: Metabolic Effects of the pksCT Gene in M. aurantiacus Li As3.4384
14
*
To whom correspondence should be addressed. (Telephone: +86-0791-88305177-8109; Fax: +86-0791-88333708;
E-mail:
[email protected]).
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
15
ABSTRACT
16
Monascus spp. have been used for the production of natural pigments and
17
bioactive compounds in China for several centuries. Monascus also can produce
18
the mycotoxin citrinin, restricting its use. Disruption of the pksCT gene in Monascus
19
aurantiacus Li AS3.4384 reduces citrinin production capacity of this strain
20
(Monascus PHDS26) by over 98%. However, it is unclear how other metabolites of
21
M. aurantiacus Li AS3.4384 (the wild-type strain) are affected by the pksCT gene.
22
Here, we used metabolomics analyses to compare red yeast rice (RYR) metabolite
23
profiles of the wild-type strain and Monascus PHDS26 at different stages of
24
solid-state fermentation. Eighteen metabolites forming components within the
25
glycolysis, acetyl-CoA, amino acid, and TCA cycle metabolic processes were found
26
to be altered between the wild-type strain and Monascus PHDS26 at different
27
stages of solid-state fermentation. Thus, these findings provide important insights
28
into the metabolic pathways affected by the pksCT gene in M. aurantiacus.
29 30
KEYWORDS: Monascus aurantiacus Li AS3.4384, pksCT gene disrupted strain,
31
metabolomics, GC-TOF-MS
32
ACS Paragon Plus Environment
Page 2 of 38
Page 3 of 38
Journal of Agricultural and Food Chemistry
33
INTRODUCTION
34
Monascus and its metabolites have been used in food and medicines for thousands
35
of years in China.1 Monascus can produce many bioactive compounds, including
36
Monascus pigments, monacolin K, γ-aminobutyric acid (GABA), and ergosterol.2,3
37
Monacolin K is usually used to reduce serum triglyceride and cholesterol levels.4
38
Some
39
immunomodulatory, and anti-inflammatory effects. Moreover, some components
40
increase levels of high-density lipoprotein cholesterol and have anti-atherosclerosis
41
and anticancer effects.5-8 However, Monascus also can produce the mycotoxin
42
citrinin, which is found after fermentation of Monascus.9 Citrinin has been shown to
43
have nephrotoxic, carcinogenic, and hepatotoxic effects.10 Notably, the citrinin
44
production capacity of Monascus is decreased following disruption of the pksCT
45
gene in M. purpureus; in this strain, which was constructed by homologous
46
recombination, citrinin is not produced.11 Additionally, in our previous study,
47
disruption of the pksCT gene in strain Monascus PHDS26 decreased citrinin
48
production capacity by over 98%.12,13 However, it is unclear how other metabolites
49
of Monascus are affected by the pksCT gene.
Monascus
pigments
have
been
shown
to
have
anti-obesity,
50
Metabolomics can be used to compare different metabolites and changes in
51
different biological samples with changes in their environment, including analysis of
52
global small molecules in biofluids.14 Metabolite analysis is usually carried out using
53
multiple experimental platforms, such as gas chromatography-mass spectrometry
54
(GC-MS),
liquid
chromatography-mass
spectrometry
ACS Paragon Plus Environment
(LC-MS),
capillary
Journal of Agricultural and Food Chemistry
55
electrophoresis-mass spectrometry (CE-MS), and nuclear magnetic resonance
56
(NMR).15-17 Unlike in genomics and proteomics studies, metabolomics cannot be
57
used to completely characterize all metabolites present in biological samples using
58
only a single analytical platform.18 Nevertheless, GC-MS-based metabolomics
59
approaches aim to target the metabolism of some intermediates, such as fatty acids,
60
nucleotides, carbohydrates, hydroxyl acids, amino acids, and related compounds.
61
GC-MS-based metabolomics is considered the most mature analytical platform;
62
more than 250,000 compounds and small molecules are included in standard
63
repositories of mass spectra, and compounds can be searched using the NIST and
64
Wiley libraries. In metabolomics approaches, large amounts of analytical data have
65
been obtained using these analytical platforms, and the data can be analyzed and
66
processed with multivariate analytical methods. Multivariate analyses usually
67
include principal component analysis (PCA), partial least-squares discriminate
68
analysis (PLS-DA), and orthogonal partial least-squares discriminate analysis
69
(OPLS-DA). These data analysis methods can cluster of metabolite profiles or be
70
used to discover important changes between different groups of samples. The
71
methods have been used to analyze metabolite profiles of traditional meju during
72
fermentation,19 fast-fermented bean paste,20 and transgenic Beauveria bassiana.21
73
To the best of our knowledge, no metabolomic studies have been conducted in
74
Monascus.
75
Therefore, the goal of this work was to determine differences in the metabolic
76
profiles of red yeast rice (RYR) from the wild-type strain and Monascus PHDS26
ACS Paragon Plus Environment
Page 4 of 38
Page 5 of 38
Journal of Agricultural and Food Chemistry
77
after cultivation for different fermentation times (6, 10, 14, 18, and 22 days),
78
followed by GC-TOF-MS analysis and subsequent multivariate statistical analysis.
79
Our results provide important insights into the metabolic profiles of Monascus spp.,
80
particularly with regard to the role of the pksCT gene.
81 82
MATERIALS AND METHODS
83
Chemicals and materials
84
Rice was purchased from local supermarkets. M. aurantiacus Li As3.4384 was
85
obtained from the Institute of Microbiology, Chinese Academy of Science (Beijing,
86
China). The Monascus PHDS26 strain was constructed by disrupting the pksCT gene
87
in the wild-type strain.12,13 Ribitol, methoxyamine hydrochloride, and pyridine were
88
purchased from Sigma-Aldrich (St. Louis, MO, USA). N,O-bis(trimethylsilyl)
89
trifluoroacetamide (BSTFA) and trimethylchlorosilane (TMCS) were obtained from
90
Regis Technologies, Inc. (Morton Grove, IL, USA). High-performance liquid
91
chromatography (HPLC)-grade methanol was purchased from Merck (Darmstadt,
92
Germany). All other chemicals (analytical grade) were obtained from Sinopharm
93
Chemical Reagent Co. Ltd. (Beijing, China). Distilled water was purified by the Milli-Q
94
Plus System (18.2 MΩ cm; Millipore Inc., Billerica, MA, USA).
95
Preparation and extraction of RYR
96
RYR was obtained according to our previously reported methods.22 Rice was
97
soaked in water for 12 h. After draining, the rice (100 g) was placed in an Erlenmeyer
98
flask (250 mL) and then sterilized using standard procedures at 121°C for 15 min.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
99
Each Erlenmeyer flask was inoculated with the seed cultures (2 mL) in an incubator at
100
30°C and fermented for 6, 10, 14, 18, or 22 days. After different periods of incubation,
101
the RYR was collected, immediately frozen in liquid nitrogen, and then stored at
102
-80°C for further freeze-drying. Each RYR (100 mg) sample was pulverized in a 2 mL
103
centrifuge tube containing a steel metal ball in a mixer-mill grinder at 70 Hz for 5 min,
104
and 60 µL of 0.2 mg/mL of ribitol (as an internal standard) was then added. The
105
samples were mixed by vortexing for 1 min, and 0.4 mL cold (-20°C)
106
methanol-chloroform (1:1, v/v) was then added, followed by mixing again with
107
vortexing for 1 min. Each sample was extracted by sonication for 30 min and then
108
shaken for 3 h at 4°C. After extraction, the sample was centrifuged at 10000 × g for
109
10 min at 4°C (Sigma-Laborzentrifugen, 2K15; Osterode, Germany), and 0.35 mL of
110
the supernatant was then collected and dried in a vacuum concentrator at 45°C for
111
1.5 h. Finally, the dry extracts were derivatized by methoxyamination using 80 µL of
112
methoxyamine-HCl (20 mg/mL in pyridine) at 37°C for 3 h. The samples were then
113
trimethylsilylated for 1 h using 100 µL of BSTFA with 1% TMCS (v/v) at 70°C with
114
shaking. After silylation, 50 µL of n-heptane was added.
115
GC-TOF-MS analysis
116
GC-TOF-MS analysis was carried out using an Agilent 7890A gas chromatograph
117
system (Agilent, Atlanta, GA, USA) coupled with a Pegasus 4D TOF mass
118
spectrometer (LECO Corp., St. Joseph, MI, USA). Each derivatized sample (1 µL)
119
was separated on a 30 m × 250 µm i.d. DB-5MS capillary column (0.25 µm film
120
thickness) coated with 5% diphenyl crosslinked with 95% dimethylpolysiloxane
ACS Paragon Plus Environment
Page 6 of 38
Page 7 of 38
Journal of Agricultural and Food Chemistry
121
(Agilent J&W Scientific, Folsom, CA, USA). The helium (purity 99.9999%) gas flow
122
rates of the front inlet purge and through the column were 3 and 1 mL/min,
123
respectively. The temperatures set were as follows: initial temperature of 90°C for
124
0.25 min, increased to 180°C at a rate of 10°C/min, then increased to 240°C at a rate
125
of 5°C/min, and finally increased to 285°C at a rate of 20°C/min, after which the
126
samples were held at 285°C for 11.5 min. The injection, ion source, and transfer line
127
temperatures were 280, 220, and 245°C, respectively. Ions were generated by 70 eV
128
electro-impact (EI) ionization at a current of 3.0 mA, and EI spectra were acquired in
129
the m/z range of 20–600 at a rate of 100 spectra/s after a solvent delay of 7 min. Each
130
sample was analyzed with six replications.
131
Data analysis
132
The GC-TOF-MS raw data were exported to the netCDF file format using LECO
133
ChromaTOF version 4.44 software (LECO Corp., St. Joseph, MI, USA) for mass
134
spectral deconvolution and peak finding, followed by retention index matching,
135
normalization, peak merging, and statistical significance tests. The signal-to-noise
136
(S/N) ratios of peaks lower than 30 were removed. All compounds were identified
137
by comparing their mass spectra and standards with the LECO/Fiehn Metabolomics
138
Library, yielding similarity values for compound identification accuracy. If the
139
similarity was more than 700, the peak was assigned as a metabolite, and the
140
metabolite identification was considered reliable. If the similarity was less than 200,
141
the library would only use “analyte” for the compound name. If the similarity was
142
between 200 and 700, the compound name was given a putative annotation.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
143
Obvious differences in the metabolites were assigned with a variable importance in
144
the projection (VIP) value of greater than 1.0. Some derivatives was assigned by
145
increasing numbers according to the retention index, e.g., serine 1, serine 2, and
146
serine 3.23
147
For multivariate analysis, the relative data matrix containing peak area
148
information as variables and sample names as observations was processed using
149
SIMCA-P version 12.0 software (Umetrics, Umeå, Sweden). PCA, PLS-DA, and
150
OPLS modelling were used to compare different metabolites from the wild-type
151
strain and Monascus PHDS26 during different fermentation stages. The
152
parameters of the models were defined as described previously.24
153 154
RESULTS AND DISCUSSION
155
GC-TOF-MS-based metabolite profiling of Monascus
156
The color value of the pigments and citrinin in the RYR obtained from the wild-type
157
strain and Monascus PHDS26 were measured by ultraviolet-visible (UV-Vis)
158
spectrophotometry and HPLC, respectively, according to previous methods.25,26 As
159
shown in Figure 1(A), no significant differences of the contents of pigments in the
160
RYR were found in the wild-type strain after different fermentation periods from 6 to
161
16 days; however, the contents of pigments increased obviously from 16 to 22 days.
162
Furthermore, as shown in Figure 1(B), the content of citrinin increased significantly
163
for M. aurantiacus Li As3.4384 during the fermentation periods from 6 to 14 days. In
164
addition, the contents of two orange pigments from M. aurantiacus Li As3.4384 also
ACS Paragon Plus Environment
Page 8 of 38
Page 9 of 38
Journal of Agricultural and Food Chemistry
165
increased from 6 to 16 days.27 Therefore, the fermentation period from 6 to 22 days
166
was selected for analysis of the metabolic effects of pksCT gene in M. aurantiacus
167
Li As3.4384.
168
Next, the metabolic profiles of the wild-type strain and Monascus PHDS26
169
during solid-state fermentation (6, 10, 14, 18, and 22 days) were investigated by
170
GC-TOF-MS. Ribitol was used as an internal standard to evaluate the performance
171
of the GC-TOF-MS system. The retention time of the ribitol shift ranged from 0.01 to
172
0.03 min, and the RSD of the peak area ranged from 2.0% to 5.1%. The results
173
showed that the GC-TOF-MS system provided stable sample analysis. The extracts
174
of RYR from the solid-state fermentation of wild-type strain and Monascus PHDS26
175
were analyzed. The GC-TOF-MS total ion current chromatographs of RYR extracts
176
obtained from fermentation after 22 days are shown in Figure 2. Some metabolites
177
exhibited obvious increases in relative abundance, whereas others showed
178
decreasing relative abundances, and some remained unchanged after different
179
fermentation times (Figure 2). In total, 696 metabolites, including organic acids,
180
saccharides, amino acids, and amines, were detected in metabolomic profiling of
181
the RYR. Moreover, 85 metabolites of the RYR from the solid-state fermentation of
182
the wild-type strain and Monascus PHDS26 were integrated by NIST mass spectral
183
similarity scores of no less than 700 (Table 1). These known metabolites included
184
amino acids and amines (e.g., alanine, glycine, valine, proline, serine, aspartic acid,
185
oxoproline, glutamic acid, phenylalanine, glutamine, and tyrosine), organic acids
186
(e.g., oxamic acid methylmalonic acid, succinic acid, fumaric acid, L-malic acid, and
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
187
citric acid), and sugars and sugar alcohols (e.g., fructose, mannose, sorbitol,
188
sucrose, trehalose, isomaltose, and galactinol; Table 1). Based on the P-values (
1.0), a total of 91
190
metabolites with important roles in different groups were selected. Of these, 10
191
metabolites were identified based on similarity scores of no less than 700, 12
192
metabolites were putatively annotated according to the similarities between 200
193
and 700, and 69 exhibited unknown metabolic profiles in the RYR after fermentation
194
for 6 days by the wild-type strain and Monascus PHDS26 (Supplemental Table S1).
195
Moreover, 124 metabolites differed significantly after fermentation for 10 days (12
196
known metabolites, 17 putatively annotated metabolites, and 95 unknown
197
metabolites; Table S2). After 14 days of fermentation, 127 metabolites differed
198
significantly (seven known metabolites, 16 putatively annotated metabolites, and
199
104 unknown metabolites; Table S3). After fermentation for 18 days, 134
200
metabolites differed significantly, including seven known metabolites, 25 putatively
201
annotated metabolites, and 101 unknown metabolites (Table S4). Finally, after
202
fermentation for 22 days, 17 metabolites differed, including eight known metabolites,
203
three putatively annotated metabolites, and six unknown metabolites (Table S5). As
204
shown in Table 2, some metabolites were identified according to similarity scores of
205
no less than 700.
206
In order to study the metabolic variations between different fermentation
207
periods (from 6 to 22 days), the metabolic profiles of the wild-type strain and
208
Monascus PHDS26 after different fermentation times, e.g., 6, 10, 14, 18, and 22
ACS Paragon Plus Environment
Page 10 of 38
Page 11 of 38
Journal of Agricultural and Food Chemistry
209
days, were compared. Figure 3 shows a region of the TIC detection of RYR
210
extracts comparing the results obtained from different fermentation periods (6 to 22
211
days). Several metabolites were measured in this region of the TIC, including
212
glutamine, trehalose, diglycerol, aminomalonic acid, fructose, L-malic acid, sorbitol,
213
threitol, citric acid, and mannose. As shown in Figure 3, the mean levels of some
214
known metabolites, such as maleimide, hydroxylamine, threitol, L-malic acid, and
215
mannose increased or decreased and then increased from 6 to 14 days, reaching
216
maximum levels in the fermentation of Monascus PHDS26 after 14 days. The
217
maximum levels of alanine, diglycerol, citric acid, and fructose were reached after
218
the fermentation of Monascus PHDS26 for 18 days. Moreover, the maximum levels
219
of
220
2,3-dihydroxybenz acid, glutamine, trehalose, and sorbitol were reached in the
221
fermentation of Monascus PHDS26 for 22 days. Additionally, the level of
222
2-hydroxybutanoic acid increased obviously from 6 to 22 days for M. aurantiacus Li
223
As3.4384 and remained unchanged from 6 to 22 days for Monascus PHDS26. The
224
contents of some metabolites, such as 2-hydroxybutanoic acid, phenylalanine,
225
2,3-dihydroxybenz acid, myristic acid, fructose, and sorbitol in the RYR produced by
226
Monascus PHDS26 were much lower than those in the RYR produced by the
227
wild-type strain. However, the levels of maleimide and trehalose were much higher
228
in the RYR produced by Monascus PHDS26, and the contents of threitol, and
229
L-malic
230
produced by the wild-type strain during the fermentation period from 10 to 22 days.
leucine,
aminomalonic
acid,
oxoproline,
phenylalanine,
asparagine,
acid were also much higher in Monascus PHDS26 than in the RYR
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
231
The contents of hydroxylamine, citric acid, and mannose in the RYR produced by
232
Monascus PHDS26 were also much higher than those in the RYR produced by the
233
wild-type strain during the fermentation period from 6 to 14 days. The levels of
234
L-malic
235
higher than those in the wild-type strain under normal fermentation conditions.
236
Analysis of data by PCA, PLS-DA, and OPLS-DA
acid, hydroxylamine, citric acid, and trehalose in Monascus PHDS26 were
237
To understand the overall differences in the metabolism of the RYR from the
238
wild-type strain and Monascus PHDS26, the acquired metabolomic data different
239
fermentation stages were analyzed using PCA, which was performed to discover
240
principal components for the evaluation the majority of differences in the data on
241
the basis of significantly altered metabolites. PCA has been proven to be useful for
242
discrimination and grouping in the analysis of chemical profiles of medicines and
243
food.28, 29 The total of PCA scores plot is shown in Figure 4(A), and the values of
244
R2X and Q2 were 0.609 and 0.477, respectively. The changes in metabolic profiles
245
in the RYR produced by M. aurantiacus Li As3.4384 and Monascus PHDS26 from 6
246
to 22 days are also shown in the 3D PCA scores plot graph (Figure 4(B)). As
247
shown in Figure 4, the metabolic profiles of the wild-type strain and Monascus
248
PHDS26 were changed dramatically during fermentation from 6 to 18 days, with
249
slight changes from 18 to 22 days. The PCA scores graph (R2X and Q2 were 0.79
250
and 0.468, respectively) of changes in the metabolic profiles of the wild-type strain
251
and Monascus PHDS26 are shown in Figure 4(C).
252
Next, we used PCA to compare the various metabolic components produced by
ACS Paragon Plus Environment
Page 12 of 38
Page 13 of 38
Journal of Agricultural and Food Chemistry
253
the wild-type strain and Monascus PHDS26 during the same fermentation period.
254
As shown in Figure 5(A) and Figure S1–4(A), the values of R2X and Q2 were 0.624
255
and 0.453, 0.675 and 0.346, 0.826 and 0.548, 0.721 and 0.443, and 0.576 and
256
0.379, respectively, for differences at 6, 10, 14, 18, and 22 days of fermentation.
257
PCA score plots showed that the different RYR samples were distributed in two
258
separate areas, indicating significantly different metabolomes between the
259
wild-type and mutant Monascus strains at different fermentation stages.
260
The PLS-DA mold is widely performed for clustering and comparing nontargeted
261
metabolic profiling data of various sample groups.20,30 Thus, to identify correlations
262
of data (X) and other variables (Y, packet information), the normalized data were
263
analyzed by the partial least squares (PLS-DA) method using SIMCA-P software.
264
As shown in Figure 5(B and C) and Figure S1–4(B and C), the PLS-DA score plot
265
for the whole data set was used to categorize RYR samples from the wild-type
266
strain and Monascus PHDS26 into different fermentation stages (6, 10, 14, 18, and
267
22 days). The R2X, R2Y, and Q2 values of the PLS-DA score plot of the wild-type
268
strain and Monascus PHDS26 at the same fermentation stage were 0.623, 0.998,
269
and 0.987 for 6 days; 0.547, 0.998, and 0.981 for 10 days; 0.57, 0.999, and 0.986
270
for 14 days; 0.563, 0.999, and 0.987 for 18 days, and 0.62, 0.983, and 0.956 for 22
271
days, respectively. These results indicated that the PLS-DA models exhibited
272
differences between the different RYR samples from the wild-type strain and
273
Monascus PHDS26 in different fermentation periods.
274
To further compare the differences in metabolites closely associated with the
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
275
pksCT gene between the wild-type strain and Monascus PHDS26 at different
276
fermentation periods (6 to 22 days), we used the OPLS-DA model, a supervised
277
pattern recognition method that can facilitate the identification of some marker
278
metabolites responsible for maximum division by removing systematic variation
279
unrelated to grouping.21 As shown in Figure 5(D) and Figure S1–4(D), the
280
OPLS-DA model showed a clear separation between the RYR from the wild-type
281
strain and Monascus PHDS26 at the same fermentation time points (6 to 22 days)
282
according to the first component (OPLS1) in the score plots. The R2X, R2Y, and Q2
283
of the OPLS-DA score plots of the wild-type strain and Monascus PHDS26 at the
284
same fermentation time points were 0.623, 0.998, and 0.972 for 6 days; 0.547,
285
0.998, and 0.968 for 10 days; 0.57, 0.999, and 0.978 for 14 days; 0.563, 0.999, and
286
0.985 for 18 days; and 0.62, 0.983, and 0.957 for 22 days, respectively. These
287
results suggested that the OPLS-DA model was a reliable predictive method,
288
yielding R2Y and Q2 values close to 1.0, and the RYR samples of the wild-type
289
strain and Monascus PHDS26 were clearly separated into two groups.
290
Pathway analysis
291
To unravel the changes in metabolites in the metabolic network, we evaluated the
292
metabolic pathways of several identified metabolites in RYR extracts (Figure 6).
293
Some of the differential metabolites formed components in the glycolysis, acetyl-CoA,
294
amino acid, and TCA cycle metabolic pathways.
295
For glycolysis, the identified compounds included trehalose, myo-inositol,
296
phenylalanine, glucose, alanine, sorbitol, fructose, mannose, and leucine. A large
ACS Paragon Plus Environment
Page 14 of 38
Page 15 of 38
Journal of Agricultural and Food Chemistry
297
amount of polyhydroxyl compounds, including cyclitols, sugars, and acyclic polyols,
298
are used by fungi to overcome stress conditions, such as cold, desiccation, heat, and
299
others.31 Trehalose is a nonreducing disaccharide and plays roles in protein and cell
300
membrane fluidity and stability and in maintenance of the lipid bilayer under stress
301
conditions.21,32 Trehalose has been shown to protect yeast pyrophosphatase under
302
exposure at 50°C much better than other sugars, including glucose, fructose, and
303
sucrose.31,33 In this study, the RYR was obtained by Monascus solid-fermentation in
304
the dry environment; the spores of these fungi are adapted to withstand desiccation.
305
As shown in Figure 3, the levels of trehalose produced by the wild-type strain and
306
pksCT-disrupted strain increased from 6 to 22 days. This could be explained by the
307
observation that trehalose was synthesized during the drying and ripening of the
308
spore and played a role in protection of membranes.
309
The metabolic flux moved to acetyl-CoA after glycolysis. Acetyl-CoA, a
310
biomass-derivable building block, can be indirectly transformed into malic acid and
311
asparagine and may be involved in fatty acid biosynthesis, followed by conversion
312
to myristic acid. Malic acid, another building block involved in metabolite synthesis,
313
can be produced by some fungal species, including M. araneosus, Schizophyllum
314
commune, Aspergillus flavus, and Zygosaccharomyces rouxii.34 The malic acid
315
biosynthesis pathway includes the following three routes: oxidation of citrate by the
316
TCA cycle; direct reduction of oxaloacetate; and production from acetyl-CoA via
317
glyoxylate metabolism.35 However, malic acid is primarily produced from
318
carboxylation of pyruvate during glucose metabolism, followed by reduction of
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 38
319
oxaloacetate to malate.34 As shown in Figure 6, malic acid was produced from
320
acetyl-CoA via glyoxylate metabolism in the wild-type strain and Monascus
321
PHDS26. Acetyl-CoA is also essential for polyketide-dependent secondary
322
metabolite biosynthesis of Monascus pigments and citrinin and is formed through
323
tetraketide as a precursor intermediate, yielding polyketides.1 During cultivation of
324
the pksCT-disrupted M. aurantiacus strain Li AS3.4384 under rice, the production of
325
polyketide pigment was increased; however, citrinin production was inhibited. This
326
could be explained by the shifting of the metabolic flux of acetyl-CoA to the TCA
327
cycle, followed by switching to the biosynthetic polyketide pathway for the formation
328
of
329
2-hydroxybutanoic acid is formed from the TCA cycle via propanoyl-CoA. Several
330
compounds, including hydroxylamine, oxoproline, glutamate, and glutamine, are
331
produced from the TCA cycle via isocitrate and 2-oxoglutarate. Additionally,
332
pyroglutamic acid (oxoproline) is a γ-aminobutyric acid intermediate produced
333
during the γ-glutamyl cycle.
pigments. Citric
acid is
an intermediate in the TCA
cycle,35 and
334
In summary, the metabolic effects of the pksCT gene in M. aurantiacus Li
335
As3.4384 were monitored by GC-TOF-MS analysis. From this analysis, 85
336
compounds were identified, and 18 metabolites were found to be altered between
337
the wild-type strain and Monascus PHDS26 at different stages of solid-state
338
fermentation.
339
2,3-dihydroxybenz acid, myristic acid, fructose, allo-inositol, myo-insitol, and
340
sorbitol in the wild-type strain were higher than those in Monascus PHDS26 under
The
levels
of
2-hydroxybutanoic
ACS Paragon Plus Environment
acid,
phenylalanine,
Page 17 of 38
Journal of Agricultural and Food Chemistry
341
normal fermentation conditions (Figure 3).These metabolites forming components
342
within the glycolysis, acetyl-CoA, amino acid, and TCA cycle metabolic processes
343
were found to be altered between the wild-type strain and Monascus PHDS26 at
344
different stages of solid-state fermentation.
345
ACKNOWLEDGMENT
346
This project was financially supported by the National Natural Science Foundation of
347
China (No. 31160308 and No. 21165012). This work was also supported by the
348
Project of the Free-exploration Research Project of State Key Laboratory of Food
349
Science and Technology (SKLF-ZZB-201518).
350
Notes
351
The authors declare no competing financial interest.
352
REFERENCES
353
(1) Lin, W. Y.; Chang, J. Y.; Hish, C. H.; Pan, T. M. Profiling the Monascus pilosus
354
proteome during nitrogen limitation. J. Agric. Food Chem. 2008, 56, 433–441
355
(2) Gao, J.M.; Yang, S.X.; Qin, J.C. Azaphilones: Chemistry and Biology.
356
Chem.Rev. 2013, 113(7), 4755–4811.
357
(3) Cheng, M. J.; Wu, M. D.; Chen, I. S.; Yuan, G. F. Secondary metabolites from
358
the mycelia of the fungus Monascus pilosus BCRC 38072. Chem. Pharma. Bull.
359
2008, 56, 394–397.
360
(4) Tobert, J.A.. Lovastatin and beyond: The history of the HMG-CoA reductase
361
inhibitors. Nat. Rev. Drug Discov. 2003, 2(7), 517–526.
362
(5) Chang, Y. Y.; Hsu, W. H.; Pan, T. M. Monascus secondary metabolites
363
monascin and ankaflavin inhibit activation of RBL-2H3 cells. J. Agric. Food Chem.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
364
2015, 63 (1), 192–199.
365
(6) Lee, C. L.; Hung, Y. P.; Hsu, Y. W.; Pan, T.M. Monascin and ankaflavin
366
have more anti-atherosclerosis effect and less side effect involving increasing
367
creatinine phosphokinase activity than monacolin K under the same dosages.
368
J. Agric. Food Chem. 2013, 61 (1), 143–150.
369
(7) Patakova, P. Monascus secondary metabolites: Production and biological
370
activity. J. Ind. Microbiol. Biotechnol. 2013, 40,169–181.
371
(8) Lee, C. L.; Wen, J. Y.; Hsu, Y. W.; Pan, T. M. Monascus-fermented yellow
372
pigments monascin and ankaflavin showed antiobesity effect via the suppression of
373
differentiation and lipogenesis in obese rats fed a high-fat diet. J. Agric. Food
374
Chem. 2013, 61 (7), 1493–1500.
375
(9) Atar, N.; Eren, T.; Yola, M. L. A molecular imprinted SPR biosensor for sensitive
376
determination of citrinin in red yeast rice. Food Chem. 2015, 184, 7-11.
377
(10) Liu, B. H.; Wu, T. S.; Su, M. C.; Yu, F. Y. Evaluation of citrinin occurrence and
378
cytotoxicity in Monascus fermentation products. J. Agric. Food Chem. 2005, 53
379
(1),170–175.
380
(11) Shimizu, T.; Kinoshita, H.; Ishihara, S.; Sakai, K.; Nagai, S.; Nihira, T.
381
Polyketide synthase gene responsible for citrinin biosynthesis in Monascus
382
purpureus. Appl. Environ. Microbio. 2005, 71(7), 3453–3457.
383
(12) Li, Y. P.; Pan, Y. F.; Zou, L. H.; Xu, Y.; Huang, Z. B.; He, Q. H. Lower citrinin
384
production by gene disruption of ctnB involved in citrinin biosynthesis in Monascus
385
aurantiacus Li AS3.4384. J. Agric. Food Chem. 2013, 61(30), 7397−7402.
ACS Paragon Plus Environment
Page 18 of 38
Page 19 of 38
Journal of Agricultural and Food Chemistry
386
(13) Fu, G. M.; Xu, Y.; Li, Y. P.; Tan, W. H. Construction of a replacement vector
387
to disrupt pksCT gene for the mycotoxin citrinin biosynthesis in Monascus
388
aurantiacus and maintain food red pigment production. Asia Pac. J. Clin. Nutr.
389
2007,16 (S1), 137−142.
390
(14) Lai, Y. S.; Chen, W. C.; Kuo, T. C.; Ho, C. T.; Kuo, C. H.; Tseng, Y. J.; Lu, K. H.;
391
Lin, S. H.; Panyod, S.; Sheen, L. Y. Mass-spectrometry-based serum metabolomics
392
of a C57BL/6J mouse model of high-fat-diet-induced non-alcoholic fatty liver
393
disease development. J. Agric. Food Chem. 2015, 63 (35), 7873–7884.
394
(15) Zheng, H.; Yde, C. C.; Clausen, M. R.; Kristensen, M.; Lorenzen, J.; Astrup, A.;
395
Bertram, H. C. Metabolomics investigation to shed light on cheese as a possible
396
piece in the French paradox puzzle. J. Agric. Food Chem. 2015, 63(10), 2830-2839.
397
(16) Freitas, D. S.; Carlos, E. F.; Gil, M. C. S. S.; Vieira, L. G. E.; Alcantara, G. B.
398
NMR-based metabolomic analysis of Huanglongbing-asymptomatic and
399
–symptomatic citrus trees. J. Agric. Food Chem. 2015, 63(34), 7582-7588.
400
(17) Arapitsas, P.; Corte, A. D.; Gika, H.; Narduzzi, L.; Mattivi, F.; Theodoridis, G.
401
Studying the effect of storage conditions on the metabolite content of red wine
402
using HILIC LC-MS based metabolomics. Food Chem. 2016, 197, 1331-1340.
403
(18) Barupal, D. K.; Haldiya, P. K.; Wohlgemuth, G.; Kind, T.; Kothari, S. L.;
404
Pinkerton, K. E.; Fiehn, O. MetaMapp: mapping and visualizing metabolomic data
405
by integrating information from biochemical pathways and chemical and mass
406
spectral similarity. BMC Bioinfor. 2012, 13, 99−113
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 38
407
(19) Lee, S. Y.; Kim, H. Y.; Lee, S.; Lee, J. M.; Muthaiya, M. J.; Kim, B. S.; Oh, J. Y.;
408
Song, C. K.; Jeon, E. S.; Ryu, H. S.; Lee, C. H. Mass spectrometry-based
409
metabolite profiling and bacterial diversity characterization Korean traditional meju
410
during fermentation. J. Microbiol. Biotechnol. 2012, 22, 1523−1531.
411
(20) Kim, J.; Choi, J. N.; John, K. M. M.; Kusano, M.; Oikawa, A.; Saito, K.; Lee, C.
412
H. GC-TOF-MS- and CE-TOF-MS-based metabolic profiling of cheonggukjang
413
(fast-fermented bean paste) during fermentation and its correlation with metabolic
414
pathways. J. Agric. Food Chem. 2012, 60, 9746−9753.
415
(21) Lou, F.; Lu, R.; Zhou, H.; Hu, F.; Bao, G.; Huang, B.; Li, Z. Metabolic effect of an
416
exogenous
417
chromatography-mass spectrometry-based metabolomics. J. Agric. Food Chem.
418
2013, 61, 7008−7017.
419
(22) Huang, Z. B.; Zhang, S. Y., Xu, Y., Li, L.S., Li, Y.P. Structural characterization of
420
two new orange pigments with strong yellow fluorescence. Phytochem. Lett. 2014,
421
10, 140–144.
422
(23) Kind, T.; Wohlgemuth, G.; Lee, D. Y.; Lu, Y.; Palazoglu, M.; Shahbaz, S.; Fiehn,
423
O. FiehnLib: mass spectral and retention index libraries for metabolomics based on
424
quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal. Chem.
425
2009, 81(24),10038−10048.
426
(24) Jin, H.; Qiao, F.; Chen, L.; Lu, C.; Xu, L.; Gao, X. Serum metabolomic
427
signatures of lymph node metastasis of esophageal squamous cell carcinoma. J.
428
Proteome Res. 2014, 13, 4091−4103.
gene
on
transgenic
Beauveria
bassiana
ACS Paragon Plus Environment
using
liquid
Page 21 of 38
Journal of Agricultural and Food Chemistry
429
(25) Wang, W.; Chen, Q.; Zhang, X.; Zhang, H.; Huang, Q.; Li, S.; Yao, J.
430
Comparison of extraction methods for analysis of citrinin in red fermented rice.
431
Food Chem. 2014, 157, 408-412.
432
(26) Yongsmith, B.; Kitprechavanich, V.; Chitradon, L.; Chaisrisook, C.; Budda, N.
433
Color mutants of Monascus sp. KB9 and their comparative glucoamylases on rice
434
solid culture. J. Mole. Catal. B : Enzym. 2000, 10, 263-272.
435
(27) Huang, Z. B.; Zhang, S. Y.; Xu, Y., Li, L. S.; Li, Y. P. Development of an
436
HPLC–UV detection method for the simultaneous determination of two
437
Monascus orange pigments in red yeast rice. Food Anal. Method 2016, 9,
438
148-155.
439
(28)Liang, X. P.; Chen, X.; Liang, Q. L.; Zhang, H. Y.; Hu, P.; Wang, Y. M.; Luo, G. A.
440
Metabonomic study of Chinese medicine Shuanglong formula as an effective
441
treatment for myocardial infarction in rats. J. Proteome Res. 2011, 10, 790−799.
442
(29) Lu, Y.; Lam, H.; Pi, E.; Zhan, Q.; Tsai, S.; Wang, C.; Kwan, Y.; Ngai, S.
443
Comparative metabolomics in Glycine max and Glycine soja under salt stress to
444
reveal the phenotypes of their offspring. J. Agric. Food Chem. 2013, 61, 8711−8721.
445
(30) Kusano, M.; Fukushima, A.; Kobayashi, M.; Hayashi, N.; Jonsson, P.; Moritz, T.;
446
Ebana, K.; Saito, K.
447
one-dimensional and two-dimensional gas chromatography-time-of –flight mass
448
spectrometry to metabolic phenotyping of natural variants in rice. J. Chromatogr. B
449
2007, 855, 71-79.
450
(31) Jain, N. K.; Roy, I. Effect of trehalose on protein structure. Protein Sci. 2009,
Application of a metabolomic method combining
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
451
18(1), 24–36.
452
(32) Elbein, A. D.; Pan, Y.; Pastuszak, I.; Carroll, D. New insights on trehalose: A
453
multifunctional molecule. Glycobiology 2003, 13, 17R−27R.
454
(33) Sola-Penn, M., Meyer-Fernandes, J. R. Stabilization against thermal inactivation
455
promoted by sugars on enzyme structure and function: why is trehalose more
456
effective than other sugars? Arch. Biochem. Biophys. 1998, 360, 10–14.
457
(34) Zell, R. M.; Hulster, E.; Winden, W. A.; Waard, A. A.; Geertman, J. M. A.; Dijken,
458
J. P.; Pronk, J. T.; Maris, A. J. Malic acid production by Saccharomyces cerevisiae:
459
engineering of pyruvate carboxylation, oxaloacetate reduction, and malate export.
460
Appl. Environ. Microbiol. 2008, 74, 2766-2777
461
(35) Clausen, M. R.; Bach, V.; Edelenbos, M.; Bertram, H. C. Metabolomics reveals
462
drastic compositional changes during overwintering of Jerusalem artichoke
463
(Helianthus tuberosus L.) tubers. J. Agric. Food Chem. 2012, 60, 9495−9501.
464
ACS Paragon Plus Environment
Page 22 of 38
Page 23 of 38
Journal of Agricultural and Food Chemistry
465
Figure captions
466
Figure 1. (A) The color value of pigments and (B) citrinin content in the RYR
467
obtained from the wild-type strain and Monascus PHDS26 after different
468
fermentation times. 4384: wild-type strain of Monascus aurantiacus Li AS3.4384. P:
469
Monascus PHDS26 All the points of values are means of all data obtained in
470
triplicate analyses, and the relative standard deviation (RSD) range was from
471
1.57% to 7.41%. Some of points of the standard deviations are very small (within
472
the size of the points on the plots)
473
Figure 2. GC-TOF-MS total ion current chromatographs (n = 6) of RYR extracts
474
obtained from (A) the wild-type strain and (B) Monascus PHDS26 after fermentation
475
for 22 days.
476
Figure 3. Trajectory plots from GC-TOF-MS experiments representing the average
477
normalized relative peak areas of metabolites of RYR from different fermentation
478
stages for the wild-type strain and Monascus PHDS26. (All the points of values are
479
means of all data obtained by at least triplicate analyses, and the relative standard
480
deviation (RSD) range was from 1.0% to 5.5%)
481
Figure 4. (A) Total PCA score plot, (B) 3D PCA score plot graph, and (C) PCA
482
score graph of changes in metabolic profiles of the RYR from different fermentation
483
stages for the wild-type strain and Monascus PHDS26. 4384: wild-type strain of
484
Monascus aurantiacus Li AS3.4384. P: Monascus PHDS26.
485
Figure 5. (A) The PCA, (B) PLS-DA, and (C) OPLS-DA score plots of metabolic
486
differences in the RYR from the wild-type strain and Monascus PHDS26 after 22
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
487
days of fermentation. 4384: wild-type strain of Monascus aurantiacus Li AS3.4384.
488
P: Monascus PHDS26.
489
Figure 6. General biosynthetic pathways of some metabolites for the wild-type strain
490
and Monascus PHDS26 according to the KEGG database. The solid lines between
491
the compounds indicate a direct relationship between the compounds, and the dotted
492
lines represent indirect relationships.
ACS Paragon Plus Environment
Page 24 of 38
Page 25 of 38
Journal of Agricultural and Food Chemistry
Figure 1 8000
Color value units (U/g)
4384 P
(A)
6000
4000
2000
0 4
6
8
10 12 14 16 18 20 22 24 26 Fermentation time (d)
300 4384 P
Content of citrinin (µg/g)
250
(B)
200 150 100 50 0 4
6
8
10 12 14 16 18 20 22 24 26 Fermentation time (d)
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 2
ACS Paragon Plus Environment
Page 26 of 38
Page 27 of 38
Journal of Agricultural and Food Chemistry
Figure 3 Alanine
0.6
Leucine
As3.4384 PHDS26
As3.4384 PHDS26
0.06
0.5 0.4
0.04
0.3
0.02 0.2
0.00
0.1 4
6
8
10
12
14
16
18
20
22
24
4
6
8
10
Time (d)
2-hydroxybutanoic acid
14
16
18
20
22
24
Time (d)
3.0 2.5
12
8
As3.4384 PHDS26
As3.4384 PHDS26
L-malic acid
6
2.0 1.5
4
1.0
2 0.5
0
0.0 4
6
8
10
12
14
16
18
20
22
24
4
6
8
10
Time (d)
12
14
16
18
20
22
24
Time (d) 1.8
0.020
hydroxylamine
Oxoproline
As3.4384 PHDS26
0.016
1.5
0.012
1.2
0.008
As3.4384 PHDS26
0.9
0.004
0.6 0.000 4
6
8
10
12
14
16
18
20
22
24
0.3
Time (d)
4
6
8
10
12
14
16
18
20
22
24
22
24
Time (d)
0.015 Phenylalanine
As3.4384 PHDS26
0.012
As3.4384 PHDS26
Asparagine
0.009
0.006 0.009
0.003
0.006 0.003
0.000 4
6
8
10
12
14
16
Time (d)
18
20
22
24
4
6
8
10
12
14
16
18
Time (d)
ACS Paragon Plus Environment
20
Journal of Agricultural and Food Chemistry
0.25
As3.4384 PHDS26
0.03 2,3-dihydroxybenzoic acid
Page 28 of 38
Glutamine
As3.4384 PHDS26
0.20
0.02
0.15 0.10
0.01
0.05 0.00
0.00 4
6
8
10
12
14
16
18
20
22
24
4
6
8
10
12
Time (d)
Citric acid
16
18
20
22
24
Myristic Acid
0.0040
As3.4384 PHDS26
1.2
14
Time (d)
As3.4384 PHDS26
0.0032
0.8 0.0024
0.4
0.0016
0.0008
0.0 4
6
8
10
12
14
16
18
20
22
4
24
6
8
10
12
14
16
18
20
22
24
Time (d)
Time (d)
10 Fructose
Mannose
0.09
As3.4384 PHD26
8
0.06
6
As3.4384 PHDS26
4 0.03
2 0.00
0 4
6
8
10
12
14
16
18
20
22
24
4
6
8
10
Time (d)
12
14
16
18
20
22
24
Time (d)
0.012 Allo-inositol
0.05 Myo-inositol
As3.4384 PHDS26
0.009
As3.4384 PHDS26
0.04 0.03
0.006
0.02
0.003
0.01 0.00
0.000 4
6
8
10
12
14
16
18
20
22
4
24
6
8
10
12
14
16
18
20
22
24
Time (d)
Time (d)
2.5 Trehalose
Sorbitol
0.0020
As3.4384 PHDS26
As3.4384 PHDS26
2.0
0.0015
1.5
0.0010
1.0
0.0005
0.5
0.0000
0.0
4
6
8
10
12
14
16
Time (d)
18
20
22
24
4
6
8
10
12
14
16
Time (d)
ACS Paragon Plus Environment
18
20
22
24
Page 29 of 38
Journal of Agricultural and Food Chemistry
Figure 4
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
ACS Paragon Plus Environment
Page 30 of 38
Page 31 of 38
Journal of Agricultural and Food Chemistry
Figure 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
ACS Paragon Plus Environment
Page 32 of 38
Page 33 of 38
Journal of Agricultural and Food Chemistry
Figure 6 D-glucose-6-phosphate Phenylalanine
Shikimate
Trehalose Glucose
Sorbitol
Myo-inositol
Pyruvate
Alanine Leucine
Fructose
Mannose
Asparagine Hydroxylamine Malic acid
Acetyl-CoA
Fatty acid biosynthesis
Ammonia
Myristic acid
Oxoproline Glutamine
Citrate Propanoyl-CoA
Glutamate 2-oxoglutarate
Isocitrate
2-hydroxybutanoic acid
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 34 of 38
Table 1 List of metabolites from RYR identified by GC-TOF-MS according to similarity score no less than 700. No.
Metabolite
Similarity
1
Oxamic acid
705
2
Methylmalonic acid
774
RT
Mass
No.
Metabolite
Similarity
RT
7.11
190
44
Glutamic acid
923
14.18
246
7.37
147
45
Phenylalanine
907
14.42
192
(min)
Mass
(min)
3
Alanine
935
7.54
116
46
Lyxose
709
14.92
217
4
Maleimide
868
7.58
154
47
Asparagine
915
14.96
231
5
Glycine
792
7.76
102
48
D-arabitol
811
15.39
217
6
Hydroxylamine
814
7.79
146
49
2,3-Dihydroxybenzoic acid
891
16.07
355
7
2-Butyne-1,4-diol
745
7.82
147
50
5-Methoxytryptamine
744
16.23
174
8
3-Hydroxypyridine
903
8.12
152
51
D-(glycerol 1-phosphate)
900
16.27
299
9
Methyl phosphate
791
8.43
241
52
Glucose-1-phosphate
806
16.38
217
712
16.48
297
4-Hydroxy-3-methoxybenz 10
2-Hydroxybutanoic acid
776
8.85
131
53
oic acid
11
Valine
902
8.92
144
54
Glutamine
931
16.57
156
12
Urea
882
9.35
189
55
Alizarin
733
16.84
369
13
Benzoic acid
923
9.48
179
56
β-Mannosylglycerate
770
16.94
204
14
Leucine
751
9.64
158
57
Citric acid
938
17.21
273
15
Phosphate
924
9.71
211
58
Ornithine
953
17.24
142
16
Isoleucine
884
9.92
158
59
Myristic acid
816
17.80
285
17
Proline
937
10.05
142
60
Fructose
901
17.96
103
18
Glycine
960
10.12
174
61
Mannose 1
889
18.43
79
19
Succinic acid
877
10.16
247
62
Mannose 2
943
18.82
319
20
2,3-Dihydroxypyridine
881
10.23
240
63
Tyrosine
802
19.31
179
21
D-glyceric acid
934
10.28
189
64
Allo-inositol
737
19.38
318
22
Oxalic acid
709
10.49
147
65
Myo-inositol 1
776
20.47
318
23
Fumaric acid
927
10.59
245
66
Palmitic acid
952
20.97
117
24
Serine
954
10.66
204
67
Myo-inositol 2
910
21.52
305
25
Mevalonic acid lactone
909
11.05
145
68
Galactose 1
839
22.09
319
26
2,3-Dihydroxypyridine
745
11.06
240
69
Galactose 2
824
22.18
103
27
β-alanine
843
11.6
248
70
Heptadecanoic acid
808
22.44
132
28
D-erythronolactone
844
11.67
247
71
Sorbitol
870
23.16
319
29
5-Methoxytryptamine
724
11.69
174
72
Oleic acid
960
23.38
339
30
5-Methylresorcinol
900
11.86
253
73
β-mannosylglycerate 1
911
23.59
132
31
Aminomalonic acid
806
12.11
218
74
β-mannosylglycerate 2
812
24.01
337
32
Citramalic acid
724
12.15
247
75
Linoleic acid
743
24.11
80
33
Glutamine
843
12.23
155
76
Gluconic lactone
752
24.48
204
34
L-malic acid
950
12.31
73
77
β-mannosylglycerate
720
24.97
204
35
Threitol 1
917
12.40
217
78
Sorbitol
833
25.33
319
36
Threitol 2
949
12.55
217
79
2-Monopalmitin
802
25.62
218
ACS Paragon Plus Environment
Page 35 of 38
Journal of Agricultural and Food Chemistry
4-Hydroxy-6-methyl-2-pyr 37
one
38
Aspartic acid
957
12.79
232
81
Sucrose
894
26.08
361
39
Oxoproline
914
12.96
156
82
Adenosine
882
26.16
236
40
4-Aminobutyric acid
923
13.01
174
83
Trehalose
930
26.90
361
3,6-Anhydro-D-galactose
704
13.30
247
84
Isomaltose
833
27.88
361
85
Galactinol
749
28.41
204
41
723
12.60
170
80
1-Monopalmitin
896
25.84
371
2-Hydroxy-3-isopropylbut 42
anedioic acid
941
13.60
275
43
Diglycerol
748
13.82
103
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 36 of 38
Table 2 List of various metabolites from Monascus aurantiacus Li AS3.4384 and Monascus PHDS26 during the different fermentation stages from 6 to 22 days. Time
Metabolite
Similarity
RT (min)
Mass
VIP
p-Value
q-Value
Fold change
Hydroxylamine
814
7.79
146
1.381
0.021
0.033
0.207
2-Butyne-1,4-diol
745
7.82
147
2.106
0.011
0.022
7.792E+06
(days)
6
10
D-Erythronolactone
844
11.67
247
1.146
0.004
0.012
5.228
Diglycerol
748
13.82
103
2.888
0.000
0.004
0.000
2,3-Dihydroxybenzoic acid
891
16.07
355
2.347
0.036
0.046
3.806E+06
Glutamine
931
16.57
156
1.987
0.037
0.047
5.629
Myristic acid
816
17.80
285
1.280
0.001
0.004
4.034
Myo-inositol
776
20.47
318
2.070
0.006
0.015
8.747
β-Mannosylglycerate
720
24.97
204
1.413
0.007
0.017
14.858
Trehalose
930
26.90
361
2.519
0.006
0.014
0.003
Maleimide
868
7.58
154
2.461
0.001
0.004
0.000
Leucine
751
9.64
158
2.093
0.007
0.012
8.241E+07
D-Erythronolactone
844
11.67
247
1.067
0.000
0.003
7.139
Aminomalonic acid
806
12.11
218
1.281
0.021
0.022
3.177
Diglycerol
748
13.82
103
2.633
0.001
0.005
0.000
Phenylalanine
907
14.42
192
1.159
0.002
0.006
3.839
Asparagine
915
14.96
231
1.742
0.000
0.003
13.524
Glutamine
931
16.57
156
1.743
0.000
0.003
6.780
Fructose
901
17.96
103
2.022
0.042
0.035
150.846
Allo-inositol
737
19.38
318
1.316
0.026
0.025
1.409E+06
β-Mannosylglycerate
812
24.01
337
2.672
0.007
0.011
0.000
Trehalose
930
26.90
361
2.809
0.029
0.027
0.000
Maleimide
868
7.58
154
2.330
0.004
0.007
0.000
D-Erythronolactone
844
11.67
247
1.066
0.000
0.000
8.000
950
12.31
73
2.996
0.003
0.006
0.000
L-Malic
14
acid
Asparagine
915
14.96
231
1.022
0.008
0.012
2.944
β-Mannosylglycerate
812
24.01
337
2.542
0.000
0.000
0.000
Sorbitol
833
25.33
319
1.713
0.000
0.000
5.157E+05
Trehalose
930
26.90
361
2.783
0.010
0.014
0.000
Maleimide
868
7.58
154
1.989
0.000
0.001
0.111
2-Hydroxybutanoic acid
776
8.85
131
1.152
0.000
0.000
364.884
ACS Paragon Plus Environment
Page 37 of 38
Journal of Agricultural and Food Chemistry
L-Malic
18
Threitol
950
12.31
73
2.807
0.000
0.002
0.000
917
12.40
217
2.439
0.001
0.002
0.000
Threitol
949
12.55
217
2.924
0.001
0.002
0.000
β-Mannosylglycerate
812
24.01
337
1.456
0.044
0.030
0.000
Trehalose
930
26.90
361
1.522
0.015
0.014
0.197
Alanine 1
935
7.54
116
1.480
0.036
0.024
0.412
2-Hydroxybutanoic acid
776
8.85
131
16.208
0.000
0.000
815.109
950
12.31
73
4.816
0.001
0.002
0.000
Threitol
949
12.55
217
5.547
0.000
0.002
0.000
Oxoproline
914
12.96
156
5.308
0.002
0.004
0.458
Glutamine
931
16.57
156
1.010
0.003
0.004
0.237
Citric acid
938
17.21
273
4.357
0.000
0.002
2.385
Mannose
943
18.82
319
3.765
0.007
0.007
2.259
L-Malic
22
acid
acid
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Graphic for manuscript
Metabolomics
ACS Paragon Plus Environment
Page 38 of 38