Application of Targeted Metabolomics to Investigate Optimum Growing

Oct 11, 2017 - Therefore, taking into account the overall results, it could be concluded that RF provide higher sensitivity compared with the PLS-DA m...
0 downloads 6 Views 1MB Size
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