Identification and characterization of cherry (Cerasus pseudocerasus G. Don) genes responding to parthenocarpy induced by GA3 through transcriptome analysis
BMC Genetics volume 20, Article number: 65 (2019)
Fruit set after successful pollination is key for the production of sweet cherries, and a low fruit-setting rate is the main problem in production of this crop. As gibberellin treatment can directly induce parthenogenesis and satisfy the hormone requirement during fruit growth and development, such treatment is an important strategy for improving the fruit-setting rate of sweet cherries. Previous studies have mainly focused on physiological aspects, such as fruit quality, fruit size, and anatomical structure, whereas the molecular mechanism remains clear.
In this study, we analyzed the transcriptome of ‘Meizao’ sweet cherry fruit treated with gibberellin during the anthesis and hard-core periods to identify genes associated with parthenocarpic fruit set. A total of 25,341 genes were identified at the anthesis and hard-core stages, 765 (681 upregulated, 84 downregulated) and 186 (141 upregulated, 45 downregulated) of which were significant differentially expressed genes (DEGs) at the anthesis and the hard-core stages after gibberellin 3 (GA3) treatment, respectively. Based on DEGs between the control and GA3 treatments, the GA3 response mainly involves parthenocarpic fruit set and cell division. Exogenous gibberellin stimulated sweet cherry fruit parthenocarpy and enlargement, as verified by qRT-PCR results of related genes as well as the parthenocarpic fruit set and fruit size. Based on our research and previous studies in Arabidopsis thaliana, we identified key genes associated with parthenocarpic fruit set and cell division. Interestingly, we observed patterns among sweet cherry fruit setting-related DEGs, especially those associated with hormone balance, cytoskeleton formation and cell wall modification.
Overall, the result provides a possible molecular mechanism regulating parthenocarpic fruit set that will be important for basic research and industrial development of sweet cherries.
Cherry genes respond to parthenocarpy and promote fruit setting induced by GA3
Fruit set is an important step in fruit growth and development. In this process, the ovary becomes enlarged with development of the embryo after fertilization, inducing fruit formation . A variety of hormonal synergies play a major role in controlling fruit set; GA, auxin, and cytokinin alone cause the fruit to grow and develop to a certain stage, but the fruit develops normally due to their synergistic effect [2, 3].
Exogenous spraying of 5 mg / L GA3 maintains the activity of citrus parietal cells and promote their division, thereby increasing the rate of fruit setting , and it was reported that treatment of grape flower spikes with 30 mg/L GA3 before full bloom significantly increased the fruit-setting rate . Moreover, silencing of SIGH3 and upregulation of SIARF7, a negative regulator of fruit set, in transgenic tomato plants decreased the IAA content; the number of middle and endocortical cells was also reduced, but their volume increased [6, 7]. In addition to auxin and gibberellin, the fruit-setting rate is affected by cytokinin, which influences fruit growth by promoting cell division . Transcriptome analysis of normally pollinated and gibberellin-treated cucumbers showed that gibberellin promotes fruit set through interaction with auxin and cytokinin . Therefore, knowledge of the transcriptome can help in an understanding of the molecular mechanism by which gibberellin regulates fruit set.
Parthenocarpy is an agronomic trait in which fruit develops without pollination and fertilization under natural conditions or with human intervention , and exogenous gibberellin can induce parthenocarpy during crop production. For example, application of 1500 mg / L GA3 during anthesis led to high-quality parthenogenetic fruit in Annona lucidum , and 500 mg / L GA4 + 7 induced parthenocarpy in sand pear . Gibberellins are also active in inducing parthenogenesis and fruit development in other crops such as cucumber , lemon , tomato , and banana . Overall, gibberellin affects plant growth and development, and there are two main genes that function in the regulation of GA signaling. The protein encoded by the DELLA gene represses gibberellin signaling; binding of gibberellins to the gibberellin receptor GID1, stimulates ubiquitination of the DELLA protein, resulting in 26S proteasome-mediated degradation. In SIDELLA-silenced transgenic tomato, gibberellin synthesis pathway-related gene expression is altered, and parthenocarpy is induced via promotion of ovarian cell division and growth [16,17,18,19,20].
Sweet cherry is a delicious and juicy fruit, the production of which has become a major industry, providing farmers with increasing income. However, the low fruit-setting rate caused by self-fruiting or its low self-setting rate has become a bottleneck restricting the development of this industry. As exogenous gibberellin can induce parthenocarpy in fruit, thereby increasing the fruit-setting rate, we sought to elucidate the molecular mechanism by which gibberellin induces parthenocarpic fruit set. To this end, we used next-generation sequencing technology to carry out transcriptome sequencing (RNA-seq) at two key developmental periods after GA3 treatment of ‘Meizao’ European sweet cherry grown under greenhouse conditions. This is the first application of RNA-seq to analyze developmental metabolic pathways and gene expression patterns in sweet cherry fruit after GA treatment. Genes that promote fruit setting, parthenogenesis, cell division and growth were among those analyzed. The aim of this study is to provide evidence for the molecular mechanism of fruit enlargement and fruit setting and to examine the effect of plant hormone regulation on the growth and development of fruit in parthenogenetic European sweet cherry treated with GA3. The differentially expressed genes identified through transcriptome sequencing provide a molecular foundation explaining the induction of parthenocarpic fruit set by exogenous gibberellin in sweet cherry.
Growth and development characteristics of fruit after exogenous GA3 treatment
Full growth and development of fruit was observed after treatment of ‘Meizao’ sweet cherry trees in the anthesis period (7 DAT) and the hard-core period (27 DAT) with 300 mg/L GA3. The size of GA3-treated sweet cherry is shown in Table 1. At 7 DAT, the vertical and transverse diameters of sweet cherry treated with GA3 were significantly higher than those of the control. Among them, the increase in the vertical diameter was higher than the increase in the transverse diameter. In the anthesis stage (7 DAT), the diameter of GA3-treated fruits increased by 123.9% compared with the control group, and the vertical diameter increased by 214.09% (Fig. 1). In the first expansion period (17 DAT), the growth rate of the transverse diameter was not significantly different from that of the control, but the vertical diameter reached 34.16%. After the hard-core stage (27 DAT), all the fruits of the control group had dropped.
Fruit setting rate and parthenogenetic rate of fruit after exogenous GA3 treatment
After a week of spraying GA3 at the anthesis stage (7 DAT), the fruit set rate was 77.33% and the control 12.55%. During the hard-core period (27 DAT), the untreated sweet cherries had completely dropped, whereas the GA3-treated sweet cherries maintained a fruit setting rate of 57.23%. As shown in Fig. 2h, spraying GA3 in the blossom drop period can both induce parthenogenesis and significantly increase the fruit-setting rate. After GA3 treatment, the rate at which fruit dropped from the first expansion to the second enlargement was 34.3%, and the parthenogenetic rate was 100% (Fig. 2i). With the growth of the fruit, the development of the seeds was inhibited, and embryo abortion occurred.
Analysis of transcriptional sequence data
Illumina platform sequencing revealed 137,658,266 (total of CK11, CK12, CK13) and 134,689,144 (total of CK21, CK22, CK23) and 140,680,232 (total of T11, T12, T13) and 138,579,108 (total of T21, T22, T23) raw reads for the control and GA-treatment groups, respectively (Table 2). After removing linkers, low-quality reads and reads with an N ratio greater than 5%, 136,440,870, 133,557,948, 139,421,700 and 137,336,514 clean reads were selected for further analysis. The proportion of clean reads to total reads obtained from the four libraries was 99.11% (average of CK11, CK12, CK13), 99.16% (average of CK21, CK22, CK23), 99.11% (average of T11, T12, T13) and 99.10% (average of T21, T22, T23). Among these reads, 83, 83.67, 83, and 79.67% were mapped to the European sweet cherry genome (https://www.rosaceae.org/species/prunus_avium/genome_v1.0.a1), including 81.65, 80.87, 80.86, and 81.22% mapped to exonic regions, 8.62, 9.45, 9.29, and 9.12% mapped to intronic regions, and 9.74, 9.69, 9.85, and 9.66% mapped to intergenic regions (Table 2).
Identification and classification of DEGs at different stages
In the two key developmental stages of European sweet cherry, changes in gene expression were determined by comparing C1 versus T1 and C2 versus T2. In these two periods, 765 and 186 significant differentially expressed genes were obtained. In the first period, 685 genes were upregulated, and 84 genes were downregulated. In the second period, 141 upregulated and 45 downregulated genes were identified by RNA-seq (Fig. 3a). The number of upregulated genes was significantly higher than that of downregulated genes in the two periods of European sweet cherry development. However, Venn diagram analysis indicated that only 32 genes showed a pattern of upregulated expression in both key periods (Fig. 3b), and specific upregulated genes were detected at every time point.
Changes in expression patterns of genes related to gibberellin metabolism and signal transduction
After treatment with gibberellin, the fruit rate, fruit size and rate of parthenocarpy changed significantly. Gibberellin 2-beta-dioxygenase (GA2ox) is a key enzyme of GA catabolism , and the level of GA2ox expression changed significantly during anthesis; it was also detected during the hard-core period. Pav_sc0000095.1_g1110.1.mk was significantly upregulated at 7 DAT (Additional file 2); the annotation is PavGA2ox, which was used to generate GA51, GA29, GA34 and GA8 in GA9, GA20, GA4 and GA1, respectively . In Arabidopsis thaliana, the proteins Scarecrow-like (SCL) and DELLA are antagonists that promote gibberellin signaling . Transcriptome sequencing revealed certain genes at the anthesis stage after GA treatment, such as pavscl1, pavscl3 and GA2ox. Encoding a GA receptor, GID1B (Pav_sc0000848.1_g080.1.mk), was downregulated after GA3 treatment in the hard-nucleus phase, whereas GA2ox and DELLA were upregulated. These results showed that exogenous GA3 is associated with parthenogenesis.
In GA signal transduction, active GA regulates the hydrolysis and inactivation of DELLA by stimulating the GID1-DELLA-GID2 complex, thus regulating downstream gene expression . In this study, we found that the gene encoding GID1 was downregulated after GA treatment and that expression of the DELLA gene was upregulated; however, the gene encoding GID2 was not detected. Pav_sc0000848.1_g080.1.mk was further annotated as PavGID1B and only noted in the hard-core period. During the hard-core stage, the DELLA protein Pav_sc0000464.1_g350.1.mk significantly increased. The above results indicate that after treatment of ‘Meizao’ sweet cherry with exogenous GA3, no significant expression of GID1-DELLA-GID2-related genes was detected at the full anthesis stage. We speculate that GID1-DELLA-GID2-related genes respond very quickly and are correlated with seed development.
Changes in expression patterns of cytoskeleton and cell wall modification genes
After GA3 treatment, significant changes in the size and shape of ‘Meizao’ cherries were noted. After 7 days, the treated fruits were significantly swollen (Fig. 2). Upregulation of plant hormone synthesis, degradation, control and signal transduction genes and downregulation of genes involved in the cytoskeleton, cell wall modifications and other related biological processes synergistically caused irreversible fruit expansion (Additional file 3). Expression of genes encoding actin and tubulin changed significantly after GA3 treatment. At 7 DAT, only one gene encoding actin, EDH1 (Pav_sc0001084.1_g180.1.mk), was downregulated. Conversely, genes encoding tubulin were not detected during the hard-core period. In addition, different genes involved in cell division, including a cyclin-dependent kinase and five histones, were downregulated at the anthesis stage after GA treatment.
Expression of cell wall modification genes was significantly altered after GA treatment. Expansion proteins catalyze the extension of cell walls . In our study, 26 cell wall-associated DEGs were detected at the anthesis stage after GA treatment, and 5 cell wall-associated DEGs were detected during the hard-core period, including those encoding expansin, xyloglucan glycosyltransferase and glucanase; most of these genes (27/31) were upregulated. Clearly, changes in expression of genes associated with cell division and cell wall metabolism occur during fruit development.
DEG analysis of transcription factors
Analysis of all DEGs at the two stages uncovered 403 transcription factors, which were further divided into 46 and 30 gene families (Fig. 4a). At anthesis, the transcription factor families with the most DEGs were the bHLH and B3 families (7.94%), followed by the NAC family (5.96%), the WRKY family (4.96%), the MYB-related family (4.71%), the C3H, C2H2, ERF and bZIP families (4.47%), the FAR1 family (4.22%), the MYB family (3.97%), the HSF family (3.23%), the TALE family (2.98%), the ARF family (2.73%), and the M-type, HB-other, and TCP families (2.23%). DEGs were strongly enriched in these 17 families of transcription factors. In the hard-core stage, the B3 family of transcription factors had the most DEGs (12.9%), followed by the NAC, bHLH and HD-ZIP families (7.52%), the C3H family (6.45%), the C2H2 family (5.38%), and the MYB-related, GRAS and MYB families (4.3%) (Fig. 4b). The abovementioned families of transcription factors contained a higher than average number of DEGs. Based on the expression patterns of the DEGs identified, most DEG transcription factors exhibited a trend of upregulation compared with the control at the two key stages after gibberellin treatment. Although some of the DEGs belonging to the Nin-like, HSF, ARF, B3, HB-other, MYB, Bzip, GeBP, C3H and NF-YC families were downregulated after the first treatment, the NAC and Trihelix families were downregulated after the second treatment. As previous studies have shown, many transcription factors play an important role in seed development in a wide variety of plants . For example, the B3 superfamily is key in seed development [27,28,29]. Moreover, all NAC genes are reportedly involved in the synthesis of cell secondary metabolites , and transcription factors (TFs) such as GRAS and HB participate in ABA and GA signaling  (Figs. 5, 6). In our study, transcription factors related to seed development were differentially expressed between seeded and parthenogenetic progeny, indicating that they are associated with parthenogenetic phenotypes.
Analysis of DEG pathway enrichment to reveal possible pathways influencing parthenogenesis
To further understand differences in the metabolic pathways of ovule development between seeded and parthenogenetic sweet cherry, we analyzed metabolic pathways showing differential gene enrichment in the two stages (Fig. 7), revealing 27 metabolic pathways enriched during the anthesis and hard-core periods. These 27 enriched metabolic pathways are involved in a variety of biological processes, including the synthesis and metabolism of secondary metabolites, plant hormone signal transduction, flavonoid biosynthesis, starch and sucrose metabolism, Cutin, and suberine and wax biosynthesis. Although no pathways were significantly enriched at the second developmental stage, two metabolic pathways, phenylalanine metabolism and phenylpropanoid biosynthesis, were significantly enriched at the first stage of cherry development, and most DEGs in these two pathways were upregulated. Collectively, these results suggest that parthenogenetic fruit set may be related to phenylalanine metabolism and phenylpropanoid biosynthesis pathways.
Functional analysis of GOs of differential genes indicates factors that may lead to parthenogenetic fruit set
To explore the key biological processes and functional classification related to cherry parthenogenetic fruit set, DEGs were further analyzed for biological function GOs in the two critical stages of cherry development (Fig. 8). The biological functional classification of GOs involves the following three categories: biological process, molecular function and cellular component. Most of the DEGs were upregulated, and those detected at the first stage were highly represented in cellular process, metabolic process, single-organism process, response to stimulus, biological regulation, and cellular component organization of biogenesis groups in the biological process category. Most of the differences in gene-rich GOs were related to plant growth and development, response to external stimuli and signal transduction, glucose and lipid metabolism, and epigenetic regulation; the abovementioned groups and parthenogenetic fruit set are closely related. The status of differential gene enrichment was generally the same in the anthesis stage after GA treatment.
Based on the classification of GO-based molecular functions, DEGs were enriched in catalytic, binding, transporter and nucleic acid-binding transcription factor categories. These enriched GOs are closely related to the regulation of biological processes and anabolism of functional substances. In each subcategory, more genes were upregulated than downregulated in the first stage of sweet cherry development.
Further analyses of the GO classification of cellular components were performed, showing most DEGs to be enriched in the following groups: intrinsic component of membrane, integral component of membrane, membrane and membrane part. Similar to the GO enrichment of molecular function and biological process, in the GO enrichment of the cell group, most DEGs were upregulated in the anthesis stage after GA treatment compared with the control.
To further explore differences in the GO functions enriched among DEGs during the two critical stages of sweet cherry development, we focused on GO enrichment in the biological process category. During the two critical stages of sweet cherry development, 709/24 upregulated and 55/0 downregulated genes were found to be significantly enriched compared to the control. The GO items enriched in these two developmental stages are summarized below. In the anthesis stage, GO items with upregulated DEGs are mainly associated with cell wall development, with further division into the groups of cell wall synthesis and cytokinesis, carbohydrate and polysaccharide biosynthesis, secondary metabolism, and lignin and phenylpropanoid anabolism. During the anthesis stage, downregulation of GO items related to cell wall development was associated with DEGs significantly enriched in the cell cycle process, and GO items with upregulated DEGs are mainly associated with photosynthesis and sucrose biosynthesis. As indicated above, no downregulated genes were enriched in the second stage.
Transcriptional dynamics analysis of DEGs
We conducted further analysis of genes associated with fruit set, cell division and parthenocarpy in the anthesis and hard-core stages of sweet cherry using k-means clustering to classify genes with relatively consistent expression. The different genes (15,167 in total) were subjected to this analysis of temporal and spatial expression patterns, with all differentially expressed genes with similar expression patterns during the process of cherry development grouped into one group, as presented in Fig. 9. Six coexpressed clusters in the CK and GA3 treatments were identified during sweet cherry development, of which 4 were consistent in expression pattern. For parthenogenetic sweet cherry, DEGs were consistently upregulated in CK at the early stage of development in cluster 1, cluster 2, cluster 4, and cluster 6, and the only difference was the degree of upregulation. These four clusters displayed lower expression in the anthesis period after GA treatment. Cluster 5 exhibited different expression patterns, and the transcriptional level of the DEGs showed continuous downregulation. After GA3 treatment, the DEG patterns of cluster 5 and clusters 1, 2, 4, and 6 were reversed and correlated over time. The level of transcription was increased or decreased during anthesis in sweet cherry, and the differential genes between the five clusters exhibited completely opposite spatiotemporal expression patterns in the CK and GA3 treatments.
Based on clustering results, 5 of the transcription factors were selected for analysis, and it was found that the greatest numbers of transcription factors in cluster 5 belong to the AP2 and MYB families, followed by the MADS family. In cluster 5, the different genes are mainly enriched in GO items of seed and embryo development, cell differentiation, cell wall formation, defense response, hormone regulation, and transcriptional regulation, among others.
RNA-seq validation of DEGs by quantitative real-time PCR (qRT-PCR)
To validate the reliability of the transcriptome data, we applied qRT-PCR to detect expression of 30 DEGs related to fruit set, cell division and parthenocarpy in samples from the treatment and control groups (Fig. 10) and further carried out correlation analysis of gene expression fold changes. As shown in Fig. 10d, the results were highly consistent and thus verified the reliability of the transcriptome sequencing data. These 30 genes are mainly involved in the following processes: fruit setting (Fig. 10a), including PavYUCCA, PavYUCCA10, PavACS, PavACO and PavCYP707A3; cell division (Fig. 10b), including PavH4, PavCDC2D, PavH3.2, and PavCDC45; parthenocarpy (Fig. 10c), including PavSAUR32, PavSCL1, PavCYCPA3, and PavDELLA. The qRT-PCR results were consistent with the RNA-seq data, which strongly increased the reliability of our data.
Transcriptome analysis in response to GA3 was performed during the growth and development of ‘Meizao’ cherry during two critical periods of development. In total, 765 and 186 differentially expressed genes were identified, which were then compared to determine their biological functions based on GO and metabolic pathway enrichment. In this study, GA3 treatment was carried out at two key stages of sweet cherry development: anthesis and the hard-core stage. A one-to-one comparison of parthenocarpic fruit development after GA3 treatment was conducted, and differences were compared. The findings might result in a relatively comprehensive follow-up analysis to identify important genes linked to sweet cherry fruit production and that participate in functional pathways. Furthermore, it will be helpful to explore the potential mechanisms regulating and controlling the production of sweet cherry fruit.
The effect of gibberellin on gene expression in the development of parthenogenetic sweet cherry
In agriculture, the production of parthenocarpic seedless grapes often involves treating flower buds with gibberellin before flowering . After GA3 application, the response regulates expression of DELLA and GID1 via feedback. GA application decreases expression of auxin-responsive proteins, leading to parthenocarpy. Moreover, other factors appear to have a regulatory role in parthenocarpy; for example, GRAS was upregulated after GA3 application (Fig. 11). Therefore, analysis of hormone-related gene expression during the development of sweet cherry in two stages can help to further determine the influence of gene expression differences on sweet cherry parthenocarpy and elucidate the mechanism of parthenocarpy. Previous studies have shown that SAUR can promote cellular expansion in Arabidopsis [33, 34], and SAUR was upregulated in the first stage in the IAA signal transduction pathway, which is consistent with the results shown in Fig. 2. Expression of DEGs induced by GA3 is shown in Additional file 2. Differentially expressed genes associated with auxin synthesis and the auxin inhibitor IAA also showed patterns of upregulated expression at the first stage. Accordingly, these auxin-related differential gene expression patterns may be a manifestation of auxin balance feedback adjustment. The GH3 gene maintains the dynamic balance of auxin in plants by encoding IAA-amido synthase, which converts active IAA into a bound form [7, 35]. In the hard-core stage of sweet cherry, expression of the GH3 gene and the auxin inhibitor AUX/IAA was downregulated, further explaining the auxin feedback regulation. Different hormones can play cooperative roles through signaling, and one hormone that acts synergistically with auxin is BR. Indeed, studies have shown that auxin promotes plant growth by synthesizing BR .
In the BR response pathway, upregulated expression of the receptor gene BRI1 indicates an association with the GA-induced BR response. In this study, genes related to BR synthesis displayed upregulated expression in the two critical stages of sweet cherry development, which is consistent with the synergistic effect of auxin and brassinolide. The study of Fu showed that external application of BR can stimulate cell division in cucumber and lead to parthenocarpy . It is speculated that expression of genes related to BR synthesis and signal transduction in sweet cherry during seed development after gibberellin treatment is upregulated relative to that in untreated sweet cherry. GA may also promote the parthenogenesis of sweet cherry through similar mechanisms.
It has been reported that a high level of ethylene can cause parthenogenesis in tomato . The results of our study showed that in the first period, expression of genes related to ethylene synthesis and signal transduction was upregulated after GA3 treatment. Abscisic acid also participates in the development of the seed coat in Arabidopsis and regulates cell wall metabolism by reducing cell wall PG, PE and XET [39, 40]. In the first critical period of development of sweet cherry, expression of a gene related to seed growth was detected by analyzing seed development-related gene patterns, consistent with the upregulated expression of the abscisic acid receptor and genes related to transmembrane transport at the two critical periods. However, no genes associated with seed coat development, which may only be expressed a short time after treatment, were detected in the hard-core period. In summary, genes associated with GA, auxin and BR were significantly upregulated after GA3 treatment, which may be associated with parthenogenesis. Regardless, the specific mechanism needs further study.
Gibberellin promotes cell wall relaxation and induces fruit enlargement
The cell wall is composed of cellulose microfilaments embedded in the hemicellulose stromal layer through noncovalent bonding . In the process of fruit enlargement, expansion proteins promote noncovalent dissociation, and glucanase and xylanase loosen the hemicellulose matrix . Expression of DEGs associated with cell division after GA3 treatment is presented in Additional file 3. According to our results, glucanase gene ncRNA was significantly upregulated in the first stage after GA3 treatment, and the endoglucanase gene ENDO17 was upregulated in the second stage. The two expansion genes EXPA2 and EXPB1 were significantly upregulated in the anthesis period; however, only one gene, EXPA10, was detected in the second period, and its expression was downregulated. Genes encoding glucanase, xylanase and expansion were significantly upregulated in the first period after GA3 treatment, and one expansion gene was downregulated in the second period, suggesting that these genes play an important role in the process of sweet cherry enlargement. Previous studies in grape have found that the XET gene was significantly upregulated during the first rapid growth stage and then decreased over time [42, 43], which is not consistent with the results. The main reason may be that different species exhibit different expression patterns. The growth and developmental period of sweet cherry is short. However, after GA3 treatment, the rapid growth period of sweet cherry was relatively long, and glucanase, xylanase and expansion genes were highly expressed. It is known cell wall modification plays a key role in exogenous GA3-induced fruit enlargement. In addition to genes associated with cell wall relaxation, cytoskeleton modification genes also play an important role in the process of GA3-induced fruit enlargement. According to the analysis of protein groups in the second and third period of grape development, expression of 3 actin and 10 tubulin genes was upregulated.
Gibberellin-induced fruit setting in sweet cherry
Pollination is a key step in the development of flowers to fruits. Egg cells develop into seeds after fertilization by pollen, and the developing zygotes stimulate the ovary to form a fruit. We identified three genes that may play an important role in sweet cherry fruit set. First, the effects of ethylene include promotion of fruit ripening, and ethephon is used as a thinning agent for many fruit trees [44,45,46]. AP2/ERF is one of the largest transcription factor families in plants , the members of which are mainly involved in the regulation of such processes as plant growth [48, 49], flower development , fruit development  and seed development . In the present study, six AP2/ERF transcription factors showed altered expression in the flowering stage, and four of them are ethylene-responsive transcription factors. We speculate that these genes may play an important role in the transition from flowers to fruit, thus promoting fruit set. Altered expression of three AP2/ERF transcription factors was found in the hard-core stage, possibly regulating ovule senescence via changes in the expression patterns of different genes. In addition, the ethylene synthesis pathway enzymes 1-aminocyclopropane-1-carboxylate oxidase (ACO) and 1-aminocyclopropane-1-carboxylate synthase (ACS) were upregulated after GA treatment. These enzymes may act by altering GA signaling to promote the occurrence of parthenocarpy and fruit set. Second, the YUCCA gene plays an important role in the biosynthesis of auxin and plant development . Expression of DEGs associated with fruit set after GA3 treatment are listed in Additional file 4. The YUCCA10 gene encoding indole-3-pyruvate monooxygenase was significantly upregulated after GA treatment. The YUCCA gene plays an important role in vegetative growth and reproductive development of plants and regulates the formation of the floral meristem and the development of lateral organs . Some studies have reported that the content of auxin in the ovary increases after GA treatment [55, 56]. Overall, GA may play a role in auxin accumulation in the ovary and in turn act as a signal to affect fruit set. Third, the gene coding for ABA 8′-hydroxylase was upregulated after GA treatment in the second stage, and this gene is thought to be the main reason for the decrease in ABA . A high ABA content causes fruit to drop. After GA treatment, expression of this gene was upregulated, which was beneficial to the fruit. These results are consistent with the results shown in Fig. 2, indicating that this period is critical for fruit set. These results indicate that auxin, ethylene and abscisic acid may regulate fruit set in sweet cherry after GA treatment.
To the best of our knowledge, this study is the first to provide a dynamic analysis of differences in differentially expressed genes related to parthenogenesis, fruit enlargement and fruit set in sweet cherry fruit due to GA3 treatment. We selected two critical periods during the development of sweet cherry, and genes involved in many metabolic processes changed significantly after a week of GA treatment, including hormone-related genes, transcription factors and cell division genes. Further analysis of these genes will provide a new approach for studying fruit set in sweet cherry. Additionally, our research results provide a relatively complete molecular platform for future research on mechanisms of parthenocarpy, fruit enlargement and fruit set in sweet cherry.
Trees of the ‘Meizao’ variety of sweet cherry (Cerasus pseudocerasus G. Don Meizao) were treated in the anthesis and hard-core stages at the commercial greenhouse in Jueyu (Tai’an, China) in 2017. All plant materials used in our study are regular species, which do not require any permissions. Experimental materials were obtained from three 10-year-old trees growing under similar conditions. We selected four major branches with similar growth status and direction per tree and divided them into two groups: one group was used for fruit set statistics and the other group for sampling. Every group included two branches, one for GA3 (300 mg/L of GA3 containing 0.1% (v/v) Tween 20) treatment and another as a control (growing under the natural state). After treatment with GA3, 30 cm*25 cm transparent paper bags were applied to cover flowers to avoid interference from pollen and to allow for normal management. Blossoming to maturity was divided into five periods: anthesis, the first enlargement period, the hard-core period, the second enlargement period and the mature period, respectively. Samples were collected from the two important stages for parthenogenetic fruit (anthesis and hard-core stages); ten fruits were collected from each branch at the anthesis (7 days after treatment, 7 DAT) and hard-core (27 DAT) stages and used for transcriptome analysis, with three branches of three trees used as three biological replicates. A workflow chart of the analytical process is shown in Additional file 5.
The two key developmental periods of sweet cherry, anthesis and hard-core were examined. In each stage, RNA was extracted from sweet cherries treated or not with GA3. Total RNA was extracted using the RNAprep Pure TIANGEN kit (TIANGEN, DP441, Beijing). The quality and quantity of the RNA were determined using a Kaiao K5500 spectrophotometer (Kaiao, Beijing) followed by 1% agarose gel electrophoresis (18S and 28S bands). The integrity and concentration of the RNA samples were detected using Agilent 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, CA, USA), and high-quality RNA was sent to Beijing Genomics for further quality assessment.
Library preparation for transcriptome analysis
After evaluating the total RNA samples, mRNA was enriched using magnetic beads containing oligo (dT) primers. Fragmentation buffer was added to the obtained mRNA to produce short fragments, which were then used as templates. First-strand cDNA was synthesized using six-base random primers; buffer, dNTPs, RNaseH and DNA polymerase I were added to synthesize the second strand. The cDNA obtained was purified using a QIAquick PCR kit and eluted with EB buffer. The purified double-stranded cDNA was subjected to end repair, base A addition and sequencing, and fragment sizes were then determined by agarose gel electrophoresis and PCR amplification. The constructed library was sequenced using the Illumina platform.
Illumina read processing and de novo assembly
Illumina high-throughput sequencing results were originally obtained as raw image data files. After CASSAVABase base recognition (base calling), the data were transformed into the original sequence (sequenced reads), and the results were stored in FASTQ (fq) file format. In this format, the ASCII code value corresponding to each base quality character minus 33 (Sanger quality system) is the sequencing quality score of the base. Different scores represent different nucleotide sequencing error rates, such as score values of 20 and 30, indicating base sequence error rates of 1 and 0.1%, respectively. To ensure the quality of the data for analysis, we filtered out sequencing linkers, low-quality sequences and reads with an N ratio greater than 5%. After filtering, the total number of bases with a mass value greater than 30 (error rate less than 0.1%) was divided by the total number of bases to obtain the Q30 value; the higher this number, the higher is the sequencing quality. A follow-up analysis based on the clean reads was performed.
Functional annotation and classification
High-quality clean reads were obtained by removing sequences of sequencing joints, low-quality sequences (Q is less than 19 bases), and reads with an N ratio of more than 5%. TopHat software  with Bowtie 2  were employed for alignment of the sequences to the European sweet cherry genome located at https://www.rosaceae.org/species/prunus_avium/genome_v1.0.a1.
No more than two base mismatches were allowed, and annotation results for genes were obtained by comparing them with the European sweet cherry genome data. In addition, NCBI, UniProt, GO and KEGG were used for annotations. For GO analysis, the NCBI NR database was first used for BLAST alignment, and Blast2 GO was utilized to obtain GO entries corresponding to each gene. Enrichment analysis was performed for each pathway in KEGG to determine the significance of differentially expressed genes (enrichment of pathway).
The clean reads from the two developmental stages were obtained through transcriptome sequencing as described above, and the complete dataset excluded sequences with more than two base mismatches . The amount of mRNA produced was calculated using FPKM  to estimate gene expression levels, which were compared between the two developmental stages using DESeq. Genes with |log2Ratio| ≥ 1 and q < 0.05 were selected as significant differentially expressed genes . These differentially expressed genes were subjected to GO analysis and KEGG metabolic pathway enrichment analysis based on the calculation of the number of genes for each item and hypergeometric test results after p-value correction with a threshold of q < 0.05. The GO items satisfying this condition were defined as significantly enriched among the differentially expressed genes. A pathway defined by a Q-value ≤0.05 (corrected p-value) is a significantly enriched KEGG pathway for differentially expressed genes, and the definition of a GO item meeting the condition for significant enrichment among differently expressed genes is a Q-value ≤0.05 (pathway corrected p-value) based on the difference in gene expression enrichment in a KEGG pathway.
Quantitative real-time PCR (qRT-PCR) validation
To verify the results of the transcriptome analysis, we selected genes involved in fruit set, cell division and parthenocarpy. Total RNA was reverse transcribed (1 μg per sample) in a 20-μl reaction using a cDNA reverse transcription synthesis kit (RR047A, TaKaRa, China); the primers used are listed in Additional file 1. The SYBR® Premix Ex Taq TM (Tli RNaseH Plus) kit (Thermo Fisher) was employed for the fluorescence quantitative PCR reaction. The reaction system included the following: 12.5 μL SYBR® Premix Ex Taq (2 ×), 1 μL each forward and reverse primers, 1 μL cDNA, and ddH2O to 25 μL. The experimental design included 3 technical replicates. The fluorescence quantitative PCR reaction conditions were as follows: 35 ~ 40 cycles of pre-denaturation at 95 °C for 30 s, denaturation at 95 °C for 5 s, and annealing at 60 °C for 30 s. After the reaction, a fluorescence curve and a melting curve were obtained, and the comparative Ct (2-ΔΔCt) method was used for data analysis .
Availability of data and materials
The dataset supporting the conclusions of this article is available in the GenBank data libraries with the accession number PRJNA529517 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA529517).
Days after treatment
Differentially genes expressed
Kyoto Encyclopedia of Genes and Genomes
Quantitative real-time PCR
Gillaspy G, Ben–David H, Gruissem W. Fruits: a developmental perspective. Plant Cell. 1993;5(10):1439–51.
Carrera E, Ruiz-Rivero O, Peres LE, Atares A, Garcia-Martinez JL. Characterization of the procera tomato mutant shows novel functions of the SIDELLA protein in the control of flower morphology, cell division and expansion, and the auxin-signaling pathway during fruit-set and development. Plant Physiol. 2012;160(3):1581–96.
Ruan YL, Patrick JW, Bouzayen M, Osorio S, Femie AR. Molecular regulation of seed and fruit set. Trends Plant Sci. 2012;17(11):656–65.
Mesejo C, Yuste R, Reig C, Martínez-Fuentes A, Lglesias DJ, Muñoz-Fambuena N, et al. Gibberellin reactivates and maintains ovary-wall cell division causing fruit set in parthenocarpic citrus species. Plant Sci. 2016;247:13–24.
Lu L, Liang J, Chang X, Yang H, Li T, Hu J. Enhanced vacuolar invertase activity and capability for carbohydrate import in GA-treated inflorescence correlate with increased fruit set in grapevine. Tree Genet Genomes. 2017;13(1):21.
Jong MD, Mariani C, Vriezen WH. The role of auxin and gibberellin in tomato fruit set. J Exp Bot. 2009;60(5):1523–32.
De JM, Wolters-Arts M, Garcia-Martinez JL, Mariani C, Vriezen WH. The Solanum lycopersicum AUXIN RESPONSE FACTOR 7 (SIARF7) mediates cross-talk between auxin and gibberellin signalling during tomato fruit set and development. J Exp Bot. 2011;62(2):617–26.
Matsuo S, Kikuchi K, Fukuda M, Honda I, Imanishi S. Roles and regulation of cytokinins in tomato fruit development. J Exp Bot. 2012;63(15):5569–79.
Li J, Wu Z, Cui L, Zhang T, Guo Q, Xu J, et al. Transcriptome comparison of global distinctive features between pollination and parthenocarpic fruit set reveals transcriptional phytohormone cross-talk in cucumber (Cucumis sativus L.). Plant Cell Physiol. 2014;55(7):1325–42.
Gorguet B, Eggink PM, Ocaña J, Tiwari A, Schipper D, et al. Mapping and characterization of novel parthenocarpy QTLs in tomato. Theor Appl Genet. 2008;116(6):755–67.
Santos RCD, Pereira MCT, Mendes DS, Sobral RRS, Nietsche S, Mizobutsi GP. Gibberellic acid induces parthenocarpy and increases fruit size in the ‘Gefner’custard apple (‘Annona cherimola’x’ Annona squamosa’). Aust J Crop Sci. 2016;10(3):314.
Niu Q, Wang T, Li J, Yang Q, Qian M, Teng Y. Effects of exogenous application of GA4+7 and N-(2-chloro-4-pyridyl)-N’-phenylurea on induced parthenocarpy and fruit quality in Pyrus pyrifolia ‘Cuiguan. Plant Growth Regul. 2015;76(3):251–8.
Wu Z, Zhang T, Li L, Xu J, Qin X, Zhang T, Chen J. Identification of a stablemajor-effect QTL (Parth 2.1) controlling parthenocarpy in cucumber and associated candidate gene analysis via whole genome re-sequencing. BMC Plant Biol. 2016;16(1):182.
George WL, Scott JW, Splittstoesser WE. Parthenocarpy in tomato. Hortic Rev. 1984;6:65–84.
Sardos J, Rouard M, Hueber Y, Cenci A, Hyma KE, Houwe IV, et al. A genome-wide association study on the seedless phenotype in banana (Musa spp.) reveals the potential of a selected panel to detect candidate genes in a vegetatively propagated crop. PLoS One. 2016;11(5):e0154448.
Gorguet B, Heusden A, Lindhout P. Parthenocarpic fruit development in tomato. Plant Biol. 2005;7(2):131–9.
Serrani JC, Carrera E, Ruiz-Rivero O, Gallego-Giraldo L, Peres LE, Garcia-Martinez JL. Inhibition of auxin transport from the ovary or from the apical shoot induces parthenocarpic fruit-set in tomato mediated by gibberellins. Plant Physiol. 2010;153(2):851–62.
Garciahurtado N, Carrera E, Ruizrivero O, Lópezgresa MP, Hedden P, Gong F, et al. The characterization of transgenic tomato overexpressing gibberellin 20-oxidase reveals induction of parthenocarpic fruit growth, higher yield, and alteration of the gibberellin biosynthetic pathway. J Exp Bot. 2012;63(16):5803–13.
Martí C, Orzáez D, Ellul P, Moreno V, Carbonell J, Granell A. Silencing of DELLA induces facultative parthenocarpy in tomato fruits. Plant J. 2007;52(5):865–76.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-seq. Bioinf. 2009;25(9):1105–11.
Yamaguchi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008;59(1):225–51.
Giacmelli L, Rota-Stabelli O, Masuero D, Acheampong AK, Moretto M. Gibberellin metabolism in Vitis vinifera L. during bloom and fruit-set: functional characterization of grapevine gibberelin oxidases. J Exp Bot. 2013;64(14):4403–19.
Zhang ZL, Orawa M, Fleet CM, Zentella R, Hu J, Heo JO, et al. Scarecrow-like 3 promotes gibberellin signaling by antagonizing master growth repressor DELLA in Arabidopsis. Proc Natl Acad Sci U S A. 2011;108(5):2160–5.
Sun TP. Gibberellin-GID1-DELLA: a pivotal regulatory module for plant growth and development. Plant Physiol. 2010;154(2):567–70.
Rose JK, Lee HH, Bennett AB. Expression of a divergent expansin gene is fruit-specific and riprning-regulated. Proc Natl Acad Sci. 1997;94(11):5955–60.
Skinner DJ, Hill TA, Gasser CS. Regulation of ovule development. Plant Cell. 2004;16(Suppl:S32–45.
Stone SL, Kwong LW, Yee KM, Pelletier J, Lepiniec L, Fischer RL, Goldberg RB, Harada JJ. LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development. Proc Natl Acad Sci U S A. 2001;98(20):11806–11.
Tsukagoshi H, Saijo T, Shibata D, Morikami A, Nakamura K. Analysis of a sugar response mutant of Arabidopsis identified a novel B3 domain protein in that functions as an active transcriptional repressor. Plant Physiol. 2005;138(2):675–85.
Suzuki M, Wang H, McCarty DR. Repression of the LEAFY COTYLEDON1/B3 regulatory network in plant embryo development by VP1/ABSCISIC ACID INSENSITIVE3-like B3 genes. Plant Physiol. 2007;143(2):902–11.
Zhong R, Lee C, Ye ZH. Evolutionary conservation of the transcriptional network regulating secondary cell wall biosynthesis. Trends Plant Sci. 2010;15(11):625–32.
Chai L, Li Y, Chen S, Perl A, Zhao F, Ma H. RNA sequencing reveals high resolution expression change of major plant hormone pathway genes after young seedless grape berries treated with gibberellin. Plant Sci. 2014;229:215–24.
Agüero C, Vigliocco A, Abdala G, Tizio R. Effect of gibberellic acid and uniconazol on embryo abortion in the stenospermocarpic grape cultivars Emperatriz and Perlon. Plant Growth Regul. 2000;30(1):9–16.
Chae K, Isaacs CG, Reeves PH, Maloney GS, Muday GK, Nagpal P, Reed JW. Arabidopsis SMALLAUXIN UP RNA63 promotes hypocotyl and stamen filament elongation. Plant J. 2012;71(4):684–97.
Spartz AK, Sang HL, Wenger JP, Gonzalez N, Itoh H, et al. The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J Cell Mol Biol. 2012;70(6):978–90.
Yang JH, Han SJ, Yoon EK, Lee WS. Evidence of an auxin signal pathway, microRNA167-ARFS-GH3, and its response to exogenous auxin in cultured rice cells. Nucleic Acids Res. 2006;34(6):1892–9.
Chung Y, Choe S. The regulation of brassinosteroid biosynthesis in Arabidopsis. Crit Rev Plant Sci. 2013;32(6):396–410.
Fu FQ, Mao WH, Shi K, Zhou YH, Asami T, Yu JQ. A role of brassinosteroids in early fruit development in cucumber. J Exp Bot. 2008;59(9):2299–308.
Lin Z, Arciga-Reyes L, Zhong S, Alexander L, Hackett R, Wilson I, Grierson D. Sl TPR1, a tomato tetratricopeptide repeat protein, interacts with the ethylene receptors NR and Le ETR1, modulating ethylene and auxin responses and development. J Exp Bot. 2008;59(15):4271–87.
Finkelstein RR, Gampala SSL, Rock CD. Abscisic acid signaling in seeds and seedlings. Plant Cell. 2002;14(Suppl:S15–45.
Sun R, He F, Lan Y, Xing R, Liu R. Transcriptome comparison of Cabrnet sauvignon grape berries from two regions distinct climate. J Plant Physiol. 2015;178:43–54.
Dixit R. Plant cytoskeleton: DELLA connects gibberellins to microtubules. Curr Biol. 2013;23(11):479–81.
Deluc LG, Grimplet J, Wheatley MD, Tillett RL, Quilici DR. Transcriptomic and metabolite analyses of cabernet sauvignon grape berry development. BMC Genomics. 2007;8(1):429.
Schloss J, Olsson N, Weis M, Reid K, Peng F, Lund S, Bowen P. Cellular expansion and gene expression in the developing grape (Vitis vinifera L.). Protoplasma. 2008;232(3–4):255–65.
Picton S, Barton SL, Bouzayen M, Hamilton AJ, Grierson D. Altered fruit ripening and leaf senescence in tomatoes expressing an antisense ethylene-forming enzyme transgene. Plant J. 1993;3(3):469–81.
Jones KM, Koen TB, Oakford MJ, Bound SA. Thinning ‘red Fuji’ apples using ethephon at two timings. J Hortic Sci. 1990;65(4):381–4.
Wertheim SJ. Developments in the chemical thinning of apple and pear. Plant Growth Regul. 2010;31(1–2):85–100.
Agarwal G, Garg V, Kudapa H, Doddamani D, Pazhamala LT, et al. Genome-wide dissection of AP2/ERF and HSP90 gene families in five legumes and expression profiles in chickpea and pigeonpea. Plant Biotechnol J. 2016;14(7):1–15.
Gao C, Li P, Song A, Wang H, Wang Y, et al. Isolation and characterization of six AP2/ERF transcription factor genes in chrysanthemum nankingense. Int J Mol Sci. 2015;16(1):2052–65.
Lee SY, Hwang EY, Seok HY, Tarte VN, Jeong MS. Arabidopsis AtERF71/HRE2 functions as transcriptional activator via cis-acting GCC box or DRE/CRT element and is involved in root development through regulation of root cell expansion. Plant Cell Rep. 2015;34(2):223–31.
Huala E, Sussex IM. LEAFY interacts with floral homeotic genes to regulate Arabidopsis floral development. Plant Cell. 1992;4(8):901–13.
Chung MY, Vrebalov J, Alba R, Lee J, Mcquin R, et al. A tomato (Solanumly copersicum) APETALA2/ERF gene, SlAP2a, is a negative regulator of fruit ripening. Plant J. 2010;64(6):936–47.
Niu X, Helentjaris T, Bate NJ. Maize ABI4 binds coupling element1 in abscisic acid and sugar response genes. Plant Cell. 2002;14(10):2565–75.
Cheng C, Xu X, Singer SD, Li J, Zhang H, Gao M, Wang L, Song J, Wang X. Effect of GA3 treatment on seed development and seed-related gene expression in grape. PLoS One. 2013;8:e80044.
Gallavotti A, Barazesh S, Malcomber S, Hall D, Jackson D, Schmidt RJ, et al. Sparse inflorescencel encodes a monocot-specific YUCCA-like gene required for vegetative and reproductive development in maize. Proc Natl Acad Sci U S A. 2008;105(39):15196–201.
Sastry KK, Muir RM. Gibberellin: effect on diffusible auxin in fruit development. Sci. 1963;140(3566):494–5.
Niu Q, Zong Y, Qian M, Yang F, Teng Y. Simultaneous quantitative determination of major plant hormones in pear flowers and fruit by uplc/esi-ms/ms. Anal Methods. 2014;6(6):1766.
Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, et al. The arabidopsis, cytochrome P450 CYP707A encodes ABA 8′-hydroxylases: key enzymes in ABA catabolism. EMBO J. 2004;23(7):1647–56.
Langmead B, Trapnell C, Pop M, Salzberg S. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment. Bioinf. 2009;25(15):1966–7.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G. Vanbaren MJ et al. transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinf. 2010;26(1):136–8.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative pcr and the 2-ΔΔct method. Methods. 2001;25(4):402–8.
We are grateful to the editor and the reviewers for their valuable suggestions and comments which have improved our manuscript. We thank Wei Xiao and Qiuping Tan for their technical assistance.
This study was funded by Shandong Province Modern Agricutural Technology System Fruit Innovation Team (Grant number: SDAIT-06-01). The funding body had no role in the design of the study, data collection and analysis, interpretation of data, or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
The author declares that he has no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The primers used in this article. (DOC 78 kb)
Expression of DEGs induced parthenocarpy by GA3. (DOC 34 kb)
Expression of DEGs associated with cell division after GA3 treatment. (DOC 20 kb)
Expression of DEGs associated with fruit setting after GA3 treatment. (DOC 22 kb)
the workflow chart of the analytical process. (JPG 105 kb)
About this article
Cite this article
Wen, B., Song, W., Sun, M. et al. Identification and characterization of cherry (Cerasus pseudocerasus G. Don) genes responding to parthenocarpy induced by GA3 through transcriptome analysis. BMC Genet 20, 65 (2019). https://doi.org/10.1186/s12863-019-0746-8