Subscriber access provided by LAURENTIAN UNIV
Article
Application of targeted metabolomics to investigate optimum growing conditions to enhance bioactive content of strawberry Ikram Akhatou, Ana Sayago, Raúl González-Domínguez, and Angeles Fernández-Recamales J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b03701 • Publication Date (Web): 11 Oct 2017 Downloaded from http://pubs.acs.org on October 12, 2017
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
Application of Targeted Metabolomics to Investigate Optimum Growing
2
Conditions to Enhance Bioactive Content of Strawberry
3 4
Ikram Akhatou†, ‡, Ana Sayago†, ‡, Raúl González-Domínguez†,
5
Fernández-Recamales†, ‡ *
‡ *
, Ángeles
6 7
†
8
21007 Huelva, Spain. ‡International Campus of Excellence CeiA3, University of
9
Huelva, 21007 Huelva, Spain.
Department of Chemistry, Faculty of Experimental Sciences, University of Huelva,
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
26
Abstract
27
A simple, sensitive and rapid assay based on liquid chromatography coupled to tandem
28
mass spectrometry was designed for simultaneous quantitation of secondary metabolites
29
in order to investigate the influence of variety and agronomic conditions on the
30
biosynthesis of bioactive compounds in strawberry. For this purpose, strawberries
31
belonging to three varieties with different sensitivity to environmental conditions
32
('Camarosa', 'Festival', 'Palomar') were grown in soilless system under multiple
33
agronomic conditions (electrical conductivity, substrate type and coverage). Targeted
34
metabolomic analysis of polyphenolic compounds, combined with advanced
35
chemometric methods based on learning machines, revealed significant differences in
36
multiple bioactives, such as chlorogenic acid, ellagic acid rhamnoside, sanguiin H10,
37
quercetin 3-O-glucuronide, catechin, procyanidin B2, pelargonidin 3-O-glucoside,
38
cyaniding 3-O-glucoside and pelargonidin 3-O-rutinoside, which play a pivotal role in
39
organoleptic properties and beneficial healthy effects of these polyphenol-rich foods.
40 41 42
Keywords: Phenolic compounds, targeted metabolomics, strawberry, variety,
43
agronomic practices, UHPLC-ESI-MS/MS
44 45 46 47 48 49 50
2 ACS Paragon Plus Environment
Page 2 of 35
Page 3 of 35
Journal of Agricultural and Food Chemistry
51
Introduction
52
Phenolic compounds, also called polyphenols, constitute a miscellaneous group of
53
secondary metabolites widely distributed in plant kingdom. Because of their structure,
54
containing aromatic rings and hydroxyl groups, and their chemical properties as
55
hydrogen donors and singlet oxygen quenchers,1 polyphenols are commonly involved in
56
defense mechanisms by acting as mediators of plant response against biotic and abiotic
57
stresses, including those related to drought, salinity, low temperatures, heavy metal
58
exposure or UV radiation. Some phenolic compounds (e.g. anthocyanins) are also
59
responsible for distinct plant characteristics such as color and smell of flowers and
60
fruits, attractants for pollinators and seed-dispersing animals. On the other hand, several
61
phenolics exhibit broad-spectrum antimicrobial activity, thus protecting plants against
62
viruses, parasites and other externally damaging agents; whereas others possess
63
phytoalexin properties, and their synthesis can be induced in response to wounding,
64
feeding by herbivores or infection by pathogens. Thereby, it has been demonstrated that
65
proanthocyanidins can contribute to defense and stress resistance, while ellagic acid and
66
their derivatives play an important role by protecting plants from predators and
67
regulating their growth.2 Since phenolic compounds are ubiquitous in all plants, they are
68
an integral part of the human diet through the consumption of edible plants and plant-
69
derived products, including fresh and cooked vegetables and fruits, legumes, spices and
70
beverages such as fruit juices, wine, tea, coffee and other infusions. In recent years, the
71
interest for these bioactives has largely increased due to their health benefits and their
72
impact on food quality. Thus, numerous efforts have been made in order to characterize
73
polyphenol-rich foods and to study the influence of various factors such as variety,
74
pedoclimatic conditions and geographical origin on the production of these bioactives,
75
thus serving as a very valuable tool for increasing nutritional quality of plant-derived
3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
76
foods in breeding programs, as well as for assessing the authenticity and traceability of
77
these foods.3,4
78
Strawberries (Fragaria x ananassa Duch.) are one of the most economically important
79
and widely cultivated fruit crops across the world. These summer fruits are very popular
80
among consumers due to their attractive red color and highly desirable taste and flavor.
81
Furthermore, epidemiological studies have revealed that consumption of berries is
82
associated with the prevention and improvement of chronic-degenerative diseases.5 In
83
this line, a significant number of in vivo and in vitro studies have demonstrated the
84
effect of strawberry consumption on the regulation of oxidative stress,6 inflammation,7
85
obesity,8 diabetes,9 cancer10 and cardiovascular diseases.11 These protective effects
86
appear to be associated with the high content of bioactive compounds, including both
87
micronutrients and phytochemicals.3 In particular, strawberries are a very rich source of
88
anthocyanins, flavan-3-ols, flavonols, as well as ellagic acid and derivatives, which are
89
well known for their antioxidant, cardioprotective and anti-inflammatory properties.12
90
However, the concentration of polyphenols in strawberry significantly vary between
91
genotypes,13-17 and many other factors related to growth conditions and agricultural
92
practices may contribute to this phenolic profile. Moreover, their synthesis can also
93
change in response to biotic and abiotic stress.18-20 Therefore, breeding programs might
94
be used to develop strawberry cultivars with better organoleptic characteristics and
95
enhanced production of health-related metabolites.
96
Based on the above-reported considerations, this study was aimed to investigate the
97
polyphenolic profile of three strawberry cultivars with different sensitivity to
98
environmental conditions, which were grown in soilless systems under different
99
agronomic conditions (electrical conductivity, substrate type and coverage). A targeted
100
metabolomic approach based on ultra-high performance liquid chromatography coupled
4 ACS Paragon Plus Environment
Page 4 of 35
Page 5 of 35
Journal of Agricultural and Food Chemistry
101
to tandem mass spectrometry (UHPLC-MS/MS) was applied to identify and quantitate
102
the main polyphenol compounds in strawberry fruits. Then, partial least squares
103
discriminant analysis and random forest methods were employed to explore the
104
differences in metabolite profiles between strawberry cultivars, and to assess the
105
influence of agronomic practices on polyphenol production.
106
Materials and Methods
107
Chemicals and standards
108
Acetonitrile, methanol and glacial acetic acid (LC/MS grade) were purchased from
109
Fisher Scientific (Leicestershire, UK), while the ultrapure water was obtained from a
110
Milli-Q water purification system (Millipore, Bedford, MA). Pure standards of gallic
111
acid, protocatechuic acid, sinapic acid, (-) epicatechin gallate, ellagic acid, chlorogenic
112
acid, (-) epicatechin, ferulic acid and quercetin were acquired from Sigma-Aldrich
113
(Steinheim, Germany). Vanillic acid, vanillin, caffeic acid, p-hydroxyphenylacetic acid,
114
cinnamic acid, p-hydroxybenzoic acid, m-coumaric acid, p-coumaric acid and tyrosol
115
were from Fluka (Steinheim, Germany). Apigenin, luteolin, pelargonidin 3-O-glucoside,
116
cyaniding 3-O-glucoside, (+) catechin, procyanidin B1, procyanidin B2, quercetin 3-O-
117
glucoside,
118
hesperidin and hyperoside were purchased from Extrasynthèse (Genay, France).
119
Experimental design and sampling
120
Strawberries (Fragaria x ananassa Duch.) were cultivated in soilless systems, in
121
plantations managed by the University of Huelva (southwest Spain, latitude 37° 14ʹ N,
122
longitude 6° 53ʹ W, altitude 23 m). In this study, three cultivars were investigated
123
('Palomar', 'Festival' and 'Camarosa'), which were grown under different electrical
124
conductivities (1, 2 and 3 dS/m), two macrotunnel types (covered and uncovered) and
125
three substrates (coconut fiber, perlite and rockwool). Detailed description about the
kaempferol
3-O-glucoside,
isorhamnetin
3-O-glucoside,
5 ACS Paragon Plus Environment
quercetrin,
Journal of Agricultural and Food Chemistry
126
experimental design can be found in our previous works.4,17 For each study group,
127
several fruits were collected and pooled, and then washed and homogenized using a
128
kitchen mixer. These samples were stored for up to 2 months at -21 °C until further
129
analysis.
130
Sample preparation
131
Homogenized fruits (5.0 g) were mixed with 10 mL of methanol, then submitted to
132
sonication for 15 min and finally centrifuged at 10000 rpm during 10 min at 4 °C.
133
Supernatants were taken to dryness at 40 °C by using a rotary evaporator, and residues
134
were redissolved in 3 mL of 50% methanol (v/v). The concentrated extracts were
135
filtered through 0.20 µm nylon filter prior to injection.
136
Targeted metabolomic analysis
137
An Agilent 1200 series ultra-high performance liquid chromatography system (Agilent
138
Technologies, Santa Clara, CA), equipped with binary pump system, vacuum degasser,
139
cooled autosampler and thermostated column compartment, was employed for targeted
140
metabolomic analysis of strawberry extracts. The column used was a 50 mm × 2.1 mm
141
i.d., 1.8 µm, Zorbax SB-C18 (Agilent Technologies, Santa Clara, CA), thermostated at
142
30 °C. Optimum separation was achieved using a binary gradient delivered at 0.4
143
mL/min flow rate, consisting of 0.2% acetic acid in water (v/v) at pH 3.10 (solvent A),
144
and pure acetonitrile (solvent B). The gradient elution program was: 0-3 min, 5% B; 3-
145
15 min, 5-40% B; 15-15.5 min, 40-100% B; 15.5-17 min, 100% B isocratic; and finally
146
returning to initial conditions in 6 min. The injection volume was 5 µL, and the
147
autosampler was maintained at 4 °C. The eluent from the column was directly
148
introduced into a model 6410 triple quadrupole mass spectrometer equipped with
149
electrospray ionization source (ESI), and controlled by MassHunter Workstation
150
Software (Agilent Technologies, Santa Clara, CA). Source working conditions were:
6 ACS Paragon Plus Environment
Page 6 of 35
Page 7 of 35
Journal of Agricultural and Food Chemistry
151
capillary voltage 4000 V, gas flow rate 10 L/min, gas temperature 300 °C and nebulizer
152
pressure 35 psi. Phenolic acids, flavonoids, ellagitannins and ellagic acid derivates were
153
detected in the negative ion mode, while anthocyanins were analyzed under positive ion
154
mode. Detection was performed in dynamic multiple reaction monitoring (dMRM)
155
mode using retention times and detection windows, which improves chromatographic
156
peaks and provides increased reproducibility and accuracy of quantitation. For
157
maximum sensitivity and specificity, instrumental parameters were empirically
158
determined to optimize MRM transitions (two for each analyte). The most sensitive
159
transition was used for quantitation purposes, and the other one as qualifier transition.
160
Identification and quantitation of phenolic compounds was performed using commercial
161
standards, when possible. On the contrary, all phenolic derivatives for which
162
commercial standards were not available were identified according to previously
163
published literature.21 Briefly, ellagic acid derivatives presented characteristic fragments
164
at m/z 301 and 207 Da; MS/MS spectra of quercetin derivatives were dominated by
165
daughter ions of m/z 301, 179 and 151 Da; while fragmentation of kaempferol and
166
pelargonidin derivatives yielded intense signals at m/z 285 and 271 Da, respectively. In
167
accordance, the estimated concentrations of ellagic acid derivatives were expressed as
168
ellagic acid equivalents, quercetin derivatives as quercetin equivalents, kaempferol
169
derivatives as kaempferol 3-O-glucoside equivalents, and pelargonidin derivatives as
170
pelargondin 3-O-glucoside equivalents.
171
Analytical quality control
172
The proposed methodology was validated according to UNE 82009-1:1998 normative.
173
The within-laboratory repeatability and reproducibility were assessed by analyzing an
174
extract six times within the same day and over a period of 1 month in duplicate,
175
respectively. Calibration curves were built at seven concentration levels (0.5-40 µg/mL
7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
176
for anthocyanins; 0.05-8.0 µg/mL for the rest of phenolic compounds) by diluting the
177
stock solutions (1000 µg/mL in methanol). Recoveries were determined in pooled
178
samples by spiking with pure standards at three concentration levels, and three
179
replicates from each sample were analyzed.
180
Statistical Analysis
181
Metabolomic data were subjected to supervised classification methods in order to
182
investigate the sample classification in predefined groups and to find the variables with
183
the highest discriminant power. As a final step, the Kruskal-Wallis test was performed
184
on previously selected metabolites in order to validate these results. Statistical treatment
185
of data was performed by using Statistica 8.0 (StatSoft, Tulsa, UK), SIMCA-P v11.5
186
(UMETRICS, Umeå, Sweden) and R statistical software 2009 (R Foundation for
187
Statistical Computing, Vienna, Austria) packages. Two data mining methods, partial
188
least square discriminant analysis (PLS-DA) and random forest (RF), were applied for
189
building predictive models, which were subsequently compared for evaluating their
190
classification and prediction abilities during cross-validation procedure. In addition,
191
sensitivity and specificity of these models were also evaluated using confusion matrices.
192
Results and Discussion
193
Metabolite Profiling Analysis
194
Phenolic profiles of three strawberry cultivars, grown using multiple agronomic
195
practices, were evaluated by liquid chromatography-mass spectrometry in MRM
196
scanning mode. Twenty-nine different compounds were detected and quantitated, which
197
can be classified into different groups of polyphenols: a) phenolic acids (gallic, vanillic,
198
protocatechuic, p-hydroxybenzoic, chlorogenic, cinnamic, caffeic, m-coumaric and p-
199
coumaric acids); b) ellagic acid and derivatives (HHDP galloyl glucose, ellagic acid
200
pentoside, ellagic acid rhamnoside, and sanguiin H10); c) flavonols (quercetin,
8 ACS Paragon Plus Environment
Page 8 of 35
Page 9 of 35
Journal of Agricultural and Food Chemistry
201
quercetin 3-O-glucoside, quercetin 3-O-galactoside, quercetin 3-O-glucuronide,
202
kaempferol 3-O-glucoside and kaempferol 3-O-acetylhexoside); d) flavones (apigenin);
203
e) flavan-3-ols (catechin, epicatechin, epicatechin gallate, procyanidins B1 and B2); f)
204
anthocyanins (pelargonidin 3-O-glucoside, cyaniding 3-O-glucoside, pelargonidin 3-O-
205
rutinoside). The structures of major phenolic compounds detected in strawberries are
206
shown in Figure 1. Relative standard deviation obtained for the within-laboratory
207
repeatability and reproducibility of the proposed methodology were below 7.3% and
208
10.4%, respectively. The higher RSD for the repeatability study corresponded to
209
chlorogenic acid, whereas the worst reproducibility was showed by quercetin 3-O-
210
glucoside. The method here developed presented a good linearity in the studied
211
concentration range (0.5-40 µg/mL for anthocyanins; 0.05-8.0 µg/mL for the rest of
212
phenolic compounds), with correlation coefficients higher than 0.99 for all the
213
compounds. The limits of detection and quantitation were calculated as the
214
concentration necessary to produce a chromatographic signal three and ten times higher
215
than the baseline noise, respectively. Results obtained were satisfactory to quantitate
216
target compounds at the concentrations found in strawberry samples. Finally, recoveries
217
ranged from 54-121% for phenolic acids, from 59-106% for flavonols, from 67-128%
218
for flavan-3-ols, and from 75-95% for anthocyanins. Recovery for apigenin was 75%.
219
Table 1 shows the mean, minimum and maximum concentration values for the
220
polyphenolic compounds detected in the three strawberry cultivars. As can be observed,
221
each cultivar showed a characteristic phenolic profile, 'Camarosa' being the strawberry
222
variety with greatest total phenolic content mainly due to the high amounts of ellagic
223
acid and its derivatives, as well as anthocyanins, which accounted for 44.5% and 35.7%
224
of the total polyphenols, respectively. It should be noted that these results agreed with
225
our previous study findings3 and with results reported by other authors.22 In this context,
9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
226
pelargonidin 3-O-glucoside was the major anthocyanin found in the three cultivars, in
227
agreement with previously published results.23 On the other hand, ellagic acid and
228
derivatives were the most concentrated phenolic acids in strawberries, in consonance
229
with other author’s findings.24-25 Particularly, conjugated ellagic acid rhamnoside was
230
the most important metabolite in terms of quantity in the three cultivars (141.43, 123.40
231
and 107.87 µg/g of fresh fruit, for 'Camarosa', 'Festival' and 'Palomar', respectively).
232
Flavonols and flavan-3-ols were also abundant in strawberries, and accounted for 9.0
233
and 8.0% of the total phenolic content, respectively. In absolute terms, the higher levels
234
of flavonols were found in cv. 'Camarosa' (46.5 µg/g of fresh fruit,), followed by
235
'Palomar' (38.9 µg/g of fresh fruit,) and 'Festival' (37.4 µg/g of fresh fruit,). Quercetin 3-
236
O-glucuronide was the most concentrated flavonol in the three surveyed strawberry
237
cultivars, with a mean value of 37.7 µg/g. In this line, glucosides and glucuronides of
238
quercetin and kaempferol have previously been identified as the major flavonols in
239
strawberry.25-26 Regarding the family of flavanols, the higher levels were found in cv.
240
'Festival' (10.23 µg/g of fresh fruit,) and the lower in cv. 'Palomar' (5.8 µg/g of fresh
241
fruit,). Procyanidin B1 was the predominant flavanol in the three studied strawberry
242
cultivars (mean value 22.6 µg/g), followed by catechin (14.03 µg/g), results that agree
243
with those reported by other authors.27-28 Finally, it is noteworthy that the total content
244
of phenolic acids in the three strawberry cultivars was very similar, except for
245
chlorogenic acid, ranging from 2.2 µg/g (cv. 'Festival') to 3.2 µg/g of fresh fruit (cv.
246
'Palomar'). On the contrary, the concentration of chlorogenic acid was highly variable,
247
'Festival' being the cultivar with greater content (12.7 µg/g) followed by 'Palomar' (5.8
248
µg/g) and 'Camarosa' (3.5 µg/g). In this sense, previous studies conducted on
249
strawberries have also reported variable composition and contents of phenolic acids.28-31
10 ACS Paragon Plus Environment
Page 10 of 35
Page 11 of 35
Journal of Agricultural and Food Chemistry
250
To conclude, the flavone apigenin was detected in the three strawberry varieties, with a
251
mean content of 0.065 µg/g.
252
Chemometric discrimination of strawberry cultivars based on PLS-DA
253
Partial least square discriminant analysis (PLS-DA) was applied to identify potential
254
markers for cultivar differentiation. This supervised multivariate statistical tool is based
255
on the search of latent variables (or components), able to discriminate between
256
previously defined categories. In order to optimize the number of these components data
257
are usually subjected to cross-validation, so that the statistical performance of the model
258
can be defined by the following parameters: RY2, the proportion of variance of the
259
response variable explained by the model; RX2, the proportion of variance in the data
260
explained by the model; and Q2, the predictive ability of the model.
261
Before performing statistical analyses, data were submitted to Pareto scaling to reduce
262
the relative importance of larger values. Then, three binary classification models were
263
built (model A, 'Camarosa' vs. 'Festival'; model B, 'Camarosa' vs. 'Palomar'; model C,
264
'Festival' vs. 'Palomar'), and their quality was evaluated using seven-fold cross
265
validation. In all cases, a two-component model was obtained with good quality of fit
266
(R2X = 0.645, 0.702 and 0.659, for model A, B and C, respectively), validity (R2Y =
267
0.655, 0.751 and 0.733, for model A, B and C, respectively) and prediction (Q2 = 0.527,
268
0.663 and 0.664, for model A, B and C, respectively), which resulted in a satisfactory
269
separation of sample classes in the corresponding scores plots (Figure 2A-C). Then,
270
discriminant metabolites were selected according to the variable importance for the
271
projection parameter (VIP), which indicates the importance of the variable in the model.
272
For this, phenolic metabolites with VIP values higher than 1.0 were picked and
273
compared by using the Kruskal-Wallis test. These results are summarized in box-plots
274
represented in Figure 3. As can be observed, 'Camarosa' presented significantly higher
11 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
275
levels of anthocyanins, ellagic acid rhamnoside, quercetin 3-O-glucuronide and
276
procyanidin B2, while cv. 'Festival' was richer in catechin and chlorogenic acid, and cv.
277
'Palomar' in sanguiin H10.
278
Numerous evidences emphasize the genotype as a pivotal player in the compositional
279
variability in primary and secondary metabolites. Their effects on the nutritional and
280
sensorial quality of plant derived foods are well known,13-17 and several breeding and
281
biotechnological programs are focused on them.32 In this context, various authors have
282
previously studied the influence of variety in the content of phytochemicals in
283
strawberry.33,34 Other works focused on investigating metabolic changes in strawberry
284
during ripening, as well as the differences between organic crops, conventional crops,
285
and breeding lines with wild genotypes.35,36 The differential accumulation pattern of
286
bioactive componds during fruit development is another important issue, as previously
287
reported by Fait et al.37 They found that metabolism in strawberry plant during
288
development is organ specific, with a shift towards the conversion of proanthocyanidins
289
into anthocyanins during maturation in the receptacle, and the accumulation of
290
ellagitannins and flavonoids in the achene during early and late development. In the
291
present work, targeted metabolomic profiling by UHPLC-MS/MS revealed significant
292
changes in phenolic compounds between three strawberry cultivars with different
293
sensitivity to the environment (i.e. 'Camarosa', 'Palomar' and 'Festival'), including
294
anthocyanins, ellagic acids and ellagitannis, flavonols, flavan-3-ols and phenolic acids.
295
Phenolic profiles showed that the most resistant cultivar (cv. 'Camarosa') presents
296
higher levels of anthocyanins, ellagic acid rhamnoside, quercetin 3-O-glucuronide and
297
procyanidin B2 compared with more sensitive strawberries, which might indicate an up-
298
regulation of the shikimate and malonate pathway. In this sense, it is well known that
299
anthocyanins, flavonol glycosides and proanthocyanidins are synthesized in response to
12 ACS Paragon Plus Environment
Page 12 of 35
Page 13 of 35
Journal of Agricultural and Food Chemistry
300
environmental and developmental signals.38 Furthermore, it has been previously
301
described that increased content of flavonoids could be associated with tolerance
302
mechanisms against abiotic stress, and that tolerance depends on the species, genotype
303
and stage of development of the plant. Thus, it has been demonstrated that strawberry
304
production of polyphenols is highly influenced by drought stress,25,39-40 salt stress,19,20
305
availability of nutrients in growing medium,41,42 extreme temperatures43,44 and
306
ultraviolet radiation.45,46
307
Predictive modeling of strawberry cultivars using Random Forest
308
Complementarily, random forest (RF) analysis was also applied in order to compare its
309
statistical performance with that provided by PLS-DA models, since this statistical tool
310
has demonstrated excellent efficiency during the last decade for the analysis of complex
311
data sets in many research areas, including food science. RF is a non-parametric and
312
non-linear classification and regression algorithm based on ensemble learning, which
313
operates by generating multiple decision trees on bootstrap samples taken from original
314
data and aggregating their results. Trees are split to many nodes using different subsets
315
of randomly selected input variables (m). Thus, the most important parameters to
316
optimized RF analysis are the value of m (attributes number made available at each
317
node during the growing of trees) and the number of decision trees. Moreover, the
318
model quality can be evaluated by the out-of-bag (OOB) error, which is the error
319
calculated by predicting data not used for growing the trees. In practice, other
320
measurements can also be used for validation purposes, such as Gini index (a measure
321
of node purity). Furthermore, the Gini index is also used to evaluate the potential
322
discriminant ability of each feature, i.e. the importance of each variable in the model.
323
Thus, the mean decrease of the Gini index (MDGI) can be used to select the most
324
important metabolites for inter-class discrimination.
13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
325
In the case of our three-class strawberry dataset, samples were perfectly classified using
326
1000 trees and two predictors for each node. During tree construction, one third of the
327
training samples were taken as the test set and the out-of-bag (OOB) data were then
328
used to estimate the classification error. Thereby, the statistical performance of the built
329
model was 97.4% accuracy rate for training set, 93.3% accuracy rate for testing set, and
330
2.56% error rate for OOB. Finally, the RF model was also used for feature selection,
331
since variables can be ranked according to their importance in the obtained
332
classification. The results indicated that the content of chlorogenic acid, ellagic acid
333
rhamnoside, cyaniding 3-O-glucoside, quercetin 3-O-glucoside, epicatechin, sanguiin
334
H10, pelargonidin 3-O-rutinoside and hyperoside were by far the eight most
335
discriminant variables, in agreement with results provided by PLS-DA. Finally, we
336
compared the statistical performance of the two pattern recognition tools here employed
337
by computing the sensitivity (SENS) and specificity (SPEC) parameters from the
338
obtained confusion matrixes. Sensitivity and specificity mean values for each cultivar,
339
as well as for the overall dataset, are listed in Table 2. The results show that the three
340
cultivars are predicted with low SENS (values range from 66-88.8%) when PLS-DA
341
modeling was applied, while the SPEC values were close to 100%. On the other hand,
342
RF modeling significantly improved SENS values for the three groups, especially for
343
'Camarosa' (100%), maintaining similar specificity values. Therefore, taking into
344
account the overall results, it could be concluded that RF provide higher sensitivity
345
compared with the PLS-DA model, by increasing about 24% (models A and B) and
346
15.6% (model C), but slightly decreasing the specificity.
347
Metabolomic alterations associated with agronomic practices
348
The nutritional value of berries is also significantly affected by agronomic and
349
environmental conditions in which plants are cultivated, including soil type, climate and
14 ACS Paragon Plus Environment
Page 14 of 35
Page 15 of 35
Journal of Agricultural and Food Chemistry
350
critical temperatures.33,47-49 For this reason, metabolomic data were re-analyzed by PLS-
351
DA considering the different agronomic conditions as categorical factors with the aim to
352
elucidate the effect of cultivation practices on the bioactive content of strawberries.
353
Each cultivar was statistically treated separately, and the performance of models
354
obtained was evaluated by using the parameters Rx2, Ry2 and Q2, listed in Table 3.
355
Furthermore, discriminant metabolites were selected following the procedure previously
356
described (Table 4).
357
Considering the macrotunnel type as the categorical factor (i.e. covered and uncovered),
358
PLS-DA models were able to explain 67.2, 68.1 and 64.9% of the variance for cv.
359
'Camarosa', 'Festival' and 'Palomar', respectively. The VIP scores indicated that
360
anthocyanins (pelargonidin 3-O-glucoside, cyaniding 3-O-glucoside and pelargonidin 3-
361
O-rutinoside) and ellagic acid derivatives (pentoside and rhamnoside) were the most
362
important variables to discriminate between the two groups. For the three cultivars,
363
anthocyanins showed higher concentrations in strawberries from covered macrotunel.
364
Particularly, higher content of pelargonidin 3-O-glucoside, cyaniding 3-O-glucoside,
365
pelargonidin 3-O-rutinoside, ellagic acid pentoside and ellagic acid rhamnoside was
366
found in 'Camarosa' strawberries grown in covered macrotunnels. Similarly,
367
concentrations of anthocyanins, procyanidin B2, catechin and quercetin 3-O-
368
glucuronide were higher in cv. 'Festival'; and pelargonidin 3-O-glucoside, pelargonidin
369
3-O-rutinoside, ellagic acid, ellagic acid rhamnoside, sanguiin H10, procyanidin B2 and
370
catechin for cv. 'Palomar'. These differences could be allocated to the controlled
371
climatic conditions that strawberries experiment when are grown in covered
372
macrotunnel (reduced light intensity, higher temperature, controlled moisture). Under
373
these conditions, previous studies have demonstrated an accumulation of anthocyanins
374
in strawberries.3,45 Josuttis et al.45 also found a higher content of flavonols and a lower
15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
375
total concentration of anthocyanins in open-field strawberry fruits. On the other hand,
376
Wang and Zheng43 observed that when strawberry fruits were submitted to relatively
377
high temperatures (22-30 °C), they increased the antioxidant concentration in the fruits
378
with respect to the plants that grew at lower temperatures (12-18 °C).
379
In relation to the electrical conductivity, a poor separation was observed in PLS-DA
380
models between samples, yielding low predictability (Q2 = -0.21, -0.08 and 0.179, for
381
'Camarosa', 'Festival' and 'Palomar', respectively). In that case, the most discriminant
382
variables were pelargonidin 3-O-glucoside, ellagic acid, ellagic acid rhamnoside,
383
quercetin 3-O-glucuronide, procyanidin B2 and catechin for cv. 'Camarosa';
384
pelargonidin 3-O-glucoside, pelargonidin 3-O-rutinoside, ellagic acid rhamnoside,
385
quercetin 3-O-glucuronide, procyanidin B2,
386
'Festival'; and pelargonidin 3-O-glucoside, pelargonidin 3-O-rutinoside, ellagic acid
387
rhamnoside, procyanidin B2 and catechin for cv. 'Palomar'. For the three cultivars,
388
increased pelargonidin 3-O-glucoside content was observed for higher irrigation
389
conductivity. However, only significant changes of pelargonidin 3-O-glucoside (33%)
390
and ellagic acid rhamnoside (12%) were found for cv. 'Palomar'. These results are in
391
accordance with those published by Keutgen and Pawelzik,19,20 who found that
392
moderate salt stress significantly increased the antioxidant capacity, total phenolic
393
compounds and anthocyanins in cv. 'Korona'.
394
Finally, PLS-DA models according to the growing substrate explained 65.4, 64.0 and
395
59.8% of the variance for 'Camarosa', 'Festival' and 'Palomar', but low predictability
396
values were observed (Table 3). In this case, discriminant compounds between study
397
groups were pelargonidin 3-O-glucoside, pelargonidin 3-O-rutinoside, cyaniding 3-O-
398
glucoside, ellagic acid pentoside, ellagic acid rhamnoside, quercetin 3-O-glucuronide,
399
procyanidin B2 and catechin for cv. 'Camarosa'; pelargonidin 3-O-glucoside,
catechin and chlorogenic acid for cv.
16 ACS Paragon Plus Environment
Page 16 of 35
Page 17 of 35
Journal of Agricultural and Food Chemistry
400
pelargonidin 3-O-rutinoside, ellagic acid rhamnoside, ellagic acid pentoside and
401
chlorogenic acid for cv. 'Festival'; and pelargonidin 3-O-glucoside, ellagic acid, ellagic
402
acid pentoside, ellagic acid rhamnoside, sanguiin H10, procyanidin B2 and catechin for
403
cv. 'Palomar'. The higher content of pelargonidin 3-O-glucoside was observed in
404
strawberries grown in perlite, regardless the cultivar. For 'Camarosa' cultivar, the use of
405
coconut fiber yielded increased pelargonidin 3-O-rutinoside, ellagic acid rhamnoside
406
and quercetin 3-O-glucuronide, while strawberries cultivated in perlite were richer in
407
procyanidin B2 and catechin. For 'Festival' cultivar, ellagic acid rhamnoside was
408
increased in strawberries cultivated in coconut fiber, while perlite induced the
409
accumulation of pelargonidin 3-O-rutinoside and ellagic acid pentoside. Finally,
410
increased levels of ellagic acid and procyanidin B2 were observed in 'Palomar'
411
strawberries grown in coconut fiber, while those grown in perlite were richer in
412
catechin, ellagic acid pentoside and sanguiin H10. It should be noted that these results
413
agree with previous findings that point to coconut fiber as the best soilless substrate for
414
strawberry cultivation.48
415
Therefore, it could be concluded that the application of targeted metabolomics and
416
subsequent multivariate statistics by means of PLS-DA and RF allows the identification
417
of changes in levels of strawberry polyphenols associated with variety and the
418
application of various agronomic practices, including different electrical conductivities
419
of irrigation, macrotunnel type and substrate. The most discriminant bioactive
420
compounds were anthocyanins (pelargonidin 3-O-glucoside, pelargonidin 3-O-
421
rutinoside, cyaniding 3-O-glucoside), ellagic acid derivatives (ellagic acid rhamnoside,
422
sanguiin H10), flavan-3-ols (catechin, procyanidin B2), chlorogenic acid, and quercetin
423
3-O-glucuronide. Thus, it should be noted that the understanding the influence of
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
424
genotype and agronomic conditions on the biosynthesis of strawberry antioxidants
425
might be crucial for the promotion of strawberry fruits as a nutraceutical food.
426 427
Abbreviations. MRM, multiple reaction monitoring; PLS-DA, partial least squares
428
discriminant analysis; RF, random forest; SENS, sensitivity; SPEC, specificity; VIP,
429
variable importance for the projection.
430
Supporting Information. This material is available free of charge via the Internet at
431
http://pubs.acs.org.
432
Table S1. Retention time and MRM parameters of phenolic compounds profiled by
433
UHPLC-MS/MS analysis.
434
Figure S1. Discriminant variables obtained from Random Forest analysis ranked
435
according to the Gini index.
436 437
Author Information
438
*Corresponding authors (Tel: +34 959219975, Fax: +34 959219942; E-mail:
439
[email protected]), (Tel: +34 959219958, Fax: +34 959219942; E-mail:
440
[email protected])
441 442
References
443
(1) Łabanowska, M.; Filek, M.; Kurdziel, M., Bidzińska; E., Miszalski, Z.;
444
Hartikainen, H. EPR spectroscopy as a tool for investigation of differences in
445
radical status in wheat plants of various tolerances to osmotic stress induced by
446
NaCl and PEG-treatment. J. Plant Physiol. 2013, 170, 136-145.
18 ACS Paragon Plus Environment
Page 18 of 35
Page 19 of 35
Journal of Agricultural and Food Chemistry
447
(2) Cheynier, V.; Comte, G.; Davies, K. M.; Lattanzio, V.; Martens, S. Plant phenolics:
448
recent advances on their biosynthesis, genetics, and ecophysiology. Plant Physiol.
449
Biochem. 2013, 72, 1-20.
450
(3) Akhatou, I.; Fernández-Recamales, A. Nutritional and nutraceutical quality of
451
strawberries in relation to harvest time and crop conditions. J. Agric. Food Chem.
452
2014, 62, 5749-5760.
453
(4) Akhatou, I.; Fernández-Recamales, A. Influence of cultivar and culture system on
454
nutritional and organoleptic quality of strawberry. J. Sci. Food Chem. 2014, 94,
455
866-875.
456
(5) Forbes-Hernandez, T. Y.; Gasparrini, M.; Afrin, S.; Bompadre, S.; Mezzetti, B.;
457
Quiles, J. L.; Giampieri, F. ; Battino, M. The healthy effects of strawberry
458
polyphenols: which strategy behind antioxidant capacity? Crit. Rev. Food Sci. Nutr.
459
2016, 56, S46-S59.
460
(6) González-Sarrías, A.; Núñez-Sánchez, M. A.; Tomás-Barberán, F. A.; Espín, J. C.
461
Neuroprotective effects of bioavailable polyphenol-derived metabolites against
462
oxidative stress-induced cytotoxicity in human neuroblastoma SH-SY5Y cells. J.
463
Agric. Food Chem. 2017, 65, 752-758.
464
(7) Gasparrini, M.; Forbes-Hernandez, T. Y.; Giampieri, F.; Afrin, S.; Alvarez-Suarez,
465
J. M.; Mazzoni, L.; Mezzetti, B.; Quiles, J. L.; Battino, M. Anti-inflammatory
466
effect of strawberry extract against LPS-induced stress in RAW 264.7
467
macrophages. Food Chem. Toxicol. 2017, 102, 1-10.
468 469
(8) McDougall, G. J.; Kulkarni, N. N.; Stewart, D. Berry polyphenols inhibit pancreatic
lipase activity in vitro. Food Chem. 2009, 115, 193-199.
470
(9) da Silva Pinto M.; de Carvalho, J. E.; Lajolo, F. M.; Genovese, M. I.; Shetty, K.
471
Evaluation of antiproliferative, anti-type 2 diabetes, and antihypertension potentials
19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
472
of ellagitannins from strawberries (Fragaria × ananassa Duch.) using in vitro
473
models. J. Med. Food 2010, 13, 1027-1035.
474
(10) McDougall, G. J.; Ross, H. A.; Ikeji, M.; Stewart, D. Berry extracts exert different
475
anti-proliferative effects against cervical and colon cancer cells grown in vitro. J.
476
Agric. Food Chem. 2008, 56, 3016-3023.
477 478
(11) Basu, A.; Rhone, M.; Lyons, T. J. Berries: emerging impact on cardiovascular
health. Nutr. Rev. 2010, 68, 168-177.
479
(12) Battino, M.; Beekwilder, J.; Denoyes-Rothan, B.; Laimer, M.; McDougall, G. J.;
480
Mezzetti, B. Bioactive compounds in berries relevant to human health. Nutr. Rev.
481
2009, 67, S145-S150.
482
(13) Tulipani, S.; Mezzetti, B.; Capocasa, F.; Bompadre, S.; Beekwilder, J.; De Vos, C.
483
H.; Capanoglu, E.; Bovy, A.; Battino; M. Antioxidants, phenolic compounds, and
484
nutritional quality of different strawberry genotypes. J. Agric. Food Chem. 2008,
485
56, 696-704.
486
(14) Capocasa, F.; Scalzo, J.; Mezzetti, B.; Battino, M. Combining quality and
487
antioxidant attributes in the strawberry: The role of genotype. Food Chem. 2008,
488
111, 872-878.
489
(15) Crespo. P.; Giné Bordonaba, J.; Terry, L. A.; Carlen, C. Characterisation of major
490
taste and health-related compounds of four strawberry genotypes grown at different
491
Swiss production sites. Food Chem. 2010, 122, 16-24.
492
(16) Padula, M. C.; Lepore, L.; Milella, L.; Ovesna, J.; Malafronte, N.; Martelli, G.; de
493
Tommasi, N. Cultivar based selection and genetic analysis of strawberry fruits with
494
high levels of health promoting compounds. Food Chem. 2013, 140, 639-646.
495
(17) Akhatou, I.; González-Domínguez, R.; Fernández-Recamales, A. Investigation of
496
the effect of genotype and agronomic conditions on metabolomic profiles of
20 ACS Paragon Plus Environment
Page 20 of 35
Page 21 of 35
Journal of Agricultural and Food Chemistry
497
selected strawberry cultivars with different sensitivity to environmental stress.
498
Plant Physiol. Biochem. 2016, 101, 14-22.
499
(18) Cevallos-Cevallos, J. M.; Futch, D. B.; Shits, T.; Folimonova, S. Y.; Reyes-de-
500
Corcuera, J. I. GC-MS metabolomic differentiation of selected citrus varieties with
501
different sensitivity to citrus huanglongbing. Plant Physiol. Biochem. 2012, 53, 69-
502
76.
503 504 505 506
(19) Keutgen, A. J.; Pawelzik, E. Modification of strawberry fruit antioxidant pool and
fruit quality under NaCl salinity. J. Agric. Food Chem. 2007, 55, 4066-4072. (20) Keutgen, A. J.; Pawelzik, E. Quality and nutritional value of strawberry fruit under
long term salt stress. Food Chem. 2008, 107, 1413-1420.
507
(21) Del Bubba, M.; Checchini, L.; Chiuminatto, U.; Doumett, S.; Fibbia, D.; Giordanic,
508
E. Liquid chromatographic/electrospray ionization tandem mass spectrometric
509
study of polyphenolic composition of four cultivars of Fragaria vesca L. berries
510
and their comparative evaluation. J. Mass Spectrom. 2012, 47, 1207-1220.
511
(22) Lopes da Silva, F.; Escribano-Bailón, M. T.; Pérez Alonso, J. J.; Rivas-Gonzalo, J.
512
C.; Santos-Buelga, C. Anthocyanin pigments in strawberry. LWT Food Sci. Tech.
513
2007, 40, 374-382.
514
(23) Zhang, J.; Wang, X.; Yu, O.; Tang, J.; Gu, X.; Wan, X.; Fang, C. Metabolic
515
profiling of strawberry (Fragaria x ananassa Duch.) during fruit development and
516
maturation. J. Exp. Bot. 2011, 62, 1103-1118.
517
(24) Häkkinen, S.; Heinonen, M.; Kärenlampi, S.; Mykkänen, H.; Ruuskanen, J.;
518
Törrónen, A. Screening of selected flavonoids and phenolic acids in 19 berries.
519
Food Res. Intern. 1999, 32, 345-353.
21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
520
(25) Weber, N.; Zupanc, V.; Jakopic, J.; Veberic, R.; Mikulic-Petkovsek, M.; Stampar,
521
F. Influence of deficit irrigation on strawberry (Fragaria × ananassa Duch.) fruit
522
quality. J. Sci. Food Agric. 2017, 9, 849-857.
523
(26) Seeram, N. P.; Lee, R.; Scheuller, H. S.; Heber, D. Identification of phenolic
524
compounds in strawberries by liquid chromatography electrospray ionization mass
525
spectroscopy. Food Chem. 2006, 97, 1-11.
526
(27) De Pascual-Teresa, S.; Santos-Buelga, C.; Rivas-Gonzalo, J. C. Quantitative
527
analysis of flavan-3-ols in Spanish foodstuffs and beverages. J. Agric. Food Chem.
528
2000, 48, 5331-5337.
529
(28) Buendía, B.; Gil, M. I.; Tudela, J. A.; Gady, A. L.; Medina, J. J.; Soria, C.;
530
López, J. M.; Tomás-Barberán, F. A. HPLC-MS analysis of proanthocyanidin
531
oligomers and other phenolics in 15 strawberry cultivars. J. Agric. Food Chem.
532
2010, 58, 3916-3926.
533 534 535 536
(29) Schuster, B.; Herrmann, K. Hydroxybenzoic and hydroxycinnamic acid derivatives
in soft fruits. Phytochem. 1985, 24, 2761-2764. (30) Mattila, P.; Hellstrom, J.; Törrönen, R. Phenolic acids in berries, fruits, and
beverages. J. Agric. Food Chem. 2006, 54, 7193-7199.
537
(31) Määttä-Riihinen, K. R.; Kamal-Eldin, A.; Törrönen, A. R. Identification and
538
quantification of phenolic compounds in berries of Fragaria and Rubus species
539
(family Rosaceae). J. Agric. Food Chem. 2004, 52, 6178-6187.
540
(32) Giampieri, F.; Alvarez-Suarez, J. M.; Battino, M. Strawberry and human health:
541
Effects beyond antioxidant activity. J. Agric. Food Chem. 2014, 62, 3867-3876.
542
(33) Cocco, C.; Magnani, S.; Maltoni, M. L.; Quacquarelli, I.; Cacchi, M.; Antunes, L.
543
E. C.; D’Antuono, L. F.; Faedi, W.; Baruzzi, G. Effects of site and genotype on
22 ACS Paragon Plus Environment
Page 22 of 35
Page 23 of 35
Journal of Agricultural and Food Chemistry
544
strawberry fruits quality traits and bioactive compounds. J. Berry Res. 2015, 5, 145-
545
155.
546
(34) Pradas, I.; Medina, J. J.; Ortiz, V.; Moreno-Rojas, J. M. ‘Fuentepina’ and ‘Amiga’,
547
two new strawberry cultivars: Evaluation of genotype, ripening and seasonal effects
548
on quality characteristics and health-promoting compound. J. Berry Res. 2015, 5,
549
157-171.
550
(35) D'Urso, G.; d'Aquino, L.; Pizza, C.; Montor, P. Integrated mass spectrometric and
551
multivariate data analysis approaches for the discrimination of organic and
552
conventional strawberry (Fragaria ananassa Duch.) crops. Food Res. Int. 2015, 77,
553
264-272.
554
(36) Diamanti, J.; Mazzoni, L.; Balducci, F.; Cappelletti, R.; Capocasa, F.; Battino, M.;
555
Dobson, G.; Stewart, D.; Mezzetti, B. Use of wild genotypes in breeding program
556
increases strawberry fruit sensorial and nutritional quality. J. Agric. Food Chem.
557
2014, 62, 3944-3953.
558
(37) Fait, A.; Hanhineva, K.; Beleggia , R.; Dai, N.; Rogachev, I.; Nikiforova, V. J.;
559
Fernie, A. R.; Aharoni, A. Reconfiguration of the achene and receptacle metabolic
560
networks during strawberry fruit development. Plant Physiol. 2008, 148, 730-750.
561
(38) Peng, H.; Yang, T.; Whitaker, B. D.; Shangguan, L.; Fang, J. Calcium/calmodulin
562
alleviates substrate inhibition in a strawberry UDP-glucosyltransferase involved in
563
fruit anthocyanin biosynthesis. BMC Plant Biol. 2016, 16, 197-207.
564
(39) Terry, L. A.; Chope, G. A.; Giné-Bordonaba, J. Effect of water deficit irrigation
565
and inoculation with Botrytis cinerea on strawberry (Fragaria x ananassa) fruit
566
quality. J. Agric. Food Chem. 2007, 55, 10812-10819.
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
567
(40) Giné-Bordonaba, J.; Terry, L. A. Effect of deficit irrigation and methyl jasmonate
568
application on the composition of strawberry (Fragaria x ananassa) fruit and
569
leaves. Sci. Horti. 2016, 199, 63-70.
570
(41) Sousa. C; Pereira, D. M.; Pereira, J. A.; Bento, A.; Rodrígues, M. A.; Dopico-
571
García, S.; Valentao, P.; Lopes, G.; Ferreres, F.; Seabra, R. M.; Andrade, P. B.
572
Multivariate analysis of Trochuda cabbage (Brassica oleracea L. var. costata DC)
573
phenolics: influence of fertilizers. J. Agr. Food Chem. 2008, 56, 2231-2239.
574
(42) Valentinuzzi, F.; Mason, M.; Scampicchio, M.; Andreotti, C.; Cesco, S.; Mimmo,
575
T. Enhancement of the bioactive compound content in strawberry fruits grown
576
under iron and phosphorus deficiency. J. Sci. Food Agric. 2015, 95, 2088-2094.
577
(43) Wang, S. Y.; Zheng, W. Effect of plant growth temperature on antioxidant capacity
578
in strawberry. J. Agric. Food Chem. 2001, 49, 4977-4982.
579
(44) Matsushita, K.; Sakayori, T.; Ikeda, T. The effect of high air temperature on
580
anthocyanin concentration and the expressions of its biosynthetic genes in
581
strawberry ‘Sachinoka’. Environ. Control Biol. 2016, 54, 101-107.
582
(45) Josuttis, M.; Dietrich, H.; Treutter, D.; Will, F.; Linnemannstöns, L.; Krüger, E.
583
Solar UVB response of bioactives in strawberry (Fragaria x ananassa Duch. L.): A
584
comparison of protected and open-field cultivation. J. Agric. Food Chem. 2010, 58,
585
12692-12702.
586
(46) Rubin de Oliveira, I.; Rodrigues Crizel, G.; Severo, J.; Renard, C.; Clasen Chaves,
587
F.; Valmor Rombsaldi, C. Preharvest UV-C radiation influences physiological,
588
biochemical, and transcriptional changes in strawberry cv. Camarosa. Plant
589
Physiol. Biochem. 2016, 108, 391-399.
24 ACS Paragon Plus Environment
Page 24 of 35
Page 25 of 35
Journal of Agricultural and Food Chemistry
590
(47) Nestby, R.; Martinussen, I.; Krogstad, T.; Uleberg, E. Effect of fertilization, tiller
591
cutting and environment on plant growth and yield of European blueberry
592
(Vaccinium myrtillus L.) in Norwegian forest fields. J. Berry Res. 2014, 4, 79-95.
593
(48) Martínez, F.; Oliveira, J. A.; Calvete, E. O.; Palencia, P. Influence of growth
594
medium on yield, quality indexes and SPAD values in strawberry plants. Sci.
595
Hortic. 2017, 217, 17-27.
596 597
(49) Nestby, R.; Takeda, F. Method to reduce low temperature stress (LTS) in red
raspberry (Rubus idaeus L.). J. Berry Res. 2015, 5, 219-230.
598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
615
FIGURE CAPTIONS
616
Figure 1. Major phenolic compounds detected in strawberry fruits.
617
Figure 2. Scores plots of PLS-DA models for two-class comparisons. (A) 'Festival' vs.
618
'Camarosa', (B) 'Palomar' vs. 'Camarosa', (C) 'Palomar' vs. 'Festival'.
619
Figure 3. Box plots with whiskers for discriminant polyphenols for differentiation of
620
strawberry cultivars.
621
26 ACS Paragon Plus Environment
Page 26 of 35
Page 27 of 35
Journal of Agricultural and Food Chemistry
Table 1.
Mean, Minimum and Maximum Concentration Values for Phenolic
Metabolites Detected in the Three Strawberry Varieties. 'Camarosa'
'Festival'
'Palomar'
(n=18)
(n=18)
(n=18)
Compounds Mean
Min
Max
Mean
Min
(µg/g)
(µg/g)
(µg/g)
(µg/g)
(µg/g) (µg/g) (µg/g) (µg/g) (µg/g)
gallic acid
0.12
0.09
0.16
0.12
0.10
0.19
0.15
0.10
0.22
protocatechuic acid
0.06
0.05
0.08
0.06
0.05
0.09
0.06
0.06
0.07
p-hydroxybenzoic acid
0.24
0.18
0.44
0.21
0.15
0.36
0.22
0.15
0.39
vanillic acid
0.02
0.01
0.04
0.04
0.02
0.10
0.03
0.01
0.07
cinnamic acid
0.71
0.36
1.25
0.64
0.24
2.53
0.84
0.40
2.00
caffeic acid
0.62
0.31
1.08
0.46
0.31
1.03
0.54
0.28
1.31
p-coumaric acid
1.22
0.48
2.47
0.67
0.30
1.97
1.35
0.47
3.97
m-coumaric acid
0.01
0.01
0.02
0.01
0.01
0.02
0.02
0.01
0.09
chlorogenic acid
3.50
2.64
4.79
12.73
3.51
16.60
5.80
3.56
8.00
(-) epicatechin
0.01
0.00
0.05
0.01
0.00
0.04
0.06
0.01
0.11
(+) catechin
14.00
4.63
24.66
16.57
1.33
24.65
11.52
4.69
17.06
(-) epicatechin gallate
0.06
0.05
0.07
0.06
0.05
0.07
0.06
0.05
0.07
procyanidin B2
0.15
0.12
0.21
0.18
0.11
0.92
0.18
0.14
0.23
procyanidin B1
23.13
6.15
50.54
30.03
4.14
49.32
14.71
2.95
29.07
quercetin 3-O-glucoside
0.30
0.20
0.39
0.19
0.12
0.31
0.29
0.13
0.39
quercetin 3-O-galactoside
0.13
0.05
0.20
0.05
0.01
0.14
0.12
0.03
0.20
quercetin
0.10
0.02
0.39
0.07
0.03
0.28
0.13
0.02
0.68
quercetin 3-O-glucuronide
43.07
23.43
59.09
34.57
24.96
47.44
35.59
28.76 43.29
27 ACS Paragon Plus Environment
Max
Mean
Min
Max
Journal of Agricultural and Food Chemistry
kaempferol 3-O-glucoside
Page 28 of 35
1.06
0.29
1.73
1.05
0.51
1.52
1.15
0.56
1.62
hexoside
1.89
0.39
3.13
1.45
0.50
1.87
1.60
0.82
2.33
apigenin
0.09
0.00
0.27
0.05
0.00
0.16
0.05
0.00
0.22
ellagic acid pentoside
25.51
18.88
33.10
22.33
17.15
28.58
22.24
17.09 26.99
HHDP galloyl hexose
9.23
7.87
17.33
9.57
7.50
13.60
9.15
7.65
11.45
sanguiin H10
2.37
1.54
5.23
1.74
1.14
3.23
3.52
1.65
4.89
ellagic acid rhamnoside
141.43 116.97 170.49 123.40
103.82 147.38 107.87 94.49 126.54
ellagic acid
25.40
18.37
44.78
24.80
18.39
30.38
22.28
17.02 30.66
cyanidin 3-O-glucoside
12.53
7.68
21.17
10.06
6.69
16.27
7.29
5.70
pelargonidin 3-O-glucoside
121.29 56.88
213.00 105.20
5.78
185.25 93.99
44.38 136.39
52.87
7.45
41.51
7.43
kaempferol
3-O-acetyl-
pelargonidin 3-O-rutinoside 29.56
17.72
18.86
28 ACS Paragon Plus Environment
15.69
8.57
34.66
Page 29 of 35
Journal of Agricultural and Food Chemistry
Table 2. Comparison of the Statistical Performance for the Two Classification Methods: Partial Least Squares Discriminant Analysis (PLS-DA) and Random Forest (RF). 'Camarosa'
'Festival'
'Palomar'
Overall
SENS
SPEC
SENS
SPEC
SENS
SENS
SPEC
Cam-Fes
66.6
94.4
88.8
100
77.7
97.2
Cam-Pal
72.2
100
Model SPEC
PLS83.3
100
77.7
100
DA Fes-Pal RF
100
94
77.7
100
88.8
94.4
83.3
97.2
94.4
100
94.4
100
96.3
96.3
29 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Table 3. Statistical Parameters of PLS-DA Models for Classification of Strawberries According to Growing Conditions (Type of Coverage, Electrical Conductivity and Substrate). Cultivar
Factor
R2 X
R2Y
Q2
Coverage
0.672
0.512
0.241
0.648
0.157
-0.21
Substrate
0.654
0.183
-0.21
Coverage
0.681
0.399
-0.21
EC
0.629
0.308
-0.08
Substrate
0.640
0.284
-0.05
Coverage
0.649
0.504
0.084
EC
0.665
0.438
0.179
Substrate
0.598
0.124
-0.21
'Camarosa' EC
'Festival'
'Palomar'
30 ACS Paragon Plus Environment
Page 30 of 35
Page 31 of 35
Journal of Agricultural and Food Chemistry
Table 4. Discriminant Polyphenolic Compounds Associated with Different Growing Conditions of Strawberries Cultivated in Soilless Systems (p-values in brackets). Coverage Compound
CAM
EC
Substrate
FES
PAL
CAM
FES
NS
NS
NS
NS
PAL
CAM
FES
PAL
NS
NS
NS
Pelargonidin 37% 3-O-
33%
(0.0405)
(0.0244)
glucoside Pelargonidin 61%
69%
3-O-
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
NS
(0.0003) (0.0189) rutinoside Cyanidin 3- 43% O-glucoside
38%
(0.0037) (0.0127)
Ellagic acid
12% NS
NS
rhamnoside
-12% NS
(0.0178)
(0.0222)
Ellagic acid
-17% NS
NS
NS
NS
NS
NS
NS
NS
pentoside
(0.0082)
Sanguiin
*
33% NS
H10
NS *
NS
NS
NS
NS
(0.0153) * perlite vs. rockwool; N.S. non significant
31 ACS Paragon Plus Environment
NS
NS
NS
Journal of Agricultural and Food Chemistry
FIGURE 1
32 ACS Paragon Plus Environment
Page 32 of 35
Page 33 of 35
Journal of Agricultural and Food Chemistry
FIGURE 2
33 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
FIGURE 3
34 ACS Paragon Plus Environment
Page 34 of 35
Page 35 of 35
Journal of Agricultural and Food Chemistry
TABLE OF CONTENTS GRAPHIC
35 ACS Paragon Plus Environment