Biochemical and Molecular Factors Governing Peel-Color

Apr 11, 2019 - To identify factors governing peel-color development in mandarins, we examined carotenoid content and composition and the expression of...
1 downloads 0 Views 943KB Size
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