Skip to main content

Advertisement

Comparative transcriptome analysis revealing the potential mechanism of seed germination stimulated by exogenous gibberellin in Fraxinus hupehensis

Article metrics

Abstract

Background

Fraxinus hupehensis is an endangered tree species that is endemic to in China; the species has very high commercial value because of its intricate shape and potential to improve and protect the environment. Its seeds show very low germination rates in natural conditions. Preliminary experiments indicated that gibberellin (GA3) effectively stimulated the seed germination of F. hupehensis. However, little is known about the physiological and molecular mechanisms underlying the effect of GA3 on F. hupehensis seed germination.

Results

We compared dormant seeds (CK group) and germinated seeds after treatment with water (W group) and GA3 (G group) in terms of seed vigor and several other physiological indicators related to germination, hormone content, and transcriptomics. Results showed that GA3 treatment increases seed vigor, energy requirements, and trans-Zetain (ZT) and GA3 contents but decreases sugar and abscisic acid (ABA) contents. A total of 116,932 unigenes were obtained from F. hupehensis transcriptome. RNA-seq analysis identified 31,856, 33,188 and 2056 differentially expressed genes (DEGs) between the W and CK groups, the G and CK groups, and the G and W groups, respectively. Up-regulation of eight selected DEGs of the glycolytic pathway accelerated the oxidative decomposition of sugar to release energy for germination. Up-regulated genes involved in ZT (two genes) and GA3 (one gene) biosynthesis, ABA degradation pathway (one gene), and ABA signal transduction (two genes) may contribute to seed germination. Two down-regulated genes associated with GA3 signal transduction were also observed in the G group. GA3-regulated genes may alter hormone levels to facilitate germination. Candidate transcription factors played important roles in GA3-promoted F. hupehensis seed germination, and Quantitative Real-time PCR (qRT-PCR) analysis verified the expression patterns of these genes.

Conclusion

Exogenous GA3 increased the germination rate, vigor, and water absorption rate of F. hupehensis seeds. Our results provide novel insights into the transcriptional regulation mechanism of effect of exogenous GA3 on F. hupehensis seed germination. The transcriptome data generated in this study may be used for further molecular research on this unique species.

Background

Fraxinus hupehensis Chu, Shang et Su. is a woody plant of the Oleaceae family that is officially listed as a national rare and endangered tree species in China [1,2,3]. The species has high commercial value due to its slow growth, interlaced roots, intricate tree shape, and easy to shape. Furthermore, F. hupehensis has great potential to improve and protect the environment on account of its strong adaptability to arid and cold climates and resistance to diseases and insects [4, 5]. Previous studies have focused on the resource investigation [2], seedling technology [6], tissue culture [7], chemical composition, and pharmacological value of F. hupehensis [8]. In general, the fruit of F. hupehensis produces seeds after ripening and falling off in the spring of the second year. Under conventional sowing conditions, F. hupehensis seeds take more than 1 year of dormancy to germinate [9]. Considering the rapid development of the bonsai industry of F. hupehensis in China, resulting in serious damage to wild resources and a sharp decline in population. Therefore, its resources should be protected, and its reproductive and growth cycles should be accelerated through artificial technologies.

Seed germination refers to the resumption of embryo growth after water absorption and seed expansion until the radicle breaks through the endosperm and seed coat before germination. When stimulated by germination conditions (e.g., light, low temperature, water, and hormones), immature seeds ripen to form mature seeds and undergo germination [10]. Woody-plant seeds are naturally difficult to germinate, and dormancy is disrupted by external stimulation to promote germination. Exogenous hormones could promote the germination of dormant seeds. In particular, gibberellin (GA3) could break seed dormancy and play an endogenous signaling role during seed germination [11, 12]. In recent years, great advances have been achieved in the mechanism of GA3 on seed germination. GA3 plays two important roles in plant seed germination: (1) it overcomes the mechanical constraints imposed by seed mulch by weakening the tissues around the radicle [13], and (2) it increases the growth potential of embryos [14]. GA3 induces the degradation of the plant growth inhibitor DELLA protein by binding to its receptor, thereby promoting plant germination [15, 16]. The amounts of some hormones in the seeds affect the gene expression required for germination [13], and most of the genes involved in GA3 biosynthesis are up-regulated during seed germination [17]. Transcriptome methods have been used to study the gene expression of plant seeds during dormancy [18,19,20]. Transcriptome analysis reveals that the genes of many GA3 response elements are differentially expressed in Arabidopsis seeds with different dormancy levels [21]. In our recent work has shown that the seed germination of F. hupehensis was substantially promoted by GA3 treatment, followed by cryogenic stratification [22]. However, the physiological and molecular mechanisms by which GA3 promotes F. hupehensis seed germination remain unclear.

In this study, we constructed three independent cDNA libraries of three treated F. hupehensis seeds, including a CK group (dormant seeds, CK1, CK2, and CK3), a W group (germinated seeds after treatment with water, W1, W2 and W3), and a G group (germinated seeds after treatment with GA3, G1, G2, and G3), for Illumina HiSeq sequencing. We compared the seed vigor of these groups, as well as some physiological indicators related to their germination, hormone content, and transcriptomics among these seed groups. We also identified differentially expressed genes (DEGs) and transcription factors (TFs) related to seed germination and validated the expression genes involved in seed germination by qRT-PCR. Our data revealed that the physiological and transcriptomic aspects of the promotive effect of GA3 on F. hupehensis seed germination and provide insights into the involvement of GA3 in the seed germination mechanism of woody plants.

Results

Effects of GA3 treatment on physiological indexes during F. hupehensis seed germination

In this study, the germinated seeds were obtained by low-temperature stratification after treatment with water (W group) or GA3 (G group); whereas the CK group comprised dormant seeds. After 50 days of treatment, the germination rates of water- and GA3-treated seeds were 46.67 and 26.03%, respectively. The CK seeds showed no signs of germination (Fig. 1a and b). Thus, GA3 and low temperature treatment significantly increased the seed germination rate of F. hupehensis. Triphenyltformazan (TTF) staining and content analyses of TTF revealed that the germinated seeds treated with GA3 had a significantly higher TTF contents than water-treated and dormant seeds (Fig. 1c), which indicated that GA3 treatment increases seed vigor.

Fig. 1
figure1

Comparison of physiological indexes during the germination of F. hupehensis seeds treated by three methods. CK, W and G represent dormant seeds, germinated seeds treated with water, and germinated seeds treated with GA3, respectively. a Comparison of F. hupehensis seed germination under three treatments. b Comparison of F. hupehensis seed germination rates under three treatments. c Comparison of F. hupehensis seed vigor under three treatments. d Comparison of F. hupehensis seed water absorption rates under three treatments. Data were analyzed by SPSS, followed by Duncan’s honestly significant difference test at p ≤ 0.05. All Statistical analyses of data had three biological repeats

To clarify the physiological mechanism by which GA3 promotes seed germination, we measured the water absorption rate of seeds among the CK, W, and G groups. The water absorption rate of seeds in the G group was significantly higher than those of seeds in the W and CK groups before saturation (within 8 h). The G seed group also revealed the highest water absorption rate (Fig. 1d). Thus, GA3-treated seeds germinated faster than the W and CK seed because the former absorbed more water during germination. Overall, GA3 treatment shortened the physiological dormancy period, improved seed vigor, and promoted seed germination in F. hupehensis.

RNA-seq and de novo assembly of F. hupehensis reference transcriptome

Nine cDNA libraries (CK1, CK2, CK3, W1, W2, W3, G1, G2, and G3) were sequenced by Illumina HiSeq to detect the transcriptome level of gene expression information during F. hupehensis seed germination. Sequencing readings were deposited at the National Center for Biotechnology Information (NCBI) under SRA accession number SUB4314435. High-quality, clean reads were obtained by filtering low-quality reads obtained by sequencing. The Q30 of nine samples exceeded 94%, and the base content was uniform (Additional file 1: Table S1). A total of 116,932 unigenes were obtained after de novo assembly using Trinity software. The average length of unigenes was 82,979 bp, and the length of N50 was 1346 bp. The length of 68,984 (58.83%) unigenes were over 400 bp, while those of 30,024 (25.67%) unigenes were over 1000 bp (Table 1). Thus, the assembly quality of the transcriptome was satisfactory.

Table 1 Statistics of unigene stitching results

Function annotation and classification of F. hupehensis unigenes

Trinotate was used to compare the sample unigene sequences with a common functional database (Table 2); here, 68,715 (58.765%) annotated unigenes were obtained from NR database; 41,776 (35.727%) were obtained from NT database; and 49,676 (42.483%), 33,484 (28.635%), 30,293 (25.907%), 47,270 (40.425%), and 23,771 (20.329%) were obtained from Swiss-Prot, PFAM, Cluster of Orthologous Group (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively.

Table 2 Functional annotations of unigenes in the NR, NT, Swiss-Prot, PFAM, COG, GO and KEGG databases

Nr is the non-redundant NCBI collection of nucleotide and protein sequence database. A total of 68,715 unigenes were annotated to the NR database (Fig. 2a). F. hupehensis transcripts were highly similar to those Sesamum indicum (13.85%), Vitis vinifera (1.84%), and Coffea canephora (1.75%).

Fig. 2
figure2

Functional annotations of the unigenes of the F. hupehensis seed transcriptome. a NR annotated species distribution map similar to the F. hupehensis transcriptome. b KOG function annotation of F. hupehensis seeds. c GO function annotation of F. hupehensis seeds

COG is a database of orthologous gene families. A total of 30,293 unigenes of F. hupehensis were annotated to 25 COG pathways (Fig. 2b). A total of 8157 (19.01%) unigenes were annotated to the general function prediction classification, followed by 4776 (11.13%), 3507 (8.17%), 2790 (6.50%), 2009 (4.68%), and 834 (1.94%) annotations related to gene expression, including post-translational modification; protein turnover and chaperones; signal transduction mechanisms; translation, ribosomal structure, and biogenesis; and transcription, replication, recombination, and repair, respectively.

GO is a gene function database. A total of 47,270 unigenes from F. hupehensis were annotated into 59 GO pathways (Fig. 2c). The functions of unigenes in biological process classifications contained cellular process (67.76%), metabolic process (59.43%), and biological regulation (28.76%). Cell part (79.75%), organelle (49.04%), and organelle part (33.66%) were the most abundant functions in terms of cellular component classifications. In the molecular function classification, binding (63.47%), catalytic activity (53.55%), and transporter activity (7.49%) were more abundant. The main GO entries revealed that the cells divided frequently during seed germination, and some catalytic, metabolic, and binding activities were relatively high.

Comparative analysis of DEGs

DEGs were analyzed was performed using the RPKM method to determine the degree of overlap between the three seed groups. Compared with the dormant seeds, the total numbers of up-regulated genes in the germinated seeds treated with water and GA3 were 24,371 and 24,412, respectively. The numbers of down-regulated genes in the W and G groups were 7485 and 8776, respectively. Compared with the W seeds, G seeds contained 680 up-regulated genes and 1376 down-regulated genes (Fig. 3a, Additional file 2: Table S2). Various process genes were compared using Venn diagrams. A total of 31,856 DEGs were obtained between the W and CK groups, 33,188 DEGs were obtained between the G and CK groups, and 2056 DEGs were obtained between the G and W groups. A total of 1458 DEGs (331 commonly up-regulated and 75 commonly down-regulated) were shared among the three treatments (Fig. 3b-d, Additional file 2: Table S2), thus implying that these 1458 DEGs might be responsible for F. hupehensis seed germination. Correlation heat map analysis of the expression among samples revealed that the results were reproducible (Additional file 3: Figure S1). Based on similarities in gene expression, DEGs in the samples were generated by hierarchical clustering combined with K-means clustering. Hierarchical clustering of the gene expression profiles of the CK, W, and G seed groups showed that the DEGs could be divided into eight clusters (Additional file 4: Figure S2), and the genes of the same subclass had similar expression patterns (Additional file 2: Table S2, Additional file 5: Figure S3).

Fig. 3
figure3

Statistical analysis of differentially expressed unigenes (DEGs) during F. hupehensis seed germination. a Statistical analysis of up/down-regulated unigenes in the CK, W, and G groups. b Venn diagram of all DEGs. c Venn diagram of up-regulated genes. d Venn diagram of down-regulated genes

KEGG pathway analysis of DEGs

KEGG is a signal pathway database with an extremely rich signal pathway map, and the map of interaction between genes contained in a pathway. Enrichment of DEGs in the KEGG pathway was analyzed at a significance level of p < 0.05. The KEGG annotations indicated that 136 pathways between the G–CK groups and W–CK and 108 pathways between the G and W groups were enriched. In particular, the sesquiterpenoid and triterpenoid biosynthesis (map 00909), endocytosis (map 04144), MAPK signaling pathway-plant (map 04016), pyruvate metabolism (map 00620), ascorbate and aldarate metabolism (map 00053), biosynthesis of amino acids (map 01230), and carbon fixation in photosynthetic organisms (map 00710) pathways between the G and CK groups were significantly enriched. DEGs were significantly enriched in the sesquiterpenoid and triterpenoid biosynthesis (map 00909) pathways between the G and W groups. No significant difference in pathway was observed between the W and CK groups (Fig. 4). Some terpenoids, such as gibberellin, abscisic acid, and other plant hormones, are necessary for plant growth and development, whereas some species were necessary in regulating the relationship between plant and environment [23].

Fig. 4
figure4

Analysis of KEGG enrichment of DEGs in F. hupehensis

DEGs related to sugar metabolism in seed germination

The highest content of soluble sugar was found in dormant seeds, followed by the W seeds, and finally G seeds (Fig. 5a). The KEGG annotations suggested that eleven DEGs were related to the sugar pathway in F. hupehensis seeds. Cluster analysis revealed that the expression levels of all eleven DEGs significantly differs among the three seed groups (Fig. 5b). Phosphoglycerate kinase (PGK), phosphoglucomutase (PGM), pyruvate kinase (PK), enolase (ENO), 2,3-bisphosphoglycerate-dependent phosphoglycerate mutase (PGAM), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), glucose-6-phosphate isomerase (GPI), and 6-phosphofructo-2-kinase (PFK) were more highly expressed in germinated seeds was higher than that in dormant seeds, especially in G group seeds (Additional file 6: Table S3). The expression levels of DEGs were consistent with the trends of soluble sugar content, thereby indicating that these eleven DEGs might play important roles in promoting seed germination. Hence, exogenous GA3 could promote the expression level of key genes in the sugar metabolic pathway of F. hupehensis seeds.

Fig. 5
figure5

Sugar content and expression profiles of DEGs related to sugar metabolism during F. hupehensis seed germination. a Comparison of sugar content in F. hupehensis seed under three treatments. b Expression profiles of DEGs related to sugar metabolism in F. hupehensis seed under three treatments. The sample names are indicated at the bottom of the figure. Changes in expression level are represented by a change in color; green indicates a lower expression level, whereas red indicates a higher expression level. All data shown indicate the average mean of three biological replicates (n = 3). Means with different letters in each treatment represent a significant difference at p ≤ 0.05

Hormone concentrations and hormone-related DEGs in seed germination

To investigate the roles of endogenous hormones in seed germination, we determined the ZT, GA3, and ABA contents in F. hupehensis seeds. Both ZT and GA3 contents were significantly the highest, but ABA content was significantly the lowest in the G group. By contrast, both ZT and GA3 contents were significantly the lowest whereas ABA content was significantly the highest in the CK group (Fig. 6a). Therefore, GA3 treatment could promote seed germination through endogenous hormone accumulation. Furthermore, exogenous GA3 induced seed germination by increasing ZT and GA3 concentrations and decreasing ABA concentration, thus conforming to the transcriptome data. Pairwise comparison of the treatments (W–CK, G–CK, and G–W) revealed that two DEGs were related to ZT biosynthesis in F. hupehensis. Cluster analysis of genes revealed that the expression levels significantly differs among the three seed groups (Fig. 6b). The expression levels of adenylate isopentenyltransferase (IPT) and zeatin O-glucosyltransferase (ZOG) were higher in germinated seeds than that in dormant seeds, especially in the G group seeds, (Additional file 6: Table S3). These results were in accordance with those of endogenous hormone ZT content and indicated that the identified DEGs had a vital function in accelerating F. hupehensis seed germination. These results revealed that exogenous GA3 treatment could advance ZT biosynthesis by regulating the key genes.

Fig. 6
figure6

Endogenous ZT, GA3, and ABA concentrations and expression profiles of DEGs related to hormone pathways involved in F. hupehensis seed germination. a Comparison of ZT, GA3, and ABA concentrations in F. hupehensis seeds under three treatments. b Expression profiles of DEGs related to the ZT pathway in F. hupehensis seeds under three treatments. c Expression profiles of DEGs related to the GA3 pathway in F. hupehensis seeds under three treatments. d Expression profiles of DEGs related to the ABA pathway in F. hupehensis seeds under three treatments. The sample names are shown at the bottom of the figure. Changes in expression level are indicated by a change in color; green indicates a lower expression level, whereas red indicates a higher expression level. All data shown reflect the average mean of three biological replicates (n = 3). Means with different letters in each treatment represent a significant difference of p ≤ 0.05

From the F. hupehensis transcriptome, we obtained a key gene cytochrome P450 (CYP450) in the GA biosynthetic pathway, CYP450 was significantly highly expressed in the G seeds, followed by W seeds, and least expressed in dormant seeds (Fig. 6c, Additional file 6: Table S3). Moreover, cluster analysis revealed that the expression levels of two DEGs associated with the GA3 signal transduction pathway significantly differed among the three seed groups (Fig. 6c). The expression levels of GAI1 and gibberellin receptor (GID1) were significantly lower in the G group than in the W group and lowest in the CK group (Additional file 6: Table S3). Thus, the expression levels of GAI1 and GID1 were down-regulated by GA3 treatment, thereby increasing the hormonal signal transduction of GA3 during F. hupehensis seed germination.

One DEG was related to the ABA biosynthetic pathway and three DEGs in the ABA signal transduction pathway were observed. The expression levels of four DEGs significantly differs among the three groups (Fig. 6d). (+)-Abscisic acid 8′-hydroxylase (CYP707A), ABSCISIC ACID-INSENSITIVE 5 like 5 (ABI5.5), and protein phosphatase 2C (PP2C) were highly expressed in the G group but least expressed in the CK group. By contrast, the expression level of ABSCISIC ACID-INSENSITIVE 5 like 6 (ABI5.6) in dormant seeds was significantly higher than that in germinated seeds, and its expression level in GA3-treated seeds was significantly higher than that in water- treated seeds (Additional file 6: Table S3). Thus, exogenous GA3 may be conducive to alleviate the inhibitory effect of ABA signal transduction on the seed germination of F. hupehensis.

Other important DEGs related to seed germination

We examined several impartment pathways involved in seed germination. KEGG pathway analysis revealed that one DEG (pyruvate dehydrogenase E1 component beta subunit) related to the TCA cycle was clustered, and its highest and lowest expressions were observed germinated seeds (G group) and the CK group, respectively (Additional file 6: Table S3). Moreover, one DEG (auxin-responsive GH3 gene family) related to the auxin signaling pathway was expressed the highest in the germinated seeds (G group) and the least in the CK group (Additional file 6: Table S3). Seven DEGs were related to antioxidants, and the expression levels of catalase (CAT) and glutathione S-transferase (GST) were significantly higher in the G group than in the W group but lowest in the CK group. Moreover, superoxide dismutase (SOD), peroxiredoxin (POD), and glutamate-cysteine ligase (GCL) were expressed the highest in the germinated seeds (W group) but the least in the CK group (Additional file 6: Table S3). Twelve DEGs were associated with mRNA degradation. Among these 12 DEGs, the expression levels of seven subunits of CCR4-nontranscriptional complex (CCR4-NOT), enhancer of mRNA-decapping protein (Edcp), and 5′-3′exoribonuclease (XRN2) were significantly greater in the W group than in the G group but lowest in the CK group. Meanwhile, the expression level of U6 snRNA-associated Sm-like protein LSm4 was highest in the G group (Additional file 6: Table S3).

DEGs related to TFs in seed germination

The gene expression network regulated by TFs plays an important role in the growth and development of plants [24]. A total of 155,244 TFs were annotated in the transcriptome of F. hupehensis seeds and classified into 56 families (Fig. 7a). To better understand the molecular mechanism underlying the effects of GA3 treatment on seed germination, we measured the differential expression of TFs by cluster analysis and obtained 13 TFs with significant differences from 5 TF families among the three seed groups (Fig. 7b). Pairwise comparison (W–CK, G–CK, and G–W) showed that significantly up-regulated TFs included two TF members belonging to the MYB family (MYB44 and MYB86), four TFs from the WRKY family (WRKY14, WRKY22, WRKY28, and WRKY33), three TFs from the ERF family (ERF3, ERF12, and ERF25), and three TFs from the bHLH family (bHLH112, bHLH123, and bHLH137). The significantly down-regulated TFs included one TGA1 belonging to the bZIP family (Additional file 7: Table S4). These TFs might play significant roles in expediting the seed germination of F. hupehensis. Moreover, the mechanism by which GA3 regulates F. hupehensis seed germination involved an extremely intricate and complex transcriptional network. These findings provided basic information for studying the role of TFs in the promoting effect of GA3 on the seed germination of F. hupehensis.

Fig. 7
figure7

Analysis of transcription factors (TFs) involved in F. hupehensis seed germination. a Distribution of TF families in the transcriptome of F. hupehensis seeds. b Expression profiles of differentially expressed TFs involved in F. hupehensis seed germination. The sample names are shown at the bottom of the figure. Changes in expression level are represented by a change in color; green indicates a lower expression level, whereas red indicates a higher expression level. All data shown reflect the average mean of three biological replicates (n = 3). Means with different letters in each treatment represent a significant difference at p ≤ 0.05

Validation of DEGs by qRT-PCR

To verify the accuracy and reproducibility of the RNA-seq results, we randomly selected 24 genes in the related pathways of seed germination for qRT-PCR validation (Additional file 8: Table S5). The expression levels of the selected genes were calculated using the 2-ΔΔCt method. We compared the expression data of the three groups obtained by RNA-seq and qRT-PCR (Fig. 8a). The correlation between RNA-Seq results (RPKM) and qPCR results (2-ΔΔCt) results for the 24 DEGs was calculated using log2 fold variation measurements to produce a scatter plot. The qRT-PCR results of 24 DEGs were significantly similar to the RNA-seq results (R2 = 0.38, p < 0.01; Fig. 8b), which indicated that the RNA-seq data were reliable and accurate.

Fig. 8
figure8

qRT-PCR verification diagram of DEGs during F. hupehensis seed germination. a Comparison of the expression levels determined by qRT-PCR and RNA-seq from three treated seeds. b Correlation plot of the RNA-Seq results (RPKM) and qPCR (2-ΔΔCt) results. Result were calculated using log2 fold variation measurements. The R2 value represents the correlation between the RNA-seq and qPCR results. All data indacate mean ± SE (n = 3 with 500 seeds per replicate). Means with different letters in each treatment representa significant difference at p ≤ 0.05

Discussion

F. hupehensis reference transcriptome

Seed germination is based on the prominence of the radicle, some criteria for judging seed germination are based on the fact that some biological processes during seed germination mainly include nutrient metabolism, transcription, DNA repair, cell elongation and cell division recovery [25, 26]. Therefore, detailed information on gene expression is critical to understand the molecular mechanisms of any developmental process. An increasing number of analytical tools are being used to study the mechanism of seed dormancy, germination, and development [13, 19, 27]. Although F. hupehensis has great economic value, complete gene information for this species remains unavailable. In this current study, we describe the gene expression profiling of F. hupehensis seed germination under GA3 treatment. This work represented the first attempt to use Illumina sequencing technology to further understand the transcriptome of F. hupehensis seeds. RNA-seq was performed using Illumina HiSeq sequencing, which assembled 116,932 unigenes. The unigenes were used to perform BLASTX-based searches and annotations on the NT, NR, Swiss-Port, PFAM, GO, COG, and KEGG databases. Several key pathways associated with seed germination were also obtained by KEGG annotation analysis to identify DEGs in these pathways based on RPKM values combined with qRT-PCR data. Our transcriptome data suggested that a large number of DEGs were involved in various metabolic pathways, most of which were related to the regulation of gene expression, followed by energy production and metastasis. Only a small fraction of the DEGs found were related to signal hormone transduction and reactive cell division.

Role of sugar in F. hupehensis seed germination

During seed germination, soluble sugars are mostly used to synthesize or transform other substances, thus providing energy for seed germination and seedling growth [28]. Exogenous GA3 treatment can reduce the sugar contents [29, 30], which was consistent with our measurements of the lowest sugar content in germinated seeds treated with GA3. In addition, eleven DEGs (PGM, three GPIs, PFK, GAPDH, two PGKs, PGAM, ENO, and PK) involved in the glycolytic pathway were identified from the RNA-seq data, which implied that exogenous GA3 up-regulated these key genes to convert a-D-glucose-1P into pyruvate, which is the final product of the glycolytic pathway and is used for the metabolic conversion of intermediates in other substances (Fig. 9). Germinating seeds require more nutrients, and exogenous GA3 promoted seed germination by accelerating the oxidative decomposition of sugars to affect their mobilization of sugar and up-regulating the expression of key genes in the glycolytic pathway in a consistent direction to provide energy for F. hupehensis seed germination.

Fig. 9
figure9

Regulation model of exogenous GA3 and cold treatment promoting F. hupehensis seed germination. The colored arrow pointing upward indicates the up-regulation of genes, whereas the colored arrow pointing downward indicates the down-regulation of genes

Energy requirements for F. hupehensis seed germination

A considerable amount of energy is required during seed germination, and all biochemical reactions in plants enter the TCA cycle, providing power for the mitochondrial electron transport chain to generate ATP [31]. The energy required during seed germination mainly comes from glycolysis, pentose phosphate pathway, and the TCA cycle of cellular respiration [32]. The ATP content in dry seeds is low but gradually increases after seed germination, which proves that the respiratory pathway plays an important role in seed germination [33, 34]. The degradation products of energy reserve always undergo glycolysis, after which ATP synthesis is synthesized through the TCA cycle and mitochondrial electron transport. In the present study, KEGG annotation analysis suggested that eight genes (PGM, GPI, PFK, GAPDH, PGK, PGAM, ENO, and PK) related to the glycolytic pathway were hardly expressed in the dormant seeds but highly expressed in germinated seeds, especially in the GA3-treated seeds (Fig. 9). Hence, GA3 and water treatment required large amounts of energy to germinate seeds, and the energy requirement for most seed germination events was primarily achieved by glycolysis [10]. However, when the energy provided by anaerobic respiration cannot satisfy the needs of germinating seeds, the seeds undergo the TCA cycle to provide considerable energy for germination under oxygen-rich conditions [18]. In this study, KEGG metabolic pathway analysis identified PDHB as a group of rate-limiting enzymes in the TCA cycle that catalyzed the irreversible oxidative decarboxylation of pyruvate to acetyl-CoA; those enzymes connected the aerobic oxidation of sugars with the TCA cycle, thereby playing an important role in the energy metabolism of the mitochondrial respiratory chain [35]. The expression level of PDHB significantly differed among the three groups, with the highest expression found in the germinated seeds (especially in the G group) and the lowest expression found in the dormant group. Our results were consistent with the results of Weitbrecht et al. [31], who demonstrated that the pyruvate dehydrogenase transcript accumulates during cold stratification and GA3 treatment. Thus, although both water and GA3 treatments promoted seed germination and required more energy for seed germination, exogenous GA3 accelerated the TCA cycle to release more energy for seed germination by up-regulating PDHB. These respiratory pathways are critical for providing energy for various cellular functions during seed germination, as described in our study [32, 36]. Most of this energy was produced during germination, especially after GA3 treatment. Thus, F. hupehensis seeds must consume more sugar for energy conversion during germination.

Comparison of ZT, GA3, and ABA involved in F. hupehensis seed germination

ZT is the first type of cytokinin isolated and identified in plant growth and development; it is mainly present in seeds with high metabolic metabolism and presents antagonistic effects on seed germination inhibitors [37]. ZT was discovered in the hypocotyl of chickpea and found to promote nutrient storage and metabolism in cotyledons [38]. Pea seeds have peak ZT contents during germination, especially during radicle prominence [39]. Similarly, ZT contents in F. hupehensis seeds were highest after treatment with GA3, followed by that after treatment with water, and then in dormant seeds. These results suggest that ZT promoted the protrusion and extension of seed radicle. When plant needs ZT, ZT is formed by removing the nitrogen chain of ZT and then modifying to form ZT again [31]. In the present study, the expression level of IPT, which is the first key enzyme in the biosynthetic pathway of ZT, was significantly higher in F. hupehensis seeds germinated by GA3 treatment than in those subjected to other treatments. Trans-ZT was synthesized by cytokinin trans-hydroxylase (CYP735A) [40], which was also highly expressed in the germinated seeds treated with GA3. ZOG is an important enzyme in the trans-ZT metabolic pathway [41, 42]. Our data showed that the expression level of ZOG was significantly higher in the germinated seeds treated with GA3 than in those treated with water, indicating that the anabolism of ZT was relatively strong under GA3 treatment (Fig. 9). This result might be related to ZT mobilization of reserve metabolism and promotion of radicle protrusion [39]. Taken together, GA3 treatment could accelerate the anabolism of ZT to facilitate seed germination in F. hupehensis.

GA3 regulates various developmental processes, including seed germination and seedling development, throughout the plant life cycle [14, 39]. GA3 can break through the mechanical constraints of seed coats during seed germination to promote radicle protrusion [14, 43]. In addition, GA3 promotes cell division, which is vigorous during seed germination and related to the synthesis and catabolism of GA3 [44]. Consistent with the results of previous studies, our data showed that the GA3 content was the highest in F. hupehensis seeds treated with GA3. GA3 is biosynthesized from geranylgeranyl-PP, which is converted into GA12 by ent-copalyl diphosphate synthase. GA12 is then converted into GA3 by KAO, CYP450, GA20ox, and GA3ox enzymes [45]. In the present study, GA3 treatment significantly up-regulated the expression of the CYP450 gene in the GA3 biosynthetic pathway of F. hupehensis seeds. Some key genes in the GA3 signaling pathway also determined seed germination. DELLA protein is a plant growth inhibitor, while GID1 is a receptor for GA3 that degrades DELLA protein in plants by binding to GID1 receptors [46, 47]. GA3 expedites germination by down-regulating some germinated inhibitory proteins, such as DELLA [15, 48]. In the present work, GAI1 and GID1 expression levels were significantly lower in the GA3-treated seeds than in the water-treated seeds, and highest in dormant seeds. Overall, although both water and GA3 treatment could regulate the expression of related genes in GA3 pathway, exogenous GA3 could significantly up-regulate genes in GA3 biosynthesis and down-regulate negative regulators in GA3 signal transduction pathways. These phenomena increased the accumulation levels of GA3 and contribute to seed germination; this supposition was supported by the observed increase in endogenous GA3 concentration in F. hupehensis seeds (Fig. 9).

ABA regulates the accumulation of phytochemicals in seeds and controls the dehydration of seed development at later stages [49]. ABA maintains a relatively high level in dry seeds and rapidly declined after germination [50, 51]. Moreover, the adequate reduction of endogenous ABA content is a major prerequisite for complete germination because ABA inhibits the weakening and rupture of endosperm [52, 53]. GA3 can promote seed germination and counteract ABA inhibition in seeds [54]. Similar results were obtained in the present experiments, in which ABA contents were highest in dormant seeds and lowest in fast-germinating seeds (GA3-treated germinated seeds). ABA is synthesized in the carotenoid pathway, starting with zeaxanthin. NCED rapidly catalyzes the synthesis of ABA when seeds absorber water and is a key enzyme for ABA synthesis [55], whereas CYP707A is the core enzyme of ABA degradation [56]. The expression levels of ABA metabolic genes are associated with inhibition of seed germination [20, 32]. In the current study, the expression level of NCED in the dormant seeds of F. hupehensis was significantly higher than that in germinated seeds, but no significant difference was observed in the seeds after germination. This result suggested that the amount of ABA synthesized was relatively small during the germination. By contrast, CYP707A expression in the GA3-treated seeds was significantly higher than that in other seeds, which indicated that CYP707A mainly degraded ABA content in F. hupehensis seeds during germination. Our results revealed that water and GA3 treatment accelerated the full degradation of ABA content, but GA3 treatment was effectively increased the expression of the degrading gene CYP707A and resisted the inhibitory effect of ABA on seed germination. In the signaling pathway, PP2C is negatively correlated with the regulation of the ABA signaling pathway in spinach and Arabidopsis seeds [57, 58]. In the present work, PP2C1 and PP2C2 were highly expressed in germinated seeds treated with GA3 but minimally expressed in dormant seeds. ABI5 is a member of the basic leucine zipper (bZIP) family, which mediates cell responses to ABA in seeds and vegetative tissues [59]. The ABI5 gene negatively regulates seed germination in Arabidopsis thaliana [60,61,62]. Our RNA-seq results revealed that the expression of ABI5.6 in dormant seeds was significantly higher than that in germinated seeds, thus further verifying the role of ABI5 in seed germination. Interestingly, the expression level of ABI5.5, another member of the ABI5 family, contrasted that of ABI5.6, which was also inconsistent with previous studies. This contradiction might be attributed to differences in ABI5 family members or differences among species. The role of these two members of the ABI5 family in the seed germination of F. hupehensis requires further verification. The results implied that GA3 significantly increased the expression of negative regulatory genes in the ABA signaling pathway and inhibited the effect of ABA on seed germination. In summary, exogenous GA3 reduced the ABA content of seed by degradation of gene related to ABA pathway and then significantly up-regulated the negative feedback factors in the ABA signaling pathways to promote the seed germination in F. hupehensis (Fig. 9). The decrease in endogenous ABA during seed germination also supported this view.

Candidate TFs associated with F. hupehensis seed germination

In the present study, 155,244 TFs were obtained and classified into 56 families from the F. hupehensis transcriptome. Among these TF families, MYB, bHLH, WRKY, C2H2, FAR1, C3H, bZIP, and ERF accounted for a relatively large proportion. Our results were basically similar to those of the Tibetan Sophora moorcroftiana [63], Lepidium apetalum [64] and Chrysanthemum [65] transcriptomes, in which NAC, MYB_related, WRKY, C2H2, B3, FAR1, C3H, bZIP, ERF, and HD-ZIP TFs exert certain regulatory effects on plant growth and development and signal transduction [66].

The MYB family is the largest class of TFs in plants. Members of this family are mainly involved in plant growth and development and response to primary and secondary metabolic reactions [67]. In Arabidopsis seeds, the level of MYB44 transcript was up-regulated by 4 °C treatment as a negative regulator of ABA signaling [68]. Interestingly, MYB44 was also differentially up-regulated in germinated seeds after treatment with GA3 in F. hupehensis. We identified another MYB66 with the same up-regulated expression pattern as MYB44 in the G group of F. hupehensis seeds. Therefore, MYB might play an important regulatory role in seed germination with exogenous GA3 treatment. bHLH is the second largest TF class in plant seeds, and it plays an important role in regulating plant growth and development [69]. This TF class a light-stable repressor of seed germination and mediates the germination response to temperature [70]. Phytochrome-interacting factors, which are members of the Arabidopsis bHLH TF family, participate in light and gibberellin signal transduction in Arabidopsis and rice [71]. WRKY TFs are involved in various plant activities, such as development and metabolism [72], growth and senescence [73], and response to biotic and abiotic stresses [74]. Many proteins in the ERF family have been associated with different functions in cell processes, such as hormone signaling [75], regulation of metabolism [76], and developmental processes in various plant species [77]. In the present study, the expression levels of three bHLHs, four WRKYs, and three ERFs were significantly higher in the G group than in the two other groups of F. hupehensis seeds, thereby indicating that these 10 TFs might play some positive regulatory roles in the GA3 regulatory seed germination pathway (Fig. 9).

bZIP regulates various biological processes, such as pathogen defense, light and stress signals, seed maturation, and flower development [78]. bZIP16 is mainly expressed in seeds and could activate the GA pathway and inhibit ABA action, thereby promoting seed germination [79]. TGAla has been isolated from tobacco as a member of the bZIP class of TFs [80], whereas TGA1 and TGA4 regulated SA biosynthesis [81]. Comparison of our transcriptome data revealed that a TGA1 TF was a DEG among the different groups. This TF was highly expressed in dormant seeds but showed the lowest transcription level in GA3-treated seeds. Therefore, exogenous GA3 could play a negative role in the expression of TGA1, and the expression of TGA1 could inhibit inhibited seed germination (Fig. 9).

Candidate genes related to RNA degradation associated with F. hupehensis seed germination

In eukaryotes, two main pathways of mRNA decay are triggered by the shortening of poly(A) of mRNA. In the 5′–3′ pathway, mRNA is degraded by 5′–3′ exonuclease [82], removed by an enhancer of mRNA-decapping protein [83], and required to accurately cleave the development of the U6 snRNA-associated Sm-like protein [84]. In the 3′–5′ pathway, the CCR4-NOT transcription complex subunit degrades mRNA and inhibits translation and transcription [85]. Degradation of mRNAs stored in seeds is a prerequisite for germination [86, 87], and these two pathways of RNA degradation contribute to seed germination [88]. In the present study, KEGG pathway analysis revealed that twelve DEGs were associated with mRNA degradation. As a result, CCR4-NOT was significantly highly expressed in germinated seeds and no transcripts were detected in dormant seeds, consistent with the findings above. The positive effects of GA3 treatment on seed hypocotyl growth are related to their effects on nucleic acid [89], and this treatment delays RNA degradation [90]. In the present study, similar results were found in the seeds of F. hupehensis transcriptome, that was, the expression level of seven subunits of the CCR4-NOT complex were significantly higher in water-treated seeds than that in GA3-treated seeds. Conversely, the expression level of U6 snRNA-associated Sm-like protein LSm4 was the highest in GA3-treated seeds, probably due to LSm4, which precisely cleaved the mRNA required for development. Thus, GA3 treatment leaded seeds to a more active germination state. In summary, increased in the gene expression of degraded RNA promoted the seed germination of F. hupehensis, GA3 could delay the degradation of RNA and induce F. hupehensis seeds to a favorable germination state for subsequent germination (Fig. 9).

Conclusion

Exogenous GA3 increased seed vigor and water absorption rates to promote the seed germination of F. hupehensis. Transcriptomics approaches were used to study the effects of GA3 on F. hupehensis seed germination. A total of 116,932 unigenes were obtained from F. hupehensis seeds. Many DEGs involved in seed germination were obtained by comparing the transcriptomes of the seeds. Some key genes related to sugar, energy, hormones, RNA degradation, and some important TFs were involved in the seed germination of F. hupehensis (Fig. 9). The sugar metabolism level, the role of hormones, and the expression patterns of some important genes in pathways related to seed germination were further verified by measuring several physiological indicators and endogenous hormone contents, as well as qRT-PCR analysis. Differential expression analysis showed that GA3 also accelerated the oxidative decomposition of sugar by up-regulating key genes in the glycolytic pathway to release energy for germination. The expression levels of key genes related to hormone synthesis and signal transduction were affected by GA3 treatment, as well as the contents of three endogenous hormones (ZT, GA3, and ABA) in F. hupehensis seeds. GA3 could delay the RNA degradation of germination seeds to maintain a favorable germination state. In addition, some TF genes such as MYB, WRKY, ERF, bHLH, and bZIP were up-regulated by GA3, thus suggesting that these TFs played important roles in seed germination (Fig. 9). The transcriptome data of F. hupehensis can be used by further studies on low-fecundity and endangered species and provide theoretical evidence for protecting and utilizing these valuable resources.

Methods

Plant materials

Seeds were collected from a 25-year-old F. hupehensis and grown in Jingshan Botanical Garden (31°03′N, 113°11′E) in Yongxing Town, Jingshan County, Jingmen City, Hubei Province. Sampling of the F. hupehensis seeds was approved by Jingmen Forestry Bureau before collection. After the winged perianth was removed, the seeds were divided into three experimental treatments: dormant seed (CK) and germinated seeds subjected to low-temperature lamination treated with water (W) and GA3 (G). Each treatment had three biological replicates with 500 seeds per replicate. All the biological replicates were from the same seed stock. The sample was frozen in liquid nitrogen and stored in a − 80 °C refrigerator for transcriptome sequencing and the determination of physiological indices.

Determination of seed TTF, ZT, GA3, and ABA contents in seeds

Seed vigor was determined by the 2,3,5-triphenyltetrazolium chloride (TTC) method, which measured TTF contents. Oxidized colorless TTC produces hydrogen from dehydrogenase in the living cell tissues of seed embryos, and TTC in seed embryos is reduced to red TTF by hydrogen [91]. Other physiological indicators were determined by grinding seeds after liquid nitrogen treatment. Soluble sugar content was measured by anthrone colorimetry [92]. The contents of ZT, GA3, and ABA contents were measured as described by Ding et al. [93].

Library construction and transcriptome sequencing

The total RNA of F. hupehensis seeds was extracted using a TAKARA MIniBEST Plant RNA Extraction Kit in accordance with the manufacturer’s instructions. RNA samples were tested for degradation and impurities by using 1% agarose electrophoresis. Sample purity was measured using a NanoDrop 2000 microspectrophotometer (IMPLEN, CA, USA), and the integrity and concentration of the RNA sample were detected using an Agilent 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, CA, USA). Sequencing libraries were obtained using the NEBNext1Ultra™ RNA Library Prep Kit for Illumina* (NEB, USA) in accordance with the manufacturer’s instructions. In brief, total mRNA was isolated with Oligo (dT), broken into short fragments, and then used to synthesize the first strand of cDNA. The purified double-stranded cDNA was subjected to terminal repair, addition of base A, and sequencing ligation. Finally, target-size fragments were isolated for PCR amplification to complete the construction of nine cDNA libraries. The concentration and quality of these libraries were tested using Agilent 2100 and Qubit 2.0, respectively. The nine cDNA libraries were then sequenced by Annoroad Gene Technology Corporation (Beijing, China) using Illumina HiSeq 2500.

De novo assembly and functional annotation of unigenes

The original sequence data were obtained by Illumina HiSeq sequencing. High-quality reads (clean reads) were obtained by removing low-quality reads (bases with a mass value of Q ≤ 19, accounting for over 15% of the total bases), joint-pollution reads (the number of bases contaminated by the linker in reads is greater than 5 bp), and reads with N ratios greater than 5%. All clean reads of the nine libraries were obtained using Trinity software and de Bruijn method to assemble the full-length transcript sequences [94], and the longest transcripts of each gene were considered as unigenes based on the transcripts. The transcriptome assembly sequence was annotated with Trinotate, and functional annotations were performed using the databases PFAM (protein domain identification), Nr (NCBI nonredundant protein sequences), Swiss-Prot (a manually annotated and reviewed protein sequence database), GO annotation, COG annotation, and KEGG. PlantTFDB (Plant Transcription Factor Database) was used to annotate TFs.

DEG analysis

RPKM value were used to represent the expression abundance of reads corresponding to Unigenes. In this study, we used DEseq242 to compare the treatment group with the reference group and selected | log 2 Ratio | > 1 and Q < 0.05 genes as DEGs [95]. Enrichment analysis of DEGs was analyzed using the GO and KEGG databases to obtain a detailed description of the DEGs during seed germination. The DEGs were clustered with p < 0.05, which indicated that the cluster distribution was significant.

qRT-PCR analysis

DEGs were selected for qRT-PCR analysis, and 18S was used as the internal reference gene. Specific primers were designed based on the sequence of unigenes using Primer 5.0 and are listed in Table S4. The RNA extracted from each sample (600 ng) was used to synthesize single-strand cDNA with a PrimeScript RT Reagent Kit following the manufacturer’s instructions. A LineGene 9600 Plus real-time PCR instrument (Bori, Hangzhou) was used to perform PCR. The relative quantitative expression of the genes was calculated using the 2-ΔΔCt method [96]. Each sample was prepared in triplicate (500 seeds per replicate). The normalized values of relative expression and RPKM values were calculated using log2 fold variation measurements, and the correlation between RNA-seq and qPCR results was analyzed using these values.

Statistical analysis of data

Data were analyzed using Excel and SPSS by ANOVA followed by Tukey’s significant difference test at p ≤ 0.05. All data had three biological repeats.

Abbreviations

AAO3:

Abscisic-aldehyde oxidase

ABA:

Abscisic acid

ABA2:

Xanthoxin dehydrogenase

ABF:

ABA responsive element binding factor

ABI5.5:

ABSCISIC ACID-INSENSITIVE 5 like 5

ABI5.6:

ABSCISIC ACID-INSENSITIVE 5 like 6

ALDO:

Fructose-bisphosphate aldolase, class I

B-ARR:

Two-component response regulator ARR-B family

BPGM:

Bisphosphoglycerate/phosphoglycerate mutase

CAT:

Catalase

CCR4-NOT:

CCR4-nontranscriptional complex

CKX:

Cytokinin dehydrogenase

CPS:

Copalyl pyrophosphate synthase

CRE1:

Arabidopsis histidine kinase 2/3/4 (cytokinin receptor)

CYP450:

Cytochrome P450 monooxygenase 1

CYP707A:

(+)-Abscisic acid 8′-hydroxylase

CYP735A:

Cytokinin trans-hydroxylase

DEGs:

Differentially expressed genes

Edcp:

Enhancer of mRNA-decapping protein

ENO:

Enolase

FBP:

Fructose-1,6-bisphosphatase I

GA20OX :

GA20 oxidase

GA3 :

Gibberellin

GA3ox:

GA3 oxidase

GAI1:

DELLA protein GAI1

GAP2:

Glyceraldehyde-3-phosphate dehydrogenase (NAD(P)

GAPDH:

Glyceraldehyde 3-phosphate dehydrogenase

GAPN:

Glyceraldehyde-3-phosphate dehydrogenase (NADP+)

GAPOR:

Glyceraldehyde-3-phosphate dehydrogenase (ferredoxin)

GCL:

Glutamate-cysteine ligase

GID1:

Gibberellin receptor

GPI:

Glucose-6-phosphate isomerase

GST:

Glutathione S-transferase

IPT:

Adenylate isopentenyltransferase

KAO:

Ent-kaurenoic acid monooxygenase

LSm4:

U6 snRNA-associated Sm-like protein LSm4

NCED:

9-cis-epoxycarotenoid dioxygenase

PDHB:

Pyruvate dehydrogenase E1 component beta subunit

PFK:

6-phosphofructo-2-kinase

PfkC:

ADP-dependent phosphofructokinase/glucokinase

PFP:

Diphosphate-dependent phosphofructokinase

PGAM:

2,3-bisphosphoglycerate-dependent phosphoglycerate mutase

PGK:

Phosphoglycerate kinase

PGM:

Phosphoglucomutase

PK:

Pyruvate kinase

POD:

Peroxiredoxin

PP2C:

Protein phosphatase 2C

PVP:

Polyvinylpyrrolidone

PYR/PYL:

Abscisic acid receptor PYR/PYL family

QPMI:

2,3-bisphosphoglycerate-independent phosphoglycerate mutase

SLY/DID2:

F-box protein GID2

SnRK2:

Serine/threonine-protein kinase SRK2

SOD:

Superoxide dismutase

TCA cycle:

Tricarboxylic acid cycle

TF:

Phytochrome-interacting factor 4

TFs:

Transcription factors

TTC:

2,3,5-Triphenyltetrazolium chloride

TTF:

Triphenylformazan

XRN2:

5′-3′ exoribonuclease

ZEP:

Zeaxanthin epoxidase

ZOG:

Zeatin O-glucosyltransferase

ZT:

Zeatin

References

  1. 1.

    Su PL, Ming J, Liao HR, Hu GQ, Wang JC, Wu RC, Li ZY. Distribution of Fraxinus hupehensis germplasm resources. Hubei Forest Sci. Technol. 1994;1:35–6.

  2. 2.

    Su PL, Ming J, Liao HR, Hu GQ, Wang JC, Wu RC, Li ZY. The comprehensive utilization value and prospect of the development of Fraxinus hupehensis resources. Hubei Forest Sci. Technol. 1995;4:1–5.

  3. 3.

    Peng FS, Li HJ. The second batch of rare and endangered plants in Hubei. Hubei Forest Sci Technol. 1990;8:38–9.

  4. 4.

    Zhang MH. Seasonal management of the festival ash in the Lingnan area. Chin Flower Bonsai. 2003;8:51.

  5. 5.

    Zhang J. Research overview and prospects of Fraxinus hupehensis. Hubei Forest Sci. Technol. 2015;44:27–9.

  6. 6.

    Wang CW, Liu R, Yi MB. Fraxinus hupehensis seed seedling raising technique. Hubei Forest Sci. Technol. 2013;42:87–8.

  7. 7.

    Wang CY, Zhao MG, Yao YJ, Wang CJ, Shu ZM. Study on the technique of Fraxinus hupehensis cutting. J Beijing Forest Univ. 2001;23:18–20.

  8. 8.

    Ma L, Yan H, Qu J, Zhang Y. Extraction technics of total coumarins from the leaf of Fraxinus Hupehensis. Herald Med. 2010;29:925–8.

  9. 9.

    Su PL. Reproduction of Fraxinus hupehensis. Hubei Forest Sci. Technol. 1995;3:14–20.

  10. 10.

    Yang P, Liu X, Wang X, Chen H, Chen F, Shen S. Proteomic analysis of rice (Oryza sativa) seeds during germination. Proteomics. 2010;7:3358–68.

  11. 11.

    Ogawa M, Hanada A, Yamauchi Y, Kuwahara A, Kamiya Y, Yamaguchi S. Gibberellin biosynthesis and response during Arabidopsis seed germination. Plant Cell. 2003;15:1591–604.

  12. 12.

    Shibakusa R. Effect of GA3, IAA, ethrel and BA on seed germination of Picea glehnii masters. J Jpn For Society. 2008;62:440–3.

  13. 13.

    Pawłowski TA. Proteome analysis of Norway maple (Acer platanoides L.) seeds dormancy breaking and germination: influence of abscisic and gibberellic acids. BMC Plant Biol. 2009;9:48–61.

  14. 14.

    Sun TP, Gubler F. Molecular mechanism of gibberellin signaling in plants. Annu Rev Plant Biol. 2004;55:197–223.

  15. 15.

    Hussain A, Cao D, Peng J. Identification of conserved tyrosine residues important for gibberellin sensitivity of Arabidopsis RGL2 protein. Planta. 2007;226:475–83.

  16. 16.

    Cao D, Cheng H, Wu W, Soo HM, Peng J. Gibberellin mobilizes distinct DELLA-dependent transcriptomes to regulate seed germination and floral development in Arabidopsis. Plant Physiol. 2006;142:509–25.

  17. 17.

    Penfield S, Li Y, Gilday AD, Graham S, Graham IA. Arabidopsis ABA INSENSITIVE4 regulates lipid mobilization in the embryo and reveals repression of seed germination by the endosperm. Plant Cell. 2006;18:1887–99.

  18. 18.

    Yu Y, Guo G, Lv D, Hu Y, Li J, Li X, Yan Y. Transcriptome analysis during seed germination of elite Chinese bread wheat cultivar Jimai 20. BMC Plant Biol. 2014;14:20–39.

  19. 19.

    Mangrauthia SK, Agarwal S, Sailaja B, Sarla N, Voleti SR. Transcriptome analysis of Oryza sativa (Rice) seed germination at high temperature shows dynamics of genome expression associated with hormones signalling and abiotic stress pathways. Trop Plant Biol. 2016;9:215–28.

  20. 20.

    Sun J, Wang P, Zhou T, Rong J, Jia H, Liu Z. Transcriptome analysis of the effects of shell removal and exogenous gibberellin on germination of Zanthoxylum seeds. Sci Rep. 2017;7:8521–31.

  21. 21.

    Hauvermale AL, Tuttle KM, Takebayashi Y, Seo M, Steber CM. Loss of Arabidopsis thaliana seed dormancy is associated with increased accumulation of the GID1 GA hormone receptors. Plant Cell Physiol. 2015;56:1773–85.

  22. 22.

    Liu Y, Zhang J, Gang LI, Wang Q, Xu SY. The research for the germination characteristics of Fraxinus hupehensis seeds. Seed. 2016;35:27–31.

  23. 23.

    Yamaguchi S. Gibberellin biosynthesis in Arabidopsis. Phytochem Rev. 2006;5(1):39–47.

  24. 24.

    Lee BH, Henderson DA, Zhu JK. The Arabidopsis cold-responsive transcriptome and its regulation by ICE1. Plant Cell. 2005;17:3155–75.

  25. 25.

    Barrôco RM, Van Poucke K, Bergervoet JH, De Veylder L, Groot SP, Inzé D, Engler G. The role of the cell cycle machinery in resumption of postembryonic development. Plant Physiol. 2005;137:127–40.

  26. 26.

    Masubelele NH, Dewitte W, Menges M, Maughan S, Collins C, Huntley R, Nieuwland J, Scofield S, Murray JA. D-type cyclins activate division in the root apex to promote seed germination in Arabidopsis. Proc Natl Acad Sci U S A. 2005;102:15694–9.

  27. 27.

    Catusse J, Job C, Job D. Transcriptome- and proteome-wide analyses of seed germination. C R Biol. 2008;331:815–22.

  28. 28.

    Downie B, Bewley JD. Soluble sugar content of white spruce (Picea glauca) seeds during and after germination. Physiol Plant. 2000;110:1–12.

  29. 29.

    Sun P, Jiangui LI, Wang N. Effect of foliar fertilizers and GA3 on sugar accumulation in fruit of Junzao jujube during its growth and development. Acta Agric Boreali-Occidentalis Sin. 2011;20:98–102.

  30. 30.

    Fan WL, Li TL, Feng Y. Effect of gibberellin (CA_3) on carbohydrate metabolism of Aralia continentalis seeds in the dormancy-breaking process. J Shenyang Agric Univ. 2013;44:708–11.

  31. 31.

    Weitbrecht K, Müller K, Leubnermetzger G. First off the mark: early seed germination. J Exp Bot. 2011;62:3289–309.

  32. 32.

    Bao YZ, Yao ZQ, Cao XL, Peng JF, Xu Y, Chen MX, et al. Transcriptome analysis of Phelipanche aegyptiaca seed germination mechanisms stimulated by fluridone, TIS108, and GR24. PLoS One. 2017;12:e0187539.

  33. 33.

    Botha FC, Potgieter GP, Botha AM. Respiratory metabolism and gene expression during seed germination. Plant Growth Regul. 1992;11:211–24.

  34. 34.

    Benamar A, Rolletschek H, Borisjuk L, Avelange-Macherel MH, Curien G, Mostefai HA, et al. Nitrite-nitric oxide control of mitochondrial respiration at the frontier of anoxia. Biochim Biophys Acta Bioenerg. 2008;1777:1268–75.

  35. 35.

    Li A, Zhang Y, Zhao Z, Wang M, Zan L. Molecular characterization and transcriptional regulation analysis of the bovine PDHB gene. PLoS One. 2016;11:e0157445.

  36. 36.

    Yao Z, Tian F, Cao X, Xu Y, Chen M, Xiang B, et al. Global transcriptomic analysis reveals the mechanism of Phelipanche aegyptiaca seed germination. Int J Mol Sci. 2016;17:1139–58.

  37. 37.

    Wu YQ, Leng YF, Xia C, Zhou SF, Lan H. Comparative analysis of endogenous hormones in maize seeds with different dormancy characteristics. Acta Prataculturae Sin. 2015;24:213–9.

  38. 38.

    Villalobos N, Martin L. Involvement of cytokinins in the germination of chick-pea seeds. Plant Growth Regul. 1992;11:277–91.

  39. 39.

    Stirk WA, Novák O, Zižková E, Motyka V, Strnad M, Staden JV. Changes in endogenous cytokinins during germination and seedling establishment of Tagetes minuta L. Plant Growth Regul. 2005;47:1–7.

  40. 40.

    Cai L, Zhang L, Fu Q, Xu ZF. Identification and expression analysis of cytokinin metabolic genes IPTs, CYP735A and CKXs in the biofuel plant Jatropha curcas. Peerj. 2018;6:e4812.

  41. 41.

    Miyawaki K, Tarkowski P, Miho MK, Kato T, Sato S, Tarkowska D, et al. Roles of Arabidopsis ATP/ADP isopentenyltransferases and tRNA isopentenyltransferases in cytokinin biosynthesis. Proc Natl Acad Sci U S A. 2006;103:16598–603.

  42. 42.

    Martin RC, Mok MC, Mok DW. Isolation of a cytokinin gene, ZOG1, encoding zeatin O-glucosyltransferase from Phaseolus lunatus. Proc Natl Acad Sci U S A. 1999;96:284–9.

  43. 43.

    Richards DE, King KE, Tahar A, Harberd NP. How gibberellin regulates plant growth and development: a molecular genetic analysis of gibberellin signaling. Annu Rev Plant Physiol Plant Mol Biol. 2001;52:67–88.

  44. 44.

    Fabian T, Lorbiecke R, Umeda M, Sauter M. The cell cycle genes cycA1;1 and cdc2Os-3 are coordinately regulated by gibberellin in planta. Planta. 2000;211:376–83.

  45. 45.

    Nelson DR, Schuler MA, Paquette SM, Werck-Reichhart D, Bak S. Comparative genomics of Rice and Arabidopsis. Analysis of 727 cytochrome P450 genes and pseudogenes from a monocot and a dicot. Plant Physiol. 2004;135:756–72.

  46. 46.

    Olszewski N, Sun TP, Gubler F. Gibberellin signaling: biosynthesis, catabolism, and response pathways. Plant Cell. 2002;14 Suppl:S61–80.

  47. 47.

    Ueguchitanaka M, Nakajima M, Motoyuki A, Matsuoka M. Gibberellin receptor and its role in gibberellin signaling in plants. Annu Rev Plant Biol. 2007;58:183–98.

  48. 48.

    Fu X, Richards DE, Tahar A, Hynes LW, Ougham H, Peng J, Harberd NP. Gibberellin-mediated proteasome-dependent degradation of the barley DELLA protein SLN1 repressor. Plant Cell. 2002;14:3191–200.

  49. 49.

    Tunnacliffe A, Wise MJ. The continuing conundrum of the LEA proteins. Naturwissenschaften. 2007;94:791–812.

  50. 50.

    Hermann K, Meinhard J, Dobrev P, Linkies A, Pesek B, Hess B, et al. 1-Aminocyclopropane-1-carboxylic acid and abscisic acid during the germination of sugar beet (Beta vulgaris L.): a comparative study of fruits and seeds. J Exp Bot. 2007;58:3047–60.

  51. 51.

    Zhu G, Ye N, Zhang J. Glucose-induced delay of seed germination in rice is mediated by the suppression of ABA catabolism rather than an enhancement of ABA biosynthesis. Plant Cell Physiol. 2009;50:644–51.

  52. 52.

    Müller K, Tintelnot S, Gerhard LM. Endosperm-limited Brassicaceae seed germination: abscisic acid inhibits embryo-induced endosperm weakening of Lepidium sativum (cress) and endosperm rupture of cress and Arabidopsis thaliana. Plant Cell Physiol. 2006;47:864–77.

  53. 53.

    Linkies A, Müller K, Morris K, Turecková V, Wenk M, Cadman CS, et al. Ethylene interacts with abscisic acid to regulate endosperm rupture during germination: a comparative approach using Lepidium sativum and Arabidopsis thaliana. Plant Cell. 2009;21:3803–22.

  54. 54.

    Kucera B, Cohn MA, Gerhard LM. Plant hormone interactions during seed dormancy release and germination. Seed Sci Res. 2005;15:281–307.

  55. 55.

    Hwang SG, Chen HC, Huang WY, Chu YC, Shii CT, Cheng WH. Ectopic expression of rice OsNCED3 in Arabidopsis increases ABA level and alters leaf morphology. Plant Sci. 2010;178:12–22.

  56. 56.

    Lee HG, Lee K, Seo PJ. The Arabidopsis MYB96 transcription factor plays a role in seed dormancy. Plant Mol Biol. 2015;87:371–81.

  57. 57.

    Gosti F, Beaudoin N, Serizet C, Webb AA, Vartanian N, Giraudat J. ABI1 protein phosphatase 2C is a negative regulator of abscisic acid signaling. Plant Cell. 1999;11:1897–909.

  58. 58.

    Hirayama T, Shinozaki K. Perception and transduction of abscisic acid signals: keys to the function of the versatile plant hormone ABA. Trends Plant Sci. 2007;12:343–51.

  59. 59.

    Wang Z, Cheng K, Wan L, Yan L, Jiang H, Liu S, et al. Genome-wide analysis of the basic leucine zipper (bZIP) transcription factor gene family in six legume genomes. BMC Genomics. 2015;16:1053–68.

  60. 60.

    Piskurewicz U, Jikumaru Y, Kinoshita N, Nambara E, Kamiya Y, Lopez-Molina L. The gibberellic acid signaling repressor RGL2 inhibits Arabidopsis seed germination by stimulating abscisic acid synthesis and ABI5 activity. Plant Cell. 2008;20:2729–45.

  61. 61.

    Jin D, Wu M, Li B, Bücker B, Keil P, Zhang S, et al. The COP9 signalosome regulates seed germination by facilitating protein degradation of RGL2 and ABI5. PLoS Genet. 2018;14:e1007237.

  62. 62.

    Zhao FY, Jiao CJ, Chen Q, WangTS TCF, Gao YM. PKS5 kinase is involved in ABA response through phosphorylating ABI5 in arabidopsis. Plant Physiol. J. 2015;51:1719–28.

  63. 63.

    Li H, Yao W, Fu Y, Li S, Guo Q. De novo assembly and discovery of genes that are involved in drought tolerance in Tibetan Sophora moorcroftiana. PLoS One. 2015;10:e111054.

  64. 64.

    Zhou Q, Zhao HX, Li PP, Zeng WJ, Li YH, Ge FW, et al. De novo characterization of the seed transcriptome of Lepidium apetalum Willd. China Biotechnol. 2016;36:38–46.

  65. 65.

    Xu Y, Gao S, Yang Y, Huang M, Cheng L, Wei Q, Fei Z, Gao J, Hong B. Transcriptome sequencing and whole genome expression profiling of chrysanthemum under dehydration stress. BMC Genomics. 2013;14:662–77.

  66. 66.

    Wang XC, Zhao QY, Ma CL, Zhang ZH, Cao HL, Kong YM, et al. Global transcriptome profiles of Camellia sinensis during cold acclimation. BMC Genomics. 2013;14:415–29.

  67. 67.

    Waites R, Selvadurai HR, Oliver IR, Hudson A. The PHANTASTICA gene encodes a MYB transcription factor involved in growth and dorsoventrality of lateral organs in Antirrhinum. Cell. 1998;93:779–89.

  68. 68.

    Nguyen XC, Hoang MH, Kim HS, Lee K, Liu XM, Kim SH, et al. Phosphorylation of the transcriptional regulator MYB44 by mitogen activated protein kinase regulates Arabidopsis seed germination. Biochem Biophys Res Commun. 2012;423:703–8.

  69. 69.

    Toledo-Ortiz G, Huq E, Quail PH. The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003;15:1749–70.

  70. 70.

    Penfield S, Josse EM, Kannangara R, Gilday AD, Halliday KJ, Graham IA. Cold and light control seed germination through the bHLH transcription factor SPATULA. Curr Biol. 2005;15:1998–2006.

  71. 71.

    Leivar P, Monte E. PIFs: systems integrators in plant development. Plant Cell. 2014;26:56–78.

  72. 72.

    Kato N, Dubouzet E, Kokabu Y, Yoshida S, Taniguchi Y, Dubouzet JG, et al. Identification of a WRKY protein as a transcriptional regulator of benzylisoquinoline alkaloid biosynthesis in Coptis japonica. Plant Cell Physiol. 2007;48:8–18.

  73. 73.

    Miao Y, Laun T, Zimmermann P, Zentgraf U. Targets of the WRKY53 transcription factor and its role during leaf senescence in Arabidopsis. Plant Mol Biol. 2004;55:853–67.

  74. 74.

    Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, et al. Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008;6:486–503.

  75. 75.

    Ohmetakagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995;7:173–82.

  76. 76.

    Zhang JY, Broeckling CD, Blancaflor EB, Sledge MK, Sumner LW, Wang ZY. Overexpression of WXP1, a putative Medicago truncatula AP2 domain-containing transcription factor gene, increases cuticular wax accumulation and enhances drought tolerance in transgenic alfalfa (Medicago sativa). Plant J. 2005;42:689–707.

  77. 77.

    George C, Michael M, Elizabeth K, Sarah H, Schmidt RJ. The control of spikelet meristem identity by the branched silkless1 gene in maize. Science. 2002;298:1238–41.

  78. 78.

    Marc J, Bernd W, Wolfgang DL, Jesus VC, Jens T, Thomas K, et al. bZIP transcription factors in Arabidopsis. Trend Plant Sci. 2002;7:106–11.

  79. 79.

    Hsieh WP, Hsieh HL, Wu SH. Arabidopsis bZIP16 transcription factor integrates light and hormone signaling pathways to regulate early seedling development. Plant Cell. 2012;24:3997–4011.

  80. 80.

    Katagiri F, Lam E, Chua NH. Two tobacco DNA-binding proteins with homology to the nuclear factor CREB. Nature. 1989;340:727–30.

  81. 81.

    Sun T, Busta L, Zhang Q, Ding P, Jetter R, Zhang Y. TGACG-BINDING FACTOR 1 (TGA1) and TGA4 regulate salicylic acid and pipecolic acid biosynthesis by modulating the expression of SYSTEMIC ACQUIRED RESISTANCE DEFICIENT 1 (SARD1) and CALMODULIN-BINDING PROTEIN 60g (CBP60g). New Phytol. 2017;217:344–54.

  82. 82.

    Cheng CP, Jaag HM, Jonczyk M, Serviene E, Nagy PD. Expression of the Arabidopsis Xrn4p 5′–3′ exoribonuclease facilitates degradation of tombusvirus RNA and promotes rapid emergence of viral variants in plants. Virology. 2007;368:238–48.

  83. 83.

    Martin FG, Christy F, Bodil N, Jens LA. Multiple processing body factors and the ARE binding protein TTP activate mRNA decapping. Mol Cell. 2005;20:905–15.

  84. 84.

    Carlos PR, Tamara HV, Rosa LC, María  MC, Julio S. LSM proteins provide accurate splicing and decay of selected transcripts to ensure normal Arabidopsis development. Plant Cell. 2012;24:4930–47.

  85. 85.

    Sandler H, Kreth J, Timmers HT, Stoecklin G. Not1 mediates recruitment of the deadenylase Caf1 to mRNAs targeted for degradation by tristetraprolin. Nucleic Acids Res. 2011;39:4373–86.

  86. 86.

    Rajjou L, Duval M, Gallardo K, Catusse J, Bally J, Job C, et al. Seed germination and vigor. Annu Rev Plant Biol. 2012;63:507–33.

  87. 87.

    Fleming MB, Richards CM, Walters C. Decline in RNA integrity of dry-stored soybean seeds correlates with loss of germination potential. J Exp Bot. 2017;68:2219–30.

  88. 88.

    Isabelle BS, Stéphanie P, Françoise C, Juliette L, Christophe B. 5′ to 3′ mRNA decay contributes to the regulation of Arabidopsis seed germination by dormancy. Plant Physiol. 2017;173:1709–23.

  89. 89.

    Khafagy EZ, Mousa AM. Nucleic acids and protein metabolic changes during germination of cotton seed. Z Pflanzenphysiol. 1982;107:321–8.

  90. 90.

    Liu ZY, Jiang WB. Effects of GA3 on postharvest lignification of green asparagus. Sci Agric Sin. 2005;38:383–7.

  91. 91.

    Harty RL, Paleg LG, Aspinall D. Quantitative reduction of triphenyl tetrazolium chloride as a measure of viability in cereal seeds. Aust J Exp Agric. 1972;12:517–22.

  92. 92.

    Gill PK, Sharma AD, Singh P, Bhullar SS. Changes in germination, growth and soluble sugar contents of Sorghum bicolor (L.) Moench seeds under various abiotic stresses. Plant Growth Regul. 2003;40:157–62.

  93. 93.

    Ding J, Shen ZD, Fang YX, Feng XX, Li L, Ni JS. Extraction, isolation and biological identification of plant endogenous hormones. Plant Physiol J. 1979;2:27–39.

  94. 94.

    Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

  95. 95.

    Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26:136–8.

  96. 96.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.

Download references

Acknowledgements

The authors would like to thank Jingmen Forestry Bureau reviewers for providing us permission to collect Fraxinus hupehensis seeds.

Funding

This study was funded by National Natural Science Foundation of China (Grant no. 31200528), The Special projects for technological innovation in Hubei Province (Grant Recipient: Shuiyuan Cheng), and the Foundation of Science and Technology Program of Enshi City. The funding bodies did not play any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the authors on reasonable request (Feng Xu, xufeng198@126.com; Qiling Song, song70@126.com).

Author information

FX, SC, QS, and ZC conceived and designed the experiment. JZ and MZ provide seed and treatment advice. QS, GN and JY treated seeds and measured physiological indicators. QS, WZ and YL analyzed the data. QS and FX wrote the paper. All of the authors read and approved the final manuscript.

Correspondence to Feng Xu.

Ethics declarations

Ethics approval and consent to participate

Sampling of the F. hupehensis seeds was approved by Jingmen Forestry Bureau before collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Evaluation of sample sequencing data. (XLSX 9 kb)

Additional file 2:

Table S2. All DEGs from the transcriptome of F. hupehensis seeds. (XLSX 5488 kb)

Additional file 3:

Figure S1. Correlation coefficients among the RPKM of unigenes of the samples. All data shown indicate the results of three biological replicates (n = 3). (PDF 1132 kb)

Additional file 4:

Figure S2. Clusters of DEGs obtained by K-means. DEGs were divided into eight subclasses. All data shown reflect the results of three biological replicates (n = 3). (PDF 977 kb)

Additional file 5:

Figure S3. Hierarchical clustering of DEGs between the W and CK groups, the G and CK groups, and the G and W groups. All data shown indicate the results of three biological replicates (n = 3). The sample names are shown at the bottom of the figure. Changes in expression level are indicated by a change in color; blue indicates a lower expression level, whereas red indicates a higher expression level. (PDF 5411 kb)

Additional file 6:

Table S3. Key DEGs related to F. hupehensis seed germination. (XLSX 17 kb)

Additional file 7:

Table S4. Differentially expressed TFs related to F. hupehensis seed germination. (XLSX 11 kb)

Additional file 8:

Table S5. Primers used in qRT-PCR. (XLSX 12 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Song, Q., Cheng, S., Chen, Z. et al. Comparative transcriptome analysis revealing the potential mechanism of seed germination stimulated by exogenous gibberellin in Fraxinus hupehensis. BMC Plant Biol 19, 199 (2019) doi:10.1186/s12870-019-1801-3

Download citation

Keywords

  • Fraxinus hupehensis
  • Seed germination
  • Transcriptome
  • Germination
  • Exogenous gibberellin
  • Differentially expressed genes