Subscriber access provided by FLORIDA INTL UNIV
Article
Selection of Discriminant Marker for Authentication of Asian Palm Civet Coffee (Kopi Luwak): A Metabolomics Approach Udi Jumhawan, Sastia Prama Putri, Yusianto Yusianto, Erly Marwanni, Takeshi Bamba, and Eiichiro Fukusaki J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/jf401819s • Publication Date (Web): 27 Jul 2013 Downloaded from http://pubs.acs.org on August 6, 2013
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 35
Journal of Agricultural and Food Chemistry
1
Selection of Discriminant Markers for Authentication of Asian Palm
2
Civet Coffee (Kopi Luwak): A Metabolomics Approach
3 4
Udi Jumhawan1, Sastia Prama Putri1, Yusianto2, Erly Marwani3, Takeshi
5
Bamba1, Eiichiro Fukusaki*,1
6 7
1
8
University, Osaka, Japan
9
2
Indonesian Coffee and Cocoa Research Institute, East Java, Indonesia
10
3
Bandung Institute of Technology, West Java, Indonesia
Department of Biotechnology, Graduate School of Engineering, Osaka
11 12 13 14
*Address for correspondence:
15
Prof. Dr. Eiichiro Fukusaki
16
Department of Biotechnology
17
Graduate School of Engineering
18
Osaka University
19
2-1 Yamadaoka, Suita, Osaka, 565-0871, Japan
20
Tel: +81-6-6879-7424
21
Fax: +81-6-6879-7424
22
E-mail:
[email protected] ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 35
23
Abstract
24
Kopi Luwak, an exotic Indonesian coffee, is made from coffee berries that
25
have been eaten by the Asian palm civet (Paradoxurus hermaphroditus).
26
Despite being known as world’s most expensive coffee, there is no reliable,
27
standardized
28
multimarker profiling was employed to explore significant metabolites as
29
discriminant markers for authentication. Extracts of 21 coffee beans (Coffea
30
arabica and Coffea canephora) from three cultivation areas were analyzed
31
and subjected to multivariate analyses, principal component analysis, and
32
orthogonal projection to latent structures-discriminant analysis. Citric acid,
33
malic acid, and the inositol-pyroglutamic acid ratio were selected for further
34
verification by evaluating their differentiating abilities against various
35
commercial coffee products. The markers demonstrated potential application
36
in the differentiation of original, fake Kopi Luwak, regular coffee, and coffee
37
blend samples with 50 wt% Kopi Luwak content. This is the first report to
38
address the selection and successful validation of discriminant markers for the
39
authentication of Kopi Luwak.
method
for
determining
its
authenticity.
GC/MS-based
40 41
Keywords: Kopi Luwak, Asian palm civet, GC/MS, authentication, discriminant
42
markers
ACS Paragon Plus Environment
Page 3 of 35
Journal of Agricultural and Food Chemistry
43
Introduction
44
Coffee is one of the most popular beverages in the world. Kopi Luwak is
45
considered to be the world’s most expensive coffee, with a price tag of 300–
46
400 USD per kg.1 Kopi Luwak, the Indonesian words for coffee and civet,
47
respectively, is made from coffee berries that have been eaten by the Asian
48
palm civet (Paradoxurus hermaphroditus), a small mammal native to Southern
49
and Northern Asia.2 The civet climbs coffee trees and instinctively selects
50
coffee cherries. During digestion, the coffee pericarp is completely digested
51
and the beans are excreted. The intact beans are then collected, cleaned,
52
wet-fermented, sun-dried, and further processed by roasting. Kopi Luwak’s
53
high selling price is mainly attributed to its exotic and unexpected production
54
process.3
55
Despite its profitable prospects, there is no reliable, standardized method for
56
determining the authenticity of Kopi Luwak. Moreover, there is limited
57
scientific information on this exotic coffee. Recently, coffee adulterated to
58
resemble Kopi Luwak was reported in the coffee market.4 This poses serious
59
concern among consumers over the authenticity and quality of the products
60
currently available in the market. Discrimination between Kopi Luwak and
61
regular coffee has been achieved using electronic nose data.3 However, the
62
selection of a discriminant marker for authentication was not addressed. The
63
methods currently employed by Kopi Luwak producers is sensory analysis
64
include visual and olfactory testing, both of which are inadequate. For
65
example, visual examination is only possible for green coffee beans prior to
66
roasting, and very few trained experts can perform the highly subjective
67
sensory analysis to discriminate Kopi Luwak.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 35
68
Information flows in metabolic pathways are highly dynamic and represent the
69
current biological states of individual cells. Hence, the metabolome has been
70
considered as the best descriptor of physiological phenomena.5 Metabolomics
71
techniques can be powerful tools to elucidate variations in phenotypes
72
imposed by perturbations such as gene modification, environmental factors,
73
or physical stress. The “black box” process during animal digestion can be
74
translated as physical and enzymatic consequences to the coffee bean, which
75
presents a smoother surface and color changes after digestion.3 Thus, a
76
metabolomics technique was chosen to screen and select discriminant
77
markers for the authenticity assessment of Kopi Luwak. Metabolomics
78
techniques have been effectively applied to distinguish the phytochemical
79
compositions of agricultural products of different origins,6 varieties,6,7 and
80
cultivars8 for quality control and breeding.
81
In
82
spectrometry (GC-Q/MS)-based multimarker profiling was employed to
83
identify discriminant markers for the differentiation of Kopi Luwak and regular
84
coffees. A combination of gas chromatography and mass spectrometry
85
(GC/MS) provides high sensitivity, reproducibility, and the quantitation of a
86
large number of metabolites with a single-step extraction.9,10 Sample
87
classification by means of chemometrics was performed using principal
88
component analysis (PCA). Subsequently, orthogonal projection to latent
89
structures combined with discriminant analysis (OPLS-DA)11 and significance
90
analysis of microarrays/metabolites (SAM)12 were used to isolate statistically
91
significant compounds as discriminant marker candidates. The applicability of
our
study,
gas
chromatography
coupled
with
ACS Paragon Plus Environment
quadrupole
mass
Page 5 of 35
Journal of Agricultural and Food Chemistry
92
these candidates as discriminant markers was verified to differentiate various
93
commercial coffee products.
94 95
Materials and Methods
96
Materials
97
Samples were divided into experimental and validation coffee groups.
98
Experimental coffees were utilized to construct the discrimination model and
99
to select significant compounds. Verification of the applicability of the
100
discriminant markers was carried out using the validation coffee set. In this
101
report, coffee that had been digested by the animal is referred to as Kopi
102
Luwak, and the other beans as regular coffee. Kopi Luwak and regular coffee
103
samples of two species, Coffea arabica (Arabica) and Coffea canephora
104
(Robusta), were used. Coffee samples were obtained from 21 sampling points
105
in three cultivation areas in Indonesia (Java, Sumatra, and Bali). Samples
106
were roasted in a Probat-Werke von Gimborn Maschinenfabrik GmbH model
107
BRZ 2 (Probat, Rhein, Germany) at 205°C for 10 min to obtain a medium
108
degree of roasting, and then were air-cooled for 5 min. Coffee beans were
109
ground and stored in sealed 50 mL BD Falcon tubes at −30°C with light-
110
shielding prior to analysis. The sample descriptions are shown in Tables S1
111
and S2 (Supporting Information). Wet-fermentation was applied to both the
112
Kopi Luwak and regular coffees. For regular coffees, conventional protocols
113
were applied after harvesting, including de-pulping, wet fermentation, and
114
drying.
115
The validation set consisted of authentic Kopi Luwak, commercial Kopi
116
Luwak, commercial regular coffee, fake coffee, and coffee blend. Authentic
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 35
117
Kopi Luwak was produced via controlled processing to ensure the quality of
118
the beans pre- and post-digestion. The remaining samples were purchased
119
commercially. The coffee blend was utilized to examine the feasibility of the
120
method for differentiating mixed and pure coffees.
121 122
Reagents
123
Methanol, chloroform, distilled water, ribitol, and pyridine were purchased
124
from Wako Pure Chemical Industries, Ltd. (Osaka, Japan). Methoxyamine
125
hydrochloride and alkane standard solution were purchased from Sigma
126
Aldrich
127
(trimethylsilyl)trifluoroacetamide (MSTFA) was purchased from GL Science,
128
Inc. (Tokyo, Japan). The six authentic standards of the discriminant markers,
129
their providers, and purities were as follows: citric acid (Nacalai-Tesque,
130
Kyoto, Japan, 99.5%), malic acid (Nacalai-Tesque, 99%), pyroglutamic acid
131
(ICN Biomedicals, Ohio, USA, 99.5%), caffeine (Sigma Aldrich, 98.5%),
132
inositol (Wako, 99%), and glycolic acid (Sigma Aldrich, 99%).
(Milwaukee,
Wisconsin,
USA).
N-Methyl-N-
133 134
Extraction
135
Coffee beans were put into a grinding mill container, cooled for 3 min on water
136
ice cubes, and then ground with a Retsch ball mill (20 Hz, 3 min). Coffee bean
137
powder (15 mg) was transferred into a 2 mL Eppendorf tube. In addition to
138
pure samples, a 50:50 (wt%) blend of Kopi Luwak and regular coffee was
139
used. One milliliter of a single-phase extraction solvent consisting of 2.5/1/1
140
(v/v/v) methanol, distilled water, and chloroform, respectively, was added to
141
extract a wide range of metabolites. A non-specific extraction procedure was
ACS Paragon Plus Environment
Page 7 of 35
Journal of Agricultural and Food Chemistry
142
applied to avoid limiting the target analysis to specific compounds and to
143
comprehensively screen the components of Kopi Luwak. As an internal
144
standard, ribitol (60 µL, diluted with deionized water to 0.2 mg/mL) was
145
utilized. The mixture was shaken for 1 min and then centrifuged at 4°C and
146
16000 g for 3 min. The supernatant (900 µL) was transferred into a 1.5 mL
147
Eppendorf tube and diluted with 400 µL Milli-Q water (Wako). The mixture
148
was then vortexed and centrifuged for 3 min. A 400 µL portion of the aqueous
149
phase was transferred into a fresh 1.5 mL Eppendorf tube with a screw cap.
150
The solvent was removed by vacuum centrifugation for 2 h, followed by
151
freeze-drying overnight. All samples were analyzed in triplicate (n = 3).
152 153
Derivatization for GC/MS analysis
154
Oximation and trimethylsilylation were used for derivatization. Methoxyamine
155
hydrochloride (100 µL, 20 mg/mL in pyridine) was added to the dried extract
156
as the first derivatization agent. The mixture was incubated at 30°C for 90
157
min. After addition of the second derivatization agent, N-methyl-N-
158
(trimethylsilyl)trifluoroacetamide (MSTFA, 50 µL), the mixture was incubated
159
at 37°C for 30 min.
160 161
GC/MS Analysis
162
GC-Q/MS analysis was performed on a GCMS-QP 2010 Ultra (Shimadzu,
163
Kyoto, Japan) equipped with a CP-SIL 8 CB low bleed column (0.25 mm x 30
164
m, 0.25 µm, Varian Inc., Palo Alto, California, USA) and an AOC-20i/s
165
(Shimadzu) as an autosampler. The mass spectrometer was tuned and
166
calibrated prior to analysis. The derivatized sample (1 µL) was injected in split
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
167
mode, 25/1 (v/v), with an injection temperature of 230°C. The carrier gas (He)
168
flow was 1.12 mL/min with a linear velocity of 39 cm/s. The column
169
temperature was held at 80°C for 2 min, increased by 15°C/min to 330°C, and
170
then held for 6 min. The transfer line and ion source temperatures were 250
171
and 200°C, respectively. Ions were generated by electron ionization (EI) at
172
0.93 kV. Spectra were recorded at 10000 u/s over the mass range 85−500
173
m/z. A standard alkane mixture (C8–C40) was injected at the beginning and
174
end of the analysis for tentative identification.
175 176
Identification and quantitation of marker candidates
177
The discriminant marker candidates were identified and quantitated against
178
six authentic standards (malic acid, citric acid, glycolic acid, pyroglutamic acid,
179
caffeine, and inositol) at various concentrations. The final concentrations of
180
the authentic standards were adjusted to 1, 10, 50, 100, 250, 500, 750, 1000,
181
1500, and 2000 µM with the extraction solvent to construct a calibration curve.
182
For extraction, the authentic standards were processed identically to the
183
coffee bean samples. The standards were co-injected during sample analysis.
184
Two blank solutions were prepared by adding only extraction solvent and
185
distilled water, respectively. The limits of detection (LOD) and quantitation
186
(LOQ) were determined via known protocols.13,14 The construction of the
187
standard curve and quantitation were conducted using GC/MS Solution
188
software (Shimadzu). No authentic standards were detected in either of the
189
blank samples.
190 191
Data Processing
ACS Paragon Plus Environment
Page 8 of 35
Page 9 of 35
Journal of Agricultural and Food Chemistry
192
Chromatographic data were converted into ANDI files (Analytical Data
193
Interchange Protocol, *.cdf) using the GC/MS Solution software package
194
(Shimadzu). Peak detection, baseline correction, and peak alignment of
195
retention times were performed on the ANDI files using the freely available
196
software package, MetAlign.15 Spectra were normalized manually by adjusting
197
the peak intensity against the ribitol internal standard.
198
Retention indexes of the eluted compounds were calculated based on the
199
standard alkane mixture. By comparing the retention indexes and unique
200
mass spectra with our in-house reference library, tentative identifications were
201
obtained. To simplify and accelerate the tentative identifications of
202
compounds that were registered in the in-house library database, AIoutput2
203
(version 1.29) annotation software, developed in the authors’ laboratory, was
204
utilized.16 For comparison with the National Institute of Standards and
205
Technology (NIST) library, retention times were used instead.
206 207
Multivariate data analysis
208
PCA and OPLS-DA were performed. OPLS-DA with an S alphabet-like plot, or
209
S-plot, was chosen to isolate and select statistically significant and potentially
210
biochemically interesting compounds. The S-plot provides covariance and
211
correlation between metabolites and the modeled class to allow easier
212
visualization. The variables that changed significantly are plotted at the top
213
and bottom of the S-plot, and those that do not significantly contribute are
214
plotted in the middle.11 A sevenfold cross validation was carried out to assess
215
the accuracy of the discrimination model in practice. The goodness-of-fit (R2)
216
and predictability (Q2) parameters were then determined. Analysis was
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
217
performed with commercial software, SIMCA-P+ version 12 (Umetrics, Umeå,
218
Sweden). Data were Pareto-scaled to reduce the effect of noise in the
219
chromatograms.
220
To confirm the selection of significant compounds by OPLS-DA, data were
221
also subjected to significance analysis of microarrays/metabolites (SAM)
222
using MetaboAnalyst 2.0. At first, SAM was projected to microarrays analysis
223
to assign significance genes on the basis of changes in their expression.12
224
Recently, MetaboAnalyst 2.0, a web server for metabolomics analysis, was
225
developed for the application of SAM to metabolomics data.17,18
226
The “relative differences,” d(i), i.e., the differences in intensity for each peak,
227
were calculated by a formula described elsewhere.12 To generate the control
228
set, data were analyzed by n-balanced permutations. The relative differences
229
of the permutated data, dp(i), were determined. Next, the “expected relative
230
differences”, dE(i), were calculated as the average of dp(i) over n-balanced
231
permutations. Scatter plots of the expected relative differences (dE(i)) and
232
relative differences (d(i)) were constructed. Insignificant compounds were
233
identified as those which satisfied d(i)
234
assigned as a certain distance displaced from d(i) = dE(i); compounds at
235
distances larger than the threshold value were considered significant. The
236
higher the threshold value, the lower the false discovery rate (FDR, the
237
percentage of falsely significant compounds) and the number of significant
238
compounds that can be obtained. Univariate statistics, performed by the freely
239
available software, R project,19 was applied to compare the means of selected
240
significant compounds using the Student’s t-test. Box plots were also
dE(i). A threshold value (∆) was
ACS Paragon Plus Environment
Page 10 of 35
Page 11 of 35
Journal of Agricultural and Food Chemistry
241
constructed to display the differentiation among coffee samples. Differences
242
were considered significant when p < 0.05.
243 244
Results and Discussion
245
GC/MS-based metabolite profiling of Kopi Luwak
246
GC-Q/MS analysis was performed on aqueous extracts of coffee beans to
247
investigate the differences in their metabolite profiles and select discriminant
248
markers for robust authentication. In addition, this research focused on
249
increasing the scientific information about Kopi Luwak. A quadrupole mass
250
spectrometer was selected because of its availability as the most widely used
251
mass analyzer. However, a conventional Q/MS can be operated only at a
252
slow scan rate.20 With processor improvements and high-speed data
253
processing, newly developed GC-Q/MS instruments provide increased
254
sensitivity at high scan speeds of up to 10.000 u/s.21
255
Because of their broad cultivation areas and commercial profitability, Coffea
256
arabica and Coffea canephora, which represent 65% and 35% of the total
257
annual coffee trade, respectively, were utilized for metabolomics analysis.22 A
258
total of 182 peaks from 21 coffee beans were extracted using MetAlign.
259
Moreover, 26 compounds were tentatively identified by comparison with our
260
in-house library (by retention index) and the NIST library (by retention time);
261
six of these were identified by co-injection with an authentic standard.
262
Tentatively identified components consisted of organic acids, sugars, amino
263
acids, and other compounds (Table S3). Previously reported coffee bean
264
constituents, including chlorogenic, quinic, succinic, citric, and malic acids;
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
265
caffeine, one of the compounds supplying bitter taste in coffee; and sucrose,
266
the most abundant simple carbohydrate, were identified.23–27
267
In recent research, unsupervised analysis, PCA, has been employed for data
268
exploration and to visualize information based on sample variance.28,29 A PCA
269
score plot derived from the 21 coffee beans differentiated two data groups
270
based on their species, Arabica and Robusta (Figure 1), and resulted in a
271
goodness-of-fit parameter (R2) of 0.844. Caffeine and quinic acid were
272
significant for the Robusta coffee data sets, whereas the Arabica data set was
273
mainly supported by various organic acids such as malic, chlorogenic, citric,
274
and succinic acids. The data differentiation was explained by 42.9% of
275
variance along PC1. The results indicated that genetic diversity more strongly
276
influenced the data separation than animal perturbation.
277
Because of the large variance among coffee species, sample differentiation
278
based on the type of coffee, Kopi Luwak or regular coffee, could not be
279
observed. Additional analyses were carried out independently for each coffee
280
species originating from the same cultivation area. The PCA score plot
281
revealed data separation based on the type of coffee, in which Kopi Luwak
282
and regular coffee could be clearly separated (Figure S1). For the Arabica
283
coffee data set, the separation was explained by 45.5% and 27.7% variances
284
in PC1 and PC2, respectively. By PC2, Kopi Luwak was closely clustered in
285
the same region, whereas regular coffees tended to separate based on their
286
cultivation areas. In the loading plot, malic and glycolic acids contributed
287
highly to the Kopi Luwak data (data not shown). Thereby, coffee beans may
288
possess similar profiles after animal digestion. Differences in cultivation areas
289
were considered to have the least significance for data separation. In Robusta
ACS Paragon Plus Environment
Page 12 of 35
Page 13 of 35
Journal of Agricultural and Food Chemistry
290
coffee, a clear separation between Kopi Luwak and regular coffee was
291
observed, which was explained by 79.1% variance of PC1. Significant
292
compounds for separation, including inositol and pyroglutamic acid for Kopi
293
Luwak and quinic acid for regular coffee, were observed.
294 295
Discriminant analysis to select candidates for discriminant markers
296
An overview of all data samples was provided by the unsupervised analysis,
297
PCA. However, detailed information regarding compounds contributing to the
298
data differentiation between Kopi Luwak and regular coffee remained unclear.
299
Therefore, coffee bean data sets were subjected to supervised discriminant
300
analysis (OPLS-DA). For analyses having two or more classes, OPLS-DA is
301
the most suitable platform for isolating and selecting differentiation markers.
302
Compounds with reliable, high contributions to the model may possess
303
potentially biochemically interesting characteristics; thus, they can be selected
304
as biomarker candidates.11 All OPLS-DA models exhibited R2 and Q2 values
305
greater than 0.8, which would be categorized as excellent.30 In addition, all
306
models were in the range of validity after permutation tests using 200
307
variables (data not shown). The model was considered valid after permutation
308
for those that met the following criteria: R2Y-intercepts and Q2-intercepts
309
which did not exceed 0.3–0.4 and 0.05, respectively.31
310
Potential candidates for discriminant markers can be selected via S-plots by
311
setting the cut-off for covariance, p[1], and the correlation value, p[corr], to >
312
|0.2|. S-plots of the coffee data sets are shown in Figures 2B and D. In
313
addition to cut-off values, candidates for discriminant markers were selected
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
314
by variable importance in projection values (VIP). Large VIP values (> 1) are
315
more relevant for model construction.
316
The OPLS-DA score plot of Arabica coffee data sets is shown in Figure 2A.
317
Discrimination between Kopi Luwak and regular coffee was obtained. The
318
model was evaluated with R2 and Q2 values of 0.965 and 0.892, respectively.
319
Interestingly, compounds that were uncorrelated with Kopi Luwak were quinic
320
acid, caffeine, and caffeic acid. These compounds have been reported as
321
contributors of bitterness as well as acidity in coffee.23–26 In contrast,
322
compounds that were predictive to Kopi Luwak, i.e., over the cut-off value,
323
included citric, malic, and glycolic acids. The OPLS-DA score plot of the
324
Robusta coffee data sets (Figure 2C) was explained by R2 and Q2 values of
325
0.957 and 0.818, respectively. Caffeine, one of the bitter principles in coffee,
326
was found to be significantly correlated with Robusta Kopi Luwak data sets.
327
Robusta coffee has been reported to contain higher amounts of caffeine than
328
Arabica. Thus, it tends to be bitter and flavorless, whereas Arabica coffee is
329
considered to be milder, contain more aromatic compounds, and is more
330
appreciated by the consumer.32
331
We employed significance analysis of microarrays/metabolites (SAM) to
332
select significant compounds as discriminant markers, as a comparison to
333
OPLS-DA. In general, we used the same data sets, applying different
334
multivariate data analysis. By assigning the lowest possible FDR value, a total
335
of 12 compounds was considered as significant in the Arabica coffee data set
336
(Figure S2). Of these, citric acid, glycolic acid, malic acid, quinic acid, and
337
other unidentified peaks were included. In the Robusta coffee data set, nine
338
significant compounds were found (Figure S2); inositol, caffeine, pyroglutamic
ACS Paragon Plus Environment
Page 14 of 35
Page 15 of 35
Journal of Agricultural and Food Chemistry
339
acid, and six unidentified peaks exhibited high significance by SAM.
340
Candidates for discriminant markers for the authentication assessment of
341
Arabica and Robusta coffees are listed in Table 1. The selected marker
342
candidates met significant criteria in both OPLS-DA and SAM. Discriminant
343
markers were chosen independently for the Arabica and Robusta coffee.
344
To confirm whether these selected markers were generated as a result of
345
animal digestion, we investigated cause-effect relationships by quantitating
346
the discriminant marker candidates in green and roasted coffee beans from
347
controlled processing, pre- and post- animal digestion (Samples No. 5 and
348
No. 11 experimental sets). The results are displayed in Figure S3. In both the
349
raw and roasted beans, citric acid was present in higher concentration after
350
animal digestion, exhibiting a significant value difference (p < 0.05) between
351
Kopi Luwak and regular coffee. The concentration of caffeine was also
352
increased after digestion, but the difference was insignificant (p > 0.05). As a
353
result of roasting, the glycolic acid concentration increased dramatically (p
|0.2|. Triangles indicate peaks detected by GC/MS. Figure 3. PCA score plot of validation coffee set. Data plotted inside dashed circle indicates authentic Kopi Luwak. Figure 4. Box plots of peak intensity and concentrations of selected discriminant markers for validation.
ACS Paragon Plus Environment
Page 28 of 35
Page 29 of 35
Journal of Agricultural and Food Chemistry
Tables
Table 1. Candidates for discriminant markers from OPLS-DA and SAM, and analytical parameters for quantitation Discriminant marker
Glycolic acid Malic acid Pyroglutamic acid Citric acid Caffeine Inositol
RT (min)
4.96 9.05 9.43 11.61 12.18 13.45
VIP
3.93 5.53 1.7 5.6 2.28 4.47
RSD [%] (n = 3) RT Areaa 0.12 0.05 0.05 0.04 0.04 0.03
1.87 2.29 3.36 3.29 3.51 5.09
Linearity R2 0.9999 0.9996 0.9992 0.9997 0.9961 0.9974
Range (µ µM) 1-1000 1-1000 1-750 1-1000 100-2000 1-1000
a
at 100 µM
ACS Paragon Plus Environment
LOD (mg/kg)
LOQ (mg/kg)
0.021 0.043 0.054 0.504 1.531 0.082
0.066 0.132 0.164 1.526 4.638 0.247
Journal of Agricultural and Food Chemistry
Figure graphics
Figure 1
ACS Paragon Plus Environment
Page 30 of 35
Page 31 of 35
Journal of Agricultural and Food Chemistry
B
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
D
Figure 2
ACS Paragon Plus Environment
Page 32 of 35
Page 33 of 35
Journal of Agricultural and Food Chemistry
Figure 3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 4
ACS Paragon Plus Environment
Page 34 of 35
Page 35 of 35
Journal of Agricultural and Food Chemistry
Graphic for table of contents (TOC)
ACS Paragon Plus Environment