Physiological and Transcriptomic Analysis Validates Previous

Jul 31, 2017 - Primers for qRT-PCR, summary of RNA-seq data, and Venn diagrams of differentially expressed genes(PDF). View: ACS ActiveView PDF | PDF ...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIVERSITY OF THE SUNSHINE COAST

Article

Findings of Changes in Primary Metabolism for the Production of Phenolic Antioxidants in Wounded Carrots Cong Han, Peng Jin, Meilin Li, Lei Wang, and Yonghua Zheng J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b01137 • Publication Date (Web): 31 Jul 2017 Downloaded from http://pubs.acs.org on August 1, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36

Journal of Agricultural and Food Chemistry

1

Physiological and Transcriptomic Analysis Validates Previous

2

Findings of Changes in Primary Metabolism for the

3

Production of Phenolic Antioxidants in Wounded Carrots

4

5

† † † ‡ *,† Cong Han, Peng Jin, Meilin Li, Lei Wang, Yonghua Zheng

6 7 †

8 9

10 11

College of Food Science and Technology, Nanjing Agricultural University,

Nanjing, 210095, People’s Republic of China ‡

College of Agriculture, Liaocheng University, Liaocheng, 252000, People’s

Republic of China

12 13 14 15

16 17

*

Corresponding author. Tel.:+86 25 8439 9080; Fax: +86 25 8439 5618

E-mail address: [email protected] (Y. H. Zheng)

18 19 20 21 22

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

23

ABSTRACT: Wounding induces the accumulation of phenolic compounds in

24

carrot. This study uses physiological and transcriptomic analysis to validate

25

previous findings relating primary metabolism and secondary metabolites in

26

wounded carrots. Our data confirmed that increased wounding intensity

27

strengthened the accumulation of phenolics accompanied by enhancing respiration

28

and showed the loss of fructose and glucose and the increase of energy status in

29

carrots. Besides, transcriptomic evaluation of shredded carrots indicated that the

30

respiratory metabolism, sugar metabolism, energy metabolism and phenolics

31

biosynthesis related pathways, such as “Citrate cycle (TCA cycle)”, “Oxidative

32

phosphorylation” and “Phenylpropanoid biosynthesis”, were activated by wounding.

33

Also, the differentially expressed genes (DEGs) involved in the conversion of

34

sugars to phenolics were extensively up-regulated after wounding. Thus, the

35

physiological and transcriptomic data validate previous findings that wounding

36

accelerates the primary metabolisms of carrot including respiratory metabolism,

37

sugar metabolism and energy metabolism to meet the demand for the production of

38

phenolic antioxidants.

39 40

KEYWORDS: carrot, wounding, phenolic antioxidants, RNA-Seq, primary

41

metabolism

42 43 44

2 ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

Journal of Agricultural and Food Chemistry

45

■ INTRODUCTION

46

In the last few years, postharvest abiotic stresses have been used as an effective

47

strategy for enhancing the health benefits and medical properties of fresh fruit and

48

vegetables. Many research reports have shown that wounding can obviously induce

49

the production of phenolic compounds in a large variety of horticultural crops, such

50

as carrot,1 potato,2 onion3 and pitaya.4 Phenolic compounds, as the most important,

51

numerous and ubiquitous groups of compounds in plant kingdom, have been

52

considered the major contributors to the total antioxidant capacity of fruit,

53

vegetables and grains.5 Wound-induced accumulation of phenolic compounds

54

usually occurs concomitantly with enhanced antioxidant activity. In this sense, the

55

use of wounding provides us a simple and feasible way to gain more phenolic

56

antioxidants in horticultural crops.

57

Among these stressed crops, carrot has been used as a model system to study the

58

synergetic effects of wounding in combination with other abiotic stresses and the

59

potential mechanisms of wound-induced accumulation of phenolic compounds.6–12

60

To date, some of the physiological and signaling mechanisms on biosynthesis of

61

phenolic compounds have been elucidated by carrot tissue. As interpreted by Reyes

62

et al (2007),8 the increase in total soluble phenolic (TSP) content can be attributed

63

to the faster phenolic synthesis rate, as compared with its decrease or utilization rate.

64

Moreover, reactive oxygen species (ROS) have been demonstrated playing a crucial

65

role as signaling molecules on the wound-induced accumulation of phenolic

66

compounds in carrot.10,12 In addition, most recent studies have shown that the 3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

67

primary metabolism supports the secondary metabolism in horticultural crops to

68

produce secondary metabolites when subjected to wounding stress. It was reported

69

that wounding stress could result in the accumulation of shikimic acid and

70

L-phenylalanine in carrot, both of which are required as carbon sources for the

71

production of phenolic compounds.11,13 In plants, the aromatic amino acids are

72

drived from sugars.14 The increase of chlorogenic acid isomers contents in wounded

73

potato tubers was associated with a pronounced increase in reducing sugar levels,

74

indicating a tight correlation between primary carbohydrate metabolism and

75

secondary phenylpropanoid metabolism.15

76

In general, the burst in respiration has been recognized as a common

77

phenomenon in wounded tissues.16 Respiration of sugars via glycolysis, the

78

oxidative pentose phosphate pathway (OPPP), and the tricarboxylic acid (TCA)

79

cycle can provide precursors for the synthesis of organic acids, amino acids, and

80

many other secondary metabolites, including phenolic compounds.17 Accompanied

81

by mitochondrial respiration, energy is released and chiefly kept in terms of ATP.18

82

In light of these information we wanted to validate previous results that wounded

83

products may accelerate their primary metabolisms, such as the respiratory

84

metabolism, the sugar metabolism and the energy metabolism to provide sufficient

85

precursors and energy for the production of phenolic compounds, thereby repairing

86

and healing the damaged tissues.

87

To validate these previous findings, we investigated the TSP content, the

88

respiration rate, the soluble sugar concentrations and the energy status of carrot 4 ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

Journal of Agricultural and Food Chemistry

89

tissues under different wounding intensities. Particularly, transcriptomic analysis of

90

shredded carrots at different storage period was performed and differentially

91

expressed genes (DEGs) involved in the conversion of sugars to phenolics were

92

screened. The information generated by RNA-Seq would be useful since it validates

93

previous findings on the molecular mechanisms of wound-induced accumulation of

94

phenolic compounds in carrot.

95 96

■ MATERIALS AND METHODS

97

Plant Material, Processing and Storage. Carrots (Daucus carota L. cv.

98

Sanhongliucun) were obtained from a local wholesale market (Nanjing, China) and

99

transported to the laboratory within 1 h. Uniform size, color and unblemished fresh

100

tissues were selected, washed, air-dried and conditioned overnight at 20 °C before

101

processing.

102

In the first experiment, whole carrots were cut into three different cutting styles-

103

slices, pies and shreds, respectively. Based on the method described by Surjadinata

104

and Cisneros-Zevallos,1 the calculated wounding intensities (surface area/weight)

105

were 4.5, 6.2 and 19.1 cm2 g-1, respectively, for slices, pies, and shreds. After

106

wounding treatment, all wounded samples were placed in polypropylene containers

107

and kept at 20 °C for 2 days. Non-wounded carrots were used as control. The TSP

108

content and the antioxidant capacity were determined every 24 h, while the

109

respiration rate, the soluble sugar contents and the energy status were measured at 0

110

(before wounding), 1, 3, 6, 12, 24, 36 and 48 h. At each sampling time, fresh 5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

111

samples were taken, frozen in liquid nitrogen and kept in polyethylene bags at

112

–80 °C.

113

In the second experiment, shredded carrots with the same storage condition as

114

described in the first experiment were used for the transcript profiling analysis.

115

Fresh tissues were sampled individually at the following four stages of storage: (1)

116

Before wounding (Control); (2) 6 h after wounding (W6); (3) 24 h after wounding

117

(W24); (4) 48 h after wounding (W48). The collected samples from each individual

118

were pooled together, frozen immediately in liquid nitrogen and kept in sterile tubes

119

at –80 °C for RNA separation.

120

Total Soluble Phenolics (TSP) Quantification. Quantification of TSP was

121

conducted as earlier described by Slinkard and Slingleton.19 Frozen tissues (5 g)

122

were extracted with cold methanol (25 mL). The tissue homogenate was preserved

123

in covered plastic tubes at 4 °C for 12 h, and centrifuged at 13,000g for 20 min. The

124

reaction mixture consisted of 500 µL clear supernatant, 1.5 mL of distilled water

125

and 1 mL of Folin−Ciocalteu reagent. The reaction was started by adding 1 mL of

126

7.5% Na2CO3 solution (1 mL). After incubating for 2 h at 25 °C, the absorbance at

127

765 nm was measured and the TSP content was expressed as mg gallic acid (GAE)

128

kg-1 (fresh weight).

129

Antioxidant Capacity Analysis. Antioxidant capacity was determined using

130

2,2-diphenyl-1-picrylhydrazyl (DPPH) free radical-scavenging activity. An aliquot

131

of 0.2 ml of the same extracts for the determination of TSP content was added to 2.8

132

mL of 0.12 mM DPPH solution (dissolved in ethanol) and kept in the dark at 25 °C. 6 ACS Paragon Plus Environment

Page 6 of 36

Page 7 of 36

Journal of Agricultural and Food Chemistry

133

After 30 min, the absorbance was measured at 517 nm (A1) and methanol was used

134

as control (A0). The results of DPPH inhibition were calculated with the following

135

equation:

136

% DPPH inhibition=[(A0−A1)/A0]×100

137

Measurement of Respiration Rate. The respiration rate was determined based on

138

Gao et al with slight modifications.20 Briefly, about 1 Kg of slices, pies, shreds or

139

whole carrot tissues were enclosed independently in a 4 L glass chamber with a

140

small hole on the lid. Then, the hole was covered by adhesive tape and the container

141

was kept at 20 °C for 1 h. Finally, the respiration rate of carrot tissues was assayed

142

as CO2 production by using the gas analyzer (CheckMate II, PBI Dansensor,

143

Denmark). The respiration rate was expressed as mg CO2 kg−1 h−1.

144

Soluble Sugars Analysis. The extraction of soluble sugars from carrot was

145

modified from the experimental data described by Macedo and Peres.21 Frozen

146

tissues (2 g) were extracted using of 80% ethanol solution (5 mL). A high

147

performance liquid chromatography (HPLC) system (Shimadzu LC-20A, Japan)

148

equipped with a quaternary pump system (Shimadzu DGU-20A, Japan), an

149

auto-injector (Shimadzu SIL-20A, Japan), an evaporative light scattering detector

150

(Shimadzu ELSD-LT II, Japan) and a column heater (Shimadzu CTO-20A, Japan)

151

was used for quantification. The column used to separate the soluble sugars was a

152

4.6 × 250 mm, 5 µm, hydrophilic interaction chromatography (HILIC) column

153

(Shodex Asahipak NH2P-50 4E, Japan), which was maintained at 40 °C. Elution

154

was performed using the mobile phase of acetonitrile (A) and water (B). The soluble 7 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 36

155

sugar content was identified according to the pure standards and quantified using

156

standard calibration curves.

157

ATP, ADP, and AMP Contents and Energy Charge Assays. Frozen tissues (2 g)

158

were finely ground in perchloric acid (5 mL). ATP, ADP, and AMP were quantified

159

using a HPLC system (Agilent 1100) based on a method previously reported by Jin

160

et al.22 Energy charge was calculated with the following equation:

161

[ATP + 1/2ADP]/ [ATP + ADP + AMP].

162

RNA-sequencing (RNA-seq). To obtain the transcriptome of the aforementioned

163

four groups, RNAs were isolated from frozen tissue samples using Plant Total RNA

164

Extraction Kit (TaKaRa 9769, Shanghai, China). For each group, we performed two

165

replications. Sequencing libraries were generated following manufacturer’s

166

recommendations (NEB, USA), which were named Control, W6, W24 and W48,

167

respectively. After cluster generation, the library preparations were sequenced using

168

Illumina HiseqTM 4000. After preprocessing, the data were in a cleaned and

169

polished form. Clean data with high quality were used for the subsequent analyses.

170

Clean

reads

were

mapped

to

the

reference

genome

of

carrot

171

(https://www.ncbi.nlm.nih.gov/genome/?term=Daucus+carota) using Tophat2 (v

172

2.1.0)23. The reads numbers mapped to each gene were counted using HTSeq (v

173

0.6.0).24 FPKM (reads per kilobase per million reads) were used to quantify the

174

mapped whole gene expression levels. Genes with an adjusted P-value 99.9

274

(Q30), and the summary of RNA-seq data is shown in Supplementary Table S2.

275

Read count data were used to identify the DEGs among samples by using DESeq R

276

packages.26 As shown in the Venn diagrams (Supplementary Figure S1), altogether

277

13407, 14448, and 14509 DEGs were identified between W6 and Control, W24 and

278

Control, and W48 and Control, respectively. There were 4865 up-regulated and

279

4837 down-regulated DEGs that were commonly shared among the three pairwise

280

comparisons.

281

KEGG Analysis of DEGs. To further look into the biological functions of genes in

282

carrot, the DEGs were searched against the reference canonical pathways in KEGG.

283

This study focused mainly on these KEGG pathways related to respiratory

284

metabolism,

285

metabolism. Here, a total of 11 pathways were listed in Table 1. It was found that

286

wounding resulted in extensive up-regulation of DEGs assigned to the pathways of

sugar

metabolism,

energy

metabolism

13 ACS Paragon Plus Environment

and

phenylpropanoid

Journal of Agricultural and Food Chemistry

287

“Pyruvate metabolism”, “Pentose phosphate pathway”, “Citrate cycle (TCA cycle)”

288

and “Oxidative phosphorylation”, especially in W24 and W48. As these pathways

289

are known closely associated with respiratory metabolism,36 this observation may

290

account for the remarkable increase in respiration rate observed through time after

291

wounding. In general, since metabolism is an integrated network, one pathway may

292

have a close link between different metabolic processes. Glycolysis, the oxidization

293

of glucose to pyruvate, is relevant to both the respiratory metabolism and the sugar

294

metabolism. As shown in the pathway of “Glycolysis / Gluconeogenesis”, the

295

number of up-regulated DEGs increased while down-regulated DEGs declined

296

through time. Similar variation trend was observed in the pathways of “Fructose

297

and mannose metabolism” and “Starch and sucrose metabolism”, indicating an

298

enhancement of the degradation of carbohydrates in wounded carrots.

299

The metabolic pathways, such as glycolysis, TCA cycle and oxidative

300

phosphorylation, are essential for the production of ATP. Pyruvate produced from

301

glycolysis is oxidised by the mitochondrial TCA cycle, and electrons from the

302

resulting reductant are transferred through the electron transport chain and

303

ultimately generating ATP by the oxidative phosphorylation.37 Accordingly, the

304

large proportion of up-regulated DEGs, especially those in “Citrate cycle (TCA

305

cycle)” and “Oxidative phosphorylation”, may further confirm the elevated

306

synthesis of ATP and the increased demand for energy consumption in wounded

307

carrots. In addition, approximately 85% of all DGEs belong to the pathway of

308

“Phenylpropanoid biosynthesis”, were up-regulated within a short time after 14 ACS Paragon Plus Environment

Page 14 of 36

Page 15 of 36

Journal of Agricultural and Food Chemistry

309

wounding. And the pathways of “Biosynthesis of amino acids”, “Phenylalanine,

310

tyrosine and tryptophan” and “Phenylalanine metabolism”, which involved in the

311

biosynthesis and degradation of some specific amino acids, were also markedly

312

stimulated by wounding.

313

To visualize the metabolic process in relation to the conversion of sugars to

314

phenolics, different KEGG pathways were linked together and depicted in Figure 4.

315

The FPKM values of DEGs encoding specific enzymes are shown in Table 2. Our

316

data revealed that the expression of genes encoding those sugar catabolic enzymes

317

such as invertase, fructokinase and hexokinase seemed not fully stimulated at 6 h,

318

but exhibited obvious upward trends at 24 h and 48 h. This observation matched

319

with the changes of soluble sugars contents in shredded carrots, indicating the

320

growing demands of carbohydrates metabolic intermediates for the biosynthesis of

321

amino acids. However, as for the genes in regulating the shikimic acid pathway or

322

subsequent conversion of chorismate to phenylalanine, their expressions were

323

extensively and dramatically up-regulated at 6 h and kept high levels during

324

prolonged storage. This observation conformed previous findings on the variations

325

of shikimic acid and L-phenylalanine contents in wounded carrots,11,13 validating

326

the coordinated switches in primary metabolism for the production of phenolic

327

compounds. Furthermore, the DEGs related with conversion of phenylalanine to

328

phenolic profiles was significantly up-regulated after wounding. This result

329

confirmed the previous work by Becerra-Moreno et al and Jacobo-Velázquez et

330

al,11,12 which validates that wounding has a significant role in activating the 15 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 36

331

expression of genes involved in phenolics biosynthesis, and TSP content was

332

increased consequently.

333

Validation of Selected DEGs by qRT-PCR Analysis. Phenylalanine ammonium

334

lyase (PAL), cinnamate-4-hydroxylase (C4H) and 4-coumarate coenzyme A ligase

335

(4CL) are three pivotal upstream enzymes in the phenylpropanoid pathway, which

336

dominate phenolics production in plants.38 In this study, 9 DEGs that encode for

337

these three enzymes, including three for PAL (LOC108223317, LOC108227736,

338

LOC108217754), two for C4H (LOC108223289, LOC108197548), and four for

339

4CL (LOC108227923, LOC108193321, LOC108206420, LOC108223931), were

340

selected to validate the accuracy and reproducibility of the RNA-Seq data by

341

qRT-PCR (Table 3). Generally, except for the DEGs of LOC108227736,

342

LOC108197548 and LOC108223931, the other 6 assayed genes were significantly

343

up-regulated after wounding and kept high expression levels during storage.

344

Correlation between the two methods was measured by scatter plotting log 2 fold

345

changes (Figure 5). The expression profiles determined by qRT-PCR are generally

346

consistent with the data derived from RNA-Seq, reflecting the reliability of results

347

from RNA-Seq.

348 349

■AUTHOR INFORMATION

350

Corresponding Author

351

*Telephone:

352

[email protected]

+86

25

8439

9080.

Fax:

+86

25

16 ACS Paragon Plus Environment

8439

5618.

E-mail:

Page 17 of 36

Journal of Agricultural and Food Chemistry

353

Funding

354

This study was financially supported by the National Natural Science Foundation of

355

China (No. 31471632).

356

Notes

357

The authors declare no competing financial interest.

358 359

■References

360

(1) Surjadinata, B. B.; Cisneros-Zevallos, L. Biosynthesis of phenolic

361

antioxidants in carrot tissue increases with wounding intensity. Food Chem. 2012,

362

134, 615–624.

363

(2) Torres-Contreras, A. M.; Nair, V.; Cisneros-Zevallos, L.; Jacobo-Velázquez, D.

364

A. Plants as biofactories: stress-induced production of chlorogenic acid isomers in

365

potato tubers as affected by wounding intensity and storage time. Ind. Crops Prod.

366

2014, 62, 61–66.

367

(3) Berno, N. D.; Tezotto-Uliana, J. V.; dos Santos Dias, C. T.; Kluge, R. A.

368

Storage temperature and type of cut affect the biochemical and physiological

369

characteristics of fresh-cut purple onions. Postharvest Biol. Technol. 2014, 93,

370

91–96.

371

(4) Li, X. A.; Long, Q. H.; Gao, F.; Han, C.; Jin, P.; Zheng, Y. H. Effect of cutting

372

styles on quality and antioxidant activity in fresh-cut pitaya fruit. Postharvest Biol.

373

Technol. 2017, 124, 1–7.

374

(5) Haminiuk, C. W. I.; Maciel, G. M.; Plata-Oviedo, M. S. V.; Peralta, R. M. 17 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

375

Phenolic compounds in fruits–an overview. Int. J. Food Sci. 2012, 47, 2023–2044.

376

(6) Lafuente, M. T.; López-Gálvez, G.; Cantwell, M.; Yang, S. F. Factors

377

influencing ethylene-induced isocoumarin formation and increased respiration in

378

carrots. J. Amer. Soc. Hort. Sci. 1996, 121, 537–542.

379

(7) Cisneros-Zevallos, L. The use of controlled postharvest abiotic stresses as a

380

tool for enhancing the nutraceutical content and adding-value of fresh fruits and

381

fegetables. J. Food Sci. 2003, 68, 1560–1565.

382

(8) Reyes, L. F.; Villarreal, J. E.; Cisneros-Zevallos, L. The increase in

383

antioxidant capacity after wounding depends on the type of fruit or vegetable tissue.

384

Food Chem. 2007, 101, 1254–1262.

385

(9) Heredia, J. B.; Cisneros-Zevallos, L. The effect of exogenous ethylene and

386

methyl jasmonate on pal activity, phenolic profiles and antioxidant capacity of

387

carrots (Daucus carota) under different wounding intensities. Postharvest Biol.

388

Technol. 2009, 51, 242–249.

389

(10) Jacobo-Velázquez, D. A.; Martinez-Hernandez, G. B.; Rodriguez, S.; Cao, C.

390

M.; Cisneros-Zevallos, L. Plants as biofactories: physiological role of reactive

391

oxygen species on the accumulation of phenolic antioxidants in carrot tissue under

392

wounding and hyperoxia stress. J. Agric. Food Chem. 2011, 59, 6583–6593.

393

(11) Becerra-Moreno, A.; Redondo-Gil, M.; Benavides, J.; Nair, V.;

394

Cisneros-Zevallos, L.; Jacobo-Velázquez, D. A. Combined effect of water loss and

395

wounding stress on gene activation of metabolic pathways associated with phenolic

396

biosynthesis in carrot. Front. Plant Sci. 2015, 6, 837. 18 ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36

Journal of Agricultural and Food Chemistry

397

(12) Jacobo-Velázquez, D. A.; González-Agüero, M.; Cisneros-Zevallos, L.

398

Cross-talk between signaling pathways: the link between plant secondary

399

metabolite production and wounding stress response. Sci. Rep. 2015, 5, 8608.

400

(13) Becerra-Moreno, A.; Benavides, J.; Cisneros-Zevallos, L.; Jacobo-Velázquez,

401

D. A. Plants as biofactories: glyphosate-induced production of shikimic acid and

402

phenolic antioxidants in wounded carrot tissue. J. Agric. Food Chem. 2012, 60,

403

11378−11386.

404 405

(14) Yoshida, S. Biosynthesis and conversion of aromatic amino acids in plants. Annu. Rev. Plant Physiol. 1969, 20, 41−62.

406

(15) Torres-Contreras, A. M.; Nair, V.; Cisneros-Zevallos, L.; Jacobo-Velázquez,

407

D. A. Effect of exogenous amylolytic enzymes on the accumulation of chlorogenic

408

acid isomers in wounded potato tubers. Annu. Rev. Plant Physiol. 2014, 62,

409

7671−7675.

410

(16) Hodges, D. M., Toivonen, P. M. Quality of fresh-cut fruits and vegetables as

411

affected by exposure to abiotic stress. Postharvest Biol. Technol. 2008, 48, 155–162.

412

(17) Dai, Z. W.; Léon, C.; Feil, R.; Lunn, J. E.; Delrot, S.; Gomès, E. Metabolic

413

profiling reveals coordinated switches in primary carbohydrate metabolism in grape

414

berry (vitis vinifera L.), a non-climacteric fleshy fruit. J. Exp. Bot. 2013, 64, 1345.

415 416 417 418

(18) Melser, S.; Lavie, J.; Bénard, G. Mitochondrial degradation and energy metabolism. Biochim. Biophys. Acta 2015, 1853, 2812–2821. (19) Slinkard, K. S.; Singleton, V. L. Total phenol analysis: automation and comparison with manual methods. Am. J. Enol. Vitic. 1977, 28, 49–55. 19 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 20 of 36

419

(20) Gao, P. S.; Zhu, Z. Q.; Zhang, P. Effects of chitosan–glucose complex

420

coating on postharvest quality and shelf life of table grapes. Carbohydr. Polym.

421

2013, 95, 371–378.

422

(21) Macedo, E. A.; Peres, A. M. Thermodynamics of ternary mixtures

423

containing sugars. SLE of D-fructose in pure and mixed solvents. Comparison

424

between modified UNIQUAC and modified UNIFAC. Ind. Eng. Chem. Res., 2001,

425

40, 4633–4640.

426

(22) Jin, P.; Zhu, H.; Wang, J.; Chen, J. J.; Wang, X. L.; Zheng, Y. H. Effect of

427

methyl jasmonate on energy metabolism in peach fruit during chilling stress. J. Sci.

428

Food Agric. 2013, 93, 1827−1832.

429 430 431

(23) Trapnell, C.; Pachter, L.; Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25, 1105−1111. (24) Anders, S. HTSeq: Analysing high-throughput sequencing data with

432

Python. EMBL

433

http://www-huber. embl. de/users/anders/HTSeq/doc/overview. html.

434 435

Heidelberg

(Genome

Biology

Unit).

2010.

Available:

(25) Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106.

436

(26) Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.;

437

Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; Yamanishi, Y. KEGG for

438

linking genomes to life and the environment. Nucleic Acids Res. 2008, 36,

439

D480−D484.

440

(27) Han, C.; Li, J.; Jin, P.; Li, X. A.; Wang, L.; Zheng, Y. H. The effect of 20 ACS Paragon Plus Environment

Page 21 of 36

Journal of Agricultural and Food Chemistry

441

temperature on phenolic content in wounded carrots. Food Chem. 2017, 215,

442

116–123.

443 444

(28) Viña, S. Z.; Chaves, A. R. Respiratory activity and phenolic compounds in pre-cut celery. Food Chem. 2007, 100, 1654–1660.

445

(29) Surjadinata, B. B.; Cisneros-Zevallos, L. Modeling wound-induced

446

respiration of fresh-cut carrots (Daucus carota L.). J. Food Sci. 2003, 68,

447

2735–2740.

448

(30) Jacobo-Velázquez, D. A.; Cisneros-Zevallos, L. An alternative use of

449

horticultural crops: stressed plants as biofactories of bioactive phenolic

450

compounds. Agriculture 2012, 2, 259–271.

451 452 453 454

(31)

Dixon,

R.

A.;

Paiva,

N.

L.

Stress-induced

phenylpropanoid

metabolism. Plant Cell, 1995, 7, 1085. (32) Suojala, T. Variation in sugar content and composition of carrot storage roots at harvest and during storage. Sci. Hortic. 2000, 85, 1–19.

455

(33) Aguiló-Aguayo, I.; Hossain, M. B.; Brunton, N.; Lyng, J.; Valverde, J.; Rai,

456

D. K. Pulsed electric fields pre-treatment of carrot purees to enhance their

457

polyacetylene and sugar contents. Innov. Food Sci. Emerg. Technol. 2014, 23,

458

79–86.

459

(34) Klotz, K. L.; Finger, F. L.; Anderson, M. D. Wounding increases glycolytic

460

but not soluble sucrolytic activities in stored sugarbeet root. Postharvest Biol.

461

Technol. 2006, 41, 48–55.

462

(35) Yi, C.; Qu, H. X.; Jiang, Y. M.; Shi, J.; Duan, X. W.; Joyce, D. C.; Li, Y. B. 21 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

463

ATP-induced changes in energy status and membrane integrity of harvested litchi

464

fruit and its relation to pathogen resistance. J. Phytopathol. 2008, 156, 365–371.

465 466

(36) Wang, J. Y.; Zhu, S. G.; Xu, C. F. Biochemistry, 3nd ed. Higher Education Press: Beijing, China, 2002, pp. 63–175.

467

(37) Alisdair, R. F.; Fernando, C.; Sweetlove, J. Respiratory metabolism:

468

glycolysis, the TCA cycle and mitochondrial electron transport. Plant Biol. 2004, 7,

469

254–261.

470 471

(38) Weisshaar, B.; Jenkins, G. I. Phenylpropanoid biosynthesis and its regulation. Curr. Opin. Plant Biol. 1998, 1, 251–257.

22 ACS Paragon Plus Environment

Page 22 of 36

Page 23 of 36

Journal of Agricultural and Food Chemistry

472

■Legends to Figures

473

Figure 1. Effect of wounding intensity on TSP content (A) and DPPH inhibition (B)

474

of wounded carrots during 2 days storage at 20 °C. Each data point is expressed as

475

the mean ± SD (n=3). Different letters at each storage time represent significant

476

differences at p < 0.05.

477

Figure 2. Effect of wounding intensity on respiration rate (A), fructose content (B),

478

glucose content (C) and sucrose content (D) of wounded carrots during 2 days

479

storage at 20 °C. Each data point is expressed as the mean ± SD (n=3).

480

Figure 3. Effect of wounding intensity on ATP content (A), ADP content (B), AMP

481

content (C) and energy charge (D) of wounded carrots during 2 days storage at

482

20 °C. Each data point is expressed as the mean ± SD (n=3).

483

Figure 4. Metabolic pathways involved in conversion of sugars to phenolics. FK,

484

fructokinase;

485

glucose-6-phosphate 1-epimerase; G6PI, glucose-6-phosphate isomerase; TK,

486

transketolase; EL, enolase; DAHP, 3-deoxy-D-arabino-heptolosonate-7-phosphate;

487

3DQ, 3-dehydroquinate; 3DD, SD, 3-dehydroquinate dehydratase/shikimate

488

dehydrogenase;

489

synthase; CM, chorismate mutase; ADT, PDT, arogenate dehydratase/prephenate 23

HK,

EPSP,

hexokinase;

A1E,

aldose

1-epimerase;

5-enolpyruvylshikimate-3-phosphate;

ACS Paragon Plus Environment

CS,

G6P1E,

chorismate

Journal of Agricultural and Food Chemistry

AAT,

PAT,

bifunctional

490

dehydratase;

491

glutamate/aspartate-prephenate amino transferase; AA, aspartate aminotransferase;

492

TA, tyrosine aminotransferase; HPT, histidinol-phosphate aminotransferase; PAL,

493

phenylalanine ammonium lyase; C4H, cinnamate-4-hydroxylase; 4CL, 4-coumarate

494

coenzyme A ligase; HCT, shikimate O-hydroxycinnamoyltransferase; C3’M,

495

coumaroylquinate

496

coumaroylshikimic acid; CSE, caffeoylshikimate esterase; COMT, caffeate

497

O-methyl transferase. Dotted line represents that some steps were omitted.

498

Figure 5. Correlation of fold change analyzed by RNA-Seq (x axis) and the data

499

obtained using qRT-PCR (y axis). The values were derived from Table 3.

(coumaroylshikimate)

aspartate

aminotransferase

Page 24 of 36

3'-monooxygenase;

500

24 ACS Paragon Plus Environment

and

CSA,

Page 25 of 36

Journal of Agricultural and Food Chemistry

Table 1 KEGG analysis of DEGs related to respiratory metabolism, sugar metabolism, energy metabolism and phenylpropanoid metabolisma W6 vs Control W24 vs Control W48 vs Control KEGG pathways Background genes Up

Down

Up

Down

Up

Down

Pyruvate metabolism

43

37

48

28

54

29

90

Biosynthesis of amino acids

103

67

129

58

117

68

232

Pentose phosphate pathway

18

22

21

16

26

13

46

Citrate cycle (TCA cycle)

35

10

45

5

39

3

54

Glycolysis / Gluconeogenesis Fructose and mannose metabolism

40

53

57

43

62

39

126

10

27

19

23

18

17

50

Phenylalanine, tyrosine and tryptophan biosynthesis Oxidative phosphorylation

27

9

28

7

28

9

54

41

44

67

24

73

18

161

Starch and sucrose metabolism

60

43

76

34

86

34

194

Phenylalanine metabolism

41

16

52

16

54

16

147

Phenylpropanoid biosynthesis

66

13

82

13

85

15

207

a

The number of up- and down-regulated DEGs is determined by comparing with control. Background genes indicate the number of all genes in each pathway.

25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 26 of 36

Table 2. Summary of DEGs involved in regulating the conversion of sugars to phenolicsa Description Invertase

HK; hexokinase

FK; fructokinase

A1E; aldose 1-epimerase G6P1E; glucose-6-phosphate 1-epimerase

Gene symbol LOC108218454 LOC108218455 LOC108209529 LOC108222817 LOC108216856 LOC108222818 LOC108206477 LOC108218494 LOC108199548 LOC108209997 LOC108203143 LOC108225510 LOC108199984 LOC108228173 LOC108195110 LOC108197533 LOC108209124 LOC108213261 LOC108220493 LOC108210430 LOC108192895 LOC108215076 LOC108206071

Control 0.00 0.19 0.42 0.06 0.00 0.01 0.00 0.63 71.50 22.91 5.82 0.10 166.89 310.58 9.50 2.29 0.75 0.06 9.07 0.49 25.94 17.96 2.18

26

ACS Paragon Plus Environment

FPKM value W6 W24 2.43 6.04 23.70 51.88 154.98 272.77 1.03 1.73 0.01 0.01 0.11 0.16 0.46 0.81 5.83 11.78 64.25 58.80 19.14 14.96 8.60 9.86 0.05 0.07 210.12 264.53 186.35 42.39 4.97 1.61 1.39 0.82 6.10 11.57 0.28 0.46 9.09 9.87 1.38 2.27 37.19 53.29 21.67 25.33 1.41 0.31

W48 3.62 70.36 168.60 4.45 1.07 0.76 3.45 24.87 60.06 17.87 9.74 0.46 451.12 52.78 3.50 1.20 28.59 2.38 16.58 4.95 77.27 29.52 0.51

Page 27 of 36

Journal of Agricultural and Food Chemistry

LOC108199656 LOC108196280 G6PI; glucose-6-phosphate isomerase LOC108227561 LOC108214821 TK; transketolase LOC108227965 LOC108209071 EL; enolase LOC108210896 LOC108217564 LOC108192892 LOC108227937 LOC108223971 DAHP synthase; LOC108214600 LOC108223353 3-deoxy-D-arabino-heptolosonate-7-phosphate LOC108218443 synthetase 3DQ synthase; 3-dehydroquinate synthase LOC108228166 3DD, SD; 3-dehydroquinate dehydratase/shikimate LOC108202931 dehydrogenase LOC108204867 LOC108206321 LOC108197089 SK; shikimate kinase LOC108204238 LOC108205161 EPSE synthase; 5-enolpyruvylshikimate-3-phosphate LOC108209808 synthase LOC108209237 CS; chorismate synthase LOC108219299 CM; chorismate mutase LOC108224600 LOC108196652

6.01 6.80 285.61 26.02 23.29 157.20 7.71 6.10 214.94 1664.68 34.20 93.41 2.60 44.70 192.46 0.11 0.25 10.86 4.95 7.19 0.19 5.93 0.35 62.40 4.91 0.98

27

ACS Paragon Plus Environment

2.80 4.89 185.25 27.04 76.77 178.48 14.20 13.79 151.94 1011.98 29.44 73.75 325.88 480.20 185.59 1.73 0.50 8.17 63.38 26.82 0.19 179.15 1.65 158.43 25.95 22.73

0.47 2.11 98.44 25.13 125.44 226.27 20.70 20.15 89.26 319.27 24.69 50.49 760.89 914.29 188.27 2.84 0.88 3.55 125.44 48.57 0.15 366.28 2.85 267.30 44.09 51.56

0.47 2.65 97.84 32.89 182.91 415.21 26.31 31.18 122.22 430.27 25.15 42.46 793.47 846.12 196.92 9.37 3.44 3.00 119.20 63.98 0.23 437.34 1.89 265.24 29.58 44.51

Journal of Agricultural and Food Chemistry

LOC108210290 ADT, PDT; arogenate dehydratase/prephenate LOC108206009 dehydratase LOC108223049 LOC108205816 LOC108205849 LOC108210565 AAT, PAT; bifunctional aspartate aminotransferase and LOC108204769 glutamate/aspartate-prephenate amino transferase AA; aspartate aminotransferase LOC108205030 LOC108204769 LOC108225092 LOC108214974 LOC108204943 TA; tyrosine aminotransferase LOC108224044 LOC108224339 HPT; histidinol-phosphate aminotransferase LOC108197048 PAL; phenylalanine ammonium lyase LOC108227736 LOC108217754 LOC108223317 C4H; cinnamate-4-hydroxylase LOC108223289 LOC108197548 4CL; 4-coumarate coenzyme A ligase LOC108223931 LOC108227923 LOC108193321 LOC108206420 HCT; shikimate O-hydroxycinnamoyltransferase LOC108223132

Page 28 of 36

6.53 18.92 2.42 3.87 0.20 2.66 5.13

18.19 12.93 5.56 9.42 79.94 208.15 71.56

33.97 6.96 9.32 16.04 173.34 466.99 145.45

37.65 6.95 8.93 15.03 141.92 371.24 216.58

0.01 5.13 1.04 327.40 83.93 0.23 1.42 5.57 1110.10 4.54 0.67 44.14 2.98 6.07 3.97 0.05 0.22 28.22

0.03 71.56 3.80 192.75 64.01 3.57 4.04 8.56 884.30 255.77 143.05 1052.64 2.22 6.48 297.41 19.26 99.63 16.60

0.11 145.45 6.82 73.06 43.71 7.26 6.07 12.82 921.58 540.68 275.75 2196.52 1.21 7.50 617.90 40.21 202.80 1.33

0.25 216.58 12.19 67.69 47.46 5.50 4.84 19.05 781.44 531.85 380.58 1729.01 2.67 10.67 535.07 32.00 145.31 0.67

28

ACS Paragon Plus Environment

Page 29 of 36

Journal of Agricultural and Food Chemistry

LOC108220879 0.41 28.22 LOC108196831 19.72 335.17 LOC108223171 1.04 34.33 C3'M; coumaroylquinate(coumaroylshikimate) LOC108227431 0.24 21.85 3'-monooxygenase LOC108209804 1.15 141.97 CSE; caffeoylshikimate esterase LOC108224003 53.65 229.35 LOC108196175 0.05 1.44 LOC108194156 0.09 37.07 LOC108227433 14.84 81.64 COMT; caffeate O-methyl transferase LOC108214501 0.46 19.98 LOC108210405 17.03 88.10 LOC108214345 0.11 0.25 LOC108215346 0.00 0.04 LOC108223084 0.47 0.63 LOC108223086 0.00 0.00 LOC108221040 0.70 6.32 LOC108215577 1.00 5.09 a The expression level of each DEG is represented by FPKM value. Control, before wounding; W6, 6 h after wounding; W48, 48 h after wounding.

29

ACS Paragon Plus Environment

62.70 38.72 663.45 576.58 67.24 146.76 45.05 85.07 282.29 218.27 420.84 296.54 2.52 1.58 77.15 52.79 148.09 136.21 37.77 205.67 177.69 243.40 0.25 0.21 0.04 0.51 0.91 381.91 0.03 28.04 15.84 723.81 9.19 143.40 wounding; W24, 24 h after

Journal of Agricultural and Food Chemistry

Page 30 of 36

Table 3 qRT-PCR verification of 9 DEGs related with the phenylpropanoid biosynthesis in carrota W6 vs Control (log2FC) W24 vs Control (log2FC) Gene description Gene ID

a

W48 vs Control (log2FC)

FPKM

qRT-PCR

FPKM

qRT-PCR

FPKM

qRT-PCR 9.45 -2.25

LOC108223317

phenylalanineammonia-lyase 1

8.43

9.21

9.49

10.33

LOC108227736

phenylalanineammonia-lyase 1









8.31 -2.58

LOC108217754

phenylalanineammonia-lyase 3

6.63

8.78

6.70

8.49

5.66

7.10

LOC108223289

trans-cinnamate 4-monooxygenase

6.83

4.61

6.07

3.27

4.23

LOC108197548

trans-cinnamate 4-monooxygenase-like

5.36 -1.59

1.20









LOC108227923

4-coumarate-CoA ligase 1

7.02

7.94

6.68

7.79

5.60

5.71

LOC108193321

4-coumarate-CoA ligase 1-like

9.24

8.39

8.73

7.85

7.70

6.65

LOC108206420

4-coumarate-CoA ligase-like 6

9.54

14.33

8.41

13.43

7.15

11.94

LOC108223931

4-coumarate-CoA ligase-like 7





1.11

0.47

1.25

-0.16

―, not identified as DEG.

30

ACS Paragon Plus Environment

Page 31 of 36

Journal of Agricultural and Food Chemistry

■FIGURE GRAPHICS

Figure 1

31 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 2

32 ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36

Journal of Agricultural and Food Chemistry

Figure 3

33 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 4

34 ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36

Journal of Agricultural and Food Chemistry

Figure 5

35 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

■TABLE OF CONTENTS GRAPHICS

36 ACS Paragon Plus Environment

Page 36 of 36