Next Article in Journal
Selective Extraction and Antioxidant Properties of Thiol-Containing Peptides in Soy Glycinine Hydrolysates
Previous Article in Journal
Isolation and Purification of Potent Growth Inhibitors from Piper methysticum Root
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 23
Issue 8
10.3390/molecules23081908
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptomic Profiling to Understand Pre- and Post-Ripening Hormonal Regulations and Anthocyanin Biosynthesis in Early Ripening Apple Fruit

by
Jakaria Chowdhury Onik
,
Xiaojia Hu
,
Qiong Lin
* and
Zhidong Wang
Institute of Food Science and Technology, Chinese Academy of Agricultural Sciences/Key Laboratory of Agro-products Quality and Safety Control in Storage and Transport Process, Ministry of Agriculture and Rural Affairs, Beijing 100193, China
*
Author to whom correspondence should be addressed.
Molecules 2018, 23(8), 1908; https://doi.org/10.3390/molecules23081908
Submission received: 11 June 2018 / Revised: 27 July 2018 / Accepted: 28 July 2018 / Published: 31 July 2018
Browse Figures
Versions Notes

Abstract

:
The ‘Hongyu’ apple is an early ripening apple cultivar and usually used for fresh marketing. Due to the short ripening period, most of the fruit are harvested at the commercial maturity stage for proper marketing distribution and a longer shelf life. Fruit ripening involves delicate changes to its metabolic and physiological traits through well-organized synchronization of several hormones and regulatory steps. A clear understanding of these hormonal alterations is crucial for extending the period from commercial to physiological ripening. This study was intended to clarify the hormonal alterations and anthocyanin biosynthesis process prior to and immediate after, the harvesting of apple fruit considering the commercial maturity stage. Fruits harvested at 120 Days after flowering (DAF) (HY_4th) was considered as commercially ripened, 110 DAF (HY_3rd) as pre-ripening and 120 DAF followed by five days storage at 20 °C (HY_20 °C_5) as post-ripening samples. Three different stages of fruit were used for transcriptome assembly using RNA-Seq. Results revealed 9187 differentially expressed genes (DEGs) in the post-ripening samples, which was comparatively lower (922 DEGs) in the pre-ripening fruits. DEGs were subjected to Gene Ontology analysis and 31 categories were significantly enriched in the groups ‘biological process,’ ‘molecular function’ and ‘cellular component.’ The DEGs were involved in hormonal signaling pathways like ethylene, abscisic acid (ABA), auxin, gibberellin (GA), brassinosteroid (BR) and anthocyanin biosynthesis pathways such as PAL, 4CL, CHI, DFR, F3H, UFGT. Several transcription factors like the MADS-box gene, MYB, bHLH, NAC, WRKY and HSF were differentially expressed between the pre- and post-ripening fruits. Selected DEGs were subjected to gene expression analysis using quantitative RT-PCR (qRT-PCR) and the results were consistent with those of RNA-Seq. Our data suggested that in addition to ethylene, ABA and other hormones also play key roles in regulating apple fruit ripening and may interact with the ethylene signaling process. Additionally, our data provided an exhibition of the expression pattern of genes in the anthocyanin biosynthesis pathway.

1. Introduction

The apple is one of the important fruit crops all over the world, receiving attention for its delicious taste and considerable nutritive value. Apple fruit endures rapid ethylene production which triggers a dramatic change in many important quality attributes, including colors, sugars-acids ratio and others, during development and ripening [1]. Previous studies have confirmed a critical regulatory role of ethylene in fruit softening by cell wall degradation [2], biosynthesis of volatiles compounds providing aroma [3] and loss of green color by the accumulation of anthocyanin [4].
Plant hormones are reportedly considered to be closely linked with fruit development as well as ripening [5]. By the response of different hormones, the significant ripening regulations seem to be controlled primarily by ethylene and ABA [6,7]. Considering the crucial role of ethylene in climacteric fruit ripening, it remains the most researched hormone [8]. Previous studies on tomatoes provide a background of ethylene signal transduction and response [9]. Involvement of several genes in ethylene biosynthesis process is well noted. Rate limiting genes such as 1-aminocyclopropane-1-carboxylic acid synthase (ACS) and 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) were primarily regulated by RIN and MADS-box gene family [10]. The silencing of MdMADS8 resulted in a significant decrease of ACS1 and ACO1 and has been observed in the transgenic apple line, which indicated a strong correlation of MADS-box genes with ethylene biosynthesis [2].
However, understanding the roles of other hormones other than ethylene during fruit ripening is limited. A transcriptome study with different ripening stages of tomato fruits has revealed that next to ethylene, auxin-related genes were the most represented in hormone response category [11]. Interactive relations of ethylene and indole-3-acetic acid (IAA) in climacteric fruits ripening has been noted in several studies [11]. A declination of free auxin level (indole-3-acetic acid, IAA) before the ripening has been observed in both climacteric (tomato) and non-climacteric (pepper, strawberry) fruits and presumed this reduction caused by the increase of IAA-Asp (IAA-aspartic acid), a conjugated form of free auxin [12,13]. In tomato and peach fruits, ethylene biosynthesis genes (ACS2, ACS4 and ACO1) and ethylene signaling genes (ETRs and ERFs etc.) were up-regulated by auxin [14,15]. In the tomato, two IAA-amino synthase genes (GH3) exhibited a pattern of expression associated with ripening, which catalyzed the reaction of IAA-Asp (IAA-aspartic acid) conjugation [16]. Overexpression of pepper GH3 gene in tomato fruits directed towards anticipatory ripening process [17] which indicated the participation of conjugated form of IAA rather than only IAA, may contribute to the temporal regulation of ripening [13].
The profound roles of auxin in the ripening of some non-climacteric fruits have been well-studied; whereas little effort given to investigate the possible roles of other plant hormones as gibberellin (GA). However, an external application of GA3 on ripening strawberry fruits has significantly retarded red color development [18] and altered the expression of genes participate in cell enlargement process [19]. Abscisic acid (ABA) is a well-known phytohormone that takes part in several means of growth and development and stress response in plants [20]. Studies have reported the role of ABA in promoting sugar accumulation in fleshy fruits [21] and regulating the ripening of climacteric fruits. The down-regulation of cell wall degrading genes such as polygalacturonase (PG) was observed in tomato fruit due to the suppression of NCED1 (9-cis-epoxycarotenoid dioxygenase), which catalyzed the key steps in ABA biosynthesis resulting in an increase in firmness and to extend shelf-life [22]. Although several studies on brassinosteroid (BR) have been conducted, its role in fruit ripening is still indistinct.
Anthocyanins are secondary metabolites, which play a substantial part in pigmentation with many health benefits as an antioxidant and with anti-tumor properties [23]. The biosynthesis and accumulation of anthocyanins are determined by metabolic networks correlated with the expression of several genes and regulatory factors [24]. Characterization of key genes involved in the anthocyanin metabolic process has been studied by several researchers and has revealed structure genes in biosynthesis. Direct involvement of the genes encoding chalcone synthase (CHS) [25], chalcone isomerase (CHI) [26], dihydroflavonol-4-reductase (DFR) [27] and glucosyltransferase (UFGT) [28] has been found in the anthocyanin biosynthesis process. The expression patterns of these genes are differentially exhibited with different fruit developmental stages and ripening in different species [29]. Moreover, several transcription factors have been identified which are functional for regulating the anthocyanin biosynthesis process through interaction with structure genes characterized by genetic analysis in different fruits as strawberry and grape [30,31].
Overall, fruit ripening initiation is mainly characterized by ethylene synthesis, color changes and cell wall dynamics [32,33], whereas other ripening characteristics have received less attention and are yet to be investigated. McAtee et al. [6] mentioned that individual ripening processes themselves might be under specific hormonal control. Reports on the interaction between hormones at several levels above biological processes are limited and unclear. However, due to recent advances in using transcriptomic, proteomic and metabolomic analytical tools, significant progress has been achieved in the characterization of hormonal responses [34,35,36,37,38] which have been directed to the identification of the distinct expression of genes associated with hormonal response across plant species [22,39,40]. Total transcriptome analysis of fruits before the onset of ripening and immediately after ripening considering commercial maturity would be helpful to understand the changes and respective regulations more clearly. Our study involved a comparative transcriptome analysis in facilitating the understanding of ripening hormonal control in early ripening apple fruit. Several genes associated with plant hormones, the anthocyanin biosynthesis process and a few transcription factors were differentially expressed during fruit ripening. Also, the expression of selected differentially expressed genes was confirmed by quantitative RT-PCR (qRT-PCR).

2. Results

2.1. Changes in Organic Acid and Sugar Content in Pre- and Post-Ripening Apple

The main organic acids and sugars were measured at different stages of the ‘Hongyu’ apple. Two organic acids, including malate and oxalate, were measured in this analysis (Figure 1). From the results, it can be seen that malate was the most prominent organic acid in apple fruit. The oxalate contents were below 0.5 g kg−1, while the malate contents varied from 10 to 13 g kg−1. The highest malate content was observed in HY_3rd samples (12.92 g kg−1) which showed a decreasing trend in the HY_4th (11.41 g kg−1) and HY_20 °C_5 samples (10.83 g kg−1). Glucose, fructose and sucrose were the main sugars in apple fruit, in which fructose content was prominent. According to the results, glucose and sucrose showed an increasing pattern of ripening (HY_4th) and post-ripening (HY_20 °C_5) samples. The changes observed in this study agreed well with previous studies that also used the date of full bloom as a reference point [41,42,43].

2.2. RNA-seq and Deferential Gene Expression Analysis

Pre- and post-ripening apple samples were subjected to RNA-seq analysis. Fruits of 110 DAF denoted as HY_3rd, 120 DAF denoted as HY_4th and stored for 5 d at 20 °C indicated as HY_20 °C_5 groups in transcriptome analysis. cDNA libraries were constructed from the total RNA of pre-ripening and post-ripening apple fruits. With Illumina sequencing technology through pair-end reading, we got more than 6G bases of clean reads with a GC percentage above 47.19% and Q20 percentage above 95.16% after removing low quality reads and trimming the adapter sequences (Table 1).
To evaluate the significance of gene expression variation between pre- and post-ripening stages, the transcript abundance of each gene was normalized to the FPKM (Fragments Per Kilobase of transcript sequence per Millions of base pairs sequenced) value. The differentially expressed genes were identified and filtered for Q value <0.005 and log2Ratio >1. Total 922 unigenes were differentially expressed in pre-ripening (HY_4th vs. HY_3rd) samples where 578 genes were up-regulated and 344 genes were down-regulated between pre-ripening and ripening stages of apple fruit. The number of up-regulated and down-regulated unigenes is shown in the (Figure 2A, Table S1). Comparatively higher DEGs of 9187 unigenes were found in the post-ripening (HY_20 °C_5 vs. HY_4th) sample where 20 °C storage for 5 d up-regulated 4546 DEGs and downregulated 4641 genes (Figure 2A, Table S1). Overlapped DEGs between pre-and post-ripening groups are shown on a Venn diagram (Figure 2B).

2.3. GO Enrichment Classification and KEGG Pathway Analysis

Gene Ontology, an international standard classification system for gene function, was used to illustrate the role of DEGs in different groups. The DEGs were clustered into three main categories of the GO classification: Biological Process, Cellular Function and Molecular Function. In the pre-ripening group, 466 unigenes were assigned to the Biological Process category, 143 unigenes were assigned to the Cellular Component category and 351 unigenes were assigned to the Molecular Function category. These genes were further classified into 31 functional subcategories by mapped homology (Figure 3A, Table S2). The most common assignments in the Biological Process category were metabolic processes, cellular processes and organic substance metabolic processes. In the Cellular Component category, the majority of unigenes were grouped into the membrane, cell and cell part subcategories. Genes in the Molecular Function category were primarily sorted into the binding, catalytic activity and ion binding subcategories (Figure 3A, Table S2). In the post-ripening group, the most abundant GO terms in the Biological process include ‘metabolic process,’ ‘cellular process’ and ‘organic substance metabolic process.’ In the molecular function category, most of the DEGs were mapped onto the ‘binding,’ ‘catalytic activity’ and ‘organic cyclic compound binding’ groups. The cellular component category comprises the greatest numbers of genes in ‘cell,’ ‘cell part’ and ‘membrane’ groups (Figure 3B, Table S2).
To identify the biological pathways activated in selected apple fruits, we mapped the annotated sequences to the reference pathways in the KEGG database [44]. In the pre-ripening samples, significant matches were found for 574 unigenes, which were assigned to 100 KEGG pathways (Figure 4A; Table S3). The KEGG pathway analysis has revealed the top 20 most enriched metabolic pathway subcategories, including the most abundant group represented as the ‘Biosynthesis of other secondary metabolites,’ ‘plant hormone signal transduction’ and ‘flavonoid biosynthesis.’ On the other hand, post-ripening samples exhibited the top 20 pathways where the largest groups fell into the ‘carbon metabolism,’ ‘biosynthesis of amino acids,’ ‘pyruvate metabolism,’ ‘glycolysis/gluconeogenesis’ and ‘arginine and proline metabolism’ categories (Figure 4B; Table S3).

2.4. Confirmation of RNA-seq Results by qRT-PCR

Genes showing differential expression involved in ethylene signaling pathway and anthocyanin biosynthesis in the transcriptome analysis were selected for qRT-PCR comparison of their expression levels in HY_3rd, HY_4th and HY_20 °C_5 samples, to confirm the accuracy and reproducibility of the transcriptome analysis. The scatter plot indicated higher integrity and correlation between transcriptome analysis through RNA-seq and qRT-PCR results (Figure 5).

2.5. Identification of DEGs Related to Ripening Hormone

Plant hormones with diversified functions are obligatory for fruit development, maturation and ripening [6]. Several differentially expressed genes were identified as plant hormone regulators in our transcriptome study between pre- and post-ripening samples, which included ethylene, GA, auxin and ABA as well as BR (Table 2 and Table 3; Table S4). The most important hormone regulating fruit ripening is ethylene, which has been studied extensively. In the present study, the DEG analysis has revealed numerous regulatory genes which were directly associated with the ethylene signaling pathway. The post-ripening samples exhibited higher DEGs involved in the ethylene signaling pathway than pre-ripening samples (Table S4). Most of these genes, including ACO, EIN2, EIN3, ERS1, EIL and several ERFs showed higher expression in the post-ripening (Hy_20 °C_5) samples. Most of these were not expressed, or were expressed at a shallow level in the pre-ripening (HY_3rd) samples. This result showed consistency with the previous reports [7,45].
In addition to ethylene, other important plant hormones and the genes characterizing these hormones were also specified in the DEG analysis results. In the ABA biosynthesis pathway, the key catabolizing gene encoding NCED1 (9-cis-epoxycarotenoid dioxygenase) was highly expressed in the post-ripening samples, suggesting possible regulation of ABA in the ripening of apple. Moreover, another gene involving in the ABA signaling, encoding Protein Phosphatase 2C, showed higher expression in the post-ripening samples and several genes encoding serine/threonine protein kinases were downregulated (Table 3; Table S4).
Previous studies have stated that conjugation of IAA-Asp was primarily regulated by the association of GH3 (indole-3-acetic acid amino synthase) gene in the auxin signaling pathway [10]. The up-regulation of this gene in post-ripening samples indicates its role in decreasing free IAA during fruit ripening (Table 3; Table S4). Furthermore, up- or down-regulation of several auxin response factors (ARF) were observed in the post-ripening samples, which may have functional regulation and play different roles in fruit ripening (Table 3; Table S4). The inactivation of GA controlled by a primary gene encoding gibberellin 2-beta-dioxygenase (GA2ox1) was observed to be highly expressed in the post-ripening samples (Table 3; Table S4). In the BR pathway, several genes encoding cytochrome p450 were downregulated in the post-ripening samples and a lower expression of other genes in the BR pathway indicates inferior transcript abundance of these genes, representing its negative regulation of the fruit ripening process (Table 3; Table S4). A heat map including all the regulating genes involved in the ethylene signaling pathway are shown in Figure 6.

2.6. DEGs Involved in Anthocyanin Biosynthesis

‘Hongyu’ apple fruits undergo a progressive color alteration to red at the later stages of ripening. Anthocyanin most prominently changed the red color of apple fruit, comprising mostly of cyanidin 3-galactoside [46]. Studies with numerous crops, including pear and grape, have specified an imperative correlation of anthocyanin content with the expression of its biosynthetic genes [47,48]. Several DEGs related to flavonoid biosynthesis were identified in the KEGG pathway assessment, which has exposed the differential expression of the genes such as PAL, 4CL, CHI, F3H, DFR and UFGT involved in the anthocyanin biosynthesis pathway in post-ripening samples. In contrast, most were not expressed or were expressed at a low level in pre-ripening samples. In this study, 20 candidate genes that are responsible for anthocyanin synthesis in fruit were identified from the DEG data including PAL (2 unigene), 4CL (3 unigenes), CHI (3 unigene), F3H (5 unigene), DFR (4 unigene) and UFGT (3 unigenes) in post-ripening samples (Table S5). On the other hand, genes encoding PAL (2 unigenes), F3H (4 unigenes) and CHI (4 unigenes) were differentially expressed in pre-ripening samples. All of the characterized DEGs regarding anthocyanin biosynthesis showed higher expression in the post-ripening samples (Table S5). A previous qRT-PCR study of the genes associated with the anthocyanin biosynthesis process also showed similar results [49]. The genes regulating the last steps of anthocyanin biosynthesis encoding UFGT was also highly expressed in the post-ripening samples. This result indicates that UFGT may play a key role in apple fruit coloration during ripening.
A heat map involving the DEGs expressed at pre-and post-ripening apple fruit in anthocyanin biosynthesis process shown in Figure 7. The biosynthesis and accumulation of anthocyanin are primarily responsible for fruit coloration which might be controlled by the expression of these genes.

2.7. Transcriptional Regulation during Fruit Ripening

To date, many transcription factors have been specified to regulate fruit ripening processes [50]. A well-known transcription factor is the RIN gene, a MADS-box gene family. In our study, numerous transcription factors including the MADS-box gene, MYB, bHLH, NAC, WRKY, HSF were identified in the DEG analysis data between pre- and post-ripening samples. In the post-ripening samples, one MADS-box gene AGL8 was downregulated whereas the other two genes AGL15 and AGL80 were upregulated. A RIN-homologue of apple MADS8/9 has been found to regulate fruit ripening by directly controlling the auxin levels [51]. Mature MADS8/9-suppressed apples demonstrated reduced expression of GH3 and exhibited a higher concentration of free IAA [51]. In our results, the higher expression of two MADS-box genes and consequently up-regulation of GH3 genes were observed in post-ripening samples. Moreover, six MYB genes and five NAC genes were upregulated in the post-ripening samples. On the other hand, stress-related protein HSF also showed differential expression (up or down) in the post-ripening samples (Table S6). According to the results, several other bHLH transcription factors were expressed to be up or down-regulated. The differential expression of these regulatory transcription factors may control hormonal regulation and anthocyanin biosynthesis during the ripening process.

3. Discussion

Fruit ripening is considered to be a crucial parameter that influences the shelf life of fruit and its overall market value. The ripening process involves the initiation of multiple gene expression and biochemical pathways. However, the underlying molecular mechanism of their regulators remains to be uncovered. A deeper understanding of hormonal regulation and anthocyanin biosynthesis in the fruit ripening process provides the potential to investigate post-harvest management and fruit quality. In this study, we compared the fruit transcriptomes in pre- and post-ripening fruits of an early ripening apple cultivar and found that 9187 genes were differentially expressed after post-harvest storage for 5 d at 20 °C. Several DEGs were included in the GO terms involved in the biosynthesis and regulation of plant hormones including ethylene and the enrichment of the flavonoid biosynthesis pathway in the KEGG analysis.
Ethylene is supposed to be critical for controlling the ripening of climacteric fruits [52]. As an important ripening hormone, various studies on model climacterics fruits, such as the tomato and the banana, have focused on genes involved in the biosynthesis and signal transduction of ethylene in ripening [45,53]. In the present study, transcriptome analysis with the fruit prior to, and immediately after, ripening exposed large numbers of DEGs related to ripening. Using microarray analysis, 106 genes were found to be highly coordinated with fruit ripening in the tree ripening of the ‘Royal Gala’ apple [54]. DEGs specified in the ethylene and another hormone signaling process were mostly expressed in the post-ripening (HY_20 °C_5) samples. ABA is an important plant hormone, which has been anticipated to play vital roles in the ripening processes of both climacteric and non-climacteric fruits [46,55]. ABA promotes ripening by promoting ethylene biosynthesis through the up-regulation of ethylene biosynthesis genes [22]. It has been considered a ripening-inducer in strawberry and grapefruits [55,56]. The expression of the key gene for ABA biosynthesis encoding NCED1 increases at the breaker stage of tomato fruits and suggests an important role for ABA in triggering ethylene biosynthesis [39]. Furthermore, a longer shelf life of tomato fruits was observed due to the down-regulation of some ripening-related cell wall genes mediated by the suppression of LeNCED1 [22]. Similarly, the reduction of NCED expression was correlated with retardation of ripening in the strawberry [55]. In our results, genes encoding NCED1 were up-regulated in the post-ripening samples and showed consistency with previous studies. In another study, the application of 1-MCP on ‘nanguo’ pear fruits extended the shelf life by subsequently suppressing the expression of the NCED1 gene [57]. The expression pattern of NCED1, which directly correlated with ABA levels, was observed in tomato (SlNCED1), persimmon (DkNCED1), sweet cherry (PacNCED1) [58] and cucumber (CsNCED1), which declined rapidly during the postharvest stage [22]. These findings also provide substantial evidence of ABA as critical and affects the regulation of fruit ripening. Another most shared plant hormone is auxin, found as a free form of IAA (indole-3-acetic acid IAA) and reported to decline in concentration prior to fruit ripening [12,59]. However, conjugation of IAA and the formation of IAA-Asp (IAA-aspartic acid), mediated by GH3 (IAA-amido synthetase,) increased in ripened fruit [13]. Similarly, the up-regulation of GH3.1, GH3.5, GH3.9 and GH3.17 in the strawberry seed tissues with high auxin levels indicated that even IAA biosynthesis might induce IAA conjugation by promoting the expression of members of this gene family [60]. Furthermore, a higher expression of GH3 and subsequent delays of fruit ripening with IAA treatment were observed in the onset of the ripening of grape berries [13]. Our results also supported this where GH3 was up-regulated in the post-ripening samples. These findings directed the functional roles of GH3 in reducing the IAA level and promoting fruit ripening. Studies with the regulation of auxin in transcript level indicated that an increase in firmness and prolonged shelf life was mediated by the down-regulation of the ARF4 gene in the tomato [61]. In our results, we also found that ARF genes were downregulated in the post-ripening samples, suggesting a role in fruit ripening.
From our DEG analysis results, we also found several key genes involved in the anthocyanin biosynthesis pathway, including PAL, 4CL, CHI, F3H and UFGT. The distinct functions of these genes have been characterized in previous studies. In a study with a green and red apple cultivar, it was observed that the genes involved in the anthocyanin biosynthesis pathway such as CHS, F3H, DFR and UFGT, were expressed at a very low level in green apple varieties, where the expression of these genes was remarkably higher in red apple cultivar [62]. Our results showed a higher expression of PAL, 4CL, DFR, F3H and UFGT genes in the post-ripening samples. A previous qRT-PCR study showed an up-regulation of these genes along with higher anthocyanin levels and ripening progress [49]. The expression of UFGT was investigated in different fruits like grape, strawberry and lychee and was found to be critical in controlling fruit color development [63]. Also, the involvement of transcription factors regulating the structure genes expression also demonstrated and suggested a key role of MYB and the basic-helix-loop-helix (bHLH) families’ genes on the regulation of anthocyanin biosynthesis related genes [64]. Another study also reported that MYB transcription factors might play a key role in regulating anthocyanin biosynthesis in some fruit [65,66]. In our study, several MYB and bHLH, including other genes, were differentially expressed in the post-ripening samples. Thus, these transcription factors may regulate the expression of key structure genes associated with the anthocyanin biosynthesis pathway.
However, future studies on critical observations of the changes of these regulatory factors are crucial. Furthermore, the association and crosstalk between the differentially expressed structure genes and corresponding transcription factors in hormonal regulation and anthocyanin biosynthesis during ripening are still indistinct and this suggests further investigation.

4. Materials and Methods

4.1. Plant Material

‘Hongyu’ apples (Malus domestica L. Borkh) were used in this study. Fruits of uniform size and color at two different developmental stages of 110 days after flowering (DAF) and 120 DAF were collected from an orchard in Changping District, Beijing, China. Another group of fruits harvested at 120 DAF (considered as commercial maturity) was allowed to ripen for 5 d at 20 °C. Fruit at 120 DAF was considered to be commercially ripened according to the defined harvesting date by the particular apple orchard. To investigate the sudden hormonal changes and anthocyanin biosynthesis during ripening, a sampling date was selected before commercial ripening (110 DAF) and immediately after harvesting by allowing the fruits to be stored for 5 d at 20 °C. The environmental conditions and cultural practices were maintained according to standard apple orchard practices. Fruits at each sampling point were transported to the laboratory as soon as possible, the flesh was taken and immediately frozen in liquid nitrogen and stored at −80 °C for further use. Each sample consisted of nine randomly selected fruits, separated into three replicates with three fruits in each. The fruits of 110 DAF, 120 DAF and stored at 20 °C for 5 d were considered for RNA-seq analysis and assigned to a pre-ripening and a post-ripening group.

4.2. Organic Acids and Sugars Measurement

Organic acid and sugar content was measured by following a modified method previously described by Lin et al. [67]. Briefly, a total of 1 g sample powder was homogenized with 4 mL of 80% ethanol. The mixture was extracted with ultra-sonication for 30 min at room temperature and centrifuged at 10,000 g for 15 min. Aliquots of 1 mL of the upper phase were dried with pure nitrogen. The residue was dissolved in 1 mL ddH2O and filtered through a membrane with 0.22 μm pore size. For sugar analysis, a volume of 10 μL for each sample was injected into the ion chromatograph (ICS-3000, Dionex, Sunnyvale, CA, USA) fitted with Carbo PacTMPA20 column (3 mm × 150 mm). The column temperature was 35 °C and the flow rate was 0.5 mL min−1. The gradient elution buffer was used as follows: A, ddH2O; B, 250 mmol L−1 NaOH; equal gradients of 92.5% A and 7.5% B were used for elution. A pulsed amperometric detector with gold electrode was used. For organic acid analysis, a volume of 25 μL for each sample was injected into the ion chromatograph (ICS-3000, Dionex, Sunnyvale, CA, USA) fitted with IonPac AS11-HC column (4 mm × 250 mm). The column temperature was 30 °C and the flow rate was 1 mL min−1. The gradient elution procedure was 0.8 mmol L−1 KOH, 0–12 min; 0.8–34 mmol L−1 KOH, 12–40 min; 34 mmol L−1 KOH, 40–50 min. The electrical conductivity detector was used for detection. Substances were identified according to the retention time (RT) of standard compounds. The contents of glucose, fructose, sucrose, malate and oxalate were calculated by comparison with standard curves.

4.3. Total RNA Extraction

Total RNA was extracted using a modified CTAB method [68]. RNA quantity and quality (purity) were analyzed using a UV-Vis spectrophotometer (Merinton Instrument, Inc., Ann Arbor, MI, USA).

4.4. Library Construction

Total RNA from 2 different fruit developmental stages (pre-ripening and ripening) and post-harvest stored apples at 20 °C (post-ripening) were pooled before library preparation and an equimolar quantity of total RNA from samples at each stage were combined into one pool. Prior to cDNA library construction, poly-T oligo-attached magnetic beads were used to purify the mRNA, which was added to a fragmentation buffer to break mRNA into short segments of approximately 200 bp followed by cDNA synthesis using random hexamers and reverse transcriptase. Then buffer, dNTPs, RNase H and DNA polymerase I were added to synthesize the second strand cDNA. Both ends of sequencing adapters were ligated after terminal repairing of double-stranded cDNA fragments. The AMPure XP system (Beckman Coulter, Beverly, MA, USA) was used to purify the final cDNA library and selectively enriched by PCR enrichment. The library preparations were sequenced by Novogene Bioinformatics Technology Co., Ltd. (Beijing, China) on an Illumina HiSeqTM PE150 platform. Three replications were used in each group for the precision of results. Fruit of 110 DAF denoted as HY_3rd, 120 DAF indicated as HY_4th and ripened for 5 d at 20 °C were designated as HY_20 °C_5 groups in transcriptome analysis.

4.5. Functional Annotation and Analysis of Differential Gene Expression

Raw reads were cleaned by removing low-quality reads and adapter sequence reads and clean reads were mapped to the apple reference genome (https://www.rosaceae.org/organism/Malus/x-domestica) using TopHat software [69]. The analysis of differentially expressed genes (DEGs) was done using a method previously reported in Reference [70] and p-value thresholds were measured by the false discovery rate (FDR) determination via multiple testing [71]. For gene expression quantification, the read numbers mapped to each gene were counted using HTSeq v0.6.1 and then normalized to FPKM (Fragments per Kilobase of transcript sequence per Millions of base pairs sequenced) [70]. For analysis of differential expression of the unigenes between two samples, the DESeq package (v1.10.1) was used and a corrected P-value threshold of 0.005 (using the Benjamini & Hochberg method) and log2 (fold change) of 1 were considered for filtering the DEGs value. For Gene ontology (GO) enrichment analysis of differentially expressed genes, GOseq (v2.12) was used [72]. KEGG enrichment analysis was performed using the KOBAS (v 2.0) software (Center for Bioinformatics, Peking University, Beijing, China) [73].

4.6. Quantitative Real-Time PCR (qRT-PCR)

First strand cDNA was synthesized from 1 μg RNA, after digestion of genomic DNA, using the iScriptTM cDNA Synthesis Kit (Bio-Rad, Hercules, California, USA). qRT-PCR was performed using the Power SYBR® Green PCR Master Mix kit (Applied Biosystems, Foster City, CA, USA) on an ABI 7500 instrument (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA), initiated by 5 min at 95 °C then followed by 45 cycles of 95 °C for 10 s, 60 °C for 10 s and 72 °C for 15 s and completed with a melting curve analysis program. No-template controls and melting curve analyses were included in every reaction. The actin gene was included as an internal control. The comparative CT method (2−ΔΔCT method) was used to analyze the expression levels of the different genes [74]. Primers used in this analysis are listed in Table S7.

4.7. Statistical Analysis

Figures were drawn using Origin 8.6 (Microcal Software Inc., Northampton, MA, USA). The MultiExperiment Viewer software (MeV v4.8.1, J. Craig Venter Institute, La Jolla, CA, USA) was used to represent the relative values of gene expressions.

5. Conclusions

The present study investigated the transcriptome profiles of an early ripening apple cultivar prior to, and immediately after, the commercial ripening stage using DGE deep-sequencing technologies through Illumina RNA-seq. This transcriptome analysis revealed that, in addition to ethylene biosynthesis genes, key genes of ABA biosynthesis (NCED1) and several auxin response factors were differentially expressed in the post-ripening samples. This result suggested that investigating the interplay of different hormones and respective genes that are differentially expressed in post-harvest fruit at the commercial maturity stage will be helpful for discovering the possible hormonal crosstalk in ripening regulation and for extending the ripening period of early ripening apples. This study also focused on apple fruits upon ripening, which also suggests that the up-regulation of DEGs involved in anthocyanin biosynthesis plays an important role in the accumulation of anthocyanin, similar to the results of previous studies on many fruits. The genes encoding F3H and UFGT were significantly upregulated in the post-harvest fruit, indicating key roles in anthocyanin biosynthesis. Several transcription factors were also specified including MADS-box genes, MYB, WRKY, NAC, bHLH and HSF, which may regulate the ripening process by the regulation of ethylene and other hormones and anthocyanin biosynthesis.

Supplementary Materials

The following are available online, Table S1: List of annotated DEGS, Table S2: List of GO enrichment analysis results, Table S3: List of KEGG pathway enrichment analysis results, Table S4: List of DEGs related to plant hormone signaling pathway between pre and post-ripening samples, Table S5: List of DEGs involved in anthocyanin biosynthesis pathway, Table S6 List of DEGs related to different transcription factors, Table S7: List of Primers used for gene expression analysis.

Author Contributions

Conceptualization, Q.L., X.H. and Z.W.; Methodology, J.C.O.; Software, J.C.O.; Validation, J.C.O., X.H., Q.L. and Z.W.; Formal Analysis, Q.L.; Investigation, Q.L., J.C.O.; Resources, Z.W.; Data Curation, J.C.O.; Writing-Original Draft Preparation, J.C.O., Q.L.; Writing-Review & Editing, J.C.O., Q.L.; Visualization, J.C.O., Q.L.; Supervision, Q.L.; Project Administration, Q.L.; Funding Acquisition, Z.W.

Funding

This research was funded by National Key Research and Development Program of China (2016YFD0400100), the Young Scientist’s Fund of National Natural Science Foundation of China (31601527) and the Agricultural Science and Technology Innovation Program (ASTIP) from the Chinese Central Government.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript and in the decision to publish the results.

References

  1. Guo, H.; Ecker, J.R. The Ethylene Signaling Pathway: New Insights. Curr. Opin. Plant Biol. 2004, 7, 40–49. [Google Scholar] [CrossRef] [PubMed]
  2. Ireland, H.S.; Yao, J.L.; Tomes, S.; Sutherland, P.W.; Nieuwenhuizen, N.; Gunaseelan, K.; Winz, R.A.; David, K.M.; Schaffer, R.J. Apple SEPALLATA1/2-like Genes Control Fruit Flesh Development and Ripening. Plant J. 2013, 73, 1044–1056. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, X.; Song, J.; Du, L.; Forney, C.; Campbell-Palmer, L.; Fillmore, S.; Wismer, P.; Zhang, Z. Ethylene and 1-MCP Regulate Major Volatile Biosynthetic Pathways in Apple Fruit. Food Chem. 2016, 194, 325–336. [Google Scholar] [CrossRef] [PubMed]
  4. Whale, S.K.; Singh, Z. Endogenous Ethylene and Color Development in the Skin of Pink Lady Apple. J. Am. Soc. Hortic. Sci. 2007, 132, 20–28. [Google Scholar]
  5. Seymour, G.B.; Østergaard, L.; Chapman, N.H.; Knapp, S.; Martin, C. Fruit Ripening and Development. Annu. Rev. Plant Biol. 2013, 64, 219–241. [Google Scholar] [CrossRef] [PubMed]
  6. McAtee, P.; Karim, S.; Schaffer, R.; David, K. A Dynamic Interplay between Phytohormones Is Required for Fruit Development, Maturation, and Ripening. Front. Plant Sci. 2013, 4, 79. [Google Scholar] [CrossRef] [PubMed]
  7. Giovannoni, J.J. Fruit Ripening Mutants Yield Insights into Ripening Control. Curr. Opin. Plant Biol. 2007, 10, 283–289. [Google Scholar] [CrossRef] [PubMed]
  8. Bapat, V.A.; Trivedi, P.K.; Ghosh, A.; Sane, V.A.; Ganapathi, T.R.; Nath, P. Ripening of Fleshy Fruit: Molecular Insight and the Role of Ethylene. Biotechnol. Adv. 2010, 28, 94–107. [Google Scholar] [CrossRef] [PubMed]
  9. Osorio, S.; Scossa, F.; Fernie, A.R. Molecular Regulation of Fruit Ripening. Front. Plant Sci. 2013, 4, 198. [Google Scholar] [CrossRef] [PubMed]
  10. Martel, C.; Vrebalov, J.; Tafelmeyer, P.; Giovannoni, J.J. The Tomato MADS-Box Transcription Factor Ripening Inhibitor Interacts with Promoters Involved in Numerous Ripening Processes in a Colorless Nonripening-Dependent Manner. Plant Physiol. 2011, 157, 1568–1579. [Google Scholar] [CrossRef] [PubMed]
  11. Kumar, R.; Khurana, A.; Sharma, A.K. Role of Plant Hormones and Their Interplay in Development and Ripening of Fleshy Fruits. J. Exp. Bot. 2014, 65, 4561–4575. [Google Scholar] [CrossRef] [PubMed]
  12. Srivastava, A.; Handa, A.K. Hormonal Regulation of Tomato Fruit Development: A Molecular Perspective. J. Plant Growth Regul. 2005, 24, 67–82. [Google Scholar] [CrossRef]
  13. Böttcher, C.; Keyzers, R.A.; Boss, P.K.; Davies, C. Sequestration of Auxin by the Indole-3-Acetic Acid-Amido Synthetase GH3-1 in Grape Berry (Vitis vinifera L.) and the Proposed Role of Auxin Conjugation during Ripening. J. Exp. Bot. 2010, 61, 3615–3625. [Google Scholar] [CrossRef] [PubMed]
  14. Jones, B.; Frasse, P.; Olmos, E.; Zegzouti, H.; Li, Z.G.; Latché, A.; Pech, J.C.; Bouzayen, M. Down-Regulation of DR12, an Auxin-Response-Factor Homolog, in the Tomato Results in a Pleiotropic Phenotype Including Dark Green and Blotchy Ripening Fruit. Plant J. 2002, 32, 603–613. [Google Scholar] [CrossRef] [PubMed]
  15. Trainotti, L.; Tadiello, A.; Casadoro, G. The Involvement of Auxin in the Ripening of Climacteric Fruits Comes of Age: The Hormone Plays a Role of Its Own and Has an Intense Interplay with Ethylene in Ripening Peaches. J. Exp. Bot. 2007, 58, 3299–3308. [Google Scholar] [CrossRef] [PubMed]
  16. Kumar, R.; Agarwal, P.; Tyagi, A.K.; Sharma, A.K. Genome-Wide Investigation and Expression Analysis Suggest Diverse Roles of Auxin-Responsive GH3 Genes during Development and Response to Different Stimuli in Tomato (Solanum lycopersicum). Mol. Genet. Genom. 2012, 287, 221–235. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, K.; Kang, B.C.; Jiang, H.; Moore, S.L.; Li, H.; Watkins, C.B.; Setter, T.L.; Jahn, M.M. A GH3-like Gene, CcGH3, Isolated from Capsicum Chinense L. Fruit Is Regulated by Auxin and Ethylene. Plant Mol. Biol. 2005, 58, 447–464. [Google Scholar] [CrossRef] [PubMed]
  18. Martinez, G.A.; Chaves, A.R.; Anon, M.C. Effect of Exogenous Application of Gibberellic Acid on Color Change and Phenylalanine Ammonia-Lyase, Chlorophyllase, and Peroxidase Activities during Ripening of Strawberry Fruit (Fragaria x Ananassa Duch). J. Plant Growth Regul. 1996, 15, 139–146. [Google Scholar] [CrossRef]
  19. De La Fuente, J.I.; Amaya, I.; Castillejo, C.; Sánchez-Sevilla, J.F.; Quesada, M.A.; Botella, M.A.; Valpuesta, V. The Strawberry Gene FaGAST Affects Plant Growth through Inhibition of Cell Elongation. J. Exp. Bot. 2006, 57, 2401–2411. [Google Scholar] [CrossRef] [PubMed]
  20. Hirayama, T.; Shinozaki, K. Perception and Transduction of Abscisic Acid Signals: Keys to the Function of the Versatile Plant Hormone ABA. Trends Plant Sci. 2007, 12, 343–351. [Google Scholar] [CrossRef] [PubMed]
  21. Pan, Q.H.; Li, M.J.; Peng, C.C.; Zhang, N.; Zou, X.; Zou, K.Q.; Wang, X.L.; Yu, X.C.; Wang, X.F.; Zhang, D.P. Abscisic Acid Activates Acid Invertases in Developing Grape Berry. Physiol. Plant. 2005, 125, 157–170. [Google Scholar] [CrossRef]
  22. Sun, L.; Sun, Y.; Zhang, M.; Wang, L.; Ren, J.; Cui, M.; Wang, Y.; Ji, K.; Li, P.; Li, Q.; et al. Suppression of 9-Cis-Epoxycarotenoid Dioxygenase, Which Encodes a Key Enzyme in Abscisic Acid Biosynthesis, Alters Fruit Texture in Transgenic Tomato. Plant Physiol. 2012, 158, 283–298. [Google Scholar] [CrossRef] [PubMed]
  23. Kelebek, H.; Selli, S. Evaluation of Chemical Constituents and Antioxidant Activity of Sweet Cherry (Prunus avium L.) Cultivars. Int. J. Food Sci. Technol. 2011, 46, 2530–2537. [Google Scholar] [CrossRef]
  24. Jaakola, L. New Insights into the Regulation of Anthocyanin Biosynthesis in Fruits. Trends Plant Sci. 2013, 18, 477–483. [Google Scholar] [CrossRef] [PubMed]
  25. Montefiori, M.; Espley, R.V.; Stevenson, D.; Cooney, J.; Datson, P.M.; Saiz, A.; Atkinson, R.G.; Hellens, R.P.; Allan, A.C. Identification and Characterisation of F3GT1 and F3GGT1, Two Glycosyltransferases Responsible for Anthocyanin Biosynthesis in Red-Fleshed Kiwifruit (Actinidia chinensis). Plant J. 2011, 65, 106–118. [Google Scholar] [CrossRef] [PubMed]
  26. Goto-Yamamoto, N.; Wan, G.H.; Masaki, K.; Kobayashi, S. Structure and Transcription of Three Chalcone Synthase Genes of Grapevine (Vitis vinifera). Plant Sci. 2002, 162, 867–872. [Google Scholar] [CrossRef]
  27. Tsuda, T.; Yamaguchi, M.; Honda, C.; Moriguchi, T. Expression of Anthocyanin Biosynthesis Genes in the Skin of Peach and Nectarine Fruit. J. Am. Soc. Hortic. Sci. 2004, 129, 857–862. [Google Scholar]
  28. Kobayashi, S.; Ishimaru, M.; Ding, C.K.; Yakushiji, H.; Goto, N. Comparison of UDP-Glucose:Flavonoid 3-O-Glucosyltransferase (UFGT) Gene Sequences between White Grapes (Vitis vinifera) and Their Sports with Red Skin. Plant Sci. 2001, 160, 543–550. [Google Scholar] [CrossRef]
  29. Yamazaki, M.; Shibata, M.; Nishiyama, Y.; Springob, K.; Kitayama, M.; Shimada, N.; Aoki, T.; Ayabe, S.I.; Saito, K. Differential Gene Expression Profiles of Red and Green Forms of Perilla Frutescens Leading to Comprehensive Identification of Anthocyanin Biosynthetic Genes. FEBS J. 2008, 275, 3494–3502. [Google Scholar] [CrossRef] [PubMed]
  30. Carbone, F.; Preuss, A.; De Vos, R.C.H.; D’Amico, E.; Perrotta, G.; Bovy, A.G.; Martens, S.; Rosati, C. Developmental, Genetic and Environmental Factors Affect the Expression of Flavonoid Genes, Enzymes and Metabolites in Strawberry Fruits. Plant Cell Environ. 2009, 32, 1117–1131. [Google Scholar] [CrossRef] [PubMed]
  31. Hichri, I.; Heppel, S.C.; Pillet, J.; Léon, C.; Czemmel, S.; Delrot, S.; Lauvergeat, V.; Bogs, J. The Basic Helix-Loop-Helix Transcription Factor MYC1 Is Involved in the Regulation of the Flavonoid Biosynthesis Pathway in Grapevine. Mol. Plant 2010, 3, 509–523. [Google Scholar] [CrossRef] [PubMed]
  32. Giovannoni, J.J.J. Genetic Regulation of Fruit Development and Ripening. Plant Cell Online 2004, 16, S170–S180. [Google Scholar] [CrossRef] [PubMed]
  33. Saladie, M.; Matas, A.J.; Isaacson, T.; Jenks, M.A.; Goodwin, S.M.; Niklas, K.J.; Xiaolin, R.; Labavitch, J.M.; Shackel, K.A.; Fernie, A.R.; et al. A Reevaluation of the Key Factors That Influence Tomato Fruit Softening and Integrity. Plant Physiol. 2007, 144, 1012–1028. [Google Scholar] [CrossRef] [PubMed]
  34. Alba, R.; Payton, P.; Fei, Z.; McQuinn, R.; Debbie, P.; Martin, G.B.; Tanksley, S.D.; Giovannoni, J.J. Transcriptome and Selected Metabolite Analyses Reveal Multiple Points of Ethylene Control during Tomato Fruit Development. Plant Cell 2005, 17, 2954–2965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Vriezen, W.H.; Feron, R.; Maretto, F.; Keijman, J.; Mariani, C. Changes in Tomato Ovary Transcriptome Demonstrate Complex Hormonal Regulation of Fruit Set. New Phytol. 2008, 177, 60–76. [Google Scholar] [CrossRef] [PubMed]
  36. Matas, A.J.; Yeats, T.H.; Buda, G.J.; Zheng, Y.; Chatterjee, S.; Tohge, T.; Ponnala, L.; Adato, A.; Aharoni, A.; Stark, R.; et al. Tissue- and Cell-Type Specific Transcriptome Profiling of Expanding Tomato Fruit Provides Insights into Metabolic and Regulatory Specialization and Cuticle Formation. Plant Cell 2011, 23, 3893–3910. [Google Scholar] [CrossRef] [PubMed]
  37. Osorio, S.; Alba, R.; Damasceno, C.M.B.; Lopez-Casado, G.; Lohse, M.; Zanor, M.I.; Tohge, T.; Usadel, B.; Rose, J.K.C.; Fei, Z.; et al. Systems Biology of Tomato Fruit Development: Combined Transcript, Protein, and Metabolite Analysis of Tomato Transcription Factor (nor, Rin) and Ethylene Receptor (Nr) Mutants Reveals Novel Regulatory Interactions. Plant Physiol. 2011, 157, 405–425. [Google Scholar] [CrossRef] [PubMed]
  38. Rohrmann, J.; Tohge, T.; Alba, R.; Osorio, S.; Caldana, C.; McQuinn, R.; Arvidsson, S.; Van Der Merwe, M.J.; Riaño-Pachón, D.M.; Mueller-Roeber, B.; et al. Combined Transcription Factor Profiling, Microarray Analysis and Metabolite Profiling Reveals the Transcriptional Control of Metabolic Shifts Occurring during Tomato Fruit Development. Plant J. 2011, 68, 999–1013. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, M.; Yuan, B.; Leng, P. The Role of ABA in Triggering Ethylene Biosynthesis and Ripening of Tomato Fruit. J. Exp. Bot. 2009, 60, 1579–1588. [Google Scholar] [CrossRef] [PubMed]
  40. Gao, C.; Ju, Z.; Li, S.; Zuo, J.; Fu, D.; Tian, H.; Luo, Y.; Zhu, B. Deciphering Ascorbic Acid Regulatory Pathways in Ripening Tomato Fruit Using a Weighted Gene Correlation Network Analysis Approach. J. Integr. Plant Biol. 2013, 55, 1080–1091. [Google Scholar] [CrossRef] [PubMed]
  41. Ingle, M.; D’Souza, M.C.; Townsend, E.C. Fruit Characteristics of York Apples during Development and after Storage. HortScience 2000, 35, 95–98. [Google Scholar]
  42. Zude-Sasse, M.; Herold, B.; Geyer, M.; Huyskens-Keil, S. Influence of Maturity Stage on Physical Properties in Apple. Acta Hortic. 2001, 553, 109–110. [Google Scholar] [CrossRef]
  43. Kvikliene, N. Influence of harvest date on physiological and biochemical processes in apple fruit. Hortic. Veg. Grow. 2004, 23, 412–420. [Google Scholar]
  44. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for Linking Genomes to Life and the Environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, A.; Tan, D.; Tatsuki, M.; Kasai, A.; Li, T.; Saito, H.; Harada, T. Molecular Mechanism of Distinct Ripening Profiles in “Fuji” Apple Fruit and Its Early Maturing Sports. Postharvest Biol. Technol. 2009, 52, 38–43. [Google Scholar] [CrossRef]
  46. Jiang, Y.; Joyce, D.C.; Macnish, A.J. Effect of Abscisic Acid on Banana Fruit Ripening in Relation to the Role of Ethylene. J. Plant Growth Regul. 2000, 19, 106–111. [Google Scholar] [CrossRef] [PubMed]
  47. Feng, S.; Wang, Y.; Yang, S.; Xu, Y.; Chen, X. Anthocyanin Biosynthesis in Pears Is Regulated by a R2R3-MYB Transcription Factor PyMYB10. Planta 2010, 232, 245–255. [Google Scholar] [CrossRef] [PubMed]
  48. Soubeyrand, E.; Basteau, C.; Hilbert, G.; Van Leeuwen, C.; Delrot, S.; Gomès, E. Nitrogen Supply Affects Anthocyanin Biosynthetic and Regulatory Genes in Grapevine Cv. Cabernet-Sauvignon Berries. Phytochemistry 2014, 103, 38–49. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, Y.; Che, F.; Wang, L.; Meng, R.; Zhang, X.; Zhao, Z. Fruit Coloration and Anthocyanin Biosynthesis after Bag Removal in Non-Red and Red Apples (Malus × Domestica Borkh). Molecules 2013, 18, 1549–1563. [Google Scholar] [CrossRef] [PubMed]
  50. Shan, W.; Kuang, J.F.; Chen, L.; Xie, H.; Peng, H.H.; Xiao, Y.Y.; Li, X.P.; Chen, W.X.; He, Q.G.; Chen, J.Y.; et al. Molecular Characterization of Banana NAC Transcription Factors and Their Interactions with Ethylene Signalling Component EIL during Fruit Ripening. J. Exp. Bot. 2012, 63, 5171–5187. [Google Scholar] [CrossRef] [PubMed]
  51. Schaffer, R.J.; Ireland, H.S.; Ross, J.J.; Ling, T.J.; David, K.M. SEPALLATA1/2-suppressed mature apples have low ethylene, high auxin and reduced transcription of ripening-related genes. AoB Plants 2013, 5, pls047. [Google Scholar] [CrossRef] [PubMed]
  52. Giovannoni, J. Molecular Biology of Fruit Maturation and Ripening. Annu. Rev. Plant Physiol. Plant Mol. Biol. 2001, 52, 725–749. [Google Scholar] [CrossRef] [PubMed]
  53. Xiao, Y.Y.; Chen, J.Y.; Kuang, J.F.; Shan, W.; Xie, H.; Jiang, Y.M.; Lu, W.J. Banana Ethylene Response Factors Are Involved in Fruit Ripening through Their Interactions with Ethylene Biosynthesis Genes. J. Exp. Bot. 2013, 64, 2499–2510. [Google Scholar] [CrossRef] [PubMed]
  54. Janssen, B.J.; Thodey, K.; Schaffer, R.J.; Alba, R.; Balakrishnan, L.; Bishop, R.; Bowen, J.H.; Crowhurst, R.N.; Gleave, A.P.; Ledger, S.; et al. Global Gene Expression Analysis of Apple Fruit Development from the Floral Bud to Ripe Fruit. BMC Plant Biol. 2008, 8, 16. [Google Scholar] [CrossRef] [PubMed]
  55. Jia, H.-F.; Chai, Y.-M.; Li, C.-L.; Lu, D.; Luo, J.-J.; Qin, L.; Shen, Y.-Y. Abscisic Acid Plays an Important Role in the Regulation of Strawberry Fruit Ripening. Plant Physiol. 2011, 157, 188–199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Chai, Y.M.; Zhang, Q.; Tian, L.; Li, C.L.; Xing, Y.; Qin, L.; Shen, Y.Y. Brassinosteroid Is Involved in Strawberry Fruit Ripening. Plant Growth Regul. 2013, 69, 63–69. [Google Scholar] [CrossRef]
  57. Huang, G.; Li, T.; Li, X.; Tan, D.; Jiang, Z.; Wei, Y.; Li, J.; Wang, A. Comparative Transcriptome Analysis of Climacteric Fruit of Chinese Pear (Pyrus ussuriensis) Reveals New Insights into Fruit Ripening. PLoS ONE 2014, 9, e107562. [Google Scholar] [CrossRef] [PubMed]
  58. Ren, J.; Chen, P.; Dai, S.; Li, P.; Li, Q.; Ji, K.; Wang, Y.; Leng, P. Role of Abscisic Acid and Ethylene in Sweet Cherry Fruit Maturation: Molecular Aspects. N. Z. J. Crop Hortic. Sci. 2011, 39, 161–174. [Google Scholar] [CrossRef]
  59. Symons, G.M.; Chua, Y.J.; Ross, J.J.; Quittenden, L.J.; Davies, N.W.; Reid, J.B. Hormonal Changes during Non-Climacteric Ripening in Strawberry. J. Exp. Bot. 2012, 63, 4741–4750. [Google Scholar] [CrossRef] [PubMed]
  60. Kang, C.; Darwish, O.; Geretz, A.; Shahan, R.; Alkharouf, N.; Liu, Z. Genome-Scale Transcriptomic Insights into Early-Stage Fruit Development in Woodland Strawberry Fragaria Vesca. Plant Cell 2013, 25, 1960–1978. [Google Scholar] [CrossRef] [PubMed]
  61. Sagar, M.; Chervin, C.; Roustant, J.P.; Bouzayen, M.; Zouine, M. Under-Expression of the Auxin Response Factor SL-ARF4 Improves Post-Harvest Behavior of Tomato Fruits. Plant Signal. Behav. 2013, 8, e25647. [Google Scholar] [CrossRef] [Green Version]
  62. Takos, A.M.; Jaffe, F.W.; Jacob, S.R.; Bogs, J.; Robinson, S.P.; Walker, A.R. Light-Induced Expression of a MYB Gene Regulates Anthocyanin Biosynthesis in Red Apples. Plant Physiol. 2006, 142, 1216–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Zhao, Z.C.; Hu, G.B.; Hu, F.C.; Wang, H.C.; Yang, Z.Y.; Lai, B. The UDP Glucose: Flavonoid-3-O-Glucosyltransferase (UFGT) Gene Regulates Anthocyanin Biosynthesis in Litchi (Litchi chinesis Sonn.) during Fruit Coloration. Mol. Biol. Rep. 2012, 39, 6409–6415. [Google Scholar] [CrossRef] [PubMed]
  64. Wada, T.; Kunihiro, A.; Tominaga-Wada, R. Arabidopsis CAPRICE (MYB) and GLABRA3 (BHLH) Control Tomato (Solanum lycopersicum) Anthocyanin Biosynthesis. PLoS ONE 2014, 9, e109093. [Google Scholar] [CrossRef] [PubMed]
  65. Telias, A.; Lin-Wang, K.; Stevenson, D.E.; Cooney, J.M.; Hellens, R.P.; Allan, A.C.; Hoover, E.E.; Bradeen, J.M. Apple Skin Patterning Is Associated with Differential Expression of MYb10. BMC Plant Biol. 2011, 11, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ravaglia, D.; Espley, R.V.; Henry-Kirk, R.A.; Andreotti, C.; Ziosi, V.; Hellens, R.P.; Costa, G.; Allan, A.C. Transcriptional Regulation of Flavonoid Biosynthesis in Nectarine (Prunus persica) by a Set of R2R3 MYB Transcription Factors. BMC Plant Biol. 2013, 13, 68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Lin, Q.; Xie, Y.; Liu, W.; Zhang, J.; Cheng, S.; Xie, X.; Guan, W.; Wang, Z. UV-C Treatment on Physiological Response of Potato (Solanum tuberosum L.) during Low Temperature Storage. J. Food Sci. Technol. 2017, 54, 55–61. [Google Scholar] [CrossRef] [PubMed]
  68. Zong, X.; Wang, W.; Wang, J.; Wei, H.; Yan, X.; Liu, Q. The application of SYBR Green I real-time quantitative RT-PCR in quantitative analysis of sweet cherry viruses in different tissues. Acta Phytophylacica Sin. 2012, 39, 497–502. [Google Scholar]
  69. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHat: Discovering Splice Junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  70. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed]
  71. Benjamini, Y.; Yekutieli, D. The Control of the False Discovery Rate in Multiple Testing under Dependency by Yoav Benjamini 1 and Daniel Yekutieli 2. Ann. Stat. 2001, 29, 1165–1188. [Google Scholar]
  72. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene Ontology Analysis for RNA-Seq: Accounting for Selection Bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [PubMed]
  73. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated Genome Annotation and Pathway Identification Using the KEGG Orthology (KO) as a Controlled Vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef] [PubMed]
  74. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔ CT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are available from the authors.
Molecules 23 01908 g001 550
Figure 1. Variation of malate, oxalate, glucose, sucrose and fructose content in pre-ripening (HY_3rd), ripening (HY_4th) and post-ripening (HY_20 °C_5) ‘Hongyu’ Apple fruit. The error bars represent the standard errors. The significant differences were calculated at the 0.05 level.
Figure 1. Variation of malate, oxalate, glucose, sucrose and fructose content in pre-ripening (HY_3rd), ripening (HY_4th) and post-ripening (HY_20 °C_5) ‘Hongyu’ Apple fruit. The error bars represent the standard errors. The significant differences were calculated at the 0.05 level.
Molecules 23 01908 g001
Molecules 23 01908 g002 550
Figure 2. The number of differentially expressed genes (DEGs) within pre-(HY_4th vs. HY_3rd) and post-ripening (HY_20 °C_5 vs. HY_4th) groups (A). Overlapped DEGs were shown by Venn diagram between different groups (B).
Figure 2. The number of differentially expressed genes (DEGs) within pre-(HY_4th vs. HY_3rd) and post-ripening (HY_20 °C_5 vs. HY_4th) groups (A). Overlapped DEGs were shown by Venn diagram between different groups (B).
Molecules 23 01908 g002
Molecules 23 01908 g003 550
Figure 3. Classification of gene ontology (GO) terms concerning to differentially expressed genes in pre- (A) and post-ripening (B) apple fruit. The GO terms were classified into three categories, including molecular function, cellular component and biological process. Top 31 enriched GO terms were exhibited in each cluster.
Figure 3. Classification of gene ontology (GO) terms concerning to differentially expressed genes in pre- (A) and post-ripening (B) apple fruit. The GO terms were classified into three categories, including molecular function, cellular component and biological process. Top 31 enriched GO terms were exhibited in each cluster.
Molecules 23 01908 g003
Molecules 23 01908 g004 550
Figure 4. KEGG pathway enrichment concerning differently expressed genes in pre- (A) and post-ripening (B) apple. Top 20 enriched KEGG pathways were exhibited in each cluster. Rich factor is the ratio of DEGs counts to this pathway in the annotated genes counts. The more the q-value is close to zero, the more significant is the enrichment.
Figure 4. KEGG pathway enrichment concerning differently expressed genes in pre- (A) and post-ripening (B) apple. Top 20 enriched KEGG pathways were exhibited in each cluster. Rich factor is the ratio of DEGs counts to this pathway in the annotated genes counts. The more the q-value is close to zero, the more significant is the enrichment.
Molecules 23 01908 g004
Molecules 23 01908 g005 550
Figure 5. The correlation analysis of gene expression value from quantitative real-time PCR (Q-PCR) and RNA-Seq.
Figure 5. The correlation analysis of gene expression value from quantitative real-time PCR (Q-PCR) and RNA-Seq.
Molecules 23 01908 g005
Molecules 23 01908 g006 550
Figure 6. Heat map represents DEGs involved in ethylene signaling pathway in pre (HY_4th vs. HY_3rd) and post (HY_20 °C_5 vs. HY_4th) ripening groups. Genes remarked with light green color denotes expression in both groups and red color denotes expressed in the post-ripening group only. ACS (1-Aminocyclopropane-1-carboxylic acid synthase); ACO (1-Aminocyclopropane-1-carboxylic acid oxidase); ERS (Ethylene response sensor); CTR1 (CONSTITUTIVE TRIPLE RESPONSE 1); EIN2 (ETHYLENE INSENSITIVE 2); EIN3 (ETHYLENE INSENSITIVE 3); EIL3 (ETHYLENE INSENSITIVE 3-like 3); and ERFs (Ethylene Response Factors).
Figure 6. Heat map represents DEGs involved in ethylene signaling pathway in pre (HY_4th vs. HY_3rd) and post (HY_20 °C_5 vs. HY_4th) ripening groups. Genes remarked with light green color denotes expression in both groups and red color denotes expressed in the post-ripening group only. ACS (1-Aminocyclopropane-1-carboxylic acid synthase); ACO (1-Aminocyclopropane-1-carboxylic acid oxidase); ERS (Ethylene response sensor); CTR1 (CONSTITUTIVE TRIPLE RESPONSE 1); EIN2 (ETHYLENE INSENSITIVE 2); EIN3 (ETHYLENE INSENSITIVE 3); EIL3 (ETHYLENE INSENSITIVE 3-like 3); and ERFs (Ethylene Response Factors).
Molecules 23 01908 g006
Molecules 23 01908 g007 550
Figure 7. Heat map represents DEGs involved in anthocyanin biosynthesis pathway in pre (HY_4th vs. HY_3rd) and post (HY_20 °C_5 vs. HY_4th) ripening groups. Genes remarked with light green color denotes expression in both groups and red color denotes expressed in the post-ripening group only. PAL (Phenylalanine ammonia-lyase); 4CL (4-coumarate-CoA ligase); CHS (Chalcone Synthase); CHI (Chalcone isomerase); F3H (Flavanone 3-hydroxylase); F3′H (Flavanone 3′-hydroxylase); DFR (Dihydroflavonol 4-reductase); ANS (Anthocyanidin synthase) and UFGT (UDP glucose: flavonol 3-O-glucosyltransferase).
Figure 7. Heat map represents DEGs involved in anthocyanin biosynthesis pathway in pre (HY_4th vs. HY_3rd) and post (HY_20 °C_5 vs. HY_4th) ripening groups. Genes remarked with light green color denotes expression in both groups and red color denotes expressed in the post-ripening group only. PAL (Phenylalanine ammonia-lyase); 4CL (4-coumarate-CoA ligase); CHS (Chalcone Synthase); CHI (Chalcone isomerase); F3H (Flavanone 3-hydroxylase); F3′H (Flavanone 3′-hydroxylase); DFR (Dihydroflavonol 4-reductase); ANS (Anthocyanidin synthase) and UFGT (UDP glucose: flavonol 3-O-glucosyltransferase).
Molecules 23 01908 g007
Table
Table 1. Statistics of deferential gene expression (DGE) library sequencing and read mapping.
Table 1. Statistics of deferential gene expression (DGE) library sequencing and read mapping.
GroupsLibraryClean BasesError Rate (%)Q20 (%)GC Content (%)Total Mapped Reads (%)Uniquely Mapped Reads (%)
HY_3rdHY0018.03G0.0395.1647.1990.5875.78
HY0026.04G0.0297.1147.2693.6478.28
HY0037.15G0.0296.3147.3191.9276.86
HY_4thHY0047.97G0.0394.9847.4790.5475.68
HY0059.56G0.0297.0047.5892.3376.47
HY0067.16G0.0395.0547.3790.9976.13
HY_20 °C_5HY0077.42G0.0295.6147.4089.8875.53
HY0086.85G0.0295.5847.4889.8175.41
HY0097.67G0.0295.6447.5690.1775.63
Table
Table 2. DEGs related to plant hormones expressed in pre-ripening (HY_4th vs. HY_3rd) group.
Table 2. DEGs related to plant hormones expressed in pre-ripening (HY_4th vs. HY_3rd) group.
DEGsLog2Fold ChangeP (adj)RegulationAnnotation
Abscisic acid
MDP00004370331.84450.027924upARATH Protein phosphatase 2C 56 Arabidopsis thaliana
MDP0000929213−0.759360.017195down9-cis-epoxycarotenoid dioxygenase NCED1 Phaseolus vulgaris
MDP00001658670.608820.013085upSerine/threonine-protein kinase EDR1 Arabidopsis thaliana
MDP0000226571−0.70090.027422downSerine/threonine-protein kinase HT1 Arabidopsis thaliana
Auxin
MDP0000132805−0.792027.43 × 10−5downAuxin-responsive protein IAA13 Arabidopsis thaliana
MDP00002683060.542850.042346upAuxin response factor 6 Arabidopsis thaliana
MDP00001374610.695950.030746upAuxin response factor 8 Arabidopsis thaliana
MDP00006665393.29117.27 × 10−10upIndole-3-acetic acid-amido synthetase GH3.3 Arabidopsis thaliana
Gibberellin
MDP00003013681.10410.010819upScarecrow-like protein 3 Arabidopsis thaliana
MDP00009299940.715058.31 × 10−5upGibberellin receptor GID1B Arabidopsis thaliana
Brassinosteroid
MDP00002618511.42180.016674upG-type lectin S-receptor-like serine/threonine-protein kinase At4g27290 Arabidopsis thaliana
MDP0000361876−2.94272.40 × 10−5downXyloglucan endotransglucosylase/hydrolase protein 9 Arabidopsis thaliana
MDP00001200441.03823.17 × 10−5upCytochrome P450 714A1 Arabidopsis thaliana
MDP0000130467−0.633880.017284downCyclin-U4-1 Arabidopsis thaliana
Table
Table 3. DEGs related to plant hormones expressed in post-ripening (HY_20 °C_5 vs. HY_4th) group.
Table 3. DEGs related to plant hormones expressed in post-ripening (HY_20 °C_5 vs. HY_4th) group.
DEGsLog2Fold ChangeP (adj)RegulationAnnotation
Abscisic acid
MDP00009292131.97443.68 × 10−46up9-cis-epoxycarotenoid dioxygenase NCED1, Phaseolus vulgaris
MDP00001746070.846710.000358upProtein phosphatase 2C 57 Arabidopsis thaliana
MDP0000165966−0.67860.0035436downSerine/threonine-protein kinase HT1 Arabidopsis thaliana
MDP0000232165−0.72739.29 × 10−5downSerine/threonine-protein kinase ATM Arabidopsis thaliana
Auxin
MDP00002094322.28783.71 × 10−41upProbable indole-3-acetic acid-amido synthetase GH3.1 Arabidopsis thaliana
MDP0000153538−1.13540.0022471downAuxin response factor 6 Arabidopsis thaliana
MDP00002460421.87495.27 × 10−7upAuxin-responsive protein IAA31 Arabidopsis thaliana
MDP0000876321−2.06332.73 × 10−9downAuxin response factor 19 Arabidopsis thaliana
Gibberellin
MDP00001377056.59885.24 × 10−23upGibberellin 2-beta-dioxygenase 1 Pisum sativum
MDP00003195220.73800.000418upGibberellin receptor GID1B Arabidopsis thaliana
MDP0000901967−5.06423.65 × 10−14downGibberellin-regulated protein 14 Arabidopsis thaliana
MDP00001190931.45371.64 × 10−21upScarecrow-like protein 9 Arabidopsis thaliana
MDP0000256486−1.37160.0014764downScarecrow-like protein 4 Arabidopsis thaliana
Brassinosteroid
MDP0000210409−2.63981.74 × 10−7downG-type lectin S-receptor-like serine/threonine-protein kinase B120 Arabidopsis thaliana
MDP0000140341−2.16612.10 × 10−12downLRR receptor-like serine/threonine-protein kinase RCH1 Arabidopsis thaliana
MDP0000183140−1.38880.015187downXyloglucan glycosyltransferase 4 Arabidopsis thaliana
MDP0000127773−1.05927.65 × 10−6downCytochrome P450 90A1 Arabidopsis thaliana
MDP00001408030.69140.036655upCytochrome P450 78A6 Arabidopsis thaliana
MDP0000154776−2.39844.95 × 10−6downCytochrome P450 734A1 Arabidopsis thaliana
MDP0000176105−1.63880.011614downCyclin-D4-1 Arabidopsis thaliana

Share and Cite

MDPI and ACS Style

Onik, J.C.; Hu, X.; Lin, Q.; Wang, Z. Comparative Transcriptomic Profiling to Understand Pre- and Post-Ripening Hormonal Regulations and Anthocyanin Biosynthesis in Early Ripening Apple Fruit. Molecules 2018, 23, 1908. https://doi.org/10.3390/molecules23081908

AMA Style

Onik JC, Hu X, Lin Q, Wang Z. Comparative Transcriptomic Profiling to Understand Pre- and Post-Ripening Hormonal Regulations and Anthocyanin Biosynthesis in Early Ripening Apple Fruit. Molecules. 2018; 23(8):1908. https://doi.org/10.3390/molecules23081908

Chicago/Turabian Style

Onik, Jakaria Chowdhury, Xiaojia Hu, Qiong Lin, and Zhidong Wang. 2018. "Comparative Transcriptomic Profiling to Understand Pre- and Post-Ripening Hormonal Regulations and Anthocyanin Biosynthesis in Early Ripening Apple Fruit" Molecules 23, no. 8: 1908. https://doi.org/10.3390/molecules23081908

Article Metrics

Back to TopTop