Skip to main content

GA3 is superior to GA4 in promoting bud endodormancy release in tree peony (Paeonia suffruticosa) and their potential working mechanism

Abstract

Background

Sufficient low temperature accumulation is the key strategy to break bud dormancy and promote subsequent flowering in tree peony anti-season culturing production. Exogenous gibberellins (GAs) could partially replace chilling to accelerate dormancy release, and different kinds of GAs showed inconsistent effects in various plants. To understand the effects of exogenous GA3 and GA4 on dormancy release and subsequent growth, the morphological changes were observed after exogenous GAs applications, the differentially expressed genes (DEGs) were identified, and the contents of endogenous phytohormones, starch and sugar were measured, respectively.

Results

Morphological observation and photosynthesis measurements indicated that both GA3 and GA4 applications accelerated bud dormancy release, but GA3 feeding induced faster bud burst, higher shoot and more flowers per plant. Full-length transcriptome of dormant bud was used as the reference genome. Totally 124 110 459, 124 015 148 and 126 239 836 reads by illumina transcriptome sequencing were obtained in mock, GA3 and GA4 groups, respectively. Compared with the mock, there were 879 DEGs and 2 595 DEGs in GA3 and GA4 group, 1 179 DEGs in GA3 vs GA4, and 849 DEGs were common in these comparison groups. The significant enrichment KEGG pathways of 849 DEGs highlighted plant hormone signal transduction, starch and sucrose metabolism, cell cycle, DNA replication, etc. Interestingly, the contents of endogenous GA1, GA3, GA4, GA7 and IAA significantly increased, ABA decreased after GA3 and GA4 treatments by LC–MS/MS. Additionally, the soluble glucose, fructose and trehalose increased after exogenous GAs applications. Compared to GA4 treatment, GA3 induced higher GA1, GA3 and IAA level, more starch degradation to generate more monosaccharide for use, and promoted cell cycle and photosynthesis. Higher expression levels of dormancy-related genes, TFL, FT, EBB1, EBB3 and CYCD, and lower of SVP by GA3 treatment implied more efficiency of GA3.

Conclusions

Exogenous GA3 and GA4 significantly accelerated bud dormancy release and subsequent growth by increasing the contents of endogenous bioactive GAs, IAA, and soluble glucose such as fructose and trehalose, and accelerated cell cycle process, accompanied by decreasing ABA contents. GA3 was superior to GA4 in tree peony forcing culture, which might because tree peony was more sensitive to GA3 than GA4, and GA3 had a more effective ability to induce cell division and starch hydrolysis. These results provided the value data for understanding the mechanism of dormancy release in tree peony.

Peer Review reports

Background

Bud dormancy of woody perennial plants is an adaptive mechanism that allows them to survive in winter. The annual dormancy is classified into paradormancy, endo dormancy, and ecodormancy based on the internal or external repression signals [1]. In general, environmental signals, such as low temperatures and short-day control endodormancy induction and maintenance [2]. The effective accumulation of low temperature and several chemicals, such as hydrogen cyanamide (HC), mineral oil, potassium nitrate, gibberellin (GA) and 5-azacytidine (5-azaC), could promote dormancy release and induce bud break of deciduous trees in temperate areas [3,4,5].

Tree peony (Paeonia suffruticosa Andrews), originated from China, is well known as a horticultural plant in the world for its beautiful flower type and color. Tree peony has a short and concentrated florescence every year, so anti-season culture forms a new industrial chain and occupies a great part in tree peony industry, and the key problem is to break bud dormancy, which is necessary to ensure the quality of flowering. As well as other perennial plants, like sweet cherry, grape, apple, peach, and kiwifruit [6,7,8,9,10], the bud dormancy of tree peony is a biological characteristic and also a necessary process before bud sprouting and flowering [11].

The buds require sufficient chilling to break endodormancy and transit into the ecodormant state, then they will burst when the comfortable environment returns. GAs had been demonstrated to play important roles in the activity-dormancy-activity cycle [4, 12,13,14]. For example, Nell et al. found that foliar applications of GA3 or GA4,7 accelerated the flowering and increased the pedicel length and flower size in azalea [15]. In flower buds of Japanese apricot, GA4 significantly accelerated dormancy release, and the promotion of dormancy release was related to energy metabolism activation by proteomic and transcriptomic analysis [16]. When Paeonia lactiflora endured insufficient chilling accumulation, the GA3 treatment effectively broke dormancy, accelerated sprouting and promoted subsequent growth and flowering [17]. In summary, different kinds of gibberellins show the inconsistent effect on breaking dormancy in various plants. Then, which gibberellin (GA3 or GA4) usually used in production is more effective to break bud dormancy of tree peony?

Recent reports showed that several genes associated with GA synthesis and signal transduction were involved in bud dormancy regulation. Combination of transcript profiling and microarray analysis, Gai et al. found that the activation of GA pathway played a core role in the regulation of tree peony bud dormancy release [18]. It is well known that the key genes associated with gibberellin synthesis and signaling are GA 20-oxidase gene (GA20ox), GA 3-oxidase gene (GA3ox), GA20ox, GA3ox, GA2ox, GID1, DELLA and GID2 [19]. Bioactive GAs is synthesized through complex pathways, of which GA20ox and GA3ox are the key rate-limiting enzymes [20]. In tea, the differentially-expressed genes (DEGs), CsGA3ox and CsGA20ox, played an important role in regulating the bud activity-dormancy transition [14]. Conversely, PmRGL2 encoding a DELLA protein acted as a negative regulator of dormancy release through GA signaling pathway in Japanese apricot [21]. The F-box protein GID2, a subunit of the SCF E3 ubiquitin ligase complex, recognizes the GA-GID1-DELLA complex and helps to degrade the DELLA through the 26S proteasome, leads to de-repression and transcriptive activation of GA-responsive genes [20].

Abscisic acid (ABA) is generally believed to be another core hormone to maintain and regulate bud dormancy with antagonistic effect to GA [12]. The increasement of the ABA content benefits to induce dormancy in autumn, while its content decrease is the trigger for dormancy release in poplar, grape and pear [10, 22,23,24,25]. In poplar, ABA content increased after short day exposure and reached a peak after growth cessation in the apex, and the key genes involving in ABA biosynthesis and ABA signal transduction were also induced [26]. ABA bounds to members of the PYL/RCAR ABA receptor family that initiates signal transduction inhibiting type 2C protein phosphatases [27]. The ABA-responsive element (ABRE) is the major cis-element for ABA-responsive gene expression, and the ABRE-dependent gene expression is regulated by ABRE-binding factor (ABF). In peach, PpABF2 interacted with PpTCP20 to form heterodimers and regulate bud endodormancy. PpyABF3 could bind to the promoter of PpyDAM3 to activate its expression, while PpyABF2 interrupted this activation role by binding with PpyABF3 during the pear endodormancy maintenance [9]. Recent report showed that SHORT VEGETATIVE PHASE (SVP)-like (SVL) played a vital role in the dormancy of poplar. Prolonged chilling deceased ABA contents and SVL expression, resulting in the induction of FT1 expression and GA biosynthesis, which accelerated dormancy release [28].

Bud dormancy-growth transition is accompanied by many metabolic and developmental processes, including carbohydrate metabolism, photosynthesis, cell division, etc. As early as 1990, Bonicel et al. indicated that deep dormancy of poplar bud was characterized by low carbohydrate levels, and the end of dormancy was marked by major changes that starch levels declined and sugars rose [29]. During dormancy release of Japanese pear, a decrease of starch content in the shoot and an increase of soluble sugars in both the flower bud and the shoot were observed [30]. Starch, a polymer of glucose molecules, is the predominant form of carbohydrate storage in cells. In poplar, the activation of carbohydrate metabolism-related proteins mainly occurred at the stage of endodormancy release [31]. Carbohydrate metabolism activation and variation of plant hormone signal transduction were also found in the dormancy release process of Lilium pumilum [32] and tree peony [33]. In addition, the cell proliferation-related genes including MAP kinase, cyclins, and cyclin dependent-kinases, etc., were assumed to involve in the dormancy-growth transition. In grapevine, cell cycle genes encoding cyclin dependent-kinase (CDKB1 and CDKB2) and cyclins (CYCA, CYCB and CYCD) were activated during bud endodormancy release [34]. These results indicated that dormancy break was accompanied by cell cycle and cell division reactivation.

Additionally, recent studies have identified several dormancy release related genes. FLOWERING LOCUS T (FT) was believed to play vital role in dormancy release [35], while its paralog CENTRORADIALIS/ TERMINAL FLOWER 1 (CEN/TFL1) acted as a negative regulator of dormancy release in poplar. TFL was also assumed to be a marker of bud burst [35, 36]. AP2 type transcript factor EBB1 (EARLY BUD-BREAK 1) accelerated dormancy release and bud break, and EBB3 was part of the EBB1/SVL-mediated regulatory mechanism, which promoted the expression of CYCD3.1 and positively regulated bud-break with activation of cell cycle in poplar. SVL involved in the same pathway with a negative manner [22, 37].

In this study, the morphological changes indicated that both GA3 and GA4 greatly accelerated bud dormancy release, but the effect of GA3 was better than that of GA4 to accelerate bud burst and subsequent growth in tree peony. In order to understand how exogenous GA3 or GA4 affected dormancy release, PacBio full-length sequencing in the dormant bud of tree peony was performed to obtain high-quality transcriptome data as reference genome, the DEGs of tree peony buds with exogenous GA3 or GA4 application were obtained by illumina transcriptome sequencing, and the contents of sugar and endogenous phytohormones were measured, respectively. All results would provide valuable data for the mechanism of dormancy release in tree peony.

Results

Morphology changes after exogenous GA3 and GA4 treatments

Despite past study reported that GA application could effectively promote endodormancy release and accelerate bud sprouting in tree peony [38], the effects of different GAs were still unknown. Therefore, GA3 and GA4, commonly used in production, were applied to compare their effects on dormancy release and subsequent growth. Exogenous GA3 and GA4 induced the bud burst after 4 d applied, and almost all apical buds sprouted 3 d thereafter, but it was 1 d earlier in GA3 group than in GA4 group to reach the peak. In contrast, the sprouting rate in mock group was 19.24% after 16 d, 96.20% in GA3 treated group (p = 0.000), and 95.82% in GA4 treated group (p = 0.000) (Fig. 1b). After GA3 and GA4 treated for 16 d, their morphology was investigated, including the shoot height, the length of branch and leaf, and flowering rate. GA3 and GA4 feedings significantly promoted the leaf length, shoot height and flowering rate (Fig. 1a, c). Compared with the mock group, GA3 application increased the height of shoot (average 30.83 cm, p = 0.014) and the length of leaf (average 14.89 cm, p = 0.010), and its effect was better than GA4, while the difference of branch length was not significant in GA3 (p = 0.376) and GA4 (p = 0.373) treated groups (Fig. 1c). Additionally, the flowering rate induced by GA3 was significantly higher than that by GA4 treatment (p = 0.039) (Fig. 1d). Altogether, these changes of morphology indicated that both GA3 and GA4 applications effectively accelerate bud dormancy release, burst, growth and flowering, but GA3 was superior to GA4.

Fig. 1
figure 1

Effect of GA3 and GA4 applications in P. suffruticosa cv ‘Luhehong’. a Morphological observation after being transferred into greenhouse for 16 d by exogenous GA3 and GA4 applications, the samples treated by sterile distilled water was used as mock. Bar = 4 cm; b The sprouting percentage after being transferred into greenhouse for 16 d by exogenous GA3 and GA4 applications; c Branch height of shoots, the length of the first leaf and flowering rate after being transferred into greenhouse for 16 d by exogenous GA3 and GA4 applications; d Flowering rate after being transferred into greenhouse for 40 d by exogenous GA3 and GA4 applications. e Effect of GAs treatment on the photosynthesis of tree peony ‘Luhehong’ leaves. Duncan’s test at p < 0.05 after analysis of variance; data were shown as mean ± SD in three biological replicates. The lowercase letters on the histogram represented the differentiation

As known, photosynthesis is an energy resource for shoot growth and flowering. Therefore, the photosynthetic parameters in the leaves after GA3 and GA4 feedings for 16 d were measured, respectively. Compared to mock group, photosynthetic rate significantly increased after GA3 (p = 0.000) and GA4 (p = 0.002) treatments. The stomatal conductance and transpiration rate by GA3 application (p = 0.006) were greater than that in mock group, whereas there was no difference between GA4 treatment and mock (p = 0.504). In comparison with mock group, the intercellular CO2 concentrations were not different among GA3 (p = 0.204) and GA4 (p = 0.693) groups (Fig. 1e).

The PacBio full-length transcriptome in tree peony during chilling-induced dormancy release

According to previous data, the samples after 0, 7, 14, 21 and 28 d chilling durations represented the whole process from endodormancy, endodormancy release to ecodormancy. Until now, there is still no high-quality genome data of tree peony. Therefore, the samples of different chilling durations (0, 7, 14, 21 and 28 d) were mixed and used to construct a pool to obtain a wide coverage of the tree peony transcriptome. A total of 18 470 579 reads were obtained with the max length of 200 013 bp and the min length of 50 bp, and the mean length was 1 778 bp. After removing the low-quality raw reads through the three steps including CCS, the classify and cluster (Isoseq) (https://github.com/PacificBiosciences/IsoSeq_SA3nUP wiki# datapub), 38 288 high-quality sequences with the length ranged from 79 bp to 11 471 bp were finally obtained. Of them, the quantity of reads with the length between 1 001 and 2 000 bp (11 246) was the most, followed by the reads with the length between 2 001 and 3 000 bp (10 261) (Fig. 2a). The generated isoforms were then filtered with CD-HIT software based on an identify of 98%. The full-length transcript completion was examined using eukaryote_odb9 database containing single-copy gene sets of 100 species. After removing redundancy, there were 89.1% full-length transcripts matched to the complete BUSCO database, of which the number of complete single-copy unigenes accounted for 31.7%, complete duplicated unigenes for 57.4%, the fragmented unigenes for 0.7%, and the quantity of the missing unigenes for 10.2% (Fig. 2b). Finally, total 30 035 non-redundant full-length unigenes were obtained and used for the subsequent analysis.

Fig. 2
figure 2

The length distribution, BUSCO analysis of transcript completeness, GO analysis and KEGG pathway classification of the assembled unigenes in P. suffruticosa cv ‘Luhehong’. a The length distribution of the assembled unigenes; b BUSCO analysis of transcript completeness; c GO analysis of the assembled unigenes; d KEGG pathway classification of the assembled unigenes

The 30 035 non-redundant unigenes were performed to annotate in different public databases, and most of them were longer than 1 000 bp (Additional file 1). A total of 23 361 unigenes were annotated to GO classification terms, and there were 19 781, 21 120 and 19 843 unigenes classified into biological process, cellular component, and molecular function category, respectively (Fig. 2c). In the “biological process” category, the unigene numbers classified into cellular process (16 435), metabolic process (14 128) and biological regulation (7 388) were in the majority. In the “cellular component” category, the majority of the unigenes were represented by cell (19 827), cell part (19 802) and organelle (16 200). Terms of binding (14 670) and catalylic activity (11 185) occupied a clear majority in “molecular function” category. In total, 13 155 unigenes annotated in KEGG database, and the most majority of unigene were enriched in translation (1 942), signal transduction (1 351), carbohydrate metabolism (1 329) and folding, sorting and degradation (1 195) (Fig. 2d). In addition, there were 28 420 (94.62%), 27 609 (91.92%), 25 071 (83.47%), and 18 379 (61.19%) unigenes matched sequences in the Nr, eggnog, Swissprot, and KOG databases, respectively.

In addition, the annotated unigenes were blasted to the transcription factor (TF) database, and total 8 019 unigenes were found to encode potential TFs (Additional file 2). Most of the TFs were mainly distributed in the families of C3H, MADS, FAR1, PHD, bHLH, MYB-related, FHA, C2H2, NAC, AP2-EREBP, etc. The opening reading frames (ORFs) of these assembled unigenes were predicted by TransDecoder with priority order of Nr, Swissprot and KOG, and hit total 29 418 unigenes, occupied by 97.9% of the non-redundant unigenes. The results indicated that the transcriptome data was high quality and informative.

Illumina transcriptome sequencing in tree peony after exogenous GA3and GA4 treatments

As indicated above, exogenous GA3 and GA4 had a different effect on dormancy release and subsequent growth, and we wanted to know what molecular mechanism hid in this process? In order to answer this question, the illumina transcriptome sequencing of mock, GA3, and GA4 treated buds were carried out. After quality control, a total of 144.55 M, 143.21 M and 144.03 M clean reads were screened out from 148.19 M, 147.94 M, and 147.66 M raw reads, respectively. The average Q30 percentage (percentage of bases whose quality was greater than 30 in clean reads) and GC percentage were 94.32% and 45.14%, respectively (Table 1). Totally 124 110 459, 124 015 148 and 126 239 836 reads in mock, GA3 and GA4 groups were mapped to the above reference unigenes, and the average mapped percentages were 85.86%, 86.59% and 87.65%, respectively (Table 1).

Table 1 Summary of the Illumina sequencing output for flower buds after exogenous GA3 and GA4 applications

DEGs after exogenous GAand GA4 treatments

In order to investigate the repeatability among biological replicates and the correlation between different samples, PCA analysis was performed based on the transcription abundances of all the unigenes (Additional file 3). The result indicated that biological triplicates had high reliability for each of the repeats grouped together, respectively.

A large number of DEGs were obtained using DEGseq software in the comparisons of GA3 vs mock, GA4 vs mock, GA3 vs GA4 (Table 2). Compared with the mock group, the DEGs number of GA4 treated group was more than that of GA3 treated group. In GA3 vs mock, totally 879 DEGs were obtained, including 405 up-regulated and 474 down-regulated DEGs, respectively. In GA4 vs mock, there were 2 595 DEGs, including 1 187 up-regulated and 1 408 down-regulated DEGs, respectively. Between GA3 and GA4 treatment, a total of 1 179 DEGs were obtained, of which 547 up-regulated and 632 down-regulated DEGs.

Table 2 Statistics of the DEGs in different comparisons

The GO classification and KEGG enrichment analysis of DEGs were further performed. Totally 652, 870, and 1 922 DEGs were divided into 49, 50, and 49 sub-categories in GA3 vs mock, GA4 vs mock, and GA3 vs GA4, respectively. (Additional file 4). Among these three comparison groups, the major sub-categories were “cellular process” in “biological process” category with 1 260, 474 and 592 DEGs, “cell” in “cellular component” category with 1 445, 433, and 626 DEGs, and “binding” in “molecular function” category with 1 167, 427, and 526 DEGs, respectively. A total of 227, 838, and 354 DEGs were annotated by the KEGG database and classified into 21, 24 and 23 sub-categories, respectively (Fig. 3). In the top 20 KEGG enrichment terms, many of the DEGs were enriched in plant hormone signal transduction with 15, 57, and 22, starch and sugar metabolism with 11, 44 and 22, respectively (Fig. 3).

Fig. 3
figure 3

KEGG pathway classification and KEGG enrichment of DEGs (-value < 0.05 & |log2FC|> 1) after GA3 and GA4 applications. a GA3 vs mock; b GA4 vs mock; c GA3 vs GA4

In order to know the reason for exogenous GA3 and GA4 effectively accelerating bud dormancy release and burst, and their different effects on plant morphological changes, a Venn diagram was constructed to show the overlapping DEGs (Fig. 4a). There were 220 and 712 DEGs overlapping between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock, respectively, and 83 DEGs were common in these two comparison groups. Among the overlapping 849 DEGs, a total of 619 were assigned to GO terms, and the major sub-categories were “cellular process” and “metabolic process” in “biological process” category with 389 and 367 DEGs, “cell” and “cell part” in “cellular component” category with 419 and 419 DEGs, and “binding” and “catalytic activity” in “molecular function” category with 374 and 311 DEGs, respectively (Fig. 4b). KEGG pathway enrichment of these 849 overlapping DEGs was performed, and they were significantly enriched in 11 KEGG pathways (p < 0.05), including “plant hormone signal transduction” with 14 DEGs, “starch and sucrose metabolism” with 11 DEGs, “endocytosis” with 10 DEGs, “Galactose metabolism” with 9 DEGs, “alpha-Linolenic acid metabolism” with 6 DEGs, “pentose and glucuronate intercom” with 6 DEGs, “Cell cycle” with 6 DEGs, “Ascorbate and aldarate metabolism” with 5 DEGs, “Phenylpropanoid biosynthesis” with 5 DEGs, “DNA replication” with 4 DEGs etc. (Fig. 4c).

Fig. 4
figure 4

GO classification and KEGG enrichment of the 849 overlapping DEGs between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock, respectively. a Venn diagram of DEGs among different comparison groups; b GO classification of the 849 (137 + 83 + 629) overlapping DEGs between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock; c KEGG enrichment of the 849 (137 + 83 + 629) overlapping DEGs between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock

Identification of plant hormone and signaling transduction, sugar metabolism, and cell cycle-related DEGs after GAand GA4treatments

Based on the Venn diagram and KEGG enrichment of the overlapping 849 DEGs (Additional file 5) between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock, we constructed the phylogenetic trees of the related DEGs (Additional file 6), which enriched in the KEGG pathways about plant hormone signal transduction, starch and sucrose metabolism, cell division and DNA replication, and analyzed their expression patterns by the results of RNAseq data and real-time quantitative RT-PCR (qPCR) (Fig. 5, Fig. 6). Totally, 14 DEGs were related to plant hormone and signaling transduction, including GA, ABA, Auxin and ethylene pathway (KEGG: ko04075). The relative expression of the related unigenes validated by qPCR was similar to the RNAseq data (Fig. 5a and b, Additional file 7). Among these DEGs, a F-box protein gene PsGID2, encoding a subunit of ubiquitin E3 ligases, decreased after GA3 and GA4 treatments, and the down-regulated level of GA4 treatment was greater. The expression levels of two ABA responsive element binding factor gene (PsABF) and three ABA receptor PYR/PYL family members (PsPYL8, PsPYL12 and PsPYL1) were down-regulated after exogenous GAs treated. Among four auxin-related DEGs, two DEGs encoding indole-3-acetic acid-amido synthetase GH3 family members including PsGH3.6 and PsGH3.1 significantly increased after GAs applications, of which the expression levels of PsGH3.1 were increased 2.37-fold by GA3 and 24.22-fold by GA4 treatment, respectively. RNAseq showed that the expressions of two ethylene receptor 2 (PsETR2) were up-regulated after GA3 and GA4 treatments, with more than 2.28-fold increasing by GA4, and the ethylene-insensitive protein gene (PsEIN2) was significantly inhibited by GA4 treatment. One EIN3-binding F-box protein gene (PsEBF1) was significantly down-regulated after GA4 application. In addition, a transcription factor encoding gene, PsMYC2, was extremely induced by GAs applications (Fig. 5a, b).

Fig. 5
figure 5

The relative expression levels of plant hormone signal transduction-related DEGs and the contents of endogenous GAs, ABA, and IAA after GA3 and GA4 applications in P. suffruticosa cv ‘Luhehong’. a Heatmap of plant hormone signal transduction-related DEGs after GA3 and GA4 applications based on the RNAseq data. The colors from blue to red in heatmap represented the relative expression level of DEGs from low to high, and the blue-red colors corresponded to the 0–1 normalization (minimum–maximum normalization) values of log2(FPKM). b The relative expression levels of plant hormone signal transduction-related DEGs after GA3 and GA4 applications by qPCR. c The contents of endogenous GAs, ABA, and IAA after 48 h GA3 and GA4 feedings, and noted GA3-48 h, GA4-48 h. The mean ± SD in three biological replicates was shown. *, ** indicated significant differences of one-way ANOVA at p < 0.05, p < 0.01, respectively

Fig. 6
figure 6

The expression levels of starch and sucrose metabolism and cell cycle-related DEGs after GA3 and GA4 applications in P. suffruticosa cv ‘Luhehong’. a Heatmap of starch and sucrose metabolism-related DEGs after GA3 and GA4 applications based on the RNAseq data. b The expression levels of starch and sucrose metabolism -related DEGs after GAs applications by qPCR. c Heatmap of cell cycle-related DEGs enriched in KEGG pathway map based on the RNAseq data. d The expression levels of cell cycle-related DEGs after GAs applications by qPCR. e The expression level of PsFNR after GAs applications by qPCR. Data were represented as mean of three different determinations ± SD. Asterisks indicate statistically significant differences (one-way ANOVA, *p < 0.05, **p < 0.01)

The relative expression of starch and sugar metabolism and cell cycle-related DEGs were also validated by qPCR, and the results indicated that the expression profiling by qPCR was consistent with that of RNAseq data (Fig. 6, Additional file 7). Among the 11 starch and sugar metabolism-related DEGs (KEGG: ko00500), starch hydrolase beta-amylase gene PsBMY was only induced by GA3 treatment, and hardly detected in GA4-treated group. Two beta-glucosidase genes (PsBGLs) were down-regulated after GA4 feeding, but up-regulated by GA3 application. PsTPS5 (trehalose 6-phosphate synthase gene) showed up-regulation by GA3 treatment, but down-regulation by GA4 treatment, while PsTPS6 was inhibited by both GA3 and GA4 treatments. Two PsUGDs (PsUGD1 and PsUGD3) encoding UDP glucose 6-dehydrogenase were obviously induced by GA4 application. One glucose-1-phosphate adenylyltransferase gene PsglgC and two pectinesterase genes (PsPEs) were down-regulated both by GA3 and GA4, respectively. PsGAUT9, encoding alpha-1,4-galacturonosyltransferase involving in pectin biosynthesis was dramatically induced by GA4 (Fig. 6a, b).

In addition, among 6 DEGs enriched in cell cycle pathway (KEGG: ko04111), structural maintenance of chromosome 4 gene (PsSMC4) significantly increased after GA3 and GA4 treatments, and its expression level increased 286.13-fold after GA4 feeding compared to mock. Two DNA replication licensing factor genes PsMCM2 showed significantly decreased after GA3 treatment but increased by GA4 feeding (Fig. 6c, d). The expression level of Ferredoxin-NADP+ oxidoreductase gene (PsFNR), involving in photosynthesis pathway (KEGG: ko00195), significantly up-regulated after GA3 and GA4 treatments, and its expression level increased 4.15-fold after GA3 treatment compared with mock group (Fig. 6e).

Quantification of endogenous phytohormone, sugar and starch after GAand GAapplications

Because plant hormone and signal transduction-related DEGs were enriched in KEGG analysis, the contents of endogenous phytohormone, including GAs, ABA, and IAA, were measured. The result showed that the contents of GA1, GA3, GA4 and GA7 increased after GA3 and GA4 treatments, of which the contents of endogenous GA1 and GA3 increased about 57.19-fold and 48.75-fold after exogenous GA3 treated for 48 h, and GA4 and GA7 increased about 59.37-fold and 68.66-fold after 48 h exogenous GA4 feeding. The results indicated that intake of exogenous GAs enhanced the conversion of bioactive GAs. Additionally, the transcript of PsGA20ox, the key enzyme of gibberellin synthesis, was analyzed by qPCR (Fig. 5b). The expression levels of PsGA20ox were down-regulated by both GA3 and GA4 applications, and which decreased 5.92-fold by GA4, suggesting the exogenous application of GAs depressed their endogenous synthesis and feedback regulation. The contents of endogenous ABA significantly decreased after GA3 and GA4 treatments, which were consistent with the down-regulated expression of ABA receptor (PsPYL1, PsPYL8, PsPYL12) (Fig. 5b, c). In addition, the contents of endogenous IAA were up-regulated by GA3, but there was no significant difference between GA3 and GA4 application.

The results of RNAseq and qPCR showed that PsBMY was only induced by GA3, and the expressions of PsBGLs and PsAMY significantly increased by GA3 (Fig. 6a, b). Starch contents significantly decreased after GA3 and GA4 feedings evaluated by anthrone colorimetry method, and more dramatically in GA3 group, which suggested that GA3 application induced more polysaccharides to be hydrolyzed. To verify these, the soluble sugar contents of different components were determined by GC–MS/MS (Table 3). Totally, 13 different soluble sugar components besides lactose were determined. Among all these sugar components, sucrose was the most abundant among disaccharides, and its content decreased after both GA3 and GA4 applications. The contents of inositol, fructose and glucose were affluent among 9 monosaccharides, and glucose and fructose significantly increased in GA3 and GA4 treated groups. Interestingly, the content of inositol was significantly increased by GA3 but decreased by GA4 application. Additionally, the amounts of galactose, rhamnose and fucose significantly increased only after GA3 application, while maltose was significantly induced only by GA4. The content of sorbose rose, but no significant difference between GA3 and GA4.

Table 3 Contents of different soluble sugar components and starch in tree peony bud after GA3 and GA4 applications

Expression patterns of the known genes related to dormancy release after GAand GAtreatments

Based on the recent studies, the relative expression patterns of six genes including PsFT, PsFTL, PsSVP, PsEBB1, PsEBB3, PsCYCD were identified by qPCR (Fig. 7). PsFT showed significant increasing after both GA3 and GA4 feedings, and there was no significant difference between these two treatments, and PsTFL was upregulated by GA3 and GA4 treatments.PsSVP, a MADS-box gene with homology to SVP in Arabidopsis, were significantly inhibited by GA3 and GA4 feedings, of which its expression level was down-regulated about 1.96-fold by GA4. PsEBB1 and PsEBB3, encoding transcription factor of the AP2/ERF family, among them PsEBB1 was significantly induced by GA3 and GA4 treatments, and PsEBB3 showed significant upregulation only by GA3 feed. PsCYCD, encoding an important protein at cell cycle check point from G1 phase to S phase, quickly increased its expression level after GA3 and GA4 feedings, of which it was upregulated about 2.65-fold by GA3 treatment.

Fig. 7
figure 7

The relative expression levels of known dormancy-associated genes after GA3 and GA4 applications in P. suffruticosa cv ‘Luhehong’. Data were represented as mean of three different determinations ± SD. Asterisks indicate statistically significant differences (one-way ANOVA, *p < 0.05, **p < 0.01)

Discussion

Bud dormancy transition is a vital developmental process for perennial plant survival and is precisely regulated by complex endogenous signaling networks, such as hormone signaling, sink/ source energy conversion, and multiple environmental signals, including low temperature and photoperiod. Tree peony is a well-known ornamental and medicinal plant with endodormancy characteristics. In recent years, the study on endodormancy regulation mechanisms became more and more urgent for the industrial development of anti-season culture. Practically, exogenous GA3 is a common method to assist dormancy release and bud burst in the anti-season production of tree peony [4, 39]. In some cases, such as Japanese apricot, GA4 treatment iswidely utilized to induce earlier bud break and sub-sequential growth [16]. Until now, the physiological effects of exogenous GA3 and GA4 on the dormancy release in tree peony were unknown. In this study, the morphological results showed that both GA3 and GA4 markedly accelerated dormancy release. However, the effect of GA3 was superior to GA4 in promoting bud dormancy release, following bud burst and subsequent growth in tree peony. In order to unlock the possible molecular mechanism of exogenous GAs accelerating bud endodormancy release, we performed a comprehensive transcriptome analysis of GA3 and GA4 treatments. Our results indicated that endogenous signals such as sugar metabolism, plant hormone signal transduction, including GA synthesis and signaling, ethylene, ABA, IAA, etc., cell cycle and photosynthesis played important roles in promoting bud endodormancy release and bud burst.

Exogenous GAs applications regulate phytohormones metabolism and signaling

GAs are fundamental phytohormones that extensively regulate plant growth and development, especially bud dormancy and sprouting transition in perennial plants [16, 38, 40]. In tree peony, Mornya et al. considered that the accumulated GA was the key regulator to bud dormancy and growth [38]. Meanwhile, the expression profile of GA synthesis and signalling-releated genes and the up-regulated contents endogenous GAs during chilling-induced dormancy release, together suggested that GA pathway participated in the dormancy release in tree peony [41]. In this study, exogenous GA applications increased the contents of endogenous bioactive GAs including GA1, GA3, GA4 and GA7, of which GA1 and GA3 were more abundant in GA3-treated group than in GA4-treated group, while GA4 and GA7 were more in GA4 group (Fig. 5). The results suggested tree peony buds were more sensitive to GA3 considering the morphological effects of the two exogenous GAs, and exogenous GAs promoted dormancy release and subsequent bud development by the enhancement of endogenous bioactive GAs.

At the same time, the high levels of bioactive GAs feedback-regulate the expression of key GA synthesis enzyme genes, including GA20ox, GA3ox, GA2ox, etc. In tea, the transcripts of CsGA3ox-2 and CsGA20ox-1 were down-regulated, whereas CsGA2ox was up-regulated by exogenous GA3 treatment [14]. In addition, GA signaling transduction is also feedback-regulated by exogenous GA treatment [42]. Plant responds to GA by inducing the degradation of DELLAs, which depends on SCF E3 ubiquitin ligase complex recognizing and combining the GA-GID1-DELLA, and GID2 is an important subunit of E3 ubiquitin ligase [43]. In our results, DEGs were mainly enriched in plant hormone signal transduction after exogenous GAs applications, among which, PsGA20ox and PsGID2 were down-regulated after GA3 and GA4 applications (Fig. 5). These results suggested exogenous GA3 and GA4 treatments promoted the interconversion of endogenous GAs, and the increasing contents of endogenous GAs feedback regulated the GA biosynthesis and signal transduction in tree peony.

ABA is another key hormone affecting dormancy release and bud break. The contents of ABA were down-regulated under spring temperature in tree peony [38]. ABA is responsible for the initiation and maintenance of endodormancy, as the increasing content of endogenous ABA presented at the onset of endodormancy, while the bud dormancy release is linked to declining ABA levels in several plants [44, 45]. Meanwhile, ABA catabolism by up-regulation of VvA8H-CYP707A4, an ABA 8’-hydroxylase, accelerated dormancy release of grapevine buds [10]. Therefore, ABA regulates endodormancy throughout its synthesis, metabolism and signalling process. In this study, several genes involving in ABA metabolic and signal transduction pathway were down-regulated after exogenous GA3 and GA4 treatments. For example, three ABA receptor genes PsPYLs and two ABA responsive element binding factor genes PsABFs decreased significantly after GAs applications, consistent with the decreasing contents of endogenous ABA (Fig. 5). It was speculated that exogenous GAs accelerated dormancy release partially by inhibiting endogenous ABA synthesis and signaling.

It is believed that Dormancy-associated MADS-box (DAM) functions downstream of the ABA signalling pathway to regulate bud dormancy. Recent study found that pear ABRE-BINDING FACTOR3 (PpyABF3) directly bound to the PpyDAM3 promoter and positively enhance its expression [9]. In poplar, SVL, a DAM homolog, negatively regulates GA pathway but positively accelerates ABA to induce and maintain dormancy, and prolonged chilling depresses SVL expression and induces FT1 and GA biosynthesis [23, 28]. In our study, PsABFs and PsSVP significant decreased after GAs feedings, and the down-regulation of PsABFs in GAs treated buds suggested that GAs might also adjust dormancy through a DAM pathway, and there might be a loop circuit between GA and ABA in dormancy regulation of tree peony.

Auxin is an important hormone involving in many plant developmental processes. In Pinus Sylvestris, an increase of IAA level was induced by exogenous GA4 in the cambial region, resulting in stimulation of tracheid production in the terminal shoot of intact plants. Similarly, the IAA content and the related genes were up-regulated from the dormancy to the reactive stage in Prunus mume [46]. Besides the increase of auxin content, auxin-signaling pathways were also enriched in the vascular cambium during the transition from dormancy to active growth in Chinese fir [47]. Here, two DEGs encoding indole-3-acetic acid-amido synthetase GH3 family members, PsGH3s, were significantly up-regulated after GA3 and GA4 treatments, which was consistent with the previous results in Prunus mume [46]. The contents of endogenous IAA significantly increased only by exogenous GA3 application, which was inconsistent with the expression level of PsGH3 (Fig. 5). We suspected that GA4 could also enhance the increasement of endogenous IAA, but later than GA3 feeding in chronological order, which partially contributed to the morphological effect variation of different GAs.

Activation of carbohydrate metabolism involves in GAs-induced dormancy release

Carbohydrates not only act as an energy source but also as precursors for metabolic processes, and some of hexoses as molecular signals regulating many vital steps. Dormant buds of perennial woody plants seem to be in an inactive state with limited metabolic activities, and the reconfiguration of carbohydrate metabolism is related to the transition process from dormancy to regrowth [26, 46]. Our recent study indicated that carbohydrate metabolism, particularly the activation of PPP pathway, might play an important role in chilling-induced dormancy release of tree peony [33]. In this study, the significant degradation of starch and dramatic changes in soluble sugars emerged during GA-induced dormancy release in tree peony. Accompanied by the decreasing of starch level, the contents of soluble glucose and fructose showed up-regulated pattern after GA3 and GA4 applications, which were consistent with the expression levels of key gene involved in starch metabolism, sucrose transporter and glucose biosynthesis, including PsBMY, PsAMY, PsTPS and PsBGL (Fig. 6 and Table 3). Therefore, vast carbohydrates undoubtedly served as a strong energy sink for bud burst and subsequent growth after exogenous GAs applications in tree peony. These results agreed with those of leafy spurge, grape and Prunus mume crown buds during the dormancy-activity transition process [46, 48, 49]. Together, that exogenous GAs accelerated bud dormancy release and the following growth was associated with the induction of soluble sugar increase to satisfy the energy and substance requirement for bud regrowth.

The crosstalk between sugar and hormones including GA, ABA and IAA had been reported. The α-amylases act as key enzymes in starch degradation to generate soluble sugars, which play critical roles in the dormancy process and are tightly regulated by GA. For example, during the bud dormancy process in leafy spurge, GA promoted the synthesis and increased activity levels of α-amylases, suggesting GA might involve in the breakdown of starch. In grape, exogenous GA3 induced the expression of Vva-AMY3 and Vva-AMY4 during the GA3-induced dormancy release of grapevine bud [49]. In-depth, GA promoted the production of α-amylases mediating by GAMYB [50], while the induction of α-amylase triggered by GA was repressed by sugars, and this repression was transduced by a pathway independent of gibberellin signaling or that repression occurred at the level of GAMYB translation [51]. In general, ABA repressed the production of α-amylases by an ABA-induced protein kinase (PKABA1) and two ABA-induced WRKY [52, 53], and GA induction effect of the α-amylase gene was repressed by ectopic expression of wild-type ABA responsive element binding factor (ABF) 1 and 2 [54]. In this study, the down-regulation ABA might result in the increasing expression of α-amylase after exogenous GAs appication, which indicated the important role of interaction between sugar and hormones during GA-induced dormancy release in tree peony.

On the other hand, sugars also function as signaling molecules in plant growth regulations and stress responses. Strong evidences show that sugar signaling molecules, including sucrose, glucose, fructose and trehalose-6-phosphate, regulate development either directly or through interactions with other signaling pathways such as hormone-mediated processes [55]. For example, Glucose increase can activate target of rapamycin (TOR) kinase signaling and modulate indole-3-acetic acid (IAA) biosynthesis and transport, which will activate the meristem, promote cell division and expansion, and hence promote growth and development [56]. Trehalose-6-phosphate also plays crucial signaling roles in stimulating plant growth and development through repression of SnRK1 [57]. Our results revealed that GAs application promoted the degradation of sucrose and thus produced more glucose, fructose and trehalose-6-phosphate, accompanied by higher levels of IAA, and cell division was accelerated according to the expressions of related genes and the morphological changes (Fig. 1, 5 and 6). The results indicated that sucrose, glucose, fructose and trehalose-6-phosphate might function as sugar signaling molecules in GAs induced dormancy release, and GAs application might boost sugar signaling transduction.

As well known that photosynthesis ensures the continuous energy supply in the higher plant, which is the prerequisite for further growth and development. Our result indicated that exogenous GAs treatments increased the photosynthetic rate, while stomatal conductance and transpiration rate significantly enhanced only by GA3 in leaves, which well explained the morphological divergences after GAs treatments, similar to the effect of GA3 in the tree peony bracts [58]. FNR proteins existing in photosynthetic tissues is mainly responsible for catalyzing the final reaction of linear electron transport in photosynthesis, catalyzing electron from reduced ferredoxin (Fd) to NADP+, and the NADPH is mainly used for CO2 fixation and other metabolic processes of chloroplasts in the Calvin cycle [59]. In our study, the transcript of PsFNR significantly increased after GAs treatments, and its increased expression level by GA3 was more than by GA4, which suggested that GA3 promoting the growth of shoots better than GA4 might because that GA3 feeding accelerated photosynthesis to generate more energy substances.

GAs treatments accelerate cell cycle and promote the bud growth

Both GA3 and GA4 significantly promoted dormancy release and branch elongation. The increase of cell number is a prerequisite for the elongation of branches, which means more cells are in the Mitosis phase, and chromatin condensing is an important characteristic of Mitotic phase cell. The condensing complex promotes chromosome condensation during mitosis, and Structural Maintenance of Chromosomes 4 (SMC4) and SMC2 are the core components [60]. The abundance of SMC4 fluctuated with cell cycle and peaked at Mitotic stage and fall at interphase in budding yeast. In our study, PsSMC4 significantly increased after GAs treatments, which implied that more cell entered into Mitotic phase after GAs applications. The cyclin A/cdc2-kinase regulates the transition from G2 to Mitotic phase, and cyclin A (CYCA) is an important regulatory subunit. Here, the abundance of PsCYCA significantly increased by GA3, which might promote mores cells to enter Mitotic phase after GA3 treatment, which will accelerate bud burst and shoot growth. Whereas there was no difference between GA4 and mock. Additionally, the heterohexameric minichromosome maintenance complex (MCM2-7) serves as the central DNA replicative helicase in eukaryotes [61]. Two PsMCM2 and one PsMCM3 significantly decreased after GA3 treatment but increased by GA4 feedings, which were mainly consistent with the expression pattern of PsCYCA (Fig. 6), suggesting the condensed structure of chromosome in Mitotic phase after GA3 treatment.

The putative reasons for GA3 superior to GA4 promoting bud dormancy release

There are four major endogenous bioactive GAs, i.e. GA1, GA3, GA4 and GA7 in most plant species. In practice, exogenous GA3 is widely applied in dormancy breaking, such as in Paeonia lactiflora [17], Rhododendron simsii Planch [15], Camellia sinensis [62], Prunus armeniaca L. [63], etc. In some cases, GA4 is applied and more effectively to induce dormancy release and bud burst in hybrid aspen and Japanese apricot [35, 64]. In hybrid aspen, exogenous GA4 induces canonical bud burst and shoot elongation, but GA3 application only generates callus-like tissue. In contrast, GA3 and GA4 treatments exhibited similar morphological changes in tree peony, but GA3 induced faster bud burst, longer shoot and higher flowering rate in tree peony (Fig. 1), indicating GA3 was superior to GA4 in tree peony forcing culture.

On mechanism, Rinne et al. found GA3 and GA4 up-regulated different 1,3-β-glucanase members (glucan hudrolase family 17, GH17) in Populus that could hydrolyze 1,3-β-glucan deposited at pores and plasmodesmata in dormant stage. However, GA4 application could thoroughly reopen the signal conduit but GA3 was less efficient [35]. We found that both GA3 and GA4 induced same GH17 members but in diverse degrees in tree peony, which was in depth with the morphological variations (unpublished data). Therefore, plants do not response to phytohormone in a common manner.

Our results indicated that exogenous GA3 was superior to GA4 in accelerating dormancy release and subsequent growth (Fig. 1), despite more DEGs were screened in GA4 vs mock group. Collectively, we concluded several possible reasons. Firstly, the contents of endogenous GAs and the expression patterns of the related-genes including PsGA20ox and PsGID2 suggested that tree peony was more sensitive to exogenous GA3 than to GA4, and GA3 treatment enhanced more endogenous bioactive GA1 and GA3, which might function more effectively in tree peony buds. Secondly, PsBMY was only induced by GA3 treatment, and the higher expressions of PsAMY and PsTPS were up-regulated after GA3 feeding, which lead to the more effective hydrolysis of starch to offer more energy for break dormancy, and photosynthesis of leaves by GA3 application ensured the continuous energy supply. The quantification of starch and sugar further proved this speculation, and higher phtosythetic ability achieved after GA3 feeding accompanied by higher PsFNR expression. Thirdly, GA3 application resulted in higher IAA levels, another growth enhancer hormone. In addition, the interactions between sugar and hormones were also an important factor during GA-induced dormancy release, and GA3 induced more effectively sugar signaling transduction to stimulate cell division. The expression patterns of the cell cycle and DNA replication related genes, including PsCYCA, PsCYCD, PsSMC4, PsMCM2 and PsMCM3, suggested that exogenous GA3 feeding more effectively accelerated cell cycle to promote dormancy release and subsequent growth (Fig. 8). At last, we evaluated the expressions of several dormancy-associated genes, such as dormancy release-related FT, EBB1, EBB3 and CYCD, and dormancy-associated TFL and SVP. Higher expression levels were detected by GA3 treatment for TFL, FT, EBB1, EBB3 and CYCD, but lower for SVP, implying more efficiency of GA3.

Fig. 8
figure 8

Model of the putative mechanism about exogenous GA3 and GA4 on dormancy release in tree peony. Exogenous GA3 and GA4 promote dormancy release and subsequent growth by increasing endogenous bioactive GAs (GA1, GA3, GA4 and GA7) and IAA, decreasing ABA contents, enhancing the contents of soluble glucose and accelerating cell division, promoting the expression levels of known four dormancy release-related genes including PsFT, PsEBB1, PsEBB3, PsCYCD and one dormancy-related gene PsTFL, inhibiting that of PsSVP. Red represents the up-regulated gene and blue represents the down-regulated gene

Conclusions

High-quality Full-length transcriptome was obtained in dormant buds of tree peony, and was used to analyze the effects of GAs on bud dormancy and growth. Morphologically, exogenous GA3 and GA4 accelerated bud dormancy release, subsequent growth and flowering, and the intrinsic mechanism might be associated with the increasing endogenous bioactive GAs (GA1, GA3, GA4 and GA7) and IAA, decreasing ABA contents, enhancing the contents of soluble glucose, fructose and trehalose to provide sufficient energy and metabolic substance, and accelerating sugar signaling and cell division. The reason of GA3 superior to GA4 treatment might be that tree peony was more sensitive to GA3 than to GA4, and GA3 application might be more efficient to accelerate cell cycle and starch hydrolysis to provide more soluble sugar as energy and carbon skeleton for other metabolites synthesis along with stepped-up sugar signaling. All results would provide valuable data for the mechanism of dormancy release and the foundation for forcing culture improvement in tree peony.

Materials and methods

Plant materials and morphological observation

Tree peony (Paeonia suffruticosa) ‘Luhehong’ was formally identified by Heze Zhaolou Tree-peony Garden in 1968. It is a traditional cultivar and not protected by plant variety rights in China. Therefore, no specimen of this material is deposited in a publicly available herbarium. We have obtained the permission of Tree Peony Research Institute of Qingdao Agricultural University (Qingdao, China) to collect tree peony ‘Luhehong’. Four-year-old ‘Luhehong’ were potted and treated in cold conditions (0–4 °C, 24 h-dark) from Nov 10, 2018. The apical buds after 0, 7, 14, 21, 28 d chilling applications were collected and immediately frozen in liquid nitrogen, which were stored at -80 °C until use for PacBio full-length transcriptome sequencing, and three replicates (3 plants/replicate) per treatment were set, respectively. On Nov 17, 2018, 27 plants were transferred to a greenhouse (25 °C, 8-h-light/16-h-dark), and then GA3 (500 mgּּּ L−1) and GA4 (500 mgּ L−1) were applied on apical buds, respectively, the samples treated by sterile distilled water was used as mock. After treated for 48 h, the apical buds were harvested, immediately frozen in liquid nitrogen and stored at -80 °C until use for Illumina transcriptome sequencing, qPCR, sugar and starch evaluation, and endogenous phytohormone measurement, respectively. Three replicates (3 plants/replicate) per group were set.

Morphological observation was conducted on another 27 plants. Sprouting rate was observed every day till GA3 and GA4 applications for 16 d [18]. Branch height, the first leaf length and flowering rate were monitored and measured as described previously [5]. Treatment effects were analyzed using ANOVA followed by Duncan's multiple range tests at a significance level of 0.05, and the p-value of the comparison with the treated group and Mock was obtained after t-test by SPSS 13.0 for Windows (SPSS, USA).

Measurement of photosynthesis

The net photosynthetic rate (Pn), intercellular CO2 concentration (Ci), stomatic conductance (Gs), and transpiration rate (Tr) were evaluated 20 d after GAs applied, by a CIRAS-3 portable photosynthesis system (PP Systems, USA) under ambient CO2 concentrations at 20 °C, 1200 μmol (photon) m−2 s−1 light, and 80% relative humidity. The first leaf under flower was measured and 10 pieces of leaves per repeat.

Transcriptome sequencing and analysis

Ten buds were sampled per biological replicate to produce an independent pool. PacBio full-length sequencing and Illumina sequencing were both performed by OE Biotech Co. Ltd. (Shanghai, China). Total RNA was extracted using RNA extracting kit (Umagen, China) according to the manual, and the DNA contaminates were excluded using TURBO DNA-free Kit (Thermo, USA). RNA samples were checked on a 1.5% (w/v) agarose gel, and the integrities were evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The samples with RNA Integrity Number (RIN) ≥ 8.5, OD260/OD280: 2.0 ~ 2.2, OD260/OD230: 1.8 ~ 2.1, qubit/nano: 0.5 ~ 2.0 were subjected to the subsequent analysis.

PacBio full-length sequencing of mixed samples

In order to generate a representatively transcriptomic reference, 1 μg of total RNA extracted from each chilling-treated sample were equally pooled together and used for cDNA synthesis using SMARTer PCR cDNA Synthesis Kit (Clontech, USA). The library construction and PacBio sequencing were performed according to the official protocol as described by Pacific Biosciences (Pacific Biosciences, USA). Briefly, a total of 12 PCR cycles of amplification were performed using PrimeSTAR GXL DNA Polymerase (Clontech, USA). After purification with AMPure PB Beads, the cDNA products were then subjected to the construction of SMRTbell template libraries using the SMRTbell Template Prep Kit 1.0 (Pacific Biosciences, USA). Finally, SMRT cells were sequenced on a PacBio Sequel instrument using sequencing kit 2.1 with 10 h movie recordings. Sequencing reads were subjected to Circular Consensus Sequences (CCS) using the SMRT Analysis Software (https://www.pacb.com/products-and-services/analytical-software/devnet/).

Illumina sequencing of GA3, GAand mock group

The libraries were constructed using TruSeq Stranded mRNA LTSample Prep Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Then these libraries were sequenced on the Illumina sequencing platform (HiSeqTM 2500). The raw reads obtained by Illumina Hiseq 2500 sequencer were originally trimmed by removing the adapter sequences, and then the low-quality sequences with quality scores of less than 20 and the short reads with lengths of less than 10 bp were also removed. Finally, the high-quality clean reads were obtained for transcriptome de novo assembly using Trinity (version 20,140,717) [65], and redundant sequences were filtered using TGICL [66].

Unigenes annotation and transcriptome analysis

PacBio reads were classified into full-length and non-full-length sequences, and then reads were corrected using LoRDEC. The redundant sequences were removed using CD-HIT. The merged unigenes were annotated based on sequence similarity using BLASTx and BLASTn algorithms against the NCBI non-redundant database, and only significant BLAST results (E-value < 1e-5) were considered to be annotated, including protein functional annotation, COG functional annotation, GO functional annotation, and pathway annotation. Then, the best alignment results were used to determine the sequence direction and protein-coding-region prediction (CDS) of the unigenes. The functional category was assigned to each gene product according to the gene ontology (GO) terms (www.geneontology.org) that were analyzed using Blast2GO software (http://www.blast2go.com/b2ghome). Pathway assignment was performed according to the Arabidopsis database classification in the KEGG database (Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg). Based on the KEGG database, the complex biological behavior of the genes was analyzed through pathway annotation. Up to now, there is still no high-quality genome data in tree peony, and the clean reads were mapped to our PacBio full-length transcriptome using the short oligonucleotide analysis package SOAPaligner/soap2 (version 2.21).

Identification of differentially expressed genes

The expression levels of unigenes were generally evaluated by the quantity of fragments per kilobase of exon model per million mapped fragments (FPKM). False discovery rate (FDR) was used to correct for p-value thresholds in multiple testing [67]. After comparing each other using EDGER software (empirical analysis of digital gene expression data in R) [68], the unigenes with |log2Ratio|≥ 1 and p < 0.05 were identified as differentially expressed genes (DEGs). Hierarchical cluster analysis of DEGs was performed to explore gene expression patterns. Then, DEGs were subjected to GO functional enrichment and KEGG pathway enrichment analyses. GO terms and KEGG pathways fulfilling the criterion of a Bonferroni-corrected p-value ≤ 0.05 were defined as significantly enriched in DEGs.

Real-time quantitative RT-PCR

In order to verify the sequencing quality and the expression patterns of DEGs, the associated primer pairs were designed based on their partial cDNA sequences (Additional file 8). Two mg total RNA was used to synthesize the first-strand cDNA using the PrimerScript™ RT reagent Kit according to the manufacturer’s instructions (TaKaRa, Dalian, China), and then the cDNAs were three times diluted as qPCR template. PCR reaction was performed in 20 mL reaction mixture (2 × SYBR Green Master mix of 10 mL, 10 mmolּ L−1 each primer of 1 mL, 3 × diluted cDNA template 2 mL and 6 mL ddH2O) using SYBR® Premix Ex Taq™ II Kit (TaKaRa, Dalian, China). The PCR reactions were run in a LightCycler® 480 II (Roche, USA), and the program was 95 °C for 2 min, and followed by 40 cycles of 95 °C for 10 s, 55 °C for 30 s and 72 °C for 30 s. All reactions were conducted in triplicates for each sample, and the relative expression levels were calculated using the method of 2Ct [69], PsActin was used as the reference gene. Significance was tested using SPSS 13.0 for Windows (SPSS, USA).

Measurement of endogenous phytohormone

The endogenous phytohormone contents of three groups (mock, GA3 and GA4) were determined by Wuhan Greensword Creation Technology Company (http://www.greenswordcreation.com) based on LC–MS/MS analysis according to a previously reported method [70].

Determination of the sugar and starch content

The freeze-dried samples (mock, GA3 and GA4) were crushed using a mixer mill (MM 400, Retsch) with a zirconia bead for 1.5 min at 30 Hz. Twenty mg of powder was diluted to a 500 μL methanol: isopropanol: water (3:3:2, V/V/V), vortexed for 3 min and ultrasound for 30 min. The extract was centrifuged at 14 000 rpm for 3 min (4 °C). The supernatants were collected for 50 μL and evaporated under nitrogen gas stream, and add internal standard. The extract was evaporated under nitrogen gas stream and transferred to the lyophilizer for freeze-drying. The residue was used for further derivatization as follows: the sample of small molecular carbohydrates was mixed with 100 μL solution of methoxyamine hydrochloride in pyridine (15 mg·mL−1). The mixture was incubated at 37 °C for 2 h. Then 100 μL of BSTFA was added into the mixture and kept at 37 °C for 30 min after vortex-mixing. Dilute the mixture to the appropriate concentration with n-hexane. The mixture was analyzed by GC–MS/MS (Agilent 7890B-7000D) according to Sun et al. [71]. The starch content was measured according to the previous method [72, 73]. Three replicates of each assay were performed.

Availability of data and materials

The transcriptome data of the overlapping 849 DEGs between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock have been appended as Additional file 5. The full-length transcriptome data have been submitted to the public repository (SRA), and the BioProject ID is PRJNA720276. The RNAseq data used during the current study available from the corresponding author (Gai Shupeng: spgai@qau.edu.cn) on reasonable request.

Abbreviations

ANOVA:

Analysis of variance

cDNA:

Complementary DNA

COG:

Clusters of Orthologous Groups

GC–MS/MS:

Gas chromatography-mass spectrometry

GO:

Gene Ontology

IAA:

Indole-3-acetic acid

LC–MS/MS:

Liquid chromatography electrospray ionisation tandem mass spectrometry

PCA:

Principal component analysis

References

  1. Lang GA, Early JD, Martin GC, Darnell RL. Endo-, para-, and ecodormancy: physiological terminology and classification for dormancy research. Hortscience. 1987;22:371–7.

    Article  Google Scholar 

  2. Heide OM, Prestrud AK. Low temperature, but not photoperiod, controls growth cessation and dormancy induction and release in apple and pear. Tree Physiol. 2005;1(25):109–14.

    Article  Google Scholar 

  3. Carvalho RIND, Biasi LA, Zanette F, Rendoke JC, Santos JM, Pereira GP. Dormancy of “Imperial Gala” apple and “Hosui” pear tree buds in a region of low chill occurrence. Acta Sci Agron. 2014;36(4):429–34.

  4. Zhang T, Si F, Zhang Y, Gai S. Effect of exogenous gibberellin on DNA methylation level and expression of related enzyme genes in tree peony floral buds. Scientia Agricultura Sinica. 2018;51(18):3561–9.

    Google Scholar 

  5. Zhang Y, Si F, Wang Y, Liu C, Zhang T, Yuan Y, Gai S. Application of 5-azacytidine induces DNA hypomethylation and accelerates dormancy release in buds of tree peony. Plant Physiol Bioch. 2020;147:91–100.

    Article  CAS  Google Scholar 

  6. Hoeberichts FA, Povero G, Ibañez M, Strijker A, Pezzolato D, Mills R, Piaggesi A. Next Generation Sequencing to characterise the breaking of bud dormancy using a natural biostimulant in kiwifruit (Actinidia deliciosa). Sci Hortic-Amsterdam. 2017;225:252–63.

    Article  CAS  Google Scholar 

  7. Michailidis M, Karagiannis E, Tanou G, Sarrou E, Adamakis ID, Karamanoli K, Martens S, Molassiotis A. Metabolic mechanisms underpinning vegetative bud dormancy release and shoot development in sweet cherry. Environ Exp Bot. 2018;155:1–11.

    Article  CAS  Google Scholar 

  8. Wu R, Tomes S, Karunairetnam S, Tustin SD, Hellens RP, Allan AC, Macknight RC, Varkonyigasic E. SVP-like MADS box genes control dormancy and budbreak in Apple. Front Plant Sci. 2017;8(55):477.

    PubMed  PubMed Central  Google Scholar 

  9. Yang Q, Yang B, Li J, Wang Y, Tao R, Yang F, Wu X, Yan X, Ahmad M, Shen J, Bai S, Teng Y. ABA-responsive ABRE-BINDING FACTOR3 activates DAM3 expression to promote bud dormancy in Asian pear. Plant Cell Environ. 2020;43(6):1360–75.

    Article  CAS  PubMed  Google Scholar 

  10. Zheng C, Acheampong AK, Shi Z, Mugzech A, Halalybasha T, Shaya F, Sun Y, Colova V, Mosquna A, Ophir R. Abscisic acid catabolism enhances dormancy release of grapevine buds. Plant Cell Environ. 2018;41(10):2490–503.

    Article  CAS  PubMed  Google Scholar 

  11. Huang X, Xue T, Dai S, Gai S, Zheng C, Zheng G. Genes associated with the release of dormant buds in tree peonies (Paeonia suffruticosa). Acta Physiol Plant. 2008;30(6):797–806.

    Article  CAS  Google Scholar 

  12. Cooke JEK, Eriksson ME, Junttila O. The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms. Plant Cell Environ. 2012;35(10):1707–28.

    Article  CAS  PubMed  Google Scholar 

  13. Olsen JE. Light and temperature sensing and signaling in induction of bud dormancy in woody plants. Plant Mol Biol. 2010;73(s1–2):37–47.

    Article  CAS  PubMed  Google Scholar 

  14. Yue C, Cao H, Hao X, Zeng J, Qian W, Guo Y, Ye N, Yang Y, Wang X. Differential expression of gibberellin- and abscisic acid-related genes implies their roles in the bud activity-dormancy transition of tea plants. Plant Cell Rep. 2018;37(3):425–41.

    Article  CAS  PubMed  Google Scholar 

  15. Nell TA, Larson RA. The influence of foliar applications of GA3 GA4,7 and PBA on breaking flower bud dormancy on azalea CVS redwing and dogwood. J Horticult Sci. 1974;49(4):323–8.

    Article  CAS  Google Scholar 

  16. Zhuang W, Gao Z, Wang L, Zhong W, Ni Z, Zhang Z. Comparative proteomic and transcriptomic approaches to address the active role of GA4 in Japanese apricot flower bud dormancy release. J Exp Bot. 2013;64(16):4953–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Evans MR, Anderson NO, Wilkins HF. Temperature and GA3 effects on emergence and flowering of potted Paeonia lactiflora. Hortsci Publ Am Soc Horticult Sci. 1990;25(8):923–4.

    CAS  Google Scholar 

  18. Gai S, Zhang Y, Liu C, Zhang Y, Zheng G. Transcript profiling of Paoenia ostii during artificial chilling induced dormancy release identifies activation of GA pathway and carbohydrate metabolism. Plos One. 2013;8(2):e55297.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Middleton AM, Úbeda-Tomás S, Griffiths J, Holman T, Hedden P, Thomas SG, Phillips AL, Holdsworth MJ, Bennett MJ, King JR, Owen MR. Mathematical modeling elucidates the role of transcriptional feedback in gibberellin signaling. Proc Natl Acad Sci USA. 2012;109(19):7571–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Hedden P. Gibberellin metabolism and its regulation. J Plant Growth Regul. 2002;20(4):317–8.

    Article  CAS  Google Scholar 

  21. Lin L, Huo X, Wen L, Gao Z, Muhammad KUR. Isolation and role of PmRGL2 in GA-mediated floral bud dormancy release in Japanese apricot (Prunus mume Siebold et Zucc.). Front Plant Sci. 2018;9:27.

    Article  CAS  Google Scholar 

  22. Azeez A, Zhao YC, Singh RK, Yordano VYS, Dash M, Miskolczi P, Stojkovi K, Strauss SH, Bhalerao RP, Busov VB. EARLY BUD-BREAK 1 and EARLY BUD-BREAK 3 control resumption of poplar growth after winter dormancy. Nat Commun. 2021;12(1):1123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Singh RK, Miskolczi P, Maurya JP, Bhalerao RP. A tree ortholog of SHORT VEGETATIVE PHASE floral repressor mediates photoperiodic control of bud dormancy. Curr Biol. 2019;1(29):128–33.

    Article  CAS  Google Scholar 

  24. Tylewicz S, Petterle A, Marttila S, Miskolczi P, Azeez A, Singh RK, Immanen J, Mahler N, Hvidsten TR, Eklund DM. Photoperiodic control of seasonal growth is mediated by ABA acting on cell-cell communication. Science. 2018;360(6385):212–5.

    Article  CAS  PubMed  Google Scholar 

  25. Yang Q, Niu Q, Tang Y, Ma Y, Yan X, Li J, Tian J, Bai S, Teng Y. PpyGAST1 is potentially involved in bud dormancy release by integrating the GA biosynthesis and ABA signaling in “Suli” pear (Pyrus pyrifolia White Pear Group). Environ Exp Bot. 2019;162:302–12.

  26. Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, Boerjan W, Rohde A. A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007;19(8):2370–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Aleman F, Yazaki J, Lee M, Takahashi Y, Kim AY, Li Z, Kinoshita T, Ecker JR, Schroeder JI. An ABA-increased interaction of the PYL6 ABA receptor with MYC2 transcription factor: A putative link of ABA and JA signaling. Sci Rep-Uk. 2016;6(1):28941.

    Article  CAS  Google Scholar 

  28. Singh RK, Maurya JP, Azeez A, Miskolczi P, Tylewicz S, Stojkovič K, Delhomme N, Busov V, Bhalerao RP. A genetic network mediating the control of bud break in hybrid aspen. Nat Commun. 2018;9(1):4173.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Bonicel A, Verçosa DMRN. Variation of starch and soluble sugars in selected sections of poplar buds during dormancy and post-dormancy. Plant Physiol Bioch. 1990;28(5):577–86.

    CAS  Google Scholar 

  30. Ito A, Sakamoto D, Moriguchi T. Carbohydrate metabolism and its possible roles in endodormancy transition in Japanese pear. Sci Hortic-Amsterdam. 2012;144(1):187–94.

    Article  CAS  Google Scholar 

  31. Jin F, Li J, Ding Q, Wang Q, He X. Proteomic analysis provides insights into changes in the central metabolism of the cambium during dormancy release in poplar. J Plant Physiol. 2017;208:26–39.

    Article  CAS  PubMed  Google Scholar 

  32. Wang W, Su X, Tian Z, Liu Y, Zhou Y, He M. Transcriptome profiling provides insights into dormancy release during cold storage of Lilium pumilum. BMC Genomics. 2018;19(S2):196.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Zhang Y, Yu D, Liu C, Gai S. Dynamic of carbohydrate metabolism and the related genes highlights PPP pathway activation during chilling induced bud dormancy release in tree peony (Paeonia suffruticosa). Sci Hortic-Amsterdam. 2018;242:36–43.

    Article  CAS  Google Scholar 

  34. Noriega X, Pérez FJ. Cell cycle genes are activated earlier than respiratory genes during release of grapevine buds from endodormancy. Plant Signal Behavior. 2017;12(10):e1321189.

    Article  CAS  Google Scholar 

  35. Rinne PLH, Schoot CVD. Chilling of dormant buds hyperinduces FLOWERING LOCUS T and recruits GA-Inducible 1,3-β-glucanases to reopen signal conduits and release dormancy in Populus. Plant Cell. 2011;23(1):130.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Rozi M, Chieh-Ting W, Cathleen M, Olga S, Dye SJ, Puzey JR, Elizabeth E, Xiaoyan S, Richard M, Strauss SH. Populus CEN/TFL1 regulates first onset of flowering, axillary meristem identity and dormancy release in Populus. Plant J. 2010;62(4):674–88.

    Article  CAS  Google Scholar 

  37. Yordanov YS, Ma C, Strauss SH, Busov VB. EARLY BUD-BREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees. Proc Natl Acad Sci USA. 2014;111(27):10001–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Mornya PMP, Cheng F. Seasonal changes in endogenous hormone and sugar contents during bud dormancy in tree peony. J Appl Hortic. 2013;15(3):159–65.

    Article  Google Scholar 

  39. Guan Y, Xue J, Xue Y, Yang R, Wang S, Zhang X. Effect of exogenous GA3 on flowering quality, endogenous hormones, and hormone- and flowering-associated gene expression in forcing-cultured tree peony (Paeonia suffruticosa). J Integr Agr. 2019;018(006):1295–311.

    Article  CAS  Google Scholar 

  40. Schrader J, Moyle R, Bhalerao R, Hertzberg M, Lundeberg J, Nilsson P, Bhalerao RP. Cambial meristem dormancy in trees involves extensive remodelling of the transcriptome: Cambial meristem dormancy. Plant J. 2004;40(2):173–87.

    Article  CAS  PubMed  Google Scholar 

  41. Zhang Y, Zhang W, Li Y, Liu C, Zheng G, Gai S. The study of PsGA20ox gene participating in endo-dormancy release of flower buds by chilling treatment in tree peony. Acta Agriculturae Boreali-Sinica. 2014;003:22–6.

    Google Scholar 

  42. Li M, An F, Li W, Ma M, Feng Y, Zhang X, Guo H. DELLA proteins interact with FLC to repress flowering transition. J Integr Plant Biol. 2016;7(58):642–55.

    Article  CAS  Google Scholar 

  43. O’Neill DP, Davidson SE, Clarke VC, Yamauchi Y, Yamaguchi S, Kamiya Y, Reid JB, Ross JJ. Regulation of the gibberellin pathway by auxin and DELLA proteins. Planta. 2010;232(5):1141–9.

  44. Berruti A, Christiaens A, Scariot V, Van LMC. Changes in ABA levels in vegetative and flower buds during dormancy in Camellia. Acta Hort. 2012;953(953):247–54.

    Article  Google Scholar 

  45. Xu R, Niimi Y, Han D. Changes in endogenous abscisic acid and soluble sugars levels during dormancy-release in bulbs of Lilium rubellum. Sci Hortic-Amsterdam. 2006;111(1):68–72.

    Article  CAS  Google Scholar 

  46. Zhang Z, Zhuo XK, Zhao K, Zheng T, Zhang Q. Transcriptome profiles reveal the crucial roles of hormone and sugar in the bud dormancy of Prunus mume. Sci Rep-Uk. 2018;8(1):5090.

    Article  CAS  Google Scholar 

  47. Qiu Z, Wan L, Chen T, Wan Y, He X, Lu S, Wang Y, Lin J. The regulation of cambial activity in Chinese fir (Cunninghamia lanceolata) involves extensive transcriptome remodeling. New Phytol. 2013;199(3):708–19.

    Article  CAS  PubMed  Google Scholar 

  48. Anderson JV, Gesch RW, Jia Y, Chao WS, Horvath DP. Seasonal shifts in dormancy status, carbohydrate metabolism, and related gene expression in crown buds of leafy spurge. Plant, Cell Environ. 2005;28(12):1567–78.

    Article  CAS  Google Scholar 

  49. Rubio S, Donoso A, Pérez FJ. The dormancy-breaking stimuli “chilling, hypoxia and cyanamide exposure” up-regulate the expression of α-amylase genes in grapevine buds. J Plant Physiol. 2014;171(6):373–81.

  50. Gubler F, Raventos D, Keys M, Watts R, Jacobsen JV. Target genes and regulatory domains of the GAMYB transcription activator in cereal aleurone. Plant J. 1999;17(1):1–9.

    Article  CAS  PubMed  Google Scholar 

  51. Loreti E, Matsukura CA, Gubler F, Alpi A, Perata P. Glucose repression of α-amylase in barley embryos is independent of GAMYB transcription. Plant Mol Biol. 2000;44(1):85–90.

    Article  CAS  PubMed  Google Scholar 

  52. Lu C, Ho TD, Ho S, Yu S. Three novel MYB proteins with one DNA binding repeat mediate sugar and hormone regulation of α-Amylase gene expression. Plant Cell. 2002;14(8):1963–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Xie Z, Zhang ZL, Zou X, Yang G, Komatsu S. Interactions of two abscisic-acid induced WRKY genes in repressing gibberellin signaling in aleurone cells. Plant J. 2006;46(2):231–42.

    Article  CAS  PubMed  Google Scholar 

  54. Schoonheim PJ, Pereira DDAC, Boer AHD. Dual role for 14-3-3 proteins and ABF transcription factors in gibberellic acid and abscisic acid signalling in barley (Hordeum vulgare) aleurone cells. Plant Cell Environ. 2009;32(5):439–47.

    Article  CAS  PubMed  Google Scholar 

  55. Ruan YL. Sucrose metabolism: gateway to diverse carbon use and sugar signaling. Annu Rev Plant Biol. 2014;65:33–67.

    Article  CAS  PubMed  Google Scholar 

  56. Xiong Y, Sheen J. Novel links in the plant TOR kinase signaling network. Curr Opin Plant Biol. 2015;28:83–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Yadav UP, Ivakov A, Feil R, Duan GY, Walther D, Giavalisco P, Piques M, Carillo P, Hubberten HM, Stitt M, lunn JE. The sucrose-trehalose 6-phosphate (Tre6P) nexus: specificity and mechanisms of sucrose signalling by Tre6P. J Exp Bot. 2014;65:1051–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Liu Z, Shi Y, Xue Y, Wang X, Huang Z, Xue J, Zhang X. Non-structural carbohydrates coordinate tree peony flowering both as energy substrates and as sugar signaling triggers, with the bracts playing an essential role. Plant Physiol Bioch. 2021;159:80–8.

    Article  CAS  Google Scholar 

  59. Carrillo N, Ceccarelli EA. Open questions in ferredoxin-NADP+ reductase catalytic mechanism. FEBS J. 2010;270(9):1900–15.

    Google Scholar 

  60. Wang H, Liu Y, Yuan J, Zhang J, Han F. The condensin subunits SMC2 and SMC4 interact for correct condensation and segregation of mitotic maize chromosomes. Plant J. 2020;102(3):467–79.

    Article  CAS  PubMed  Google Scholar 

  61. Lyubimov AY, Costa A, Bleichert F, Botchan MR, Berger JM. ATP-dependent conformational dynamics underlie the functional asymmetry of the replicative helicase from a minimalist eukaryote. Proc Natl Acad Sci USA. 2012;109(30):11999–2004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Thirugnanasambantham K, Prabu G, Mandal A. Synergistic effect of cytokinin and gibberellins stimulates release of dormancy in tea (Camellia sinensis (L.) O. Kuntze) bud. Physiol Mol Biol Pla. 2020;26(5):1035–45.

    Article  CAS  Google Scholar 

  63. Kojori ZK, Rezaei M, Sarkhosh A, Gharangik S. The effect of bud-scale removal and gibberellin (GA3) on dormancy break of apricot (P. Armeniaca L.) vegetative buds. J Appl Horticult. 2018;20(1):52–6.

    Google Scholar 

  64. Zhuang W, Gao Z, Wen L, Huo X, Cai B, Zhang Z. Metabolic changes upon flower bud break in Japanese apricot are enhanced by exogenous GA4. Hortic Res-England. 2015;2:15046.

    Article  CAS  Google Scholar 

  65. Grabherr MG, Haas BJ, Moran Y, Levin JZ, Thompson DA, Ido A, Xian A, Lin F, Raktima R, Qiandong Z. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19(5):651–2.

    Article  CAS  PubMed  Google Scholar 

  67. Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29(4):1165–88.

    Article  Google Scholar 

  68. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;10(11):e576.

    Google Scholar 

  69. 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.

    Article  CAS  PubMed  Google Scholar 

  70. Chen ML, Fu XM, Liu JQ, Ye TT, Hou SY, Huang YQ, Yuan BF, Wu Y, Feng YQ. Highly sensitive and quantitative profiling of acidic phytohormones using derivatization approach coupled with nano-LC-ESI-Q-TOF-MS analysis. J Chromatography B. 2012;905:67–74.

    Article  CAS  Google Scholar 

  71. Sun S, Wang H, Xie J, Su Y. Simultaneous determination of rhamnose, xylitol, arabitol, fructose, glucose, inositol, sucrose, maltose in jujube (Zizyphus jujube Mill.) extract: comparison of HPLC-ELSD, LC-ESI-MS/MS and GC-MS. BMC Chemistry. 2016;10(1):1–9.

    CAS  Google Scholar 

  72. Hansen J, Moler I. Percolation of starch and soluble carbohydrates from plant tissue for quantitative determination with anthrone. Anal Biochem. 1975;68(1):87–94.

    Article  CAS  PubMed  Google Scholar 

  73. Zhang T, Yuan Y, Zhan Y, Cao X, Liu C, Zhang Y, Gai S. Metabolomics analysis reveals Embden Meyerhof Parnas pathway activation and flavonoids accumulation during dormancy transition in tree peony. BMC Plant Biol. 2020;20(1):484.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors greatly acknowledge Xin Hua at the college of Life Science, Qingdao Agricultural University, for providing useful suggestions about morphological observation and measurement of photosynthesis.

Funding

This work was supported by grants from the National Key R&D Program of China (2018YFD1000403) and National Natural Science Foundation of China (31872145 and 31972452). The funding bodies had no role in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

GSP and ZYX conceived and designed the experimental plan. YYC, LZJ, LF, ZT and LCY conducted the experiments. ZYX and YYC analyzed the data. ZYX and GSP prepared and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gai Shupeng.

Ethics declarations

Ethics approval and consent to participate

We have obtained the permission of Tree Peony Research Institute of Qingdao Agricultural University to collect tree peony ‘Luhehong’. The collection and usage of plant specimens in current study complied with relevant institutional, national, and international guidelines and legislation. Ethical approval was not applicable for this study.

Consent for publication

No applicable.

Competing interests

The authors declare there are no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

The statistical data of unigenes annotated in public database.

Additional file 2:

The number distribution of the unigenes encoding transcription factor (TF).

Additional file 3:

PCA analysis of the samples.

Additional file 4:

GO analysis of DEGs (p value<0.05 & |log2FC|>1) after GA3 and GA4 applications.

Additional file 5:

The transcriptomic data of the overlapping 849 DEGs between GA3 vs GA4 and GA3 vs mock, GA3 vs GA4 and GA4 vs mock, respectively.

Additional file 6:

Based on Neighbor-Joining (NJ) model by MEGA 6.0, the phylogenetic trees of the related DEGs enriched in the KEGG pathways about plant hormone signal transduction, starch and sucrose metabolism, cell division and DNA replication. Genetic distance was calculated based on nucleotide difference (p-distance) with complete deletion of gaps. Scale bar = 0.05. The number at each node indicates the percentage of bootstrapping of 1000 replications.

Additional file 7:

The correlation between the DEGs fold change of expression levels by qRT-PCR and the corresponding fold change with RNAseq.

Additional file 8:

Primer sequences used for qPCR analysis in this study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yuxi, Z., Yanchao, Y., Zejun, L. et al. GA3 is superior to GA4 in promoting bud endodormancy release in tree peony (Paeonia suffruticosa) and their potential working mechanism. BMC Plant Biol 21, 323 (2021). https://doi.org/10.1186/s12870-021-03106-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-021-03106-2

Keywords