Endogenous circadian rhythms in polyphenolic composition induce

6 days ago - There is increasing evidence that the circadian clock is a significant driver of plant phytochemicals. However, little is known about the...
3 downloads 0 Views 3MB Size
Subscriber access provided by Kaohsiung Medical University

Agricultural and Environmental Chemistry

Endogenous circadian rhythms in polyphenolic composition induce changes in antioxidant properties in Brassica cultivars Pilar Soengas, María Elena Cartea, Pablo Velasco, and Marta Francisco J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01732 • Publication Date (Web): 31 May 2018 Downloaded from http://pubs.acs.org on May 31, 2018

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 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 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.

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 26

Journal of Agricultural and Food Chemistry

1

Endogenous circadian rhythms in polyphenolic composition

2

induce changes in antioxidant properties in Brassica cultivars

3

Pilar Soengas, M. Elena Cartea, Pablo Velasco, Marta Francisco*

4 5

Group of Genetics, Breeding and Biochemistry of Brassicas, Misión Biológica de

6

Galicia, Spanish Council for Scientific Research (CSIC), Pontevedra, Spain.

7 8

Running title: Circadian rhythms in polyphenols affect antioxidant properties

9

*Corresponding author: Marta Francisco

10

Tlf: 0034 986 85 48 00

11 12

Pilar Soengas: [email protected]

13

M. Elena Cartea: [email protected]

14

Pablo Velasco: [email protected]

15

Marta Francisco: [email protected]

16 17 18

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 26

19

ABSTRACT

20

There is increasing evidence that the circadian clock is a significant driver of plant

21

phytochemicals. However, little is known about the clock effect on antioxidant

22

metabolites in edible crops. Thus, the aim of the present investigation was to study

23

whether the antioxidant potential of Brassica cultivars is under circadian regulation and

24

its relationship with polyphenol content. To accomplish that we entrain plants of four

25

Brassica cultivars to light-dark cycles prior to release into continuous light. The

26

antioxidant activity and phenolic content was monitored at four time points of the day

27

during four consecutive days: two days under light-dark conditions followed by two

28

days under continuous light. Results showed daily oscillation of antioxidant activity. In

29

addition, those variations were related with endogenous circadian rhythms in

30

polyphenolics and exhibit a species-specific pattern. Considering together, we

31

determined that Brassica cultivars have an optimal time during a single day with

32

increased levels of health-phytochemicals.

33

Keywords:

34

polyphenolics, circadian clock.

Brassica

oleracea,

Brassica

rapa,

antioxidants,

35

2 ACS Paragon Plus Environment

ABTS,

FRAP,

Page 3 of 26

Journal of Agricultural and Food Chemistry

36

1. Introduction

37

Cruciferous vegetables, in particular those included into the Brassica genus, are

38

good sources of a variety of nutrients and health-promoting phytochemicals 1. The

39

principal vegetable species of this genus are Brassica oleracea (i.e. broccoli, cabbage,

40

cauliflower, kale and Brussels sprouts) and Brassica rapa (i.e. turnip, Chinese cabbage

41

and pak choi). Diets rich in those vegetables have been reported to possess many useful

42

properties for human health including anti-inflammatory, enzyme inhibition,

43

antimicrobial, antiallergic or cytotoxic antitumor activity 1, 2. Part of these properties can

44

be related to the high concentration of dietary antioxidants presents in these crops

45

including polyphenols, vitamins E and C, and carotenoids which can reduce/scavenge

46

reactive oxygene species (ROS). Various ROS, such as singlet oxygen, peroxynitrite

47

and hydrogen peroxide, must be continually removed from cells to maintain healthy

48

metabolic function. Diminishing the concentrations of ROS can have several benefits

49

possibly associated with ion transport systems and so may affect redox signaling 3.

50

Antioxidant compounds can be extracted and analyzed separately, but a more general

51

approach to know the capacity of vegetable extracts to neutralize ROS is by using

52

antioxidant potential assays. The antioxidant potential of Brassica vegetables is high

53

compared to other vegetable crops. In fact, broccoli and kale are among the ones having

54

the highest potential in the group of vegetable foods, including spinach, carrot, potato,

55

purple onion, green pepper, rhubarb or green bean 4.

56

The antioxidant properties in food supplies have been the focus of intense

57

research. However, during the growth period, plants are frequently subjected to various

58

biotic and abiotic factors that can affect the regulation of the biosynthetic pathways

59

involved in production of plant phytochemicals, including polyphenols, and therefore,

3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

60

change their antioxidant potential. Light is an important environmental factor that

61

regulate not only plant growth and development, but also the biosynthesis of both

62

primary and secondary metabolites

63

absolutely light dependent and its biosynthetic rate is related to light intensity and cycle

64

6-8

65

consequence in the resulting total flavonoids and total hydroxicinnamates production.

66

According to previous studies, a positive and significant relationship has been observed

67

between total phenolics, total flavonoids and antioxidant activities in Brassica crops 7, 9.

68

It seems that different light conditions would have a direct effect on antioxidant

69

activities in plants with changing total polyphenolic content. This may be related with

70

the role of those compounds in protecting plants against excess light damages 10.

5

. Polyphenolic biosynthesis in Brassica is

. However, different cultivars had a different response to light conditions and in

71

Plants, like other living organisms, have evolved circadian clocks, an

72

endogenous timing system, that allow them to anticipate and prepare for cyclical

73

changes in their environment, such as those associated with the transitions from night to

74

day and between seasons 11. In this way, the plant can coordinate accurately the timing

75

at which each metabolite is required to synchronize multiple physiological and

76

developmental responses. Many genes encoding enzymes of polyphenol biosynthesis in

77

the model plant Arabidopsis thaliana increase their transcripts before dawn 11. Thus, it

78

can be hypothesized that daily oscillations in the light-dark cycle could influence

79

polyphenolic composition in plants during a single day. To our knowledge, there are

80

currently no reported studies under controlled conditions on the influence of the

81

circadian clock on the contents of those bioactive compounds in Brassica cultivars.

82

New agriculture requires developing new ecologically based approaches for

83

increase plant phytochemicals with a role in plant protection and in human health. Thus, 4 ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

Journal of Agricultural and Food Chemistry

84

identifying specific growing conditions and determining the optimal harvest times in

85

order to obtain cultivars with increased concentrations of bioactive phytochemicals

86

could offer value-added commercial opportunities to the food industry. The aim of the

87

present investigation was to study whether the antioxidant potential of Brassica

88

cultivars is under circadian regulation and its relationship with polyphenol content. To

89

investigate for species-specific patterns we studied four different phenolic-containing

90

cultivars belonging to B. oleracea (broccoli and cabbage) and B. rapa (Chinese cabbage

91

and turnip greens).

92

93

2. Material and methods

94

2.1. Plant material and growing conditions

95

We studied different phenolic-containing cultivars from two Brassica species: cabbage

96

(B. oleracea var. capitata), broccoli (B. oleracea var. botrytis), Chinese cabbage (B.

97

rapa spp. pekinensis) and turnip greens (B. rapa ssp. rapa). Seeds of commercial

98

hybrids of each cultivar were purchased from vegetable seed market (Batlle Seeds,

99

Spain). All the genotypes were evaluated at the same time within a single plant growth

100

chamber. If hybridization is accompanied by consistent genomic changes, then we

101

expect synthetic commercial hybrids to exhibit consistent changes in phenotypic

102

characteristics when they grow under controlled environmental conditions. All seeds

103

were planted in potting soil, stratified at 4°C in the dark for three days to optimize

104

germination. Plants were sown in a randomized block design experiment consisting in

105

192 pots (4 x 4 cm), each containing one plant, distributed in three blocks. Within each 5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

106

block, the complete set of genotypes was randomly organized (16 plants × 4 genotypes).

107

The plants were entrained to cycles of 16 h of light and 8 h of dark at 22°C and 60% of

108

humidity for one month. At the stage of five leaves per plant, they were release into

109

continuous light and temperature (22°C). The soil was regularly watered once per week

110

until the day before of the beginning of the experiment.

111

Leave samples for metabolite analysis and antioxidant evaluation were collected

112

at two different light cycles during four consecutive days. Two days during the light-

113

dark conditions (16 h of light/8 h of dark) and other two days during continuous light.

114

Three replicate samples (each containing all the leaves from one plant) per cultivar were

115

harvested at four distinct daily time points defined by the light-dark conditions: -1 h

116

before and after the lights turned off (end of the day light period and beginning of the

117

night period, respectively) and 1h before and after the lights turned on (end of the night

118

period and beginning of the day light period, respectively).

119

All samples were transferred to the laboratory and conserved at –80°C prior to

120

lyophilization (BETA 2–8 LD plus, Christ) during 72 h. The dried material was

121

powdered by using an IKA-A10 (IKA-Werke GmbH & Co.KG) mill, and the fine

122

powder was used for further analysis.

123

2.2. Evaluation of antioxidant activity

124

Freeze-dried and ground samples (10mg) were extracted with 1 mL of 80% aqueous

125

methanol in dark maceration for 24h. After centrifugation (3700 rpm, 5 min),

126

methanolic extracts were employed in order to determine the antioxidant activity by

127

[2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid)] cation assay (ABTS) and by 6 ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26

Journal of Agricultural and Food Chemistry

128

ferric reducing antioxidant activity (FRAP). Three technical replications were analysed

129

for each sample Standard prepared with different concentrations of Trolox® (0, 0.008,

130

0.016, 0.024, 0.032, 0.040 mM) was also measured. Antioxidant activity was

131

normalized to Trolox® equivalents per gram (g) of dry weight (DW).

132

2.2.1. ABTS assay

133

The method of decolorization of free radicals ABTS employed was a modified version

134

of that used by Samarth et al

135

generated by oxidation of ABTS 7 mM potassium persulphate 2.45 mM in water, at

136

room temperature for 16 h. For each analysis, the ABTS solution was freshly diluted

137

with water in order to obtain an initial absorbance around 0.8 at 734 nm. An aliquot of

138

10 µL methanolic extract for sample was added to 250 µL of ABTS solution.

139

Absorbances were measured at 734 nm after 30 min of incubation in the dark at room

140

temperature.

141

2.2.2. FRAP assay

142

For FRAP assay, the procedure followed the method of Benzie and Strain14. The FRAP

143

reagent was prepared by mixing 10 volumes of 300 mM acetate buffer (pH 3.6), one

144

volume of 10 mM TPTZ (2,4,6-tripyridyl-s-triazine, Sigma-Aldrich, Germany) in 40

145

mM hydrochloric acid and one volume of 20 mM ferric chloride and then incubating at

146

37 °C for 5 min. Fifty microliters of extract were added to 250 µl of freshly prepared

147

FRAP reagent and mixed thoroughly. Readings were taken at 593 nm after 20 min in a

148

microplate spectrophotometer (Spectra MR).

12

and initially reported by Re et al

7 ACS Paragon Plus Environment

13

. ABTS was

Journal of Agricultural and Food Chemistry

Page 8 of 26

149

2.3. Phenolic determination and quantification

150

The same leaf samples that were used for antioxidant activity were also used for the

151

determination of phenolic compounds. Fifty milligrams of each sample were extracted

152

in 500 µL 70% methanol and sonicated during 1 h (model: 3510E-MTH, Bransonic®,

153

Mexico) to facilitate the extraction. The suspensions were allowed to stand overnight at

154

4 ºC and afterwards they were sonicated again for 1h. Then, samples were centrifuged at

155

3700 rpm for 15 min. Supernatants were recovered and were filtered through 0.22 µm

156

PTFE filters (Multi Screen®, Ireland). For phenolics identification, three samples from

157

each cultivar were analyzed by HPLC-DAD-ESI/MSn in an Agilent HPLC 1100 series

158

equipped with a diode array detector and mass detector in series (Agilent Technologies,

159

Waldbronn, Germany) at the same conditions described by Francisco et al.

160

spectrometry data were acquired in the negative ionization mode. Individual phenolic

161

compounds were identified based on the data obtained from the standard substances or

162

published literature based on the mass spectra together with the ([M–H]−) 15-17. Table S1

163

sumarizes the identification data for the identified phenolic compounds. Quantification

164

of phenolic compounds was carried out in an Ultra-High-Performance Liquid-

165

Chromatograph (UHPLC Nexera LC-30AD; Shimadzu) equipped with a Nexera SIL-

166

30AC injector and one SPD-M20A UV/VIS photodiode array detector was used for the

167

quantification. The UHPLC column was a Kinetex™ 2.6 µm C18 82-102 Å, LC

168

Column 100 × 4.6 mm, protected with a C18 guard cartridge. The flow rate was 0.4 mL

169

/min-1 and the oven temperature was set at 30 ºC. The mobile phase consisted of two

170

solvents: water-acetic acid (1%) (A) and methanol (B), starting with 10% B and using a

171

gradient to obtain 40% B at 15 min and 60 % B at 24 min. The injection volume was 3

172

µL. Chromatograms were recorded at 330 nm and data was processed on a computer

173

with the LabSolutions software (Shimadzu). Quinic acids derivatives were quantified as 8 ACS Paragon Plus Environment

15

. Mass

Page 9 of 26

Journal of Agricultural and Food Chemistry

174

chlorogenic acid (Sigma–Aldrich Chemie GmbH, Steinheim, Germany), flavonoids of

175

kamepferol as kaempferol-3-O-glucoside, flavonoids of quercetin as quercetin-3-O-

176

glucoside, flavonoids of isorhmanetin as isorhamnetin-3-O-glucoside (Merck,

177

Darmstadt, Germany) and sinapic acid and derivatives as sinapic acid (Sigma–Aldrich

178

Chemie GmbH, Steinheim, Germany). The calibration curves were made, at least with

179

six data points, from 0.01 to 1 mM. The polyphenols were reported as µmol g-1 DW. In

180

total, twenty seven compounds were identified (Table S1): Sums of each kind of

181

compound were expressed as: total caffeoyl quinic acids, total flavonoids, total sinapoyl

182

derivatives and total phenolics.

183

2.4. Statistical analysis

184

We examined temporal patterns of 33 variables, including individual phenolic

185

compounds, sums of compounds by category (total quinic acids, total sinapoyl

186

derivatives, total flavonoids and total phenolics) as well as the antioxidant variables

187

(ABTS and FRAP). To test how a daily light cycle affect antioxidant properties and

188

changes in polyphenolic composition under light-dark conditions, we conducted a

189

mixed model ANOVA where the main effects of plant cultivars, time point (beginning

190

of the night period, end of the night period, beginning of the light period and end of the

191

light period) and their interactions were considered as fixed factors. Block and their

192

interactions with plant cultivars and time point were considered as random factors. As

193

we found a significant interaction between plant cultivar and time point, we also

194

performed these analyses individually for each plant cultivar considering the time point

195

as a fixed effect and block and the interaction between block and time point as random

196

factors. Means comparisons were done with Fisher’s protected least significant

197

difference (LSD) at the 0.05 level of probability. 9 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

198

The Pearson correlation analysis were performed between antioxidant activity

199

methods (FRAP and ABTS) and total quinic acids, total sinapics, total flavonoids and

200

total phenolics with data from all samples collected at different times of the day in the

201

plants grown under light-dark cycles. To test the contribution of each individual

202

phenolic compound to variability of antioxidant activity, multiple regression analysis

203

was performed using a stepwise linear regression model. In this model, single phenolic

204

compound contents were independent variables and antioxidant activity measured as

205

ABTS and FRAP assays were dependent variables. .

206

For comparative analysis on metabolite abundance across different samples and

207

conditions, the data was normalized. The normalization procedure consisted of mean-

208

centering and division by the standard deviation of each variable. The resulting data

209

matrix was subjected to hierarchical clustering using the Euclidean distance method

210

associated with complete-linkage based on each phenolic trait variation among the four

211

times of the day evaluated for each cultivar which grown under light-dark conditions.

212

All statistical analyses including clustering and heatmap representations were done

213

using R 18.

214

For phytochemical rhythmicity determination, we used data from plants

215

transferred into continuous light conditions after the light-dark entrainment cycles (16 h

216

of light/8 h of dark). The circadian parameters, period and phase were studied on the

217

content (µmol g-1 DW) of phenolic compounds, selected as significantly (P ≤ 0.05)

218

correlated with the antioxidant activity. The period was stated as the necessary time for

219

a cycle to be completed, measured as the time between two consecutive maximum

220

(peaks) periods. The phase was determined as the time of the day of maximal phenolic

10 ACS Paragon Plus Environment

Page 10 of 26

Page 11 of 26

Journal of Agricultural and Food Chemistry

221

trait accumulation. Period and phase analysis were carried out using the online platform

222

BioDare2 (biodare2.ed.ac.uk) 19.

223

224

3. Results and discussion

225

3.1. Phenolic analysis: relationship with antioxidant activity

226

It is well established that antioxidant activity is related with the content and composition

227

of several antioxidant metabolites in plant extracts, including polyphenols. We aim to

228

study whether the antioxidant potential of different Brassica cultivars is under circadian

229

regulation and its relationship with polyphenol content. To accomplish that, first we

230

studied the polyphenolic profile of each Brassica cultivar and then we addressed the

231

phenolic traits with a major influence on the antioxidant activity. Our results confirmed

232

the polyphenolic composition according to previously reported data on flavonoids and

233

hydroxycinnamic acids in B. oleracea and B. rapa cultivars

234

compounds were identified and quantified among cultivars (Table S1). Regard to

235

flavonoids, flavonols of kaempferol glycosilated and/or acylated by different

236

hydroxycinnamic acids were the most abundant in these Brassica cultivars.

237

Hydroxycinnamic acids were mainly represented by quinic acids of caffeic, ferulic, p-

238

coumaric and sinapic acids in conjugation with sugar moieties or other

239

hydroxycinnamics. They were divided into two major categories, sum of total quinic

240

acids and sinapoyl derivatives.

11 ACS Paragon Plus Environment

20

. In total, 27 individual

Journal of Agricultural and Food Chemistry

Page 12 of 26

241

As it was suggested by our previous results, we corroborated that the

242

polyphenolic profile was similar among cultivars belonging to the same species (Table

243

S2). The B. rapa cultivars, were rich in hydrocycinnamic acids, mainly represented by

244

the sinapic acid and its derivatives, accounting for 60-75% of total phenolics. In

245

addition,

246

(methoxycaffeoyl)diglucoside-7-O-glucoside and Kaempferol-3,7-di-O-glucoside were

247

the flavonoids present in higher contents, each corresponding to 5-9% of total

248

phenolics. The B. oleracea cultivars showed similar contents of hydroxycinnamic acids

249

and flavonoids being the major compounds the sinapoyl derivatives, 1,2,2′-

250

trisinapoylgentiobioside and 1,2′-disinapoyl-2-feruloylgentiobioside, followed by the 3-

251

caffeoyl quinic acid and the flavonoid kaempferol-3-O-(sinapoyl)-diglucoside. In these

252

cultivars, each compound corresponded to 9-17% of total phenolics.

Kaempferol-3-O-(caffeoyl)diglucoside-7-O-glucoside,

Kaempferol-3-O-

253

In order to evaluate the relationship between the content and composition of

254

phenolics and the antioxidant activity expressed by the two different antioxidant assays

255

(ABTS and FRAP), Pearson’s correlation coefficient were performed with data from all

256

samples collected at different times of the day in the plants grown under light-dark

257

cycles. Since cultivars of each species shared the same phenolic profile we presented the

258

results by specie. High and significant (P ⩽ 0.01) correlations were found between total

259

quinic acids, total flavonoids, total sinapoyl derivatives and total phenolics for both

260

antioxidants assays in B. oleracea cultivars, ABTS (0.90, 0.79, 0.88 and 0.85

261

respectively) and FRAP (0.83, 0.74, 0.79 and 0.79 respectively). For B. rapa cultivars,

262

we only found significant correlations of ABTS and FRAP with total flavonoid content

263

(0.81 and 0.63, respectively).

12 ACS Paragon Plus Environment

Page 13 of 26

Journal of Agricultural and Food Chemistry

264

Several authors have found significant and high correlations between antioxidant

265

activities measured with electron-transfer based assays and total phenolic content,

266

employing the Folin-Ciocalteu method, in different samples of Brassica cultivars

267

However, from the overwhelming diversity of plant phenolics, only a limited number of

268

compounds appear to promote protection via antioxidant action. Thus, to create a better

269

model linking individual phenolic compounds to the antioxidant properties of Brassica

270

extracts, we used a stepwise regression analysis on which single phenolic compound

271

contents were independent variables and antioxidant activity measured as ABTS and

272

FRAP assays were dependent variables. The coefficient of determination (R2) is the

273

proportion of variability in a data set that is accounted by the statistical model and

274

provides a measure of how well future outcomes are likely to be predicted by the model.

275

The analysis generated a model where variation of p-coumaroyl quinic acid in B. rapa

276

cultivars and variation of 1,2,2′-trisinapoylgentiobioside in B. oleracea cultivars

277

explained 81% and 66% of total variation of the variability in ABTS, respectively. For

278

FRAP, this analysis generated a model where variation of few phenolic traits explained

279

97% and 94% of the antioxidant variability in B. rapa and B. oleracea, respectively. For

280

the B. rapa extracts, the flavonoid kaempferol-3,7-di-O-glucoside was the compound

281

with the strongest influence on FRAP (59%) followed by kaempferol-3-O-

282

(methoxycaffeoyl)diglucoside-7-O-glucoside (25%) and p-coumaroyl quinic acid

283

(13%). For the B. oleracea extracts, the sinapoyl derivative 1,2′-disinapoyl-2-

284

feruloylgentiobioside showed the highest influence on FRAP variation (63%) followed

285

by kaempferol-3,7-di-O-glucoside (31%). In concordance with our correlation analysis,

286

the relationships between the individual compounds and antioxidant traits were positive.

287

Our results corroborate the importance of hydroxycinnamic acids together flavonoid

13 ACS Paragon Plus Environment

9, 21

.

Journal of Agricultural and Food Chemistry

288

compounds in the antioxidant behavior of the extracts and also show the significant

289

contribution of a single phenolic to the total antioxidant capacity of Brassica plants.

290

3.2. Daily polyphenolic variation in Brassica cultivars induces changes in

291

antioxidants properties

292

To study whether temporal accumulation pattern of phenolic contents can influence the

293

antioxidant activity of Brassica cultivars, we entrain plants of each Brassica cultivar to

294

light-dark cycles and then monitored the phenolic content and antioxidant activity at

295

four different time points of the day during two consecutive days. Analysis of variance

296

showed that phenolic traits and antioxidant activity significantly varied among time

297

points and cultivars (Table 1). In addition, the variation of antioxidant assays as well as

298

of most of the phenolic traits was contingent upon the time as demonstrated by the

299

significant interaction between plant cultivars and time points for those traits (Table 1).

300

The phenolic variation among the four times of the day evaluated for each cultivar by

301

each individual trait was further summarized (Table S2) and clustered (Figure 1).

302

Samples were grouped into two clades. In clade I, all samples collected during the light

303

hours as well as most of the samples collected in early night were clustered. This clade

304

also represents low content of phenolic traits. In contrast with that, in clade II samples

305

harvested at the end of the night period exhibiting high accumulation of phenolic traits

306

were clustered. Therefore, in general, all cultivars showed highest total phenolic content

307

at the end of the dark period. However, changes of phenolic traits related with

308

antioxidant potential exhibit a species-specific pattern. The B. oleracea cultivars

309

(broccoli and cabbage) accumulated higher levels of antioxidant traits during the dark

310

period. The B. rapa cultivars (Chinese cabbage and turnip greens) accumulated the 14 ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

Journal of Agricultural and Food Chemistry

311

highest quantities of phenolic traits related with antioxidant activity at the beginning of

312

the day (Figure 1). Therefore, in B. rapa cultivars not all the phenolic traits reached the

313

maximum values at the same time of the day.

314

Since the interaction between plant cultivars and time points for both antioxidant assays

315

it was also significant (Table 1) we performed individual analysis by each cultivar

316

(Figure 2). In particular, we found that broccoli and cabbage showed the lowest levels

317

of antioxidant activity at the end of the light period. The antioxidant activity in those

318

cultivars goes increasing thought darkness reaching the maximum levels at the end of

319

dark period. For Chinese cabbage and turnip greens the highest antioxidant potential

320

was found at the beginning of the light period while remain almost constant at the others

321

time points evaluated (Figure 2). Broccoli and cabbage showed the highest variation (2-

322

fold) of antioxidant activity between the end of the light period and the end of the dark

323

period. Thus, it has been evidenced the existence of daily oscillations of antioxidant

324

activity in the studied Brassica cultivars when they grow under light-dark cycles. In

325

addition, both antioxidant assays consistently showed that those fluctuations were

326

similar among cultivars belonging to the same species. This suggests that temporal

327

regulation of antioxidant traits is under genetic control, may be because antioxidant

328

compounds are essential for sustaining important plant physiological processes.

329

Phenolic compounds contribute to the color, sensory characteristics and

330

antioxidant potential of vegetables. However, those metabolites also perform a variety

331

of functions in the plant. They generally act also as potent antioxidants protecting the

332

plant in response to UV light. Once light reaches a plant, provides the energy for

333

photosynthesis, the principal source of substrates for all other biosynthetic pathways.

334

However, an increased photorespiration during the light period, induce the production 15 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

335

of ROS, which are important in signaling, but become harmful for the plants if

336

accumulated 3. The production of highly reactive ROS can damage the ultra structure

337

and function of chloroplasts, affecting PSII activity and photosynthetic pigments 22. In

338

plants acclimated to sunlight, ROS production is counterbalanced by inner defenses,

339

including antioxidant enzymes and antioxidant metabolites, such as phenolic

340

compounds. Some phenolics provide stress protection, for example, acting as

341

scavengers of ROS, as well as chelating metals that generate ROS via the Fenton

342

reaction 23. Similar mechanism of action that have been described about the beneficial

343

effects of those compounds as antioxidants in humans can be applied to plant

344

protections against excess light damages 10. In the present work, we found significantly

345

higher production of antioxidant polyphenols in late night and/or early morning.

346

Therefore, our results evidenced that Brassica plants do not simply respond to light but,

347

rather, anticipate the dawn and adjust their biology accordingly. Although it is unclear

348

whether the circadian clock is involved in the antioxidant plant defense system, the

349

existence of a master clock regulation behind this process it is suggested.

350

3.3. Endogenous circadian rhythms in antioxidant phenolic traits

351

Plants have adapted their growth and development to use the daily cycling of

352

light and dark. This is manifested at both the physiological level (i.e hypocotyl growth,

353

flowering time, stomatal opening and photosynthesis) and the molecular level with the

354

expression of some genes occurring only at certain times of the day24 . The day-night

355

cycling of these processes is called a diurnal rhythm and is achieved primarily by two

356

mechanisms: first, by light, and second, by a free-running internal circadian clock.

16 ACS Paragon Plus Environment

Page 16 of 26

Page 17 of 26

Journal of Agricultural and Food Chemistry

357

An attribute of circadian rhythms is that they are endogenously generated and

358

self-sustaining, so they persist under constant environmental conditions, typically

359

constant light and constant temperature. To date, most that is known about the effect of

360

the endogenous plant circadian clock in secondary metabolites is based on Arabidopsis,

361

and it is important to determine how the clock regulates plant phytochemicals,

362

especially in edible crop species. To address these questions, after an entrainment phase

363

of light-dark conditions we transferred the plants under constant light conditions “free

364

running period”. Samples from each cultivar were collected at the same time points

365

defined by the light-dark conditions during two consecutive days.

366

We found that after a period of continuous light, the antioxidant activity of both

367

assays in all studied crops showed a robust circadian rhythm (Figure 3). Thus, the

368

circadian parameters (period and phase) of compounds correlated with antioxidant

369

properties were studied. We also corroborated clear daily patterns in those compounds.

370

As expected for traits regulated by the circadian clock, the circadian period of the

371

studied individual phenolics was not exactly but close to 24h even in the absence of

372

external timing cues (Figure 4). It has generally been assumed that circadian clocks

373

have periods close to 24h in order to maintain a stable phase relationship to the Earth’s

374

rotational cycle. Proper matching of the internal circadian time with the environment

375

not only confers fitness advantages but also allows the clock to temporally gate the

376

responses to environmental stresses 25.

377

The phase of a circadian rhythm reflects where the maximum peak of

378

accumulation occurs within a period of 24h. Interestingly, in absence of environmental

379

time cues, all cultivars keep the rhythm and maintain the circadian phase (Figure 4).

380

Similarly to our previous results, Chinese cabbage and turnip greens showed a peak of 17 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

381

antioxidants accumulation during the subjective early morning. Broccoli and cabbage

382

still accumulated significantly more levels of antioxidant compounds when it would

383

have been night during the entrainment period (subjective night) even in absence of a

384

dark period. Thus, this phenolic phytochemicals maintain a time-of-day-specific phase

385

of accumulation even under free running conditions, implying a role of the circadian

386

clock in the regulation of these metabolites. These results are in concordance with the

387

endogenous circadian rhythmicity of antioxidative enzymes and low-molecular

388

antioxidants in different organisms, including plants

389

Chalcone Synthase (CHS), an important control step in the biosynthesis of flavonoids

390

showed circadian rhythm regulation

391

activity occurring in the late subjective night. It may be advantageous for the plant to

392

accumulate photoprotective pigments in advance of the daily photoperiod 11; the timing

393

of higher polyphenol production in our Brassica cultivars before and during dawn is

394

consistent with this notion.

395

In the present work, we reported for first time that the antioxidant activity of extracts

396

from different Brassica cultivars fluctuates rhythmically though the day. Interestingly,

397

those variations were related with endogenous circadian rhythms in polyphenol

398

composition and exhibit a species-specific pattern. This precise timing and amplitude of

399

daily phenolic variation patterns suggest a tight genetic control. Presumably this is

400

related with the fact that the antioxidant defense system's role in balancing free radicals

401

is highly essential for sustaining important plant physiological processes. The fact that

402

each Brassica cultivar grown under light-dark cycles have an optimal time during a

403

single day with increased levels of polyphenols, could be applied as a simple

404

mechanism to promote the nutritional value of Brassica-derived food. More detailed

26

. In Arabidopsis, the enzyme

27

. These authors found phase of peak enzyme

18 ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26

Journal of Agricultural and Food Chemistry

405

understanding of how precisely these rhythms vary under field conditions and how this

406

variation influences other phenotypes has broad implications for plant biology.

407

Funding

408

This research was supported by The Spanish Ministry of Economy and Competitiveness

409

through a ‘Juan de la Cierva-incorporación’ program (IJCI-2014-19653) to MF and by

410

the project AGL2015-66256-C2-R; and by the European Regional Development Fund

411

(ERDF) and by La Xunta de Galicia (IN607A 2016/13)

412

Conflict of interest

413

The authors declare that there is no conflict of interest regarding the publication of this

414

paper.

415

Supporting Information

416 417 418 419 420 421 422 423 424

Table S1. Abbreviations and identification data of phenolic compounds found in the studied Brassica cultivars (broccoli, cabbage, Chinese cabbage and turnip greens). Table S2. Phenolic variation (µmol g-1 DW) among the four times of the day (BD: beginning of the day light; ED: end of the day light; BN: beginning of the night; EN: end of the night) evaluated for each cultivar (broccoli, cabbage, Chinese cabbage and turnip greens). Samples were collected during two consecutive days in plants entrained to light-dark cycles (16h light/8h dark). Each cultivar within each time point is represented by the mean and standard deviation (SD) of six independent measurements.

425

426

19 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

427

References

428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475

1. Jahangir, M.; Kim, H. K.; Choi, Y. H.; Verpoorte, R., Health-Affecting Compounds in Brassicaceae. Comprehensive Reviews in Food Science and Food Safety 2009, 8 (2), 31-43. 2. Dinkova-Kostova, A. T.; Kostov, R. V., Glucosinolates and isothiocyanates in health and disease. Trends in Molecular Medicine 2012, 18 (6), 337-347. 3. Tripathy, B. C.; Oelmüller, R., Reactive oxygen species generation and signaling in plants. Plant Signaling & Behavior 2012, 7 (12), 1621-1633. 4. Cao, G.; Sofic, E.; Prior, R. L., Antioxidant Capacity of Tea and Common Vegetables. Journal of Agricultural and Food Chemistry 1996, 44 (11), 3426-3431. 5. Akula, R.; Ravishankar, G. A., Influence of abiotic stress signals on secondary metabolites in plants. Plant Signaling & Behavior 2011, 6 (11), 1720-1731. 6. Gerhardt, K. E.; Lampi, M. A.; Greenberg, B. M., The Effects of Far-red Light on Plant Growth and Flavonoid Accumulation in Brassica napus in the Presence of Ultraviolet B Radiation. Photochemistry and Photobiology 2008, 84 (6), 1445-1454. 7. Vale, A. P.; Cidade, H.; Pinto, M.; Oliveira, M. B. P. P., Effect of sprouting and light cycle on antioxidant activity of Brassica oleracea varieties. Food Chemistry 2014, 165 (Supplement C), 379-387. 8. Pérez-Balibrea, S.; Moreno, D. A.; García-Viguera, C., Influence of light on health-promoting phytochemicals of broccoli sprouts. Journal of the Science of Food and Agriculture 2008, 88 (5), 904-910. 9. Soengas, P.; Cartea, M. E.; Francisco, M.; Sotelo, T.; Velasco, P., New insights into antioxidant activity of Brassica crops. Food Chemistry 2012, 134 (2), 725-733. 10. Agati, G.; Azzarello, E.; Pollastri, S.; Tattini, M., Flavonoids as antioxidants in plants: Location and functional significance. Plant Science 2012, 196 (Supplement C), 67-76. 11. Harmer, S. L.; Hogenesch, J. B.; Straume, M.; Chang, H.-S.; Han, B.; Zhu, T.; Wang, X.; Kreps, J. A.; Kay, S. A., Orchestrated Transcription of Key Pathways in Arabidopsis by the Circadian Clock. Science 2000, 290 (5499), 2110. 12. Samarth, R. M.; Panwar, M.; Kumar, M.; Soni, A.; Kumar, M.; Kumar, A., Evaluation of antioxidant and radical-scavenging activities of certain radioprotective plant extracts. Food Chemistry 2008, 106 (2), 868-873. 13. Re, R.; Pellegrini, N.; Proteggente, A.; Pannala, A.; Yang, M.; Rice-Evans, C., Antioxidant activity applying an improved ABTS radical cation decolorization assay. Free Radical Biology and Medicine 1999, 26 (9), 1231-1237. 14. Benzie, I. F. F.; Strain, J. J., The Ferric Reducing Ability of Plasma (FRAP) as a Measure of “Antioxidant Power”: The FRAP Assay. Analytical Biochemistry 1996, 239 (1), 70-76. 15. Francisco, M.; Ali, M.; Ferreres, F.; Moreno, D. A.; Velasco, P.; Soengas, P., Organ-Specific Quantitative Genetics and Candidate Genes of Phenylpropanoid Metabolism in Brassica oleracea. Frontiers in Plant Science 2016, 6. 16. Francisco, M.; Moreno, D. A.; Elena Cartea, M.; Ferreres, F.; Garcia-Viguera, C.; Velasco, P., Simultaneous identification of glucosinolates and phenolic compounds in a representative collection of vegetable Brassica rapa. Journal of Chromatography a 2009, 1216 (38), 6611-6619. 17. Velasco, P.; Francisco, M.; Moreno, D. A.; Ferreres, F.; Garcia-Viguera, C.; Elena Cartea, M., Phytochemical Fingerprinting of Vegetable Brassica oleracea and Brassica napus by Simultaneous Identification of Glucosinolates and Phenolics. Phytochemical Analysis 2011, 22 (2), 144-152. 20 ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

Journal of Agricultural and Food Chemistry

476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504

18. Team, R. C., R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2014; Vol. URL http://www.Rproject.org/. 19. Moore, A.; Zielinski, T.; Millar, A. J., Online Period Estimation and Determination of Rhythmicity in Circadian Data, Using the BioDare Data Infrastructure. In Plant Circadian Networks: Methods and Protocols, Staiger, D., Ed. Springer New York: New York, NY, 2014; pp 13-44. 20. Cartea, M. E.; Francisco, M.; Soengas, P.; Velasco, P., Phenolic Compounds in Brassica Vegetables. Molecules 2011, 16 (1), 251-280. 21. Kusznierewicz, B.; Bartoszek, A.; Wolska, L.; Drzewiecki, J.; Gorinstein, S.; Namieśnik, J., Partial characterization of white cabbages (Brassica oleracea var. capitata f. alba) from different regions by glucosinolates, bioactive compounds, total antioxidant activities and proteins. LWT - Food Science and Technology 2008, 41 (1), 1-9. 22. Aro, E.-M.; Virgin, I.; Andersson, B., Photoinhibition of Photosystem II. Inactivation, protein damage and turnover. Biochimica et Biophysica Acta (BBA) Bioenergetics 1993, 1143 (2), 113-134. 23. Williams, R. J.; Spencer, J. P. E.; Rice-Evans, C., Flavonoids: antioxidants or signalling molecules? Free Radical Biology and Medicine 2004, 36 (7), 838-849. 24. Greenham, K.; McClung, C. R., Integrating circadian dynamics with physiological processes in plants. Nat Rev Genet 2015, 16 (10), 598-610. 25. Seo, P. J.; Mas, P., STRESSing the role of the plant circadian clock. Trends in Plant Science 2015, 20 (4), 230-237. 26. Lai, A. G.; Doherty, C. J.; Mueller-Roeber, B.; Kay, S. A.; Schippers, J. H. M.; Dijkwel, P. P., CIRCADIAN CLOCK-ASSOCIATED 1 regulates ROS homeostasis and oxidative stress responses. Proceedings of the National Academy of Sciences of the United States of America 2012, 109 (42), 17129-17134. 27. Thain, S. C.; Murtas, G.; Lynn, J. R.; McGrath, R. B.; Millar, A. J., The Circadian Clock That Controls Gene Expression in Arabidopsis Is Tissue Specific. Plant Physiology 2002, 130 (1), 102.

505

506

21 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 22 of 26

507

Table 1. ANOVA results comparing the antioxidant and phenolics response of four

508

Brassica cultivars (broccoli, cabbage, Chinese cabbage and turnip greens) entrained for

509

one month to light-dark cycles (16h light/8h dark) at four different time points of the

510

day. F-values, degrees of freedom (df), and associated significance levels (P) are shown.

511

Significant effects (P < 0.05) are in bold.

df

Cultivar 3 F P-value

Antioxidant assays 159.31 ABTS 162.36 FRAP Phenolic compounds 57.08 QA1 18.09 QA2 483.89 QA3 21.27 QA4 500.14 QA5 35.04 F1 22.19 F2 12.69 F3 7.22 F4 37.67 F5 57.81 F6 37.39 F7 52.92 F8 19.65 F9 28.61 F10 16.92 F11 30.88 F12 18.22 F13 136.82 F14 121.37 F15 7.69 Sinap1 75.16 Sinap2 387.71 Sinap3 176.83 Sinap4 10.25 Sinap5 21.41 Sinap6 133.79 Sinap7 113.40 Total_QA 104.57 Total_Flav 5.49 Total_Sinap Total_Phenolics 64.98

Time point Cultivar × Time point 3 9 F P-value F P-value

0.001 0.001

9.50 4.30

0.001 0.008

3.36 2.31

0.002 0.028

0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001

4.90 2.70 6.03 3.06 66.79 7.44 2.36 8.88 4.14 0.18 6.30 5.50 2.20 0.15 0.68 1.51 3.92 2.51 7.62 3.05 5.80 0.64 8.43 15.68 2.93 4.05 6.15 11.96 8.45 4.92 3.09

0.001 0.055 0.001 0.036 0.001 0.001 0.082 0.001 0.010 0.910 0.001 0.002 0.099 0.926 0.568 0.223 0.013 0.069 0.001 0.036 0.001 0.590 0.001 0.001 0.042 0.011 0.001 0.001 0.001 0.001 0.035

1.91 4.51 6.71 3.91 25.07 3.10 4.00 5.86 1.58 3.74 3.99 4.59 3.94 0.62 3.38 0.59 1.22 0.62 1.80 1.54 2.71 1.50 4.04 3.47 3.77 3.85 2.27 2.67 3.60 3.79 1.45

0.071 0.001 0.001 0.001 0.001 0.004 0.001 0.001 0.146 0.001 0.001 0.001 0.001 0.772 0.002 0.795 0.302 0.770 0.090 0.157 0.011 0.173 0.001 0.002 0.001 0.001 0.031 0.012 0.001 0.001 0.191

512 22 ACS Paragon Plus Environment

Page 23 of 26

Journal of Agricultural and Food Chemistry

513

Figure captions

514 515

Figure 1. Clustering dendrogram and heatmap of phenolic traits variation at different times of the day in Brassica cultivars.

516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538

Clustering of samples was based on phenolic variation among the four times of the day (BD: beginning of the day light (1 h after the lights turned on); ED: end of the day light (1h before lights turned off); BN: beginning of the night (1h after lights turned off) ; EN: end of the night (1h before lights turned on)). Samples were collected during two consecutive days in plants entrained to light-dark cycles (16h light/8h dark). Each cultivar within each time point is represented by the mean of six independent measurements. Samples were grouped into two clades (highlighted by a dotted line). Colors of the heatmap can be interpreted using the scale bar where red represents the lowest and yellow represents the highest relative level of phenolic content. Black squares highlight the significant highest levels of traits related with antioxidant potential by cultivar. See Table S1 for abbreviations and Table S2 for quantitative data.

539 540

Figure 3. Time series of ABTS and FRAP antioxidant activities in Brassica cultivars.

541 542 543 544 545 546 547 548

Plants of broccoli, cabbage, Chinese cabbage and turnip greens were grown in light-dark cycles (16h light/8h dark) for one month and then transferred into constant light conditions (starting at Time 0). Samples for antioxidant activity (µmol of Trolox g-1 of DW) were collected during two consecutive days at four time points per day, two times at subjective dark and two times at subjective light as were defined by the entrainment period of light-dark cycles. The error bars represent standard deviation across replicates. Each cultivar within each time point is represented by the mean of three independent measurements.

549 550

Figure 4. Circadian period and phase of the antioxidant phenolics presents in Brassica cultivars.

551 552 553 554 555 556 557 558

Circadian period (time to complete one cycle) and phase of maximal antioxidants accumulation (µmol g-1 DW) in plants of Brassica cultivars subjected to continuous light for 48h after an entrainment period of light-dark conditions (16h light/8h dark). Shown are the individual phenolic traits selected as significantly (P < 0.05) related with the antioxidant activity in the stepwise regression analysis. For the phase diagrams the circle denotes the period length of 24h. From 0 to 8h means subjective night and from 8 to 24h means subjective day. The length of the ellipse from the origin of coordinates is proportional to the circadian amplitude and its orientation gives the acrophase of the rhythm.

Figure 2. Daily oscillation of antioxidant activity in Brassica cultivars. Quantification of antioxidant activity (ABTS and FRAP) (µmol of Trolox g-1 of DW) from four Brassica cultivars (broccoli, cabbage, Chinese cabbage and turnip greens) entrained for one month to light-dark cycles (16h light/8h dark). Samples were collected during two consecutive days at four different time points defined by the light-dark conditions. BD: beginning of the day light (1 h after the lights turned on); ED: end of the day light (1h before lights turned off); BN: beginning of the night (1h after lights turned off) ; EN: end of the night (1h before lights turned on). Different letters indicate significant differences between time points (ANOVA, P < 0.05). The error bars represent standard deviation across replicates and days. Each cultivar within each time point is represented by the mean of six independent measurements.

559 23 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

560

Figure 1.

561 562 563

Figure 2.

564 565

24 ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

Journal of Agricultural and Food Chemistry

566

Figure 3

567 568 569 570

Figure 4

571 572

25 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

573

Graphical abstract

574

26 ACS Paragon Plus Environment

Page 26 of 26