Ultrafast Polyphenol Metabolomics of Red Wines ... - ACS Publications

Dec 23, 2015 - and Oliver Fiehn*,†,#. †. UC Davis Genome Center − Metabolomics, University of California, Davis, California 95616, United States...
0 downloads 0 Views 1MB Size
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