Metabolite and Transcriptome analyses indicate the involvement of

Oct 11, 2018 - ABSTRACT Texture is an important component of peach fruit quality. In the present study, an analysis of metabolite and transcriptome pr...
0 downloads 0 Views 4MB Size
Subscriber access provided by UNIV OF LOUISIANA

Agricultural and Environmental Chemistry

Metabolite and Transcriptome analyses indicate the involvement of lignin in programmed changes in peach fruit texture Ying Wang, XinFu Zhang, Shaolan Yang, and Yongbing Yuan J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b04284 • Publication Date (Web): 11 Oct 2018 Downloaded from http://pubs.acs.org on October 15, 2018

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

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

Page 1 of 49

Journal of Agricultural and Food Chemistry

1 2 3

Metabolite and Transcriptome analyses indicate the involvement of lignin in programmed changes in peach fruit texture

4 5

6

Ying Wang • Xinfu Zhang • Shaolan Yang∗• Yongbing Yuan∗

7 8

Qingdao Key Laboratory of Genetic Improvement and Breeding in Horticultural Plants,

9

College of Horticulture, Qingdao Agricultural University, No. 700 Changcheng Road,

10

Chengyang, Qingdao City 266109, Shandong Province, China

11

12

∗Authors

to whom correspondence should be addressed; E-mails: [email protected]

13

(Shaolan Yang); [email protected] (Yongbing Yuan)

14

Tel: 0086-0532-86080018

15

Fax: 0086-0532-86080018

16

17

18

19

20

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

21

ABSTRACT

22

Texture is an important component of peach fruit quality. In the present study, an analysis of

23

metabolite and transcriptome profiles of a non-melting flesh cultivar, ‘Baili’, and a melting flesh

24

cultivar, ‘Hongli’, during storage was conducted to explore the molecular mechanisms underlying

25

different fruit textures in peach. Results indicated that higher levels of anthocyanins were present

26

in ‘Hongli’ peach, while lignin monomers and ethylene precursors were higher in ‘Baili’. A

27

transcriptome analysis indicated that genes associated with lignin synthesis were more highly

28

expressed in ‘Baili’ than in ‘Hongli’, especially Pp4CL2, Pp4CL3, and PpCOMT2. Texture

29

differences between the two varieties may be the result of differential expression of two branches

30

of the phenylpropanoid metabolic pathway. One branch regulates flavonoid metabolism and was

31

highly active in ‘Hongli’ fruit, while the other branch regulates lignin synthesis and was more

32

highly active in ‘Baili’ fruit.

33

Keywords: Peach, Transcriptome, Metabolomics, Lignin

34 35

ACS Paragon Plus Environment

Page 2 of 49

Page 3 of 49

Journal of Agricultural and Food Chemistry

37

INTRODUCTION

38

Texture, an important component used to determine peach fruit quality1, not only affects fruit

39

taste and consumer choice but also affects the fruit storage and transportation; and thus receives a

40

great deal of attention by the peach fruit industry2,3,4. Changes in the texture that occur in peach

41

fruit after harvest can be divided into processes involved in softening and lignification5.

42

Softening is the most noticeable texture change that occurs during ripening and is a component

43

of the maturation of all fruits6. Previous research has indicated that changes in the cell wall and

44

the regulation of plant hormones are the most important aspects controlling fruit softening.

45

Enzymes involved in cell wall degradation are closely associated with softening7.

46

Polygalacturonase (PG) activity becomes greatly elevated during the postharvest ripening of

47

peach fruit and is considered a dominant aspect of peach fruit softening. Qian et al. 8 reported that

48

the difference in the softening characteristics of ‘Qin Wang’ vs ‘Qian Jian Bai’ peach fruit was

49

related to the level of activity of various PG enzymes, the higher the level of PG a fruit contained,

50

the more rapidly it would soften. As with PG, pectinesterase (PE) activity has been detected in a

51

very wide range of fruit. Abu-Sarra et al.

52

activity during the softening of mango, and found that PE activity increased continuously as

53

mango fruit softened, indicating that PE plays a major role in determining the rate of mango fruit

54

softening. Hayama et al.11 reported on the expression pattern of three expansin genes in

55

‘Akatsuki’ (melting) and ‘Manami’ (non-melting) peach during softening and found that PpExp3,

56

rather than PpExp1 or PpExp2, played a role in fruit softening of ‘Akatsuki’ peach. The

57

regulation of plant hormones is also considered to be an important factor of the ripening process12.

58

In particular, ethylene is considered to be the most important hormone in the regulation of fruit

59

softening. In general, fruit softens faster with the application of increasing concentrations of

60

ethylene13, and the expression of Pp-ACS1 and Pp-ACO1 genes was demonstrated to play a key

61

role in the regulation of ethylene biosynthesis during the promotion of ripening in ‘Hakuho’

9,10

studied the changes in cell wall pectinesterase (PE)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 49

62

peach14. And ABA also could accelerate the fruit ripening and softening, as Luo et al15 reported

63

that endogenous ABA initiated cherry fruit ripening by regulating maturation-related metabolic

64

pathways such as anthocyanin synthesis. Moreover, Soto et al.16 indicated that ABA regulated

65

ripening and softening progress of peach fruit by regulating ethylene biosynthesis and signaling

66

genes (ACS1, ACO1, ETR2, ERF2)15.

67

To date, most research studies have focused on peach fruit softening and less on the

68

lignification of peach fruit, even though lignification is also an important aspect affecting the

69

texture of fruit. A report on loquat fruit indicated that lignification was induced by chilling injury

70

and during senescence17. Cai et al analyzed loquat fruits stored at room temperature and found

71

that loquat fruit firmness increased during storage and was positively correlated with the

72

accumulation of lignin in the fruit flesh; in a process that was termed senescence lignification18.

73

In a subsequent study, the accumulation of lignin and an increase in firmness in the flesh of

74

loquat fruit occurred whether fruit was stored at room temperature (senescence lignification) or

75

subjected to chilling injury by storage at 0◦C19. Phenylalanine ammonialyase (PAL), caffeic acid

76

5-hydroxyferulic acid O-methyltransferase(COMT),4-Coumarate: coenzyme A ligase (4CL),

77

cinnamyl alcohol dehydrogenase (CAD), and peroxidase (POD) are all enzymes that are involved

78

in lignin synthesis, and increased during storage of loquat fruit, and were positively correlated

79

with the accumulation of lignin20. This phenomenon has also been reported to occur in pear21 and

80

mangosteen22. PpCAD1 and PpCAD2 all exhibited higher levels of expression in a hard pear fruit

81

cultivar, ‘Whangkeumbae’ relative to normal pear fruit. Additionally, transgenic tobacco

82

overexpressing PpCAD2 increased the accumulation of lignin23. Therefore, lignin accumulation

83

in fruit changes fruit texture and increases its firmness. Although lignification of peaches occurs

84

in response to chilling injury24, less is known about senescence lignification in peach fruit.

85 86

The combination of transcriptomic and metabolomics was analyses is an effective way to identify

metabolic-related

functional

genes25.

Yu

et

ACS Paragon Plus Environment

al.26

used

a

combined

Page 5 of 49

Journal of Agricultural and Food Chemistry

87

transcriptomic/metabolomic approach to examine the response of Lemna aequinoctialis during a

88

time course (0, 3 and 7 days) study of nitrogen starvation and exogenously administered sucrose.

89

They reported finding that genes involved in nitrogen metabolism were the earliest transcriptomic

90

response to nitrogen starvation. The results of the transcriptomic and metabolomic analyses were

91

in agreement with each other, thus supporting the reliability and accuracy of the obtained data25.

92

Wang et al27 performed a combined transcriptomic/metabolomic analysis of ‘Hujingmilu’ peach

93

stored at either a constant 0◦C or LTC (pre-storage at 8◦C for 5 days before storage at 0◦C)

94

conditions, and found that LTC can alleviate chilling injury in peach fruit. Metabolomic data can

95

provide information on a large number of different metabolites, thus providing a good reference

96

for the analysis of transcriptome data28.

97

The peach varieties, ‘Baili’ and ‘Hongli’, used in present manuscript are called by a joint name

98

as ‘Feitao’ produced in Feicheng county, Shandong Province, China. ‘Feitao’ has been famous

99

for more than 1000 years in China in term of big, sweet and juicy. The eating quality of ‘Hongli’

100

is usually considered better than that of ‘Baili’ because ‘Hongli’ flesh is more juicy and soluble

101

while ‘Baili’ flesh is firmer. Now days it is this traditional good quality character of ‘Hongli’ that

102

is not good for handling during postharvest. In this manuscript, the mechanism of different

103

texture changes during storage of two peach varieties was worked out by metabolite and

104

transcriptome analysis. Our objective was to identify important factors that impact the texture of

105

these two varieties. The metabolomic and transcriptomic analyses were conducted on fruit stored

106

over a 10-day period at 20 ◦C. Differentially expressed genes (DEGs) and metabolites (DEMs)

107

were identified and categorized by functional classifications and metabolic pathways. These

108

classifications were further validated by examining the expression of specific enzyme genes.

109

MATERIALS AND METHODS

110

Plant Material and sampling

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

111

Two varieties of peach ‘Baili’ and ‘Hongli’ were harvested on Aug 26, 2016 at about 120 days

112

after anthesis (Total soluble solid reaches 12 %) from a forty-year-old orchard in Feicheng,

113

Shandong, China. Standard and identical irrigation and fertilization management programs were

114

used throughout the orchard during the growing season. Fruit that was uniform in appearance,

115

without visual defects, were harvested and stored at 20◦C. Fruit firmness, rate of ethylene

116

production, respiratory rate, cellulose content, and lignin content were measured after harvest,

117

Three biological replicates, consisting of ten pooled fruit for each biorep, were sampled for each

118

variety at each time point. The mesocarp was sliced into approximately 1 cm3 cubes, immediately

119

frozen in liquid nitrogen, and then stored at -80◦C until further analysis. According to the

120

difference in the firmness decline of the two varieties, we selected the 0, 2 and 6 DAH of the

121

‘Hongli’ fruit (firmness rapid decline period 0-2 d and decline in stationary period 6 d after

122

harvest), and the 0, 4 and 6 DAH of the ‘Baili’ fruit (change of firmness in the ‘Baili’ is steady,

123

we have chosen representative days of 0, 4 and 6 d after harvest), which were 18 analyzed

124

samples (two varieties x three storage time points x three biological replicates)for metabolomic

125

and transcriptomic analysis. The samples of ‘Baili’ collected at 0, 4, 6 DAH designated as B0, B4

126

and B6, and ‘Hongli’ at 0, 2, 6 DAH designated asH0, H2 and H6, respectively.

127

Measurement of ethylene, respiration, and firmness

128

The protocol used for measuring ethylene levels was performed as described by Yang et al29.

129

Nine fruits were placed in three 6 L sealed containers (three fruits one container) for 1h and the

130

extracted gas was measured for ethylene content, with measurements repeated three times for

131

each sampling time point. The rate of ethylene production (μL · kg−1· h−1) was calculated as the

132

released amount of ethylene per unit time per unit fresh weight (FW) of fruits. Ethylene was

133

measured with a GC-2010 gas chromatograph (Shimadzu) and chromatographic conditions

134

included the use of a FID detector, GDX-502 capillary column (10 m × 0.53 mm × 1 μm), a split

135

ratio of 10, hydrogen carrier gas, 120 ◦C column temperature, 200 ◦C detector temperature, and a

ACS Paragon Plus Environment

Page 6 of 49

Page 7 of 49

136

Journal of Agricultural and Food Chemistry

1 mL injection volume. .

137

To assess the rate of respiration, nine fruit were placed in three 6 L sealed containers (three

138

fruits one container) for 1 h and the rate of respiration was measured with a COMBO 580

139

portable analyzer (Italy). The oxygen concentration should be adjusted to 20.8% when used the

140

machine, recording the carbon dioxide concentration at this time (C1), and then inserted the

141

needle connected to the instrument into the sealed container for measurement. After a few

142

minutes of stabilization, recording the concentration of carbon dioxide (C2) and oxygen at this

143

time, using the formula: Respiration rate (mg / (kg·h)) = carbon dioxide release (C2-C1) (mg) /

144

(fruit weight (kg) · time (h)). Firmness of the flesh on both sides of the fruit equator was

145

measured using a CT3-4500 Texture Analyzer (Brookfield, USA) (probe diameter 2.0 mm). Six

146

fruits were measured for each variety at each of the sampled time points. The firmness, ethylene

147

release rate and respiration rate were measured every two days after harvest.

148

Lignin content

149

Lignin content was assessed using an ELISA kit according to a previously published method30,

150

which was measured by Shanghai BIO-TECH CO, LTD (Shanghai, China). The supernatant was

151

centrifuged at 1000 x g for 20 minutes to obtain a supernatant for detection. First, all Standards

152

and Samples were added in duplicate to the Microelisa Strip plate and set standard wells, testing

153

sample wells. Added 50μl standard to standard well, 10μl testing sample and 40μl Sample Diluent

154

to testing sample well, Blank well didn’t add anything. Added 100μl of HRP-conjugate reagent to

155

each well, covered with an adhesive strip and incubated for 60 minutes at 37 ◦C. Aspirated each

156

well and wash, repeating the process four times for a total of five washes. Washed by filling each

157

well with Wash Solution (400μl) using a squirt bottle, manifold dispenser or autowasher. After

158

the last wash, removed any remaining Wash Solution by aspirating or decanting. Inverted the

159

plate and blotted it against clean paper towels. Added chromogen solution A 50μl and chromogen

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

160

solution B 50μl to each well. Gently mixed and incubated for 15 minutes at 37 ◦C (Protect from

161

light). Added 50μl Stop Solution to each well then the color in the wells should change from blue

162

to yellow. Read the Optical Density (O.D.) at 450 nm using a microtiter plate reader within 15

163

minutes. Lignin content was expressed as l/mL and three biological replicates were used for

164

each sample (variety) at each time point. Three biological replicates were used for each sample at

165

each time point.

166

Lignin staining

167

Lignin was stained with Weisner reagent (phloroglucinol/concentrated hydrochloric acid)

168

according to the method described by Blanco-Portaleset al31. The flesh was cut to a thickness of

169

approximately 2 mm, longitudinally from the fruit pedicel, and stained with the Weisner reagent

170

for 10 min. Three biological replicates were used for each sample at each time point.

171

RNA Isolation and cDNA library preparation

172

The ‘Baili’ fruit flesh at 0, 4 and 6 DAH and ‘Hongli’ fruit flesh at 0, 2 and 6 DAH were

173

selected for total RNA extract, which using the Trizol reagent kit (Invitrogen, China) according

174

to the manufacturer’s instructions. And three biological replicates were used for each sample at

175

each time point. The extracted total RNA was then treated with RNase-free DNase I (Takara,

176

Dalian, China) for 15 min at 25◦C to remove residual DNA. RNA quality was verified by RNase-

177

free agarose gel electrophoresis and total RNA concentration of each sample was measured on an

178

Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, USA) at 260 and 280 nm. Total

179

RNA was then enriched for mRNA utilizing Oligo (dT) beads. The enriched mRNA was

180

subsequently fragmented into short fragments using a fragmentation buffer and reverse-

181

transcribed into cDNA with random primers. Second-strand cDNA was synthesized using DNA

182

polymerase I, RNase H, dNTP and buffer. The cDNA libraries were prepared for sequencing

183

using a NEBNext Ultra RNA Library Prep Kit. Prepared libraries were sequenced using an

ACS Paragon Plus Environment

Page 8 of 49

Page 9 of 49

Journal of Agricultural and Food Chemistry

184

Illumina HiSeq TM 2500 platform at Gene Denovo Biotechnology Co. (Guangzhou, China).

185

RNA-seq data processing and analysis

186

Reads containing adapters, sequences with more than 10% unknown nucleotides (N), and a

187

quality rating less than 50% (Q-value ≤20) were first removed from the data set. Bowtie232 was

188

then used to map the reads to a ribosome RNA (rRNA) database. All sequences with rRNA

189

mapped reads were then removed from the data set and the remaining reads were used in the

190

transcriptome analysis. The reads remaining after filtering were mapped to the peach reference

191

genome (“Lovell”, https://www.ncbi.nlm.nih.gov/genome/?term=prunus+persica)by TopHat233

192

(version 2.0.3.12). Novel transcripts were identified with Cufflinks34 and TopHat2 software.

193

Identification of differentially expressed genes (DEGs)

194

Transcript abundance was quantified using RSEM software35 and gene expression levels were

195

normalized using the FPKM (Fragments Per Kilobase of transcript per Million mapped reads)

196

method. The EdgeR package (http://www.r-project.org/ ) was used to identify DEGs across

197

samples or groups36. Genes with a fold change ≥2 and a false discovery rate (FDR) 0.7 or < -0.7, P-value < 0.05) are

511

presented in Figure 9 and Supplementary Table S4.In total, 20 lignin synthesis-related genes, 25

512

ethylene metabolism genes, and 22 cell wall remodeling-related genes were analyzed. The

513

relationship between the majority of lignin synthesis-related and cell-wall remodeling genes was

514

negatively correlated. In contrast, the correlation between ethylene metabolism-related and cell

515

wall remodeling-related genes was variable.

516

We inferred that the high lignin content in ‘Baili’ fruit may be induced by the high level of

517

ethylene production. Ethylene has been extensively reported to play a dominant role in promoting

518

fruit softening and senescence. Ma, however, also provided evidence that suggested the

519

APETALA2/ethylene response factor (AP2/ERF) family may be involved in the biosynthesis of

520

lignin in Isatis indigotica Fortune65,66. Additionally, exogenous treatment with ethylene has been

521

reported to promote ethylene release in bamboo shoots, increases PAL, CAD, and POD activity,

522

and significantly increases lignin content67. Wang et al.68 isolated three ethylene receptor genes

523

EjETR1, EjCTR1, and EjEIL1 from ‘Luoyangqing’ loquat fruit that exhibited increased

524

expression in response to low temperature. The expression of EjETR1, and EjEIL1 were

525

especially linked to chilling injury in ‘Luoyangqing’ fruit, which confirmed a role for ethylene in

526

the regulation of lignin synthesis in loquat fruit, as increased lignification is a direct consequence

ACS Paragon Plus Environment

Page 22 of 49

Page 23 of 49

Journal of Agricultural and Food Chemistry

527

of chilling injury. In the current study of peach, the rate of ethylene production was higher in

528

‘Baili’ fruit than in ‘Hongli’ fruit (Figure 1c); the former of which also exhibited higher levels of

529

lignification. The network analysis also indicated a positive correlation between the expression

530

level of ethylene metabolism-related and lignin synthesis-related genes (Figure 9b). These data

531

indicate that ethylene may affect lignin synthesis which contributes to the higher level of firmness

532

in ‘Baili’ fruit. A systematic analysis of all these components will need to be conducted, however,

533

to further validate this premise.

534

Transcriptome and metabolite profiles were used in a combinatorial analysis to provide

535

information on the underlying basis regulating the different texture of ‘Baili’ vs ‘Hongli’ peach.

536

A comparison of the differences in the level of metabolites and gene expression in the two

537

cultivars identified 644 metabolites and 4664 pathway annotations. The analysis identified the

538

metabolites and genes that may regulate ripening and softening in peach. The results suggested

539

that lignification may be a major factor impacting the contrasting texture of ‘Baili’ vs ‘Hongli’

540

peach fruit. Results also indicated that ethylene may have the double effect on fruit texture in

541

peach. This study provides a foundation for studying the regulation of peach ripening and

542

softening. In addition, the provided transcriptomic and metabolomic data also greatly enhance the

543

molecular and metabolite information that are now available for ‘Baili’ and ‘Hongli’ peach fruit.

544

Our initial investigation has provided a solid framework of data and helped to clarify the specific

545

studies that will be needed to fully understand the underlying mechanism of the differences in

546

texture between ‘Baili’ and ‘Hongli’ peach fruit.

547

ABBREVIATIONS USED

548

DEMs, differentially expressed metabolites; DEGs,differentially expressed genes; 4CL, 4-

549

Coumarate: Coenzyme A Ligase; COMT, caffeic acid 5-hydroxyferulic acid O-methyltransferas;

550

CAD, cinnamyl alcohol dehydrogenase; POD, peroxidase; Nr, NCBI protein database;GO, Gene

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 24 of 49

551

Ontology; KEGG, Kyoto Gene and Genome Encyclopedia; UPLC-QQQ-MS/MS, ultra-high

552

performance liquid chromatography coupled to a triple quadrupole mass spectrometer; SA,

553

salicylic acid; MeJA, methyl jasmonate; CHS, chalcone synthase; F3H, flavanone 3-hydroxylase;

554

AP2/ERF,

APETALA2/ethylene

ACS Paragon Plus Environment

response

factor;

Page 25 of 49

Journal of Agricultural and Food Chemistry

555

556

CONTRIBUTION

557

Shaolan Yang and Yongbing Yuan conceived and designed the experiments, Ying Wang and

558

Xinfu Zhang performed the experiments and analyzed the data, Ying Wang and Shaolan Yang

559

wrote the manuscript. All authors read and approved the manuscript.

560

ACKNOWLEDGMENTS

561

This work was supported by the National Key Research and Development Program of China

562

(2016YFD0400100), the Project of Shandong Natural Science Foundation (ZR2017MC006), and

563

Shandong Modern Fruit Technology Industry System project (SDAIT-06-06).

564

CONFLICT OF INTEREST

565

The authors declare no conflict of interest.

566

Supporting Information.

567

Fig. S1 Principal component analysis (PCA) of peach metabolite profiles obtained from 18

568

independent samples collected from ‘Baili’ and ‘Hongli’ at 0, 4(0), and 6 d.

569

Figs. S2 and S3 Heat maps of different metabolite clusters and the clustering of all metabolites.

570

Heat maps of DEMs based on normalized FPKM values in three time points. Green indicates

571

lower expression, and red indicates higher expression.

572

Fig. S4 The postharvest decay rate in ‘Baili’ and ‘Hongli’ fruit.

573

Fig. S5 GO enrichment of DEGs between ‘Baili’ and ‘Hongli’ fruit.

574

Fig. S6 and Fig. S7 KEGG pathway analysis of DEGs.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

575

Table S1. List of gene-specific primers used in the qRT-PCR analysis.

576

TableS2. RNA-seq tabular data (number of reads, filtered reads, clean reads, etc.).

577

Table S3. Mapping results of all 18 RNA-seq libraries to the peach reference genome.

578

Table S4. Pearson correlation coefficients between lignin synthesis-related, cell wall remodeling-

579

related, and ethylene metabolism-related gene expression.

580

ACS Paragon Plus Environment

Page 26 of 49

Page 27 of 49

Journal of Agricultural and Food Chemistry

582

583

REFERENCES

584

(1) Li, J.; Li, F.; Qian, M.; Han, M.; Liu, H.; Zhang, D.; Ma, J.; Zhao, C. Characteristics and

585

regulatory pathway of the PrupeSEP1 SEPALLATA gene during ripening and softening in peach

586

fruits. Plant Science 2017, 257, 63–73. DOI: 10.1016/j.plantsci.2017.01.004

587

(2) Brookfield, P. L.; Nicoll, S.; Gunson, F. A.; Harker, F. R.; Wohlers, M. Sensory evaluation by

588

small postharvest teams and the relationship with instrumental measurements of apple texture.

589

Postharvest Biology & Technology 2011, 59 (2), 179–186. DOI: 10.1016/j.postharvbio. 2010.

590

08.021

591

(3) Corollaro, M. L.; Aprea, E.; Endrizzi, I.; Betta, E.; Demattè, M. L.; Charles, M.; Bergamaschi,

592

M.; Costa, F.; Biasioli, F.; Grappadelli, L. C. A combined sensory-instrumental tool for apple

593

quality evaluation. Postharvest Biology & Technology 2014, 96 (2), 135–144. DOI: 10.1016/j.

594

postharvbio.2014.05.016

595

(4) Toivonen, P. M. A.; Brummell, D. A. Biochemical bases of appearance and texture changes in

596

fresh-cut fruit and vegetables. Postharvest Biology & Technology 2008, 48 (1), 1–14.

597

DOI:10.1016/j.postharvbio.2007.09.004

598

(5) Li, X.; Xu, C.; Korban, S. S.; Chen, K. Regulatory mechanisms of textural changes in

599

ripening fruits. Critical Reviews in Plant Sciences 2010, 29: 222–243. DOI:10.1080/07352689

600

.2010.487776

601

(6) Contador, L.; Díaz, M.; Millanao, M.; Hernández, E.; Shinya, P.; Sáenz, C.; Infante, R. A

602

proposal for determining the flesh softening of peach and nectarine in postharvest through

603

simplified targeted modeling. Scientia Horticulturae 2016, 209, 47–52. DOI:10.1016/j.scienta.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

604

2016.06.015

605

(7) Brummell, D. A.; Cin, V. D.; Lurie, S.; Crisosto, C. H.; Labavitch, J. M. Cell wall metabolism

606

during the development of chilling injury in cold-stored peach fruit: association of mealiness with

607

arrested disassembly of cell wall pectins. Journal of experimental botany 2004, 55 (405), 2041–

608

2052. DOI: 10.1093/jxb/erh228

609

(8) Qian, M.; Zhang, Y.; Yan, X.; Han, M.; Li, J.; Li, F.; Zhang, D.; Zhao, C. Identification and

610

expression analysis of polygalacturonase family members during peach fruit softening.

611

International Journal of Molecular Sciences 2016, 17 (11), 1933–1950. DOI: 10.3390/ijms

612

17111933

613

(9) Abu-Sarra, A. F.; Abu-Goukh, A. A. Changes in pectinesterase, polygalacturonase and

614

cellulase activity during mango fruit ripening. Journal of Pomology & Horticultural Science 2016,

615

67 (4), 561–568.DOI: 10.1080/00221589.1992.11516284

616

(10) Zaharah; Sakimin, S.; Reid; James, B. Mode of action of abscisic acid in triggering ethylene

617

biosynthesis andsoftening during ripening in mango fruit. Postharvest Biology & Technology

618

2013, 75 (2), 37–44. DOI: 10.1016/j.postharvbio.2012.07.009

619

(11) Hayama, H.; Ito, A.; Moriguchi, T.; Kashimura, Y. Identification of a new expansin gene

620

closely associated with peach fruit softening. Postharvest Biology & Technology 2003, 29 (1), 1–

621

10. DOI: 10.1016/S0925-5214(02)00216-8

622

(12) Tucker, G.; Yin, X, R; Zhang, A.; Wang, M. M.; Zhu, Q.; Liu, X.; Xie, X.; Chen, K, S.;

623

Grierson, D. Ethylene and fruit softening. 2017, 1 (4), 253–267. DOI: 10.1093/fqsafe/fyx024

624

(13) Hayama, H.; Shimada, T.; Fujii, H.; Ito, A.; Kashimura, Y. Ethylene-regulation of fruit

625

softening and softening-related genes in peach. Journal of experimental botany 2006, 57 (15),

626

4071–4077. DOI: 10.1093/jxb/erl178

ACS Paragon Plus Environment

Page 28 of 49

Page 29 of 49

Journal of Agricultural and Food Chemistry

627

(14) Razzaq, K.; Singh, Z.; Khan, A. S.; Khan, S. A. K. U.; Ullah, S. Role of 1-MCP in regulating

628

‘Kensington Pride’ mango fruit softening and ripening. Plant Growth Regulation 2016, 78 (3),

629

401–411. DOI: 10.1007/s10725-015-0101-7

630

(15) Luo, H.; Dai, S. J.; Ren, J.; Zhang, C. X.; Ding, Y.; Li, Z.; Sun, Y.; Ji, K.; Wang, Y.; Li, Q.

631

The Role of ABA in the Maturation and Postharvest Life of a Nonclimacteric Sweet Cherry Fruit.

632

Journal of Plant Growth Regulation 2014, 33, 373-383. DOI: 10.1007/s00344-013-9388-7

633

(16) Soto, A; Ruiz, K, B; Ravaglia, D; Costa, G; Torrigiani, P. ABA may promote or delay peach

634

fruit ripening through modulation of ripening- and hormone-related gene expression depending

635

on the developmental stage. Plant Physiology & Biochemistry Ppb, 2013, 64(5):11-24. DOI :

636

10.1016/j.plaphy.2012.12.011

637

(17) Haji, T. Inheritance of flesh texture in peach and effects of ethylene treatment on softening

638

of the stony hard peach. Japan Agricultural Research Quarterly Jarq 2014, 48 (1), 57–61.

639

DOI:10.6090/jarq.48.57

640

(18) Cai, C.; Xu, C. J.; Li, X.; Ferguson, I.; Chen, K. S. Accumulation of lignin in relation to

641

change in activities of lignification enzymes in loquat fruit flesh after harvest. Postharvest

642

Biology & Technology 2006, 40 (2), 163–169. DOI: 10.1016/j.postharvbio.2005.12.009

643

(19) Jin, P.; Zhang, Y.; Shan, T.; Huang, Y.; Xu, J.; Zheng, Y. Low-temperature conditioning

644

alleviates chilling injury in loquat fruit and regulates glycinebetaine content and energy status. J

645

Agric Food Chem 2015,63 (14), 3654–3659. DOI: 10.1021/acs.jafc.5b00605

646

(20) Li, X.; Zang, C.; Ge, H.; Zhang, J.; Grierson, D.; Yin, X. R.; Chen, K. S. Involvement of

647

PAL, C4H, and 4CL in chilling injury-induced flesh lignification of loquat fruit. Hortscience A

648

Publication of the American Society for Horticultural Science 2017, 52 (1), 127–131. DOI:

649

10.21273/HORTSCI11304-16

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

650

(21) Cai, Y. P.; Li, G. Q.; Nie, J. Q.; Yi, L.; Fan, N.; Zhang, J. Y.; Xu, Y. L. Study of the

651

structure and biosynthetic pathway of lignin in stone cells of pear. Scientia Horticulturae 2010,

652

125 (3), 374–379. DOI: 10.1016/j.scienta.2010.04.029 

653

(22) Kamdee, C.; Imsabai, W.; Kirk, R.; Allan, A. C.; Ferguson, I. B.; Ketsa, S. Regulation of

654

lignin biosynthesis in fruit pericarp hardening of mangosteen ( Garcinia mangostana L.) after

655

impact. Postharvest Biology & Technology 2014, 97 (11), 68–76. DOI: 10.1016/j.postharvbio.

656

2014.06.004

657

(23) Lu, G. L.; Li, Z, J.; Zhang, X, F.; Wang, R.; Yang, S, L. Expression analysis of lignin-

658

associated genes in hard end pear (Pyrus pyrifolia Whangkeumbae) and its response to calcium

659

chloride treatment conditions. Journal of Plant Growth Regulation. 2015, 34, 251–262. DOI:

660

10.1007/s00344-014-9461-x

661

(24) Puig, C. P.; Dagar, A.; Ibanez, C. M.; Singh, V.; Crisosto, C. H.; Friedman, H.; Lurie, S.;

662

Granell, A. Pre-symptomatic transcriptome changes during cold storage of chilling sensitive and

663

resistant peach cultivars to elucidate chilling injury mechanisms. BMC Genomics 2015, 16, 245–

664

280. DOI: 10.1186/s12864-015-1395-6

665

(25) Lin, Q.; Wang, C, Y.; Dong, W, C; Jiang, Q.; Wang, D, L.; Li, S, J.; Chen, M.; Liu, C, R.;

666

Sun, C, D.; Chen, K, S. Transcriptome and metabolome analyses of sugar and organic acid

667

metabolism in Ponkan (Citrus reticulata) fruit during fruit maturation. Gene 2015, 554, 64–

668

74.DOI: 10.1016/j.gene.2014.10.025

669

(26) Yu, C.; Zhao, X.; Qi, G.; Bai, Z.; Wang, Y.; Wang, S.; Ma, Y.; Liu, Q.; Hu, R.; Zhou, G, K.

670

Integrated analysis of transcriptome and metabolites reveals an essential role of metabolic flux in

671

starch accumulation under nitrogen starvation in duckweed. Biotechnology for Biofuels 2017, 10,

672

167–181.DOI: 10.1186/s13068-017-0851-8

673

(27) Wang, K.; Yin, X. R.; Zhang, B.; Grierson, D.; Xu, C. J.; Chen, K. S. Transcriptomic and

ACS Paragon Plus Environment

Page 30 of 49

Page 31 of 49

Journal of Agricultural and Food Chemistry

674

metabolic analyses provide new insights into chilling injury in peach fruit. Plant Cell &

675

Environment 2017, 40 (8), 1531–1551. DOI:10.1111/pce.12951

676

(28) Chen, J, F.; Dong, X.; Li, Q.; Zhou, X.; Gao, S, H.; Chen, R, B.; Sun, L, M.; Zhang, L.;

677

Chen, W, S. Biosynthesis of the active compounds of Isatis indigotica based on transcriptome

678

sequencing and metabolites profiling. BMC Genomics 2013, 14, 857–857. DOI: 10.1186/1471-

679

2164-14-857

680

(29) Yang, S, L.; Zhang, X, F.; Lu, G, L.; Wang, C, R.; Wang R. Regulation of gibberellin on

681

gene expressions related with the lignin biosynthesis in ‘Wangkumbae’ pear ( Pyrus pyrifolia

682

Nakai) fruit. Plant Growth Regulation 2015, 76, 127–134. DOI:10.1007/s10725-014-9982-0

683

(30) Kapat, A.; Dey, S. An alternative approach to the detection of lignin: A note on the

684

application of ELISA using polyclonal antibodies. Bioprocess Engineering 2000, 22, 75–77. DOI:

685

10.1007/PL00009105

686

(31) Richardson, D, G.; Lombard, P, B. Cork spot of Anjou pear: control by calcium sprays.

687

Commun Soil Sci Plant Anal 1979, 10, 383–389. DOI: 10.1080/00103627909366902

688

(32) Langmead, B.; Alzberg, S, L. Fast gapped-read alignment with Bowtie 2. Nature methods

689

2012, 9, 357–359. DOI: 10.1038/nmeth.1923

690

(33) Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2:

691

accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.

692

Genome biology 2013, 14, R36. DOI: 10.1186/gb-2013-14-4-r36

693

(34) Trapnell, C.; Roberts, A.; Loyal, G.; Pertea, G.; Kim, D.; Kelley, D, R.; Pimentel, H.;

694

Salzberg, S, L.; Rinn, J, L.; Pachter, L. Differential gene and transcript expression analysis of

695

RNA-seq experiments with TopHat and Cufflinks. Nature protocols 2012, 7, 562–578. DOI:

696

10.1038/nprot.2012.016

697

(35) Li, B.; Dewey, C, N. RSEM: accurate transcript quantification from RNA-Seq data with or

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

698

without a reference genome. BMC bioinformatics 2011, 12, 323–323. DOI: 10.1186/1471-2105-

699

12-323

700

(36) Cosgrove, D. J. Plant cell wall extensibility: connecting plant cell growth with cell wall

701

structure, mechanics, and the action of wall-modifying enzymes. Journal of Experimental

702

Botany2016, 67, 463–476. DOI:10.1093/jxb/erv511

703

(37) Young, M. D.; Wakefield, M. J.; Smyth, G. K.; Oshlack, A. Gene ontology analysis for

704

RNA-seq: accounting for selection bias. Genome Biology 2010, 11, R14–R14. DOI: 10.1186/gb-

705

2010-11-2-r14

706

(38) Wang, A.; Li, R.; Ren, L.; Gao, X.; Zhang, Y.; Ma, Z.; Ma, D.; Luo, Y. A comparative

707

metabolomics study of flavonoids in sweet potato with different flesh colors ( Ipomoea batatas

708

(L.) Lam). Food Chemistry 2018, 260, 124–134. DOI: 10.1016/j.foodchem.2018.03.125

709

(39) Chen, W.; Gong, L.; Guo, Z.; Wang, W.; Zhang, H.; Liu, X.; Yu, S.; Xiong, L.; Luo, J. A

710

novel integrated method for large-scale detection, identification, and quantification of widely

711

targeted metabolites: application in the study of rice metabolomics. Molecular Plant 2013, 6,

712

1769–1780. DOI: 10.1093/mp/sst080

713

(40) Livak, K. J.; Schmittgen, T. D. Analysis of relative gene expression data using real-time

714

quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408.

715

DOI:10.1006/meth.2001.1262

716

(41) Haji T.; Yaegaki, H.; Yamaguchi, M. Softening of stony hard peach by ethylene and the

717

induction of endogenous ethylene by 1- aminocyclopropane-1-carboxylic acid (ACC).

718

J JPN Soc Hortic Sci, 2003, 72(3): 212-217. DOI:10.2503/jjshs.72.212

719

(42) Junqueira-Gonçalves, M.P.; Yáñez, L.; Morales, C.; Navarro, M.; Contreras, R.; Zúñiga, G.

720

E. Isolation and Characterization of Phenolic Compounds and Anthocyanins from Murta (Ugni

721

molinae Turcz.) Fruits. Assessment of Antioxidant and Antibacterial Activity. Molecules, 2015,

ACS Paragon Plus Environment

Page 32 of 49

Page 33 of 49

Journal of Agricultural and Food Chemistry

722

20(4): 257-269. DOI:10.3390/molecules20045698

723

(43) Wei, J.; Wen, X.; Tang, L. Effect of methyl jasmonic acid on peach fruit ripening progress.

724

Scientia Horticulturae 2017,220, 206–213. DOI:10.1016/j.scienta.2017.03.004

725

(44) Ji, K.; Kai, W.; Zhao, B.; Sun, Y.; Yuan, B.; Dai, S.; Li, Q.; Chen, P.; Wang, Y.; Pei, Y.

726

SlNCED1 and SlCYP707A2: key genes involved in ABA metabolism during tomato fruit

727

ripening. Journal of experimental botany 2014, 65, 5243-5255. DOI: 10.1093/jxb/eru288

728

(45) Zhang, M.; Yuan, B.; Leng, P. The role of ABA in triggering ethylene biosynthesis and

729

ripening of tomato fruit. Journal of experimental botany 2009, 60, 1579-1588. DOI:

730

10.1093/jxb/erp026

731

(46) Hentrich, M.; Sánchezparra, B.; Pérez Alonso, M. M.; Carrasco, L. V.; Carrillo, L.;

732

Vicentecarbajosa, J.; Medina, J.; Pollmann, S. YUCCA8 and YUCCA9 overexpression reveals a

733

link between auxin signaling and lignification through the induction of ethylene biosynthesis.

734

Plant Signaling & Behavior 2013, 8, e26363. DOI: 10.4161/psb.26363

735

(47) Naidoo, R.; Ferreira, L.; Berger, D. K.; Myburg, A. A.; Naidoo, S. The identification and

736

differential expression of eucalyptus grandis pathogenesis-related genes in response to salicylic

737

acid

738

10.3389/fpls2013.00043

739

(48) Costa, M. A.; Collin, R. E.; Anterola, A. M.; Cochrane, F. C.; Davin, L. B.; Lewis, N. G. An

740

in silico assessment of gene function and organization of the phenylpropanoid pathway metabolic

741

networks in Arabidopsis thaliana and limitations thereof. Phytochemistry. 2003, 64, 1097–1112.

742

DOI: 10.1016/S0031-9422(03)00517-X

743

(49) Soltani, B. M.; Ehlting, J.; Hamberger, B.; Douglas, C. J. Multiple cis -regulatory elements

744

regulate distinct and complex patterns of developmental and wound-induced expression of

and

methyl

jasmonate.

Frontiers

in

Plant

Science

ACS Paragon Plus Environment

2013,

4,

43–52.

DOI:

Journal of Agricultural and Food Chemistry

745

Arabidopsis thaliana 4CL gene family members. Planta 2006, 224, 1239–1240. DOI:

746

10.1007/s00425-006-0296-y

747

(50) Sánchez, J. P.; Ullman, C.; Moore, M.; Choo, Y.; Chua, N. H. Regulation of Arabidopsis

748

thaliana 4-coumarate: coenzyme-A ligase-1 expression by artificial zinc finger chimeras. Plant

749

Biotechnology Journal 2010, 4, 103–114. DOI: 10.1111/j.1467-7652.2005.00161.x

750

(51) Yahiaoui, N.; Marque, C.; Myton, K. E.; Negrel, J.; Boudet, A. M. Impact of different levels

751

of cinnamyl alcohol dehydrogenase down-regulation on lignins of transgenic tobacco plants.

752

Planta 1997, 204, 8–15. DOI: 10.1007/s004250050224

753

(52) Lapierre, C.; Pilate, G.; Pollet, B.; Mila, I.; Leplé, J. C.; Jouanin, L.; Kim, H.; Ralph, J.

754

Signatures of cinnamyl alcohol dehydrogenase deficiency in poplar lignins. Phytochemistry 2004,

755

65, 313–321. DOI:10.1016/j.phytochem.2003.11.007

756

(53) Gavnholt, B.; Larsen, K. Molecular biology of plant laccases in relation to lignin formation.

757

Physiologia Plantarum 2010, 116, 273–280. DOI: 10.1007/s11105-013-0660-6

758

(54) Boerjan, W.; Ralph, J.; Baucher, M. Lignin biosynthesis. Annual Review of Plant Biology

759

2003, 54, 519–546. DOI: 10.1146/annurev.arplant.54.031902.134938

760

(55) Kim, S. H.; Nam, E. Y.; Cho, K. H.; Shin, I. S.; Kim, H. R.; Hwang, H. S. Comparison of

761

transcriptome analysis between red flash peach cultivar and white flash peach cultivar using next

762

generation sequencing. Journal of Plant Biotechnology. 2012, 39, 273–280. DOI: 10.5010/JPB.

763

2012.39.4.273

764

(56) Cao, K. F.; Ding, T. Y.; Mao, D. M.; Zhu, G. R.; Fang, W. C.; Chen, C. W.; Wang, X. W.;

765

Wang, L. R. Transcriptome analysis reveals novel genes involved in anthocyanin biosynthesis in

766

the flesh of peach. Plant Physiology & Biochemistry Ppb. 2017, 123, 94–102.DOI:

767

10.1016/j.plaphy.2017.12.005

ACS Paragon Plus Environment

Page 34 of 49

Page 35 of 49

Journal of Agricultural and Food Chemistry

768

(57) Zhang, L.; Zhu, L. X.; Xu, C.; Cui, C. M.; Sheng, H. Y.; Li, R.; Wang, H. Q. The effect of

769

silencing chalcone synthase on anthocyanin metabolism in peach. Acta Horticulturae Sinica 2015,

770

42, 31–37. DOI: 10.16420/j.issn.0513-353x.2014-0554

771

(58) Rahim, M.A.; Busatto, N.; Trainotti, L. Regulation of anthocyanin biosynthesis in peach

772

fruits. Planta 2014, 240, 913–929. DOI: 10.1007/s00425-014-2078-2

773

(59) Xu, Q.; Wang, W. Q.; Zeng, J. K.; Zhang, J.; Grierson, D.; Li, X.; Yin, X. R.; Chen, K. S. A

774

NAC transcription factor, EjNAC1, affects lignification of loquat fruit by regulating lignin.

775

Postharvest Biology & Technology 2015, 102, 25–31.DOI:10.1016/j.postharvbio.2015.02.002

776

(60) Baucher, M; Monties, B; Montagu, M. V; Boerjan, W. Biosyntensis and genetic engineering

777

of lignin. Critical Reviews in Plant Sciences 1998, 12, 125–197. DOI:10.1080/07352689891304

778

203

779

(61) Gavnholt, B.; Larsen, K. Molecular biology of plant laccases in relation to lignin formation.

780

Physiol Plant 2002, 116, 273–280. DOI: 10.1034/j.1399-3054.2002.1160301.x

781

(62) Fu, T. T.; Qin, W.; Zhang, Z. Z. Effects of chilling injury on storage quality and flesh

782

lignification of postharvest "Dawuxing" loquat Fruits. Journal of Sichuan Agricultural University.

783

2011, 19, 12–14. DOI: 10.1007/s00425-014-2078-2

784

(63) Tucker, G.; Yin, X. R.; Zhang, A. D.; Wang, M. M.; Zhu, Q. G.; Liu, X. F.; Xie, X. L.; Chen,

785

K. S.; Grierso, D. Ethylene and fruit softening. Food Quality and Safety 2017, 1: 253–267. DOI:

786

10.1093/fqsafe/fyx024

787

(64) Cajuste, J. F.; Lafuente, M. T., Ethylene-induced tolerance to non-chilling peel pitting as

788

related to phenolic metabolism and lignin content in ‘Navelate’ fruit. Postharvest Biology &

789

Technology 2007, 45, 193–203. DOI: 10.1016/j.postharvbio.2007.01.019

790

(65) Zeng, J. K.; Li, X. ; Xu, Q; Chen, J. Y.; Yin, X. R.; Ferguson, I. B.; Chen, K. S. EjAP2-1, an

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 36 of 49

791

AP2/ERF gene, is a novel regulator of fruit lignification induced by chilling injury, via interaction

792

with EjMYB transcription factors. Plant Biotechnology Journal 2016, 13,1325–1334. DOI:

793

10.1111/pbi.12351

794

(66) Ma, R. F.; Xiao, Y.; Lv, Z. Y.; Tan, H. X.; Chen, R. B.; Li, Q.; Chen, J. F.; Wang, Y.; Yin, J.;

795

Zhang, L.; Chen,W. S. AP2/ERF transcription factor, Ii049, positively regulates lignan

796

biosynthesis in Isatis indigotica through activating salicylic acid signaling and lignan/lignin

797

pathway genes. Frontiers in Plant Science 2017, 8, 1361–1377. DOI: 10.3389/fpls.2017.01361

798

(67) Luo, Z. S.; Xu, X. L.; Cai,Z. Z.; Yan, B. F. Effects of ethylene and 1-methylcyclopropene (1-

799

MCP) on lignification of postharvest bamboo shoot. Food Chemistry 2007, 105, 521–527. DOI:

800

10.1016/j.foodchem.2007.04.007

801

(68) Wang, P.; Zhang, B.; Li, X.; Xu, C. J.; Yin, X. R.; Shan, L. L.; Ferguson, I. B.; Chen, K. S.

802

Ethylene signal transduction elements involved in chilling injury in non-climacteric loquat fruit.

803

Journal

of

Experimental

Botany

2010,

1,

179–190.

ACS Paragon Plus Environment

DOI:

10.2307/24038554

Page 37 of 49

Journal of Agricultural and Food Chemistry

804 805

FIGURE CAPTION

806

FIGURE 1. Physiological characteristics of ‘Baili’ and ‘Hongli’ peach fruit. A. Appearance of

807

‘Baili’ and ‘Hongli’ peach fruit at harvest (day 0). ‘Baili’ (above) and ‘Hongli’ (below). B.

808

Firmness measurements of ‘Baili’ and ‘Hongli’ peach fruit over the course of the storage period.

809

C. Respiration rate of ‘Baili’ and ‘Hongli’ peach fruit over the course of the storage period. D.

810

Ethylene generation in‘Baili’ and ‘Hongli’ peach fruit over the course of the storage period.

811

FIGURE 2. Lignin and cellulose content in ‘Baili’ and ‘Hongli’ peach over the course of the

812

storage period. A. Phloroglucinol-HCl staining of ‘Baili’ and ‘Hongli’ mesocarp fruit tissues. B.

813

Lignin (a) and cellulose (b) content in ‘Baili’ and ‘Hongli’ fruit over the course of the storage

814

period.

815

FIGURE 3. Differentially expressed metabolites (DEMs) and KEGG pathway annotations of

816

metabolites. A. DEMs in ‘Hongli’ and ‘Baili’ fruit. B0 represents ‘Baili’ fruit at 0 days after

817

harvest (DAH), H0 represents ‘Hongli’ fruit at 0 days DAH. B4 represents ‘Baili’ fruit at 4 DAH,

818

H2 represents‘Hongli’ fruit at 2 DAH. B6 represents ‘Baili’ at 6 DAH. H6 represents ‘Hongli’

819

fruit at 6 DAH. B. KEGG pathway annotation of ‘Baili’ and ‘Hongli’ peach at 0, 4(2) and 6 DAH.

820

FIGURE 4. Identification of the major DEMs (adjusted P< 0.05, and absolute log2 Fold

821

Change> 1) present in ‘Baili’ and ‘Hongli’ fruit. A. Metabolites involved in the flavonoid and

822

anthocyanin biosynthesis pathway. B. Metabolites involved in the synthesis of monolignols. C.

823

Metabolites involved in the biosynthesis of plant hormones. D. Metabolites that are precursors of

824

ethylene. The rows in the Heatmap represent metabolites, and the columns indicate samples. The

825

colours of Heatmap cells indicate scaled expression level of genes across different samples. The

826

colour gradient, ranging from green, through black, to red represents low, middle and high values

827

of the FPKM values over the course of the storage period.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

828

FIGURE 5. Summary of differentially expressed genes (DEGs) in ‘Baili’ and ‘Hongli’ over the

829

course of the storage period. A. DEGs in ‘Hongli’ vs ‘Baili’ fruit. B. Venn diagrams of DEGs. C.

830

Principal component analysis (PCA) of DEGs in ‘Hongli’ and ‘Baili’ fruit over the course of the

831

storage period. B1 represents‘Baili’ fruit at 0 days after harvest (DAH), H1 represents ‘Hongli’

832

fruit at 0 DAH, B2 represents ‘Baili’ fruit at 4 DAH, H2 represents ‘Hongli’ fruit at 2 DAH. B3

833

represents ‘Baili’ fruit at 6 DAH, H3 represents ‘Hongli’ fruit at 6 DAH.

834

FIGURE 6. KEGG pathway analysis of DEGs and overview of differentially expressed

835

metabolites and genes between pairwise comparisons of ‘Baili’ and ‘Hongli’ fruit at 0, 2(4), 6

836

days after harvest (DAH). A. KEGG pathway analysis of DEGs. Red box indicates

837

phenylpropanoid metabolism pathway. B. Overview of differences in the metabolites and genes

838

expressed in ‘Baili’and ‘Hongli’ fruit. Metabolites pathways and genes that are up- or down-

839

regulated in ‘Baili’ and ‘Hongli’ are labeled in magenta, aqua blue, jade green,or dark gray,

840

respectively. The colors represent up-regulated pathway (log>2) in ‘Baili’, up-regulated pathway

841

(2>log>1) in ‘Hongli’, down-regulated pathway (log