Subscriber access provided by Lancaster University Library
Bioactive Constituents, Metabolites, and Functions
Biochemical and Molecular Factors Governing PeelColor Development in 'Ora' and 'Shani' Mandarins Livnat Goldenberg, Matat Zohar, Lina Kirshinbaum, Yossi Yaniv, Adi Doron-Faigenboim, Ron Porat, Nir Carmi, and Tal Isaacson J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b00669 • Publication Date (Web): 11 Apr 2019 Downloaded from http://pubs.acs.org on April 15, 2019
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 32
Journal of Agricultural and Food Chemistry
1
Biochemical and Molecular Factors Governing PeelColor Development in 'Ora' and 'Shani' Mandarins Livnat Goldenberg,†,§ Matat Zohar,¥ Lina Kirshinbaum,¶ Yossi Yaniv,‡ Adi DoronFaigenboim, ⱡ Ron Porat†,* Nir Carmi‡ and Tal Isaacson¥ †
Department of Postharvest Science of Fresh Produce, ARO, The Volcani Center, P.O. Box
15159, Rishon LeZion 7505101, Israel. § Faculty
of Agricultural, Food and Environmental Quality Sciences, Hebrew University of
Jerusalem, Rehovot 76100, Israel. ¥
Newe Ya'ar Research Center, ARO, Ramat Yishay 30095, Israel.
¶
Biology and DNA Laboratory, Division of Identification and Forensic Science, Israel Police,
Jerusalem, Israel. ‡
Department of Fruit Tree Crops, ARO, The Volcani Center, P.O. Box 15159, Rishon LeZion
7505101, Israel. ⱡ
Department of Genomics and Bioinformatics, ARO, the Volcani Center, P.O. Box 15159,
Rishon LeZion 7505101, Israel.
Corresponding author:
Dr. Ron Porat Tel.:
972-3-9683617
Fax:
972-3-9683622
E-mail:
[email protected] ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 32
2 1
ABSTRACT: To identify factors governing peel-color development in mandarins, we
2
examined carotenoid content and composition and the expression of carotenoid-related
3
genes during four stages of ripening (i.e., green, breaker, yellow and orange) in two
4
varieties: ‘Ora’, which has orange fruit, and ‘Shani’, which has orange-reddish fruit. The
5
two varieties had different carotenoid compositions and ‘Shani’ had a significantly higher
6
level of total carotenoid pigments. ‘Shani’ was rich in the deep orange β-cryptoxanthin and
7
the orange-reddish β-citraurin; whereas ‘Ora’ was rich in the orange violaxanthin. RNA-Seq
8
analysis revealed significantly greater expression of the carotenoid-biosynthesis genes PSY,
9
βLCY, βCHX and CCD4b, as well as MEP-pathway genes and several ethylene-biosynthesis
10
and -signaling genes in ‘Shani’ fruit. In contrast, the expression levels of genes involved in
11
the synthesis of α-branch carotenoids (i.e., εLCY and εCHX) and ZEP, which is involved in
12
the formation of violaxanthin, were significantly higher in the ‘Ora’ fruit.
13 14 15
KEYWORDS: mandarin, peel color, carotenoids, ripening, RNA-Seq.
ACS Paragon Plus Environment
Page 3 of 32
Journal of Agricultural and Food Chemistry
3 17
INTRODUCTION
18
To enhance the visual attractiveness of mandarins, citrus breeders are interested in developing
19
new varieties with peels that have a deep orange to reddish in color.1 The achievement of this
20
objective requires a better understanding of the biochemical and molecular factors that govern
21
the development of mandarin peel color and carotenoid composition.
22
The color of mandarin peels is mainly controlled by the levels and configuration of
23
their carotenoid pigments. Most of the genes of the carotenoid-biosynthesis pathway have
24
been identified in plants2,3 and in citrus particularly (Figure 1).4,5 Earlier studies have
25
classified citrus into three different groups based on their carotenoid profiles: the carotenoid-
26
poor group includes grapefruits, pomelos, lemons and limes; the violaxanthin-abundant group
27
is comprised mainly of oranges; and the β-cryptoxanthin-abundant group is comprised mainly
28
of mandarins.6–9 The most abundant pigments identified in mandarins were violaxanthin, β-
29
cryptoxanthin, and β-citraurin.4,10,11
30
Mandarins have a profound orange to reddish color and a unique carotenoid profile.4,10 In
31
addition to having large contents of pigments, mandarins own a unique mixture of carotenoid
32
and apocarotenoid pigments, and are especially rich in β-cryptoxanthin and β-citraurin. β-
33
cryptoxanthin imparts a deep orange color,6 and β-citraurin divulges a reddish color, even at
34
low concentrations.12–14 The key gene controlling β-citraurin biosynthesis is carotenoid
35
cleavage dioxygenase 4b (CCD4b) – it was recently cloned independently by two laboratories
36
in Spain and Japan.13,14
37
The differences in carotenoid profiles between mandarins and oranges were attributed
38
to the following factors: (1) mandarins have greater gene expression levels of phytoene
39
synthase (PSY), phytoene desaturase (PDS), ζ-carotene desaturase (ZDS) and the lycopene β-
40
cyclase (βLCY) carotenoid biosynthetic genes; (2) mandarins have poorer gene expression
41
levels of the β-ring hydroxylase (βCHX) and the zeaxanthin epoxidase (ZEP) xanthophyll-
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 32
4 42
synthesis genes; (3) mandarins have greater gene expression levels of the 9-cis-
43
epoxycarotenoid dioxygenase (NCED) violaxanthin catabolism gene; and (4) mandarins have
44
greater gene expression levels of the CCD4b which is the key β-citraurin biosynthetic
45
gene.9,11
46
In commercial practice, the peel color of early-season citrus fruit is artificially
47
enhanced through exposure to the ripening hormone ethylene.15 In Satsuma mandarins,
48
ethylene treatment was reported to increase β-cryptoxanthin levels in the peel by 90% and
49
more,16,17 and to induce the expression of the PSY, PDS, βLCY, βCHX and ZEP carotenoid
50
biosynthesis genes.17,18 This indicates that ethylene plays a crucial role in the regulation of
51
carotenoid biosynthesis in mandarins.
52
We recently explored the existence of a wide genetic diversity in fruit quality traits
53
among 46 different mandarin varieties that exhibited a wide range of peel colors, from green
54
to yellow, orange and red.19 In order to elucidate the biochemical and molecular factors
55
governing peel-color development in mandarins, we evaluated carotenoid content and
56
composition and the expression patterns of carotenoid-related genes during four ripening
57
stages (i.e., green, breaker, yellow and orange) in two mandarin varieties with distinctive peel
58
colors: 'Ora', which is typically orange, and 'Shani', which has a unique orange/reddish color
59
(Figure 2A).
60 61 62
MATERIALS AND METHODS Plant Material. Fruits of two mandarin varieties, ‘Ora’ and ‘Shani’ (C. reticulata
63
Blanco), were obtained from the Israeli citrus breeding collection at the Agricultural Research
64
Organization, the Volcani Center, Rishon LeZion, Israel. Fruits of each variety were
65
harvested from the same experimental orchard, during the 2016-17 ripening season, at four
66
ripening stages: green, breaker, yellow and full orange. The fruits’ outer peel tissue (flavedo)
ACS Paragon Plus Environment
Page 5 of 32
Journal of Agricultural and Food Chemistry
5 67
was separated by hand peeler, weighed, frozen in liquid nitrogen and kept at -80°C until
68
further analysis.
69
Color Measurements. Peel color was measured with a CR-310 Chroma Meter (Minolta,
70
Tokyo, Japan). The results were expressed in terms of lightness (0–100), chroma (0–100) and
71
hue angle (H°). Data points were means SE of 10 fruits.
72
Ethylene Measurements. Ethylene production was measured at the full-orange ripening
73
stage. Mandarins were placed in 2-L sealed jars, each equipped with a septum, and the jars
74
were kept sealed for 4 h before sampling. For each variety, we used three jars, which each
75
contained five fruits. Headspace samples were taken with a 10-mL syringe and were injected
76
into a gas chromatograph equipped with a Model 3300 flame ionization detector (GC-FID;
77
Varian, Walnut Creek, CA, USA).
78
Carotenoid Extraction. Peel pigments were extracted from 200 mg of frozen flavedo
79
tissue for each replicate. The tissue was grinded and suspended in 4 mL of extraction
80
solution (hexane:ethanol:acetone 2:1:1, v/v/v) including 0.1% (w/v) butylated
81
hydroxytoluene (BHT) which was added right before extraction, and gently mixed for 5 min.
82
For saponification, 1 mL of methanolic KOH 20% (w/v) was added and the samples were
83
vortexed and then gently mixed for 2 h. Saponified carotenoids were extracted by adding 4
84
mL of diethyl ether and 4 mL of an aqueous 12% (w/v) NaCl solution to each sample. The
85
hypophases were collected, washed with water and collected again. Water residues were
86
eliminated by the addition of sodium sulfate and filtration. The hypophases were then dried
87
under a stream of N2 and the lipid extracts were re-dissolved in 200–600 μL acetone
88
(depending on the color intensity of the sample) and 50-100 μL were injected to the HPLC.
89
All procedures were conducted under dim light to avoid structural change and degradation
90
of the carotenoids.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 32
6 91
High-Pressure Liquid Chromatography (HPLC) Analysis. Carotenoid analysis was
92
performed according to Ronen et al.20 In more details: A Waters HPLC system equipped
93
with a Waters 600 pump, a Waters PDA detector 996 and a Waters 717 plus auto-sampler
94
(Milford, MA, USA) was used. A Spherisorb ODS2 C18 column (Waters; 5 µm, 4.6 × 250
95
mm) coupled with a guard cartridge system (SecurityGuard™, Phenomenex, CA, USA),
96
was used for carotenoid separation at room temperature (~24ºC). A gradient was applied at a
97
constant flow of 1.6 mL min-1 with acetonitrile:water (9:1; v/v) and ethyl acetate as
98
described in previous work.21 Spectra of eluted HPLC solvent within a wavelength range of
99
250–600 nm were recorded and absorption peaks were analyzed using Empower software
100
(Waters). The linear limit of detection was estimated to be between 10–20 ng and 1.5–2 μg,
101
depending on the carotenoid compound. Carotenoids were identified by their absorption
102
spectra and retention times, and in most cases by authentic standards. The β-carotene
103
standard was obtained from Sigma-Aldrich, the β –citraurin and phytoene standards were
104
obtained from CaroteNature (Switzerland) and β-cryptoxanthin and lutein from
105
Extrasynthase (France). All carotenoid peaks were normalized22 to correct for their specific
106
mass extinction coeffiecients23 in relation to β-carotene (= 1), using xanthophylls (1), β-
107
cryptoxanthin (1.086), ζ-carotene (1.014), phytofluene (1.920) and phytoene (2.074). Total
108
carotenoid content was determined based on β-carotene and β –citraurin calibration curves
109
prepared with authentic standards. Total carotenoid is presented as μg g-1 Fresh Weight
110
(FW).
111
RNA Isolation, cDNA Library Construction and RNA-Seq. Total RNA was extracted
112
according to the CTAB protocol24 and each combination of variety and ripening stage
113
included three biological replicates, each containing flavedo tissues collected from two
114
different fruits. The RNA samples were treated with DNase (Invitrogen; Carlsbad, CA, USA)
115
according to the supplier’s instructions. RNA concentrations were determined with a
ACS Paragon Plus Environment
Page 7 of 32
Journal of Agricultural and Food Chemistry
7 116
NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and
117
RNA purity and integrity were further verified with a Model 2100 Total RNA BioAnalyzer
118
(Agilent Technologies; Santa Clara, CA, USA). Library preparation and sequencing were
119
performed at The Center for Genomic Technologies, The Alexander Silbelman Institute of
120
Life Sciences, The Hebrew University of Jerusalem, Israel. Twenty-four single-end RNA-Seq
121
libraries were prepared using Illumina NextSeq 500 and Trueseq RNA protocols.
122
Transcriptome Analysis. The resulting raw short-reads were subjected to a filtering and
123
cleaning procedure. First, the SortMeRNA tool was used to filter out rRNA.25 Then, the
124
FASTX Toolkit was used (http://hannonlab.cshl.edu/fastx_toolkit/ index.html, version
125
0.0.13.2) to trim read-end nucleotides with quality scores < 30, using fastq_quality_trimmer,
126
and remove reads with less than 70% base pairs with a quality score ≤ 30 using the
127
fastq_quality_filter. The Bowtie 2, version 2.1 alignment tool was used to map the cleaned
128
reads on the Citrus clementina v1.0 (clementine) genes extracted from Phytozome database
129
(https://phytozome.jgi.doe.gov/pz/portal.html).
130
The transcript quantification (the number of reads per gene) from the RNA-Seq data
131
was determined using the Bowtie 2 aligner, version 2.126 and the expectation-maximization
132
method (RSEM), by estimating maximum-likelihood expression levels.27 Differential
133
expression analysis was done with the DESeq2 R package in the R environment.28 Genes that
134
were more than two-fold differentially expressed with a false discovery rate (FDR) corrected
135
statistical significance of no more than 0.05 were considered differentially expressed.29
136
Following the Trinity protocol,30 expression normalization was conducted using the trimmed
137
means of M-values (TMM), after calculating the fragments per feature kilobase per million
138
reads mapped (FPKM).
139
Statistical Analysis. One-way analysis of variance (ANOVA) was applied to the chroma-
140
meter and carotenoid data using Microsoft Excel. Tukey’s HSD and Student’s t-test were used
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 32
8 141
to identify differences in carotenoid content among the different ripening stages and varieties.
142
Those tests were performed using JMP 13.0 software (SAS Institute; Cary, NC, USA).
143 144 145
RESULTS Peel Color. The development of the peel color of ‘Ora’ and ‘Shani’ mandarin fruits
146
during ripening is illustrated in Figure 2A. As can be seen in that figure, the two varieties
147
were similar in color at the breaker stage. But later on, ‘Ora’ turned yellow and then
148
orange; whereas ‘Shani’ turned orange and then became deep orange/reddish. Lightness
149
and chroma (saturation) indices gradually increased during ripening in both varieties.
150
‘Ora’ had significantly higher lightness values during all ripening stages and significantly
151
higher chroma values at the green and orange stages (Figure 2B). In contrast, the observed
152
hue-angle values gradually decreased during ripening in both varieties, but were
153
significantly higher in ‘Ora’ at the yellow and orange stages. The hue angles of the ‘Ora'’
154
and ‘Shani’ peels at full maturation were 59° and 42°, respectively (Figure 2B). The lower
155
hue angle values of the ‘Shani’ fruits indicate that those fruits had a deeper orange/reddish
156
color than the ‘Ora’ fruits.
157
Carotenoid Content and Composition. HPLC analysis of carotenoid content and
158
composition was performed on flavedo tissues from both varieties during the four ripening
159
stages. Total carotenoid content was more or less similar in the two varieties during the
160
first and second early-ripening stages. But, later on, total carotenoid content drastically
161
increased in ‘Shani’, as compared with ‘Ora’ (Figure 3). At the full-ripening stage, the
162
flavedos of the ‘Ora’ and ‘Shani’ fruits contained 136.8 and 375.3 μg g-1 FW carotenoids,
163
respectively (Figure 3).
164 165
Detailed analysis of the carotenoid content and composition of the mandarin flavedos revealed the presence of 11 dominant carotenoids including carotenes (α-
ACS Paragon Plus Environment
Page 9 of 32
Journal of Agricultural and Food Chemistry
9 166
carotene, β- carotene, ζ-carotene, phytofluene and phytoene), xanthophylls (neoxanthin,
167
violaxanthin, antheraxanthin, lutein, zeaxanthin and β-cryptoxanthin) and the
168
apocarotenoid β-citraurin (Table 1). Examination of carotenoid profiles during ripening
169
revealed a developmental shift in carotenoid production in both varieties, resulting in a
170
gradual decrease in the levels of the α-branch carotenoids α-carotene and lutein, as well as
171
gradual increases in the levels of the β-branch carotenoids violaxanthin and β-
172
cryptoxanthin (Figure 1, Table 1).
173
Detailed examination of carotenoid composition at full ripening (the orange stage)
174
revealed significant differences between ‘Ora’ and ‘Shani’ (Figure 4). In both varieties,
175
violaxanthin was the most common carotenoid. However, it accounted for 57.6% of the
176
total carotenoid content of the ‘Ora’ flavedos, but only 43.7% of the total carotenoid
177
content of the ‘Shani’ flavedos (Figure 4). The other two main compounds contributing to
178
the difference in the carotenoid compositions of the two varieties were β-cryptoxanthin
179
and β-citraurin, which are known to provide a deep orange-reddish color (Figure 4).14 The
180
absolute levels of β-cryptoxanthin and β-citraurin were significantly higher in ‘Shani’ (as
181
compared to ‘Ora’) at all ripening stages, and at full ripening were 5 and 8 times higher,
182
respectively (Table 1). β-cryptoxanthin accounted for 5.9% and 11.2% of the total
183
carotenoid content of ‘Ora’ and ‘Shani’, respectively, and β-citraurin accounted for 5.4%
184
and 15.1% of the total carotenoid content of ‘Ora’ and ‘Shani’, respectively (Figure 4).
185
Differentially Expressed Genes. In order to elucidate the molecular factors governing
186
peel-color development in ‘Ora’ and ‘Shani’ mandarins, we performed RNA-Seq analysis of
187
flavedo from both varieties during all four ripening stages. Each of the cDNA libraries
188
yielded between 31.2 million and 46.2 million clean single-end reads, each 125 bp in length.
189
About 85% of the cleaned reads could be mapped to the Citrus Clementina genome
190
reference. In order to identify differentially expressed genes among the different time
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 32
10 191
points, we performed pair-wise comparisons and identified 5,889 and 8,698 differentially
192
expressed transcripts in ‘Ora’ and ‘Shani’, respectively, of which 4,774 transcripts were
193
common to both varieties (Figure 5). It is worth noting that we detected a much larger
194
amount of unique differentially expressed transcripts in ‘Shani’ (3,924) as compared with
195
‘Ora’ (just 1,115 transcripts) (Figure 5). Elaborated clustering analysis of the different gene
196
expression patterns of each variety during the different ripening stages is provided in
197
Supplementary Figure 1.
198
Further analysis of gene-expression patterns focused on three main pathways that
199
affect carotenoid accumulation: (i) the carotenoid-biosynthesis pathway, (ii) the
200
methylerythritol 4-phosphate (MEP) pathway, and (iii) ethylene biosynthesis and ethylene
201
signaling. Overall, the expression patterns of more than 120 different genes were examined
202
for possible differences between ‘Ora’ and ‘Shani’ during all ripening stages, with special
203
emphasis on genes that are differentially expressed during the third (yellow) and fourth
204
(orange) ripening stages.
205
With regard to genes involved in carotenoid biosynthesis, the expression levels of
206
PSY1 and CCD4b1 were significantly higher in ‘Shani’ than in ‘Ora’ at the full-maturation,
207
orange stage and the expression patterns of βLCY and βCHX were significantly higher in
208
‘Shani’ at the yellow and orange stages (Figure 6A). In contrast, the expression levels of the
209
violaxanthin biosynthetic gene ZEP and the α-branch genes lycopene ε-cyclase (εLCY) and ε-
210
carotene hydroxylase (εCHX) were significantly higher in ‘Ora’ than in ‘Shani’ at some stages
211
(Figure 1, Figure 6B).
212
Regarding the upstream MEP-pathway genes that provide the precursors for
213
carotenoid biosynthesis, it was found that expression levels of most MEP-pathway genes [i.e.,
214
DOXP synthase (DXS), DOXP reductase (DXR), CDP-ME kinase (CMK), HMB-PP reductase
215
(HDR) and DMAPP isomerase (IDI)] were significantly higher in ‘Shani’ flavedos than in
ACS Paragon Plus Environment
Page 11 of 32
Journal of Agricultural and Food Chemistry
11 216
‘Ora’ flavedos over most of the ripening stages and particularly during the later yellow and
217
orange stages (Figure 7).
218
Regarding ethylene biosynthetic and signaling, the transcript levels of the ethylene
219
biosynthetic gene ACC oxidase (ACO), the ethylene receptor 2 (ETR2) gene and the ethylene
220
response factor 1 (ERF1) gene were significantly higher in ‘Shani’ than in ‘Ora’ (Figure 8A).
221
The expression level of ACO was significantly higher at the breaker and orange stages, the
222
expression level of ETR2 was significantly higher at the yellow and orange stages and the
223
expression pattern of ERF1 was significantly higher in ‘Shani’ fruits, as compared with ‘Ora’
224
fruits, at all ripening stages (Figure 8A). Furthermore, ethylene production levels in the ripe
225
fruit (orange stage) were significantly and nearly four times higher for ‘Shani’ ripe fruits, as
226
compared with ‘Ora’ fruits (Figure 8B).
227 228
DISCUSSION
229
The color of mandarin peel is predominantly determined by its carotenoid content and
230
composition.4,10 Therefore, the development of new varieties with an attractive deep
231
orange/reddish color requires a better understanding of the unique carotenoid composition of
232
orange/reddish fruit and the molecular factors that control color development.5 In the present
233
study, we explored carotenoid content and composition, as well as the expression patterns of
234
genes related to carotenoid biosynthesis and regulation, in two mandarin varieties that differ
235
in their peel color: ‘Ora’ has an orange peel and ‘Shani’ has a deep orange/reddish peel
236
(Figure 2).
237
The key findings of the present study regarding the observed differences in carotenoid
238
content and composition between ‘Ora’ and ‘Shani’ are as follows: 1) the flavedos of ‘Shani’
239
fruits have a significantly higher level of total carotenoid pigments and 2) the flavedos of
240
‘Shani’ fruits have a different carotenoid composition than the flavedos of ‘Ora’ fruits and are
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 32
12 241
richer in the deep orange β-cryptoxanthin and the orange-reddish β-citraurin pigments than
242
the flavedo of ‘Ora’ fruits, which is richer in the orange violaxanthin pigment (Figures 3-4).
243
As so, the 'Shani' variety can be designated as β-cryptoxantin-abundant, whereas ‘Ora’ can be
244
designated as violaxanthin-abundant.6,9 Higher carotenoid content and a composition favoring
245
β-cryptoxantin and β-citraurin on the expense of violaxanthin are probably the main reasons
246
for the markedly higher intensity of peel color in ripe ‘Shani’ flavedo in comparison to ripe
247
‘Ora’ flavedo.
248
Transcriptional regulation of carotenoid biosynthesis genes was previously shown to
249
be a key factor determining carotenoid biosynthesis in many fruit such as tomato20,31,32 ,
250
pepper33 and citrus9 . Our results suggest that the different carotenoid profiles of ‘Ora’ and
251
‘Shani’ flavedos might be due to different patterns of gene expression, where elevated
252
transcription of PSY and of the β-branch carotenoid-biosynthetic genes βLCY and βCHX in
253
‘Shani’, during fruit ripening (Figure 6), leads to increased production of β-branch carotenoid
254
and higher total carotenoid in comparison to ‘Ora’ (Figure 1). Similar association between
255
gene expression and elevated levels of β-cryptoxantin were shown in the citrus cultivar
256
‘Seinannohikari’.34 The decreased expression levels of the εLCY and εCHX in ‘Shani’ (Figure
257
6) could explain the observed reduction in the α-branch carotenoid lutein during ripening
258
(Figure 1). The difference in β-cryptoxantin accumulation between the two mandarin varieties
259
is similar, to some extent, to the difference between flavedo of oranges (low β-cryptoxantin )
260
and mandarins (high- β-cryptoxantin). Kato et al.11 suggested that the high β-cryptoxantin
261
accumulation in mandarins, in contrast to oranges, is a result of up regulation of carotenoid
262
biosynthesis genes upstream to β-cryptoxantin and down regulation of genes down stream to
263
β-cryptoxantin.11 In a similar way, the lower expression of ZEP (Figure 6), responsible for
264
violaxanthin biosynthesis, in ‘Shani’ in comparison to ‘Ora’, could explain the accumulation
265
of β-cryptoxantin in ‘Shani’, while in ‘Ora’ the biosynthetic flux proceeds to violaxanthin. In
ACS Paragon Plus Environment
Page 13 of 32
Journal of Agricultural and Food Chemistry
13 266
addition, the higher expression of CCD4b, responsible for the production of β-citraurin13,14
267
can explain the higher accumulation of β-citraurin in ‘Shani’ in comparison to ‘Ora’ in the
268
flavedo of the ripe fruit (Figure 1).
269
Carotenoids are derived from isoprenoids which are synthesized by the MEP pathway
270
.35 Regulation of transcription of genes encoding the enzymes of the MEP pathway was
271
shown to influence carotenoid biosynthesis.36,37 The elevated transcription levels of the Mep-
272
pathway genes in ‘Shani’ in comparison to ‘Ora’ (Figure 7) suggest that in ‘Shani’ the
273
pathway is more active providing higher amounts of precursors to the carotenoid biosynthesis
274
pathway, and thus contributing to the higher total amount of carotenoids in flavedos of the
275
‘Shani’ variety.
276
Together, these data demonstrate that in the flavedo of ‘Shani’ variety, there seems to
277
be a high metabolic flux toward the biosynthesis of the orange/reddish pigments β-
278
cryptoxanthin and β-citraurin; whereas in ‘Ora’, there seems to be an attenuated metabolic
279
flux toward the biosynthesis of violaxanthin. It is unclear why the transcription of the
280
described sets of genes is so different between ‘Shani’ and ‘Ora’. Ethylene is known to induce
281
carotenoid biosynthesis in citrus fruit.16–18,38,39 Our finding that the transcription of ethylene
282
biosynthesis and ethylene signaling genes is higher in ‘Shani’ in comparison to ‘Ora’ (Figure
283
8A), in combination with the detection of much higher levels of ethylene in ’Shani’ fruit than
284
in ‘Ora’ fruit, suggest that in ‘Shani’ fruit the ripening process is more intensive, leading to
285
higher ethylene biosynthesis and more rapid biosynthesis of carotenoid precursors and
286
carotenoid compounds. The larger number of genes that are differentially expressed during
287
ripening in ‘Shani’ fruit as compared with ‘Ora’ (5,889 and 8,698 differentially expressed
288
transcripts respectively), further supports this hypothesis. It is interesting to note that the shift
289
in carotenoid composition toward higher levels of β-cryptoxanthin and β-citraurin in peels of
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 32
14 290
‘Shani’ in comparison to ‘Ora’ is similar to results obtained by application of ethylene to
291
peels of Gannan Newhall navel orange39 and Satsuma mandarin.17
292
Overall, our findings suggest that the development of orange/reddish peel color in
293
‘Shani’ mandarins does not result solely from the accumulation of any specific carotenoid
294
compound or from the enhanced expression of any specific gene, but rather is the cumulative
295
result of various changes in the expression patterns of genes related to carotenoid metabolism,
296
the MEP pathway and, presumably, also ethylene regulation.
297 298
FUNDING
299
This research was supported by Research Grant No. 203-1058-18 from the Chief Scientist of
300
the Israel Ministry of Agriculture & Rural Development.
301 302
DESCRIPTION OF SUPPORTING INFORMATION
303
Supplementary Table 1 – GenBank accession numbers and annotations of all of the
304
differentially expressed transcripts mentioned in this study.
305
Supplementary Figure 1 - Clustering analysis of DEGs in 'Ora' and 'Shani' varieties during the
306
green, breaker, yellow and orange stages of ripening.
307
ACS Paragon Plus Environment
Page 15 of 32
Journal of Agricultural and Food Chemistry
15 308
REFERENCES
309
(1)
fruit quality and drive consumer demand. Acta Hort. 2016, 1127, 199–202.
310 311
(2)
(3)
Yuan, H.; Zhang, J.; Nageswaran, D.; Li, L. Carotenoid metabolism and regulation in horticultural crops. Hort. Res. 2015, 2, 15036.
314 315
Nisar, N.; Li, L.; Lu, S.; Khin, N. C.; Pogson, B. J. Carotenoid metabolism in plants. Carotenoid Metab. Plants. Mol. Plant 2015, 8, 68–82.
312 313
Gmitter, F. G.; Chen, C.; Wei, X.; Yu, Y.; Yu, Q. New genetic tools to improve citrus
(4)
Alquézar, B.; Jesús Rodrigo, M.; Zacarías, L. Carotenoid biosynthesis and their regulation in citrus fruits. Tree For. Sci. Biotechnol. 2008, 2, 23–35.
316 317
(5)
Talon, M.; Gmitter, F. G. Citrus genomics. Int. J. Plant Genom. 2008, 1–17.
318
(6)
Goodner, K. L.; Rouseff, R. L.; Hofsommer, H. J. Orange, mandarin, and hybrid
319
classification using multivariate statistics based on carotenoid profiles. J. Agric. Food
320
Chem. 2001, 49, 1146–1150.
321
(7)
Fanciullino, A. L.; Dhuique-Mayer, C.; Luro, F.; Casanova, J.; Morillon, R.; Ollitrault,
322
P. Carotenoid diversity in cultivated citrus is highly influenced by genetic factors. J.
323
Agric. Food Chem. 2006, 54, 4397–4406.
324
(8)
Matsumoto, H.; Ikoma, Y.; Kato, M.; Kuniga, T.; Nakajima, N.; Yoshida, T.
325
Quantification of carotenoids in citrus fruit by LC-MS and comparison of patterns of
326
seasonal changes for carotenoids among citrus varieties. J. Agric. Food Chem. 2007,
327
55, 2356–2368.
328
(9)
Ikoma, Y.; Matsumoto, H.; Kato, M. Diversity in the carotenoid profiles and the
329
expression of genes related to carotenoid accumulation among citrus genotypes. Breed.
330
Sci. 2016, 66, 139–147.
331
(10)
Gross, J. Pigments in Fruits; Academic Press: London, UK and Orlando, FL, 1987.
332
(11)
Kato, M.; Ikoma, Y.; Matsumoto, H.; Sugiura, M.; Hyodo, H.; Yano, M. Accumulation
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 32
16 333
of carotenoids and expression of carotenoid biosynthetic genes during maturation in
334
citrus fruit. Plant Physiol. 2004, 134, 824–837.
335
(12)
Farin, D.; Ikan, R.; Gross, J. The carotenoid pigments in the juice and flavedo of
336
mandarin hybrid (Citrus reticulata) cv. Michal during ripening. Phytochemistry 1983,
337
22, 403–408.
338
(13)
Ma, G.; Zhang, L.; Matsuta, A.; Matsutani, K.; Yamawaki, K.; Yahata, M.; Wahyudi,
339
A.; Motohashi, R.; Kato, M. Enzymatic formation of β-citraurin from β-cryptoxanthin
340
and zeaxanthin by carotenoid cleavage dioxygenase4 in the flavedo of citrus fruit.
341
Plant Physiol. 2013, 163, 682–695.
342
(14)
Rodrigo, M. J.; Alquézar, B.; Alós, E.; Medina, V.; Carmona, L.; Bruno, M.; Al-Babili,
343
S.; Zacarías, L. A novel carotenoid cleavage activity involved in the biosynthesis of
344
citrus fruit-specific apocarotenoid pigments. J. Exp. Bot. 2013, 64, 4461–4478.
345
(15)
Wardowski, W. F.; Miller, W. M. Grierson W. Degreening. In Fresh Citrus Fruit;
346
Wardowski, W. F., Miller, W. M., Grierson W., Eds.; Florida Science Source:
347
Longboat Key, FL, 2006; pp. 277–298.
348
(16)
Fujii, H.; Shimada, T.; Sugiyama, A.; Nishikawa, F.; Endo, T.; Nakano, M.; Ikoma, Y.;
349
Shimizu, T.; Omura, M. Profiling ethylene-responsive genes in mature mandarin fruit
350
using a citrus 22K oligoarray. Plant Sci. 2007, 173, 340–348.
351
(17)
Ma, G.; Zhang, L.; Kato, M.; Yamawaki, K.; Kiriiwa, Y.; Yahata, M.; Ikoma, Y.;
352
Matsumoto, H. Effect of the combination of ethylene and red LED light irradiation on
353
carotenoid accumulation and carotenogenic gene expression in the flavedo of citrus
354
fruit. Postharvest Biol. Technol. 2015, 99, 99-104.
355
(18)
Matsumoto, H.; Ikoma, Y.; Kato, M.; Nakajima, N.; Hasegawa, Y. Effect of
356
postharvest temperature and ethylene on carotenoid accumulation in the flavedo and
357
juice sacs of Satsuma mandarin (Citrus unshiu Marc.) fruit. J. Agric. Food Chem.
ACS Paragon Plus Environment
Page 17 of 32
Journal of Agricultural and Food Chemistry
17
2009, 57, 4724–4732.
358 359
(19)
Goldenberg, L.; Yaniv, Y.; Kaplunov, T.; Doron-Faigenboim, A.; Porat, R.; Carmi, N.
360
Genetic diversity among mandarins in fruit-quality traits. J. Agric. Food Chem. 2014,
361
62, 4938–4946.
362
(20)
Ronen, G.; Cohen, M.; Zamir, D.; Hirschberg, J. Regulation of carotenoid biosynthesis
363
during tomato fruit development: expression of the gene for lycopene epsilon-cyclase
364
is down-regulated during ripening and is elevated in the mutantDelta. Plant J. 1999,
365
17, 341–351.
366
(21)
Isaacson, T.; Ohad, I.; Beyer, P.; Hirschberg, J. Analysis in vitro of the enzyme
367
CRTISO establishes a poly-cis-carotenoid biosynthesis pathway in plants 1. Plant
368
Physiol. 2004, 136, 4246–4255.
369
(22)
Shemesh, K.; Zohar, M.; Bar-Ya’akov, I.; Hatib, K.; Holland, D.; Isaacson, T. Analysis
370
of carotenoids in fruit of different apricot accessions reveals large variability and
371
highlights apricot as a rich source of phytoene and phytofluene. Fruits 2017, 72, 185–
372
202.
373
(23)
Handbook; Birkhäuser Verlag: Basel, Switzerland, 2004.
374 375
(24)
(25)
(26)
382
Langmead, B.; Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359.
380 381
Kopylova, E.; Noé, L.; Lè Ne Touzet, H. SortMeRNA: Fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 2012, 28, 3211–3217.
378 379
Chang, S.; Puryear, J.; Cairney, J. A simple and efficient method for isolating RNA from pine trees. Plant Mol. Biol. Report. 1993, 11, 113–116.
376 377
Britton, G.; Liaaen-Jensen, S. (Synnøve); Pfander, H. (Hanspeter). Carotenoids
(27)
Li, B.; Dewey, C. N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 2011, 12, 323.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 32
18 383
(28)
for RNA-Seq data with DESeq2. Genome Biol. 2014, 15, 550.
384 385
(29)
Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300.
386 387
Love, M. I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion
(30)
Haas, B. J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P. D.; Bowden, J.;
388
Couger, M. B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence
389
reconstruction from RNA-Seq using the Trinity platform for reference generation and
390
analysis. Nat. Protoc. 2013, 8, 1494–1512.
391
(31)
Fraser, P. D.; Truesdale, M. R.; Bird, C. R.; Schuch, W.; Bramley, P. M. Carotenoid
392
biosynthesis during tomato fruit development, evidence for tissue-specific gene
393
expression. Plant Physiol. 1994, 105, 405–413.
394
(32)
Ronen, G.; Carmel-Goren, L.; Zamir, D.; Hirschberg, J. An alternative pathway to β-
395
carotene formation in plant chromoplasts discovered by map-based cloning of beta and
396
old-gold color mutations in tomato. PNAS 2000, 97, 11102–11107.
397
(33)
Hugueney, P.; Bouvier, F.; Badillo, A.; Quennemet, J.; D’harlingue, A.; Camara, B.
398
Developmental and stress regulation of gene expression for plastid and cytosolic
399
lsoprenoid pathways in pepper fruits. Plant Physiol. 1996, 11, 61–70.
400
(34)
Ma, G.; Zhang, L.; Yungyuen, W.; Sato, Y.; Furuya, T.; Yahata, M.; Yamawaki, K.;
401
Kato, M. Accumulation of Carotenoids in a Novel Citrus Cultivar “Seinannohikari”
402
during the Fruit Maturation. Plant Physiol. Biochem. 2018, 129, 349–356.
403
(35)
Arch. Biochem. Biophys. 2010, 504, 118–122.
404 405
Rodríguez-Concepción, M. Supply of precursors for carotenoid biosynthesis in plants.
(36)
Lois, L.M.; Rodríguez‐Concepción, M.; Gallego, F.; Campos, N.; Boronat, A.
406
Carotenoid biosynthesis during tomato fruit development: regulatory role of
407
1‐deoxy‐D‐xylulose 5‐phosphate synthase. Plant J, 2000, 22, 503-513.
ACS Paragon Plus Environment
Page 19 of 32
Journal of Agricultural and Food Chemistry
19 408
(37)
Walter, M.H.; Hans, J.; Strack, D. Two distantly related genes encoding 1-deoxy-d-
409
xylulose 5-phosphate synthases: differential regulation in shoots and apocarotenoid-
410
accumulating mycorrhizal roots. Plant J. 2002, 31, 243-254.
411
(38)
Rodrigo, M. J.; Zacarias, L. Effect of postharvest ethylene treatment on carotenoid
412
accumulation and the expression of carotenoid biosynthetic genes in the flavedo of
413
orange (Citrus sinensis L. Osbeck) fruit. Postharvest Biol. Technol. 2007, 43, 14–22.
414
(39)
Hu, Y.; Wang, G.; Pan, S.; Wang, L. Influence of ethylene and ethephon treatments on
415
the peel color and carotenoids of Gannan Newhall navel orange during postharvest
416
storage. J. Food Biochem. 2018, 42, e12534.
417
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 32
1 418
FIGURE LEGENDS
419
Figure 1. The carotenoid-biosynthesis pathway. GGPP, Geranygeranyl pyrophosphate; PSY,
420
phytoene synthase; PDS, phytoene desaturase; ZISO, ζ-carotene isomerase; ZDS, ζ-carotene
421
desaturase; CRITSO, carotenoid isomerase; β-LCY, lycopene β-cyclase; ε-LCY, lycopene ε-
422
cyclase; β-CHX, β-ring carotene hydroxylase; ε-CHX, ε-ring carotene hydroxylase; ZEP,
423
zeaxanthin epoxidase; CCD4, carotenoid cleavage dioxygenase 4. The schematic diagram was
424
adapted from Alquézar et al. 4 with modifications.
425 426
Figure 2. Peel color of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange
427
stages of ripening. A) Photographs of ‘Ora’ and ‘Shani’ mandarin fruits at the green, breaker,
428
yellow and orange stages of ripening. B) Lightness (0–100), Chroma (0–100) and Hue angles
429
(0° represents a red color, 45° represents orange, 90° represents yellow and 120° represents
430
green) of ‘Ora’ and ‘Shani’ peels at the green, breaker, yellow and orange stages of ripening.
431
Data are means SE of 5 fruits and an asterisk indicates a significant difference (P 0.05)
432
between ‘Ora’ and ‘Shani’ fruits at the same stage of ripening.
433 434
Figure 3. Total carotenoid contents of flavedo of ‘Ora’ and ‘Shani’ fruits at the green,
435
breaker, yellow and orange stages of ripening. Data are means SE of 3 fruits. Different
436
capital letters indicate a significant difference among the ‘Ora’ fruits at the different ripening
437
stages. Different lower-case letters indicate a significant difference among the ‘Shani’ fruits at
438
the different ripening stages. An asterisk indicates a significant difference (P 0.05) between
439
‘Ora’ and ‘Shani’ fruits at a particular ripening stage.
440 441
Figure 4. Carotenoid composition in flavedo of ‘Ora’ and ‘Shani’ fruits at the orange (ripe)
442
stage of ripening. Numbers indicate the percent if each carotenoid of the total carotenoid
1 ACS Paragon Plus Environment
Page 21 of 32
Journal of Agricultural and Food Chemistry
2 443
content for each variety. Data points are means of 3 fruits. β-xanthopylls – neoxanthin,
444
antheraxanthin and zeaxanthin. Carotenes – α, β and ζ carotene, phytofluene and phytoene.
445 446
Figure 5. Venn diagram illustration of differentially expressed genes (DEGs) in flavedo of
447
‘Ora’ and ‘Shani’ during the green, breaker, yellow and orange stages of ripening. The
448
diagram displays the number of DEGs specific to ‘Ora’, the number of DEGs common to
449
‘Ora’ and ‘Shani’ and the number of DEGs specific to ‘Shani’.
450 451
Figure 6. Carotenoid-biosynthesis genes differentially expressed in flavedo of ‘Ora’ and
452
‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. A) Genes highly
453
expressed in ‘Shani’. B) Genes highly expressed in ‘Ora’. The results are expressed as
454
fragments per feature kilobase per million reads mapped (FPKM). An asterisk indicates a
455
significant differences (FDR 0.05) between ‘Ora’ and ‘Shani’ fruits at a particular ripening
456
stage, according to the DESeq2 analysis.
457 458
Figure 7. MEP-pathway genes differentially expressed in flavedo of ‘Ora’ and ‘Shani’ fruits
459
at the green, breaker, yellow and orange stages of ripening. The results are expressed as
460
fragments per feature kilobase per million reads mapped (FPKM). An asterisk indicates a
461
significant difference (FDR 0.05) between ‘Ora’ and ‘Shani’ fruits at the same ripening
462
stage, according to the DESeq2 analysis.
463 464
Figure 8. Ripening regulation in flavedo of ‘Ora’ and ‘Shani’ mandarins. A) Ethylene-
465
biosynthesis and ethylene-signaling genes differentially expressed in flavedo of ‘Ora’ and
466
‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. Results are
467
expressed as fragments per feature kilobase per million reads mapped (FPKM). An asterisk
468
indicates a significant difference (P 0.05) between ‘Ora’ and ‘Shani’ fruits at the same 2 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 22 of 32
3 469
ripening stage, according to the DESeq2 analysis. B) Ethylene production of ‘Ora’ and
470
‘Shani’ at the orange ripening stage. Data are means ± SE of 3 replications, each containing 5
471
fruit. An asterisk indicates a significant difference (P 0.05) between ‘Ora’ and ‘Shani’
472
fruits.
473
3 ACS Paragon Plus Environment
Page 23 of 32
Journal of Agricultural and Food Chemistry
4
Table 1. Carotenoid content of flavedo in ‘Ora’ and ‘Shani’ mandarinthe at s green, breaker, yellow and orange stages of ripening. Carotenoid content (μg g-1) Green
Breaker
Yellow
Orange
Ora
Shani
Ora
Shani
Ora
Shani
Ora
Shani
neoxanthin
5.9 ± 1.6 A
9.1 ± 2.0 a
2.4 ± 0.5 A
3.1 ± 0.8 a
3.6 ± 0.9 A
5.8 ± 1.3 a
4.3 ± 1.3 A
10.3 ± 4.6 a
violaxanthin
25.1 ± 4.6 B
35.8 ± 6.7 b
29.2 ± 6.0 B
20.3 ± 3.2 b
31.4 ± 9.4 B
78.8 ± 11.6 A
164.2 ± 42.4 a
6.2 ± 1.0 *b
1.8 ± 0.9 B
4.0 ± 0.5 b
2.4 ± 1.2 B
21.3 ± 4.5 *b
7.4 ± 1.2 A
56.5 ± 13.8 *a
5.0 ± 1.6 A
8.3 ± 0.9 a
4.0 ± 1.3 A
3.2 ± 0.5 a
7.9 ± 2.7 A
9.0 ± 1.4 a
10.3 ± 2.2 A
19.5 ± 7.2 a
15.5 ± 3.2 A
18.0 ± 3.9 a
3.6 ± 2.1 B
9.9 ± 3.6 a
6.5 ± 3.0 AB
12.0 ± 6.3 a
4.5 ± 0.6 AB
10.0 ± 1.6 *a
zeaxanthin
2.6 ± 0.7 A
2.6 ± 0.3 a
1.1 ± 0.3 A
2.3 ± 0.6 a
2.5 ± 1.3 A
3.4 ± 1.5 a
1.7 ± 0.1 A
3.7 ± 0.7 a
β-cryptoxanthin
1.6 ± 0.3 B
6.7 ± 1.0 *a
2.0 ± 0.6 B
5.6 ± 1.7 a
2.4 ± 0.3 B
27.8 ± 9.0 *a
8.1 ± 1.4 A
42.1 ± 13.7 *a
α-carotene
4.5 ± 0.7 A
4.9 ± 1.0 a
tr B
2.1 ± 0.7 ab
1.3 ± 0.2 B
nd *b
1.6 ± 0.1 B
nd *b
β-carotene
6.7 ± 1.3 A
8.5 ± 1.6 a
1.5 ± 0.5 B
2.8 ± 1.0 b
1.7 ± 0.3 B
2.6 ± 0.5 b
1.7 ± 0.0 B
3.5 ± 0.7 ab
nd C
nd b
tr C
tr ab
1.0 ± 0.0 B
2.7 ± 1.1 ab
1.6 ± 0.1 A
3.5 ± 0.8 a
phytofluene
1.8 ± 0.6 A
1.9 ± 0.5 a
1.3 ± 0.2 A
1.7 ± 0.9 a
1.6 ± 0.1 A
5.7 ± 2.0 a
2.2 ± 0.2 A
7.2 ± 1.1 *a
phytoene
2.9 ± 1.0 A
3.8 ± 0.7 a
1.9 ± 0.5 A
3.1 ± 1.8 a
2.4 ± 0.5 A
9.6 ± 4.4 a
2.1 ± 1.0 A
9.3 ± 1.3 *a
others
4.0 ± 1.5 B
7.1 ± 1.3 b
4.7 ± 1.8 AB
7.5 ± 2.8 b
3.6 ± 0.8 B
30.4 ± 3.8 *a
12.6 ± 2.7 A
45.6 ± 5.5 *a
Compound
β-citraurin
nd B
antheraxanthin lutein
ζ-carotene
81.6 ± 12.5 *ab
Data are means SE of 3 fruits. Different capital letters indicate a significant difference in ‘Ora’ among different ripening stages, different lower-case letters indicate significant 4 difference in ‘Shani’ among different ripening stages, and an asterisk indicates a significant differences (P 0.05) between ‘Ora’ and ‘Shani’ at the same ripening stage. nd – ACS Paragon Plus Environment not detected; tr – traces.
Journal of Agricultural and Food Chemistry
Page 24 of 32
5
FIGURES
Figure 1. The carotenoid-biosynthesis pathway. GGPP, Geranygeranyl pyrophosphate; PSY, phytoene synthase; PDS, phytoene desaturase; ZISO, ζ-carotene isomerase; ZDS, ζ-carotene desaturase; CRITSO, carotenoid isomerase; β-LCY, lycopene β-cyclase; ε-LCY, lycopene εcyclase; β-CHX, β-ring carotene hydroxylase; ε-CHX, ε-ring carotene hydroxylase; ZEP, zeaxanthin epoxidase; CCD4, carotenoid cleavage dioxygenase 4. The schematic diagram was adapted from Alquézar et al.4 with modifications.5 ACS Paragon Plus Environment
Page 25 of 32
Journal of Agricultural and Food Chemistry
6
Figure 2. Peel color of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. A) Photographs of ‘Ora’ and ‘Shani’ mandarin fruits at the green, breaker, yellow and orange stages of ripening. B) Lightness (0–100), Chroma (0–100) and Hue angles (0° represents a red color, 45° represents orange, 90° represents yellow and 120° represents green) of ‘Ora’ and ‘Shani’ peels at the green, breaker, yellow and orange stages of ripening. Data are means ± SE of 5 fruits and an asterisk indicates a significant difference (P ≤ 0.05) between ‘Ora’ and ‘Shani’ fruits at the same stage of 6 ripening. ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 32
7
Figure 3. Total carotenoid contents of flavedo of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. Data are means SE of 3 fruits. Different capital letters indicate a significant difference among the ‘Ora’ fruits at the different ripening stages. Different lower-case letters indicate a significant difference among the ‘Shani’ fruits at the different ripening stages. An asterisk indicates a significant difference (P 0.05) between ‘Ora’ and ‘Shani’ fruits at a particular ripening stage.
7 ACS Paragon Plus Environment
Page 27 of 32
Journal of Agricultural and Food Chemistry
8
Figure 4. Carotenoid composition in flavedo of ‘Ora’ and ‘Shani’ fruits at the orange (ripe) stage of ripening. Numbers indicate the percent if each carotenoid of the total carotenoid content for each variety. Data points are means of 3 fruits. β-xanthopylls – neoxanthin, antheraxanthin and zeaxanthin. Carotenes – α, β and ζ carotene, phytofluene and phytoene.
8 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 28 of 32
9
Figure 5. Venn diagram illustration of differentially expressed genes (DEGs) in flavedo of ‘Ora’ and ‘Shani’ during the green, breaker, yellow and orange stages of ripening. The diagram displays the number of DEGs specific to ‘Ora’, the number of DEGs common to ‘Ora’ and ‘Shani’ and the number of DEGs specific to ‘Shani’.
9 ACS Paragon Plus Environment
Page 29 of 32
Journal of Agricultural and Food Chemistry
10
Figure 6. Carotenoid-biosynthesis genes differentially expressed in flavedo of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. A) Genes highly expressed in ‘Shani’ B) Genes highly expressed in ‘Ora’. The results are expressed as fragments per feature kilobase per million reads mapped (FPKM). An asterisk indicates a significant differences (FDR 0.05) between ‘Ora’ and ‘Shani’ fruits at a particular ripening stage, according to the DESeq2 analysis.
10 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 30 of 32
11
Figure 7. MEP-pathway genes differentially expressed in flavedo of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. The results are expressed as fragments per feature kilobase per million reads mapped (FPKM). An asterisk indicates a significant difference (FDR 0.05) between ‘Ora’ and ‘Shani’ fruits at the same ripening stage, according to the DESeq2 analysis.
11 ACS Paragon Plus Environment
Page 31 of 32
Journal of Agricultural and Food Chemistry
12
Figure 8. Ripening regulation in flavedo of ‘Ora’ and ‘Shani’ mandarins. A) Ethylene-biosynthesis and ethylene-signaling genes differentially expressed in flavedo of ‘Ora’ and ‘Shani’ fruits at the green, breaker, yellow and orange stages of ripening. Results are expressed as fragments per feature kilobase per million reads mapped (FPKM). An asterisk indicates a significant difference (FDR 0.05) between ‘Ora’ and ‘Shani’ fruits at the same ripening stage, according to the DESeq2 analysis. B) Ethylene production of ‘Ora’ and ‘Shani’ at the orange ripening stage. Data are means ± SE of 3 replications, each containing 5 fruit. An asterisk indicates a significant difference (P 0.05) between ‘Ora’ and ‘Shani’ fruits.
12 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 32 of 32
13
TOC Graphics
13 ACS Paragon Plus Environment