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