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