Subscriber access provided by CMU Libraries - http://library.cmich.edu
Article
Ultrafast polyphenol metabolomics of red wines using microLC-MS/MS Yan Ma, Nobuo Tanaka, Arpana Vaniya, Tobias Kind, and Oliver Fiehn J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.5b04890 • Publication Date (Web): 23 Dec 2015 Downloaded from http://pubs.acs.org on December 26, 2015
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 28
Journal of Agricultural and Food Chemistry
Ultrafast polyphenol metabolomics of red wines using microLC-MS/MS
Yan Ma(1), Nobuo Tanaka(1,2) , Arpana Vaniya(1), Tobias Kind(1) and Oliver Fiehn(1,3) (1) UC Davis Genome Center – Metabolomics, University of California, Davis, California 95616, United States (2) GL Sciences Inc., Iruma, Saitama 358-0032, Japan (3) Department of Biochemistry, Faculty of Sciences, King Abdulaziz University, Jeddah, Saudi-Arabia
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 28
1
Abstract
2
The taste and quality of red wine is determined by its highly complex mixture of polyphenols
3
and many other metabolites. While no single method can fully cover the full metabolome, even
4
for polyphenols and related compounds, current methods proved inadequate. We optimized
5
liquid chromatography resolution and sensitivity using 1 mm I.D. columns with microLC pumps
6
and compared data-dependent to data-independent (SWATH) MS/MS acquisitions. A high-
7
throughput microLC-MS method was developed with a 4 min gradient at 0.05 mL/min flow rate
8
on a Kinetex C18 column and Sciex TripleTOF mass spectrometry. Using the novel software
9
MS-DIAL, we structurally annotated 264 compounds including 165 polyphenols in six
10
commercial red wines by accurate mass MS/MS matching. As proof of concept, multivariate
11
statistics revealed the difference in the metabolite profiles of the six red wines and regression
12
analysis linked the polyphenol contents to the taste of the red wines.
13
Keywords
14
red wines, polyphenols, metabolomics, 1 mm I.D. column, microLC-MS/MS
15 16 17 18 19 20 21 22 23
2 ACS Paragon Plus Environment
Page 3 of 28
Journal of Agricultural and Food Chemistry
24
Introduction
25
Red wine is a prime example how food quality is determined by its highly complex mixture of
26
metabolites. While all wines are characterized by their abundance of alcohol, sugar and simple
27
organic acids, red wine is specifically known for its diversity of phenolic compounds, including
28
phenolic acids, flavonoids, anthocyanins and tannins.1 For example, 111 polyphenols were
29
reported in red wines according to Phenol-Explorer,2 an online database of polyphenol
30
contents in food. All those complex components contribute together to the color, aroma, taste
31
and mouthfeel of the red wines. For example, vanillin and 4-allyl-2-methoxyphenol were
32
reported to contribute to the wood-vanillin-cinnamon aroma in aged red wines.3
33
The polyphenol components in red wines have drawn a lot of interest due to their high
34
antioxidant potential.4 One well-known example is resveratrol, which has been studied
35
intensively as an anti-aging and anti-carcinogenic agent.5 The phenolic composition also plays
36
a great role in the color stability and taste of the red wines.6 To qualitative and quantitatively
37
analyze the polyphenols and other components in red wine, reversed-phase high performance
38
liquid chromatography (HPLC) methods were developed with either ultraviolet (UV) or mass
39
spectrometry (MS) detection.7, 8 However, most of the current methods only measured a few
40
abundant components of red wines9, 10 and were inadequate for giving full coverage of these
41
complex metabolites. Hence, we set out to develop and apply a suite of advanced
42
metabolomic methods to utilize comprehensive profiling of polyphenols and related
43
compounds and test these methods by chemometrics.
44
Miniaturization of column and instrumentation is one of the leading trends in the field of liquid
45
chromatography. Ultra-high pressure liquid chromatography (UHPLC) columns with internal
46
diameters (I.D.) of 2.1 mm have long taken over conventional HPLC columns of 4.6 mm for
47
their higher sensitivity, faster separation and less consumption of organic solvents.11 For the 3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 28
48
same reasons, 1 mm I.D. columns have been developed to further improve the sensitivity
49
compared to 2.1 mm I.D. columns.12 Consequently, 1 mm I.D. columns have been applied to a
50
variety of studies of small molecules such as amino acids,13 flavonoids,14 and antibacterial
51
drugs.15
52
constraints, because both the flow rate and the extra-column dispersion source need to be
53
scaled down for optimal performance as the column volume decreases.16 Recently,
54
commercial instruments designed for micro LC performance have become available, such as
55
the Eksigent MicroLC 200 system and Dionex UltiMate 3000 Quaternary Micro LC system. A
56
previous study showed that the microLC instrument provided significantly higher efficiency for
57
1 mm I.D. columns compared to common UHPLC instrument.17 At the same time, more small-
58
sized columns have been developed.18, 19 We here evaluated different columns and optimized
59
microLC methods for comprehensive polyphenol profiling in red wines.
60
Materials and Methods
61
Comparison of 1 mm I.D. columns using polyphenol standards. An Eksigent MicroLC 200
62
was coupled to a SCIEX TripleTOF 5600 (SCIEX, Redwood City, CA, USA) for MS
63
measurements from m/z 60 to 1250 in electrospray positive mode with spray voltage 5.5 kV at
64
scan rate of 10 spectra per second. Three reversed-phase C18 columns of 1 mm I.D. and 5
65
cm length were examined, including Kinetex C18 2.6 µm (Phenomenex, Torrance, CA, USA),
66
InertSustain C18 2.0 µm (GL Sciences, Tokyo, Japan) and Acquity CSH 1.7 µm (Waters,
67
Milford, MA, USA). The Kinetex and InertSustain columns were prototype columns, while
68
Acquity was a commercially available column. The standard mixture for LC/MS experiments
69
contained (1) Umbelliferone, (2) Rutin, (3) Naringin, (4) Quercetin, (5) Naringenin, (6) 6-
70
Hydroxyflavone, (7) Genistein, (8) Genistin, (9) Daidzein, (10) Daidzin, (11) Glycitein, (12)
71
Glycitin, (13) Chrysin at 1.0 µg/mL each, in 10:90 methanol:water. Mobile phase A was water
However, practical use of 1 mm I.D. columns is still limited due to instrumental
4 ACS Paragon Plus Environment
Page 5 of 28
Journal of Agricultural and Food Chemistry
72
with 0.1% acetic acid while mobile phase B was acetonitrile with 0.1% acetic acid. The LC
73
gradient was 10% to 90% B with different flow rates (0.025, 0.05, 0.075, 0.1 mL/min) and
74
gradient time or gradient volume (4 min, 8 min, and 0.6 mL). Injection volume was 0.2 µL. All
75
experiments were carried out in duplicate. Data was processed by the SCIEX MultiQuant 2.1
76
software.
77
LC gradient optimization using a red wine example. To optimize LC gradients for total peak
78
discovery, metabolomic tests were performed with wine extracts from a 2012 Petite Sirah
79
(Concannon Vineyard, Livermore, CA). 500 µL of raw wine was added to an ice-cold pre-
80
mixed extraction solvent consisting of 500 µL water and 500 µL ethyl acetate. Samples were
81
shaken for 5 min, centrifuged for 3 min at 16,100 rcf, and then 300 µL of the ethyl acetate
82
phase was transferred to a second vial and evaporated to dryness using a speed vacuum
83
concentrator. As method blank control, instead of raw wine, 500 µL of water was added to the
84
extraction solvent and then the same protocol was applied. On the day of the experiment,
85
samples were reconstituted in 150 µL of 90:10 acetonitrile:water, briefly vortexed, sonicated for
86
5 min, centrifuged for 3 min at 16,100 rcf then transferred to HPLC vials with inserts. The
87
same gradient series and mass spectrometry measurements were applied as to the
88
polyphenol standards. The raw WIFF data files were converted to ABF files with Reifycs
89
Analysis Base File Converter (Reifycs Inc., Tokyo, Japan) and processed by the freely
90
available MS-DIAL20 software (ver. 1.58) (http://prime.psc.riken.jp/Metabolomics_Software/).
91
Metabolome composition analysis of six red wines using the optimized microLC-MS/MS
92
method. Six bottles of red wines were opened and extracted fresh, including a 2009 Syrah
93
(Putah Creek Winery, Solano county, CA), a 2012 Pinot Noir (Black Mountain Vineyard,
94
Sonoma County, CA), a 2010 Cabernet Sauvignon (Mosaic Winery, Alexander Valley, Sonoma
95
County, CA), a 2012 Zinfandel (Cline Cellars Winery, Sonoma County, CA), a 2012 Petite 5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 28
96
Sirah (Blue Fin Winery, Napa, CA), and a 2012 Merlot ( Blue Canyon Wine Company,
97
Greenfield, CA). A panel of fourteen non-trained adults tasted the remaining wines directly
98
afterwards. Overall likeability21 (averaging aroma, flavor and finish) was scored from 1 (worst)
99
to 10 (best) for each wine with blinded labels. Raw wines were extracted in the similar way as
100
the preliminary test, but instead of reconstituting the extracts with 150 µL solvent, 75 µL was
101
used to concentrate the samples. Six technical replicates were prepared for each wine as well
102
as the method blank. As result of the microLC-MS method optimization, a gradient of 4 min at
103
a flow rate of 0.05 mL/min was selected for data-dependent MS/MS acquisition for these
104
wines. To compare the data-dependent and data-independent (SWATH) method, a 12-min
105
gradient was used at a flow rate of 0.05 mL/min. MS/MS data was acquired in both positive
106
and negative mode at a collision energy of 30 ± 15 V. For data-dependent acquisition of
107
MS/MS, precursor ions from m/z 60 to 1250 were fragmented with a cycle time of 0.125 s. For
108
data-independent acquisition (SWATH), precursor ions from m/z 60 to 800 were fragmented,
109
with a cycle time of 0.4 s. A SWATH window of 50 Da was used.
110
Data processing, statistics analysis and visualization. Peak picking and alignment was
111
performed by MS-DIAL (ver. 1.58). Representative MS/MS spectra were exported from MS-
112
DIAL in MGF format for compound identification. NIST MS PepSearch was used to search the
113
MGF file against MS/MS libraries, including NIST14, MassBank,22 ReSpect23 and METLIN.24
114
For the remaining unknown compounds, the freely available MS2Analyzer25 software
115
(https://sourceforge.net/projects/ms2analyzer/) was used to search for characteristic neutral
116
losses for glycosides, including pentosides (132.063 Da), hexosides (162.053 Da),
117
deoxyhexosides (146.058 or 164.068 Da), acetylhexosides (204.063 Da), and hexuronide
118
(194.043 Da). An m/z window of 0.005 Da and a relative intensity threshold of 0.5 were
119
selected as input parameters. Each hit from MS2Analyzer output was queried by its
6 ACS Paragon Plus Environment
Page 7 of 28
Journal of Agricultural and Food Chemistry
120
neutralized molecular mass in Phenol-Explorer and by its precursor m/z value in the online
121
METLIN database to extract potential structure candidates.
122
Principle component analysis (PCA) was performed for all the annotated compounds (reverse-
123
dot score>800) in positive ionization mode with MetoboAnalyst.26 Before PCA analysis,
124
samples were normalized by sum, followed by log transformation and pareto scaling.
125
Orthogonal signal correction partial least squares (O-PLS) modeling was performed between
126
the average tasting scores and peak intensities of all the annotated compounds using Devium
127
(https://github.com/dgrapov/devium). Two latent variables (LVs) were determined based on
128
the criteria of maximum Q2 in the range of LV1 to LV10. Twenty permutation tests were used
129
to validate the model, and ten top features were selected for each sign. At last, a chemical
130
structural (Tanimoto>0.9) and mass spectral similarity (similarity>0.4) network was created for
131
all the annotated compounds in positive ionization mode with MetaMapR27 and Cytoscape.28
132
MS/MS spectra were pre-treated (m/z rounded up to 0.1, intensity filter>5% of base peak)
133
before the analysis to improve the similarity score.
134
Results and Discussion
135
Method optimization: concepts and methods. A good chromatographic method is essential
136
to LC-MS experiments to minimize ion suppression effects, to help distinguish isomers, and to
137
increase the number of features detected and annotated in untargeted metabolomics. In this
138
work, we used two parameters, peak capacity and peak height, to quantify the performance of
139
chromatographic separations. Peak capacity is defined as number of peaks that can be
140
separated with unit resolution within the gradient time with a one-channel detector. Peak
141
intensity, here defined as peak height, is a measure of method sensitivity. It is noteworthy that
142
peak intensity is not only related to sample concentration, injection volume and detector
143
sensitivity, but also strongly influenced by chromatographic performance. For example, 7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 28
144
changing the column internal diameter from 2.1 mm to 1 mm can increase peak intensity by a
145
factor of 4.41.12
146
However, chromatography is a process involving many factors, including stationary phase,
147
particle size, column size, mobile phase, flow rate, gradient, etc. To optimize the microLC-MS
148
method for phenolic compounds, we used different stationary phases, particle sizes, flow rates
149
and gradients to separate thirteen standard phenol compounds as well as a preliminary red
150
wine sample and evaluate the performance of chromatography based on peak capacity, peak
151
height and peak numbers extracted from red wine sample.
152
Method optimization: column selection with phenol standards. First we compared the
153
performance of three 1 mm I.D., 5 cm length columns in terms of peak capacity and intensity:
154
the Kinetex C18 (particle size, 2.6 µm), the InertSustain C18 (2.0 µm), and the Acquity CSH
155
(1.7 µm) columns. To test the effects of various liquid chromatography conditions on those
156
columns, experiments were performed under eleven gradient conditions at flow rates of 0.025
157
to 0.1 mL/min. Averaged peak widths and peak capacities were calculated for all the standard
158
peaks as in the eq. 1-2.
159
Peak width = ( peak width at half height / 2.35) × 4
(1)
160
Peak capacity = gradient time / peak width
(2)
161
Figure 1 shows extracted ion chromatograms (EICs) of the phenol standards separated by the
162
Kinetex column in a 4 min gradient at 0.05 mL/min, with an average peak width of around 2 s.
163
Figure 2 shows peak capacities and average peak intensities for the three columns under one
164
representative gradient condition (see all the results in Table S1). In general, the tested
165
polyphenols were eluted at later retention times using the InertSustain column in comparison
166
to the Acquity column, hence, showed increased ionization efficiency and wider peaks 8 ACS Paragon Plus Environment
Page 9 of 28
Journal of Agricultural and Food Chemistry
167
because they eluted at higher organic/water solvent ratios. Taking the Kinetex column as an
168
example to investigate the chromatography under different LC conditions, we observed the
169
following trends: 1) peak capacity increased as the gradient time and flow rate increased; 2)
170
peak height, by contrast, decreased as the gradient time and flow rate increased; 3) when the
171
gradient volume (flow rate × gradient time) was constant, peak capacity and peak height also
172
remained relatively constant. Here we estimated the optimum conditions by calculating and
173
comparing the product of the peak capacity and peak height for each gradient condition.
174
Based on the calculation (Table S1), an optimum condition for the 1 mm I.D. column was found
175
using a 4 min gradient at 0.05 mL/min for the given polyphenol standards. Overall, under the
176
conditions employed here, the Kinetex column provided the highest peak capacity and
177
intensity among the three columns. We therefore selected this column in all following wine
178
polyphenol studies. Other columns may perform better under optimized conditions for each
179
column.
180
Method optimization: gradient finalization with red wine sample. The same gradient
181
experiments were repeated with an extracted test red wine sample to evaluate the effect of the
182
LC parameters under realistic conditions. The data was processed by MS-DIAL with an
183
intensity threshold of 4,000 counts. The numbers of detected m/z-retention time (rt) features is
184
shown in Table 1. The trend was similar to the trend of peak intensities, indicating that peak
185
intensity contributed more to the overall m/z-rt count than peak capacity. With a given intensity
186
threshold for the data processing, higher peak intensities generally led to more detected peaks
187
even if peaks were a little broader, on average, than comparative methods. Hence, the peak
188
capacity parameter was less important to overall feature detections because mass
189
spectrometers are very effective multi-channel detectors, enabling selective quantification of
190
many compounds even if these are co-eluting (except for isomers). We observed that
191
separation of isomers was related to the peak capacity, but it was also affected by the 9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 28
192
selectivity of the stationary phase and mobile phase. In conclusion, the 4 min gradient at flow
193
rate of 0.05 mL/min that was selected by polyphenol standards was also confirmed to be
194
highly suitable in the comprehensive metabolomics analysis of red wine samples.
195
MicroLC-MS/MS analysis of six red wines: Compound identification. Six red wines from
196
different grape varieties and brands were analyzed to showcase the application of the
197
optimized microLC-MS/MS method. Six replicates of each red wine were extracted with ethyl
198
acetate29 and measured in positive and negative electrospray ionization mode. Data-
199
dependent MS/MS spectra were acquired to facilitate structural annotation of metabolites. MS-
200
DIAL was used for peak detection and alignment, while NIST MS PepSearch was used to
201
search the MS/MS spectra against libraries including NIST14, MassBank, ReSpect and
202
METLIN. 264 compounds were annotated with a weighted reverse-dot score30 higher than
203
800. Among the 264 compounds, 165 metabolites were classified as phenol compounds. The
204
full list of annotated compounds is presented in Table S2 including detailed metadata such as
205
retention times, precursor m/z, InChIKeys, and intensity tables. To further validate the
206
annotation results we confirmed 18 compounds by injecting authentic reference standards
207
(see Table 2).
208
For many peaks, MS/MS library searches assigned the same compound name to multiple m/z-
209
rt features. By definition, compounds that show different physicochemical properties (such as
210
retention time) but having the same exact mass are isomers. We further defined here close
211
positional isomers if peaks were well separated and MS/MS library queries yielded reverse-dot
212
score higher than 800. Among the 264 annotated compounds, 133 were found with one or
213
more isomers, due to the complexity of polyphenol structures and high possibility for positional
214
isomers. Figure 3 demonstrates the detection of a novel resveratrol isomer in addition to the
215
well-known cis- and trans-forms. Other abundant isomers include, for example, iso10 ACS Paragon Plus Environment
Page 11 of 28
Journal of Agricultural and Food Chemistry
216
dihydroquercetin,
iso-procyanidin
B2,
and
iso-caffeic
217
characterization of such isomers is impossible by MS alone, given the large number of
218
theoretically possible positional isomers, the limited availability of authentic standards and the
219
complexity of wine samples that forbids easy routes for peak purifications and nuclear
220
magnetic resonance (NMR) determinations. For example, for five dihydroquercetin isomers
221
detected in the red wines samples, 95 candidate structures were found in the structure
222
database ChemSpider, but only one compound is available to be purchased from Sigma-
223
Aldrich.
224
Although 264 compounds have been annotated using MS/MS library search, more than 1000
225
features which are related to true compounds still remain unknown. To further characterize the
226
structures or substructures of unknowns, we used MS2Analyzer in combination with online
227
databases such as Phenol-Explorer and METLIN. MS2Analyzer is a software tool using
228
characteristic neutral losses and product ions for small molecule substructure annotation.
229
Since many polyphenols are in glycoside forms, we queried the characteristic neutral losses
230
for pentosides, hexosides, deoxyhexosides, acetylhexosides, and hexuronide in all the
231
unknown MS/MS spectra with MS2Analyzer. For those MS/MS spectra containing the above
232
neutral losses, their precursors were searched in the Phenol-Explorer and METLIN databases.
233
Potential candidates with matching accurate mass and substructures for 28 features were
234
reported in Table S3. Figure 4 shows an example of MS/MS spectrum, characteristic neutral
235
loss, and potential structure of unknown m/z 521.132 eluting at 1.68 min.
236
Over the past few years, data independent MS/MS acquisition (DIA, specifically SWATH)
237
method has drawn a lot of interests from LC-MS researchers. In contrast to the conventional
238
data dependent MS/MS acquisition (DDA) method, DIA allows the acquisition of MS/MS
239
spectra for all the MS peaks simultaneously regardless of their peak intensities. However, DIA 11 ACS Paragon Plus Environment
acid.
A
complete
structural
Journal of Agricultural and Food Chemistry
Page 12 of 28
240
is not very suitable for the narrow peaks in our 4 min gradient; longer gradients need to be
241
applied to ensure there are enough data points in each chromatographic peak so MS/MS
242
deconvolution can be performed. To compare DDA and SWATH methods fairly, we extended
243
the gradient time to 12 min and analyzed the six red wines using both methods in positive
244
mode. After data processing, we observed a similar total number of m/z-rt features detected
245
for both methods (1632 for SWATH and 1667 for DDA). However, the SWATH method
246
provided more features with MS/MS spectra (1632 VS 1461) and hence more compound
247
annotations (207 vs 178 structurally annotated peaks). In summary, SWATH method expands
248
the coverage of the untargeted LC-MS profiling experiments, but may need to sacrifice the
249
chromatography performance to some extent.
250
MicroLC-MS/MS analysis of six red wines: Statistics and Visualization. Using relative
251
quantitation and principle component analysis (PCA) clustering of all annotated compounds in
252
positive ionization mode (Figure 5), clear differences were observed in the metabolite profiles
253
of the six red wines. Loadings for PC1 and PC2 were available in Table S4 in Supporting
254
Information. Specifically, Syrah and Cabernet Sauvignon were separated very well from the
255
four other wines along PC1, indicating that they have the overall most similar metabolite
256
profiles. For example, triethyl citrate and syringic acid were found to be much more abundant
257
in Cabernet Sauvignon and Syrah than in other four wines, shown in Figure 6. Similarly, the six
258
wines were separated on PC2, but not clearly clustered into two groups like on PC1. Other
259
metabolites also showed very large variance among the tested wines. For example, trans-
260
resveratrol was up to 6-fold more abundant in Pinot Noir and Merlot than in the other four
261
wines, while the novel iso-resveratrol showed a similar pattern but was less abundant in
262
Merlot. By contrast, caffeic acid and one of its isomer showed very different distribution
263
patterns in the six red wines (Figure 6). To improve the accuracy of relative quantitation, a set
264
of internal standards could be developed to normalize the peak intensities. Considering the 12 ACS Paragon Plus Environment
Page 13 of 28
Journal of Agricultural and Food Chemistry
265
complexity of compounds measured and their various physical and chemical properties, a
266
single internal standard was not adequate to address the difference in the extraction and
267
ionization processes,
268
One of the potential applications of our untargeted metabolic data was to link the taste of red
269
wines to their chemical components. As a prove of concept, we used the subjective tasting
270
scores from a panel of fourteen untrained adults (Table S5) and our metabolic data to build an
271
orthogonal signal correction partial least squares (O-PLS) model. Q2 and root mean squared
272
error of prediction (RMSEP) for LV2 was 0.990 and 0.158, respectively. Features that were
273
positively or negatively associated with the tasting scores could be selected, as shown in
274
Figure S1 and Table S6 in Supporting Information. The sample size in this sensory study may
275
not be sufficient to reveal the true association between the red wine components and overall
276
sensory experience; a larger scale study with more red wines and panelists would be needed
277
for that purpose.
278
In summary, the optimized microLC-MS/MS screening workflow combined a novel LC method
279
with chemometrics strategies from untargeted metabolomics, and was applied to the analysis
280
of red wines. It has multiple advantages compared to other published LC-MS method in the
281
red wine studies, shown in Table 3. First of all, our method is only 5.5 min long including
282
column wash and re-equilibration, allowing a cost effective high throughput operation in large
283
trials. Secondly, with a flow rate of 0.05 mL/min, each sample just needs 0.275 mL of mobile
284
phase, thus it only take 27.5 mL of mobile phase to run 100 sample which saves costs for both
285
purchase and discarding organic solvents. Last, with untargeted data acquisition and
286
chemometrics method, we are able to structurally annotate more than 200 compounds using
287
comprehensive MS/MS libraries in a single run. Compared to other methods focusing on few
13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 28
288
target known compounds, metabolomics is clearly more comprehensive and will build much
289
more robust methods for sensory tasting and quality control outputs.
290
Supporting Information
291
The Supporting Information is available free of charge on the ACS Publications website.
292
Table S1. Peak capacities, average peak heights and their products of 13 phenol standards
293
separated by three 1 mm I.D. columns on Ekigent MicroLC-MS.
294
Table S2. All annotated compound information in the six red wines.
295
Table S3. Compound class and tentative structure annotation of unknown compounds in
296
glycoside forms.
297
Table S4. Loadings for PC1 and PC2 in the PCA analysis.
298
Table S5. Range, mean and standard deviation of the tasting scores
299
Table S6. Ten most important features positively or negatively correlated to the tasting score,
300
selected by the OPLS modeling.
301
Figure S1. Orthogonal signal correction partial least squares (O-PLS) selection for the 10
302
most important features positively or negatively correlated to the average tasting score of an
303
untrained panel of 14 adults.
304
Author Information
305
Corresponding author
306
*(O.F.) Email:
[email protected]. Phone: +1-530-754-8258. Fax: +1-530-754-9658
307
Funding 14 ACS Paragon Plus Environment
Page 15 of 28
Journal of Agricultural and Food Chemistry
308
This work was financially supported by the following grants: NIH U24 DK097154, NSF MCB
309
1139644. Wines were exclusively purchased from private funds and did not use any federal
310
funds.
311
Notes
312
The authors declare no competing financial interest.
313
Acknowledgement
314
The authors thank Masayoshi Ohira of GL Sciences and Tivadar Farkas of Phenomenex for
315
providing the prototype columns.
316
317
318
319
320
321
322
323
324
325
15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 28
326
Reference
327
1.
328
The Chemistry of Wine Stabilization and Treatment. In Wiley: 2006; Vol. 2.
329
2.
330
Lobato, P.; Manach, C.; Knox, C.; Eisner, R.; Wishart, D. S., Phenol-Explorer 3.0: a major
331
update of the Phenol-Explorer database to incorporate data on the effects of food processing
332
on polyphenol content. Database 2013, 2013, bat070.
333
3.
334
Properties from Aroma Chemical Composition. Partial Least Squares Regression Models.
335
Journal of agricultural and food chemistry 2003, 51, 2700-2707.
336
4.
337
American journal of clinical nutrition 2005, 81, 215S-217S.
338
5.
339
comprehensive review of human clinical trials. Molecular nutrition & food research 2011, 55,
340
1129-1141.
341
6.
342
measurement and sensory analysis of BC commercial red wines. Food Research International
343
2007, 40, 92-100.
344
7.
345
ultraviolet determination of organic acids and phenolic compounds in red wine and must.
346
Journal of Chromatography A 2004, 1052, 211-215.
347
8.
348
MS/MS analysis of phenols for classification of red wine according to geographic origin, grape
349
variety and vintage. Food chemistry 2010, 122, 366-372.
Ribereau-Gayon, P.; Glories, Y.; Maujean, A.; Dubourdieu, D., Handbook of Enology:
Rothwell, J. A.; Perez-Jimenez, J.; Neveu, V.; Medina-Remon, A.; M'Hiri, N.; García-
Aznar, M.; López, R.; Cacho, J.; Ferreira, V., Prediction of Aged Red Wine Aroma
Scalbert, A.; Johnson, I. T.; Saltmarsh, M., Polyphenols: antioxidants and beyond. The
Smoliga, J. M.; Baur, J. A.; Hausenblas, H. A., Resveratrol and health–a
Cliff, M. A.; King, M. C.; Schlosser, J., Anthocyanin, phenolic composition, colour
Kerem, Z.; Bravdo, B.-a.; Shoseyov, O.; Tugendhaft, Y., Rapid liquid chromatography–
Jaitz, L.; Siegl, K.; Eder, R.; Rak, G.; Abranko, L.; Koellensperger, G.; Hann, S., LC–
16 ACS Paragon Plus Environment
Page 17 of 28
Journal of Agricultural and Food Chemistry
350
9.
Mirnaghi, F. S.; Mousavi, F.; Rocha, S. M.; Pawliszyn, J., Automated determination of
351
phenolic compounds in wine, berry, and grape samples using 96-blade solid phase
352
microextraction system coupled with liquid chromatography–tandem mass spectrometry.
353
Journal of Chromatography A 2013, 1276, 12-19.
354
10.
355
hydrophilic interaction LC/MS for the analysis of polar and nonpolar phenols in wine. Journal of
356
separation science 2013, 36, 1379-1388.
357
11.
358
comparative study of new trends in HPLC. International Journal of Advances in
359
Pharmaceutical Analysis 2012, 2, 53-55.
360
12.
361
Chromatography A 1999, 856, 117-143.
362
13.
363
Enantioseparation of N-derivatized amino acids by micro-liquid chromatography using
364
carbamoylated
365
Chromatography A 2014, 1363, 207-215.
366
14.
367
Confirmation of diosmetin 3-O-glucuronide as major metabolite of diosmin in humans, using
368
micro-liquid-chromatography–mass spectrometry and ion mobility mass spectrometry.
369
Analytical and Bioanalytical Chemistry 2013, 405, 8295-8310.
370
15.
371
in plasma by micro ‐ liquid chromatography – mass spectrometry: An adapted method for
372
neonates. Biomedical Chromatography 2011, 25, 827-832.
Greco, G.; Grosse, S.; Letzel, T., Serial coupling of reversed‐phase and zwitterionic
Ashokkumar, T.; Suneetha, A.; Reddy, B. V.; Appaji, S. C. V., Review article on
Vissers, J. P., Recent developments in microcolumn liquid chromatography. Journal of
Wang, Q.; Feng, J.; Han, H.; Zhu, P.; Wu, H.; Marina, M. L.; Crommen, J.; Jiang, Z.,
quinidine
functionalized
monolithic
stationary
phase.
Journal
of
Silvestro, L.; Tarcomnicu, I.; Dulea, C.; Attili, N. R. B.; Ciuca, V.; Peru, D.; Savu, S. R.,
Grondin, C.; Zhao, W.; Fakhoury, M.; Jacqz‐Aigrain, E., Determination of ciprofloxacin
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 28
373
16.
Gritti, F.; Guiochon, G., On the extra-column band-broadening contributions of modern,
374
very high pressure liquid chromatographs using 2.1 mm ID columns packed with sub-2µm
375
particles. Journal of Chromatography A 2010, 1217, 7677-7689.
376
17.
377
Furuno, M.; Nakanishi, K.; Minakuchi, H.; Mriziq, K.; Farkas, T.; Fiehn, O.; Tanaka, N.,
378
Efficiency of short, small-diameter columns for reversed-phase liquid chromatography under
379
practical operating conditions. Journal of Chromatography A 2015.
380
18.
381
Temperature Separations on a Polymer-Coated Fibrous Stationary Phase in Microcolumn
382
Liquid Chromatography. Analytical Sciences 2011, 27, 811.
383
19.
384
microcolumn liquid chromatography. Journal of separation science 2014, 37, 2633-2640.
385
20.
386
VanderGheynst, J.; Fiehn, O.; Arita, M., MS-DIAL: data-independent MS/MS deconvolution for
387
comprehensive metabolome analysis. Nature methods 2015.
388
21.
389
production practices in wines, considering asymmetric information and sensory effects.
390
Journal of Agricultural Economics 2013, 64, 483-504.
391
22.
392
Tanaka, S.; Aoshima, K., MassBank: a public repository for sharing mass spectral data for life
393
sciences. Journal of mass spectrometry 2010, 45, 703-714.
394
23.
395
K.; Sakurai, T.; Matsuda, F.; Aoki, T., RIKEN tandem mass spectral database (ReSpect) for
396
phytochemicals: a plant-specific MS/MS-based data resource and database. Phytochemistry
397
2012, 82, 38-45.
Ma, Y.; Chassy, A.; Miyazaki, S.; Motokawa, M.; Morisato, K.; Uzu, H.; Ohira, M.;
Nakane, K.; Shirai, S.; Saito, Y.; Moriwake, Y.; Ueta, I.; Inoue, M.; Jinno, K., High-
An, R.; Weng, Q.; Li, J., Silica‐particle‐supported zwitterionic polymer monolith for
Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Kanazawa, M.;
Schmit, T. M.; Rickard, B. J.; Taber, J., Consumer valuation of environmentally friendly
Horai, H.; Arita, M.; Kanaya, S.; Nihei, Y.; Ikeda, T.; Suwa, K.; Ojima, Y.; Tanaka, K.;
Sawada, Y.; Nakabayashi, R.; Yamada, Y.; Suzuki, M.; Sato, M.; Sakata, A.; Akiyama,
18 ACS Paragon Plus Environment
Page 19 of 28
Journal of Agricultural and Food Chemistry
398
24.
Smith, C. A.; O'Maille, G.; Want, E. J.; Qin, C.; Trauger, S. A.; Brandon, T. R.; Custodio,
399
D. E.; Abagyan, R.; Siuzdak, G., METLIN: a metabolite mass spectral database. Therapeutic
400
drug monitoring 2005, 27, 747-751.
401
25.
402
molecule substructure annotations from accurate mass MS/MS spectra. Analytical chemistry
403
2014, 86, 10724-10731.
404
26.
405
2.0—a comprehensive server for metabolomic data analysis. Nucleic acids research 2012, 40,
406
W127-W133.
407
27.
408
metabolomic network analysis incorporating unknowns. Bioinformatics 2015, 31, 2757-2760.
409
28.
410
features for data integration and network visualization. Bioinformatics 2011, 27, 431-432.
411
29.
412
capacity
413
spectrophotometric methods. Food Research International 2012, 45, 145-154.
414
30.
415
algorithms for compound identification. Journal of the American Society for Mass Spectrometry
416
1994, 5, 859-866.
417
31.
418
MS/MS characterization of phenolic compounds in white wine stored without added sulfite.
419
Food Research International 2014, 66, 207-215.
420
32.
421
of 20 components in red wine by LC‐MS: Application to variations of red wine components in
422
decanting. Journal of separation science 2012, 35, 2884-2891.
Ma, Y.; Kind, T.; Yang, D.; Leon, C.; Fiehn, O., MS2Analyzer – a software for small
Xia, J.; Mandal, R.; Sinelnikov, I. V.; Broadhurst, D.; Wishart, D. S., MetaboAnalyst
Grapov, D.; Wanichthanarak, K.; Fiehn, O., MetaMapR: pathway independent
Smoot, M. E.; Ono, K.; Ruscheinski, J.; Wang, P.-L.; Ideker, T., Cytoscape 2.8: new
Porgalı, E.; Büyüktuncel, E., Determination of phenolic composition and antioxidant of
native
red
wines
by
high
performance
liquid
chromatography
and
Stein, S. E.; Scott, D. R., Optimization and testing of mass spectral library search
Pati, S.; Crupi, P.; Benucci, I.; Antonacci, D.; Di Luccia, A.; Esti, M., HPLC-DAD–
Cui, Y.; Li, Q.; Liu, Z.; Geng, L.; Zhao, X.; Chen, X.; Bi, K., Simultaneous determination
19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
423 424
20 ACS Paragon Plus Environment
Page 20 of 28
Page 21 of 28
Journal of Agricultural and Food Chemistry
Figure captions Figure 1. Extracted ion chromatograms of 13 polyphenol reference standards using the 1 mm I.D. Kinetex C18 column and a 4 min acetonitrile/water gradient at 0.05 mL/min flow rate. Sensitivity and peak capacity calculations were performed on the more abundant peaks if minor impurities were detectable as well in the authentic standards chromatograms. Figure 2. Comparison of three 1 mm internal diameter C18 columns with respect to average peak intensities and peak capacities for 13 polyphenol reference standards using a 4 min acetonitrile/water gradient at 0.05 mL/min flow rate. Figure 3. Example of compound annotation using MS/MS matching. Annotation of a novel resveratrol isomer showing (A) the extracted ion chromatogram at m/z 229.085 ([M+H]+) at 1.54 min and (B) matching the accurate mass MS/MS spectra of the new isomer against cisresveratrol authentic standards. Figure 4. Example of substructure annotation for unknown compounds using MS2Analyzer. An unknown MS/MS spectrum with precursor m/z 521.132 at 1.68 min in positive mode was annotated as acetylhexoside for the characteristic neutral loss of 204.067. One potential structure candidate, petunidin 3-O-(6''-acetyl-glucoside), was extracted from Phenol-Explorer using accurate mass match. Figure 5. Principal component analysis (PCA) score plot of tested six red wines. Total quantitative variances of 177 annotated peaks in positive ionization mode were clustered to reveal the difference and relative similarities of metabolite profiles of different red wines. 95% confidence regions (Hotelling’s T2 eclipse) were displayed for each class. Figure 6. Peak intensities of selected compounds contributing to the difference of the metabolite profiles of the six red wines: (A) trans-resveratrol; (B) iso-resveratrol; (C) caffeic 21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 22 of 28
acid, negatively correlated to the tasting score; (D) iso-caffeic acid, positively correlated to the tasting score; (E) malathion, a pesticide; (F) triethyl citrate, only present in two of red wines. Figure 7. Chemical structure and mass spectral similarity network (Tanimoto>0.9, spectra similarity>0.4) for compounds annotated in positive ionization mode.
Peak intensity fold
change of Syrah (highest tasting score) versus Pinot Noir (lowest tasting score) was mapped onto the network. OPLS selected features were marked with boarders.
Tables Table 1. Number of peaks detected in one testing red wine sample using Kinetex C18 1.0 mm I.D. column under various LC conditions. The gradient was 10% to 90% mobile phase B (acetonitrile with 0.1% acetic acid) at different flow rate (FR, in mL/min) and gradient time (tG, in min).
tG FR
4
6
8
0.025
571
472
0.050
754
572
0.075
687
575
0.100
664
621
12
24 164
471
551
22 ACS Paragon Plus Environment
Page 23 of 28
Journal of Agricultural and Food Chemistry
Table 2. Example compounds annotated and confirmed by authentic reference standards in the six red wines. All compounds were detected as [M+H]+ except for gallic acid ([M-H]-).
Name Gallic acid (+)-Catechin (-)-Epicatechin Caffeic acid Syringic acid p-Coumaric acid Ferulic acid N-Acetyl-L-phenylalanine Dihydroquercetin DL-Indole-3-lactic acid Isofraxidin Myricetin 3-(2-Hydroxyethyl)indole trans-Resveratrol Quercetin cis-Resveratrol Phloretin Naringenin
Rt (min) 0.70 1.03 1.14 1.18 1.19 1.43 1.49 1.49 1.52 1.57 1.62 1.68 1.79 1.84 1.98 2.03 2.18 2.22
m/z 169.014 291.086 291.087 181.049 199.060 165.054 195.065 208.096 305.066 206.081 223.060 319.045 162.091 229.086 303.050 229.086 275.091 273.076
Error (mDa) -0.1 -0.1 -0.2 0.1 0.2 0.5 0.7 0.6 -0.1 0.6 0.6 0.0 0.4 0.2 -0.2 0.3 0.2 0.2
Table 3. Comparison of the microLC-MS method to recent alternative methods.
Study This study Pati et al. 201431 Cui et al. 201232 Greco et al. 201210 Mirnaghi et al. 20129 Jaitz et al. 20108
LC time (min) 5.5 40 40 27 20 13
flow rate (mL/min) 0.05 0.2 0.8 0.4 0.2 0.4
mobile phase 100 samples (mL) 27.5 800 3200 1080 400 520
compounds annotated 200+ 26 20 13 8 11
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure graphics Figure 1
24 ACS Paragon Plus Environment
Page 24 of 28
Page 25 of 28
Journal of Agricultural and Food Chemistry
Figure 2
Figure 3
Figure 4
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 5
Figure 6
26 ACS Paragon Plus Environment
Page 26 of 28
Page 27 of 28
Journal of Agricultural and Food Chemistry
Figure 7
27 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Table of Contents Graphics
28 ACS Paragon Plus Environment
Page 28 of 28