- Research article
- Open Access
Transcriptomic analysis of flower opening response to relatively low temperatures in Osmanthus fragrans
BMC Plant Biology volume 20, Article number: 337 (2020)
Sweet osmanthus (Osmanthus fragrans Lour.) is one of the top ten traditional ornamental flowers in China. The flowering time of once-flowering cultivars in O. fragrans is greatly affected by the relatively low temperature, but there are few reports on its molecular mechanism to date. A hypothesis had been raised that genes related with flower opening might be up-regulated in response to relatively low temperature in O. fragrans. Thus, our work was aimed to explore the underlying molecular mechanism of flower opening regulated by relatively low temperature in O. fragrans.
The cell size of adaxial and abaxial petal epidermal cells and ultrastructural morphology of petal cells at different developmental stages were observed. The cell size of adaxial and abaxial petal epidermal cells increased gradually with the process of flower opening. Then the transcriptomic sequencing was employed to analyze the differentially expressed genes (DEGs) under different number of days’ treatments with relatively low temperatures (19 °C) or 23 °C. Analysis of DEGs in Gene Ontology analysis showed that “metabolic process”, “cellular process”, “binding”, “catalytic activity”, “cell”, “cell part”, “membrane”, “membrane part”, “single-organism process”, and “organelle” were highly enriched. In KEGG analysis, “metabolic pathways”, “biosynthesis of secondary metabolites”, “plant-pathogen interaction”, “starch and sucrose metabolism”, and “plant hormone signal transduction” were the top five pathways containing the greatest number of DEGs. The DEGs involved in cell wall metabolism, phytohormone signal transduction pathways, and eight kinds of transcription factors were analyzed in depth.
Several unigenes involved in cell wall metabolism, phytohormone signal transduction pathway, and transcription factors with highly variable expression levels between different temperature treatments may be involved in petal cell expansion during flower opening process in response to the relatively low temperature. These results could improve our understanding of the molecular mechanism of relatively-low-temperature-regulated flower opening of O. fragrans, provide practical information for the prediction and regulation of flowering time in O. fragrans, and ultimately pave the way for genetic modification in O. fragrans.
Osmanthus fragrans Lour. (Oleaceae) is one of the top ten traditional ornamental plants in China with more than 2500 years’ history of cultivation . It is a small evergreen tree, grown as ornamental plants for its fragrant edible flowers. Cultivars of O. fragrans could be divided into once-flowering group and recurrent-flowering group according to different flowering habits . The once-flowering cultivars bloom in autumn and the flowering time varies greatly in different areas, such as in Hangzhou, Shanghai, Nanjing, and Suzhou, or even in different years in the same area . The researches on the flowering time of different cultivars indicated that relatively low temperature before blooming is the most important environmental factor determining the flower opening of O. fragrans [3, 4]. However, the knowledge of molecular mechanism of flower opening in O. fragrans, especially response to relatively low temperature, still remains limited.
In many higher plants, the growth of flower petals is the most remarkable process during flower opening. Flower petals are the most important component of reproductive organs and play vital roles in attracting the suitable pollinator(s). Flower color, size, shape and appearance, which are determined by flower petal, are important traits appreciated by the breeders and consumers . One of the traits, the size of flower petal is determined by cell division at early phases of petal growth and cell expansion at later stages of flower opening .
Cell expansion is accompanied by a series of process including cell wall loosening, cellulose biosynthesis, polysaccharides conversion into soluble carbohydrate, ion and water uptake, and cytoskeleton rearrangement . Cell wall loosening proteins include expansin (EXP), xyloglucan endotransglycosylase/hydrolase (XTH), endo-1,4-β-D-glucanase, and pectinase . Among them, EXP and XTH particularly participate in disrupting the non-covalent bonds between the cellulose microfibrils and the cross-linking glycans of the cell wall to increase the cell-wall creep rate, while pectinase and endo-1,4-β-D-glucanase can degrade the cell wall . In wintersweet (Chimonanthus praecox), the expressions of EXP genes increase during flower opening . The developing petals of carnation show high activities of cellulase and pectin esterase . These findings reveal that petal growth relevant to flower opening is probably attributed to cell expansion. Moreover, soluble carbohydrates depending on the degradation of polysaccharides can act as osmotically active compounds which could lower the osmotic water potential and facilitate water influx in order to allow cell expansion . The concentration of soluble carbohydrates in the petals will increase in the flower opening process of plants such as carnation , rose , chrysanthemum , Tweedia caerulea , and lisianthus .
Cell expansion is regulated by both external factors, such as temperature, humidity, and the quality and quantity of light, and internal factors, such as the circadian clock and phytohormones [18, 19]. Phytohormones are the most important mediators regulating flower opening and could be affected by circadian factors or environmental factors. As far as the present researches are concerned, most of the phytohormones, such as auxin (AUX), gibberellin (GA), ethylene (ETH), brassinosteroid (BR), jasmonic acid (JA) and abscisic acid (ABA), are proven to affect flower opening [5, 6, 20,21,22].
In this research, potted plants of O. fragrans ‘Yanhong Gui’ were employed as materials to research the effects of different temperature conditions on flower opening process under controlled relatively low temperatures. Then the transcriptomic sequencing was used to analyze the differentially expressed genes (DEGs) after different number of days’ treatment under relatively high or low temperatures to figure out the key genes involved in the regulation of flower opening in relation to relatively low temperature. This study would lay foundation on fully revealing the molecular mechanism of relatively-low-temperature-regulated flower opening of O. fragrans and provide theoretical reference for the prediction and regulation of flowering time and genetic modification in O. fragrans.
Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) analysis
Developmental stages of sweet osmanthus flowers for SEM and TEM analysis were described as follows: stage 1 (S1), the outer bud scales unfurled and the inner bud scales still furled; S2, the bud became globular-shaped and the inside bracts covering the inflorescence was visible; S3, the inflorescence burst through bracts and the florets closely crowded; S4, initial flowering stage; S5, full flowering stage; S6, pollen-scattered stage. Observation on the cell size of adaxial and abaxial petal epidermal cells with SEM revealed that cell size of adaxial petal epidermal cells increased gradually (P < 0.05, by Duncan’s multiple range test) from S1 to S6 (Fig. 1a, Table S1), while that of abaxial petal epidermal cells enlarged gradually from S1 to S5, reached the peak at S5, and then greatly decreased from S5 to S6 (Fig. 1b, Table S1). These results were coincident with results in Gaillardia grandiflora , carnation  and T. caerulea , indicating that the cell division ceased and cell expansion occurred during floral opening process. However, in rose  and E. grandiflorum  cell division and cell expansion simultaneously appeared during this process. What’s more, based on TEM observation, it was found that vacuole occupied most area of adaxial petal epidermal cells in O. fragrans (Fig. 1c, Table S1). The same situation occurred in E. grandiflorum . These results indicated that the petal cell expansion was accompanied by the enlargement of vacuole. Vacuole size increased gradually from S1 to S6. At S2, vacuole size of cells was nearly two times as large as that at S1. With the cell development, vacuole size of cells at S6 was just one time larger than that at S2 (Fig. 1c, Table S1), suggesting that the growth rate from S1 to S2 was the highest during flower development.
Bud appearance with treatments
Before 23 °C or 19 °C treatment, the buds of all plants were at S1, and they were spindle-shaped, with unfurled outer bud scales and furled inner bud scales. Bud appearance remained the same under 23 °C treatment for 2 d, 4 d, and 6 d (H2, H4, H6), as well as the 19 °C treatment for 2 d and 4 d (L2 and L4). When treated with 19 °C treatment for 6 d (L6), the bud stage reached at S2, suggesting that the relatively low temperature, like 19 °C, helps ‘Yanhong Gui’ blossom in experimental condition.
Illumina sequencing and de novo assembly of sequence reads
To reveal the effect of 19 °C on the transcriptomic change in ‘Yanhong Gui’, we collected three biological replicate samples of H2, H4, H6, L2, L4, and L6 for reference transcriptome sequencing and RNA-seq analysis. Reference transcriptome sequencing generated 164,753,984 raw reads and 110,307,366 clean reads (Table S2). A total of 152,247 transcripts with average length of 751 nt were obtained and 96,920 unigenes with average length of 873 nt were yielded after all clean reads were assembled (Table S3).
Functional annotation and classification
The assembled unigenes were aligned to Non-Redundant Protein Sequence Database (NR), Nucleotide Sequence Database (NT), Swiss-Prot Database, Kyoto Encyclopedia of Genes and Genomes (KEGG) Database, Eukaryotic Orthologous Groups of proteins (KOG) Database, Gene Ontology (GO) Database, and InterPro Database, in order to obtain the putative annotations (Table S4). Totally, 61,654 (63.61%) unigenes were successfully annotated using at least one database, while 8150 (8.41%) unigenes were annotated using all seven databases. Total 57,721 (59.56%) of unigenes were annotated using NR database, 59.26% of which have high homology to sequences from Sesamum indicum (Fig. S1). In total, 37,284 unigenes were matched in the Swiss-Prot database, which is about 38.47% of all annotated unigenes (Table S4). There were 40,764 unigenes mapped into 134 KEGG pathways which can be divided into six large pathways including “cellular processes”, “environmental information processing”, “genetic information processing”, “human diseases”, “metabolism”, and “organismal systems” (Table S5). The pathways with the highest numbers of unigenes were “metabolic pathways” (Ko01100, 8561 unigenes, 21%), followed by “biosynthesis of secondary metabolites” (Ko01110, 4440 unigenes, 10.89%), “plant-pathogen interaction” (Ko04626, 1876 unigenes, 4.6%), and “plant hormone signal transduction” (Ko04075, 1444 unigenes, 3.54%) (Table S5). Next, GO analysis was performed and a total of 16,014 unigenes (16.52% of all annotated unigenes) were categorized into 53 GO terms under three main categories: biological process, cellular component, and molecular function (Fig. S2). Proteins related to “metabolic process”, “cellular process”, and “single-organism process” were enriched in biological processes. The “single-organism process” category defines a biological process involving only one organism, and the 4806 unigenes annotated as such may be related to flowering organisms specifically. In the cellular component category, the “cell”, “cell part”, and “membrane” were the most highly presented GO terms. Under the molecular function category, the “catalytic activity”, and “binding proteins” were the most enriched (Fig. S2). We also performed KOG analysis to evaluate the functions of the assembled unigenes and 43,496 unigenes were assigned 25 KOG classification (Fig. S3). Among the 25 KOG classification, the cluster of “general function prediction only” (11,753 unigenes, 27.02%) represented the largest group, followed by “signal transduction mechanisms” (6000 unigenes, 13.79%) and “function unknown” (4761, 10.95%). “Cell motility” has the smallest proportion, only accounted for 0.20% (Fig. S3). Additionally, 46,405 (47.88%) unigenes in the transcriptome library were annotated against InterPro database (Table S4).
RNA-sequencing (RNA-seq) and mapping to the reference transcriptome database
Eighteen cDNA libraries for RNA-seq analysis were respectively sequenced, generating 23.58–24.14 Mb raw reads (Table S2). Furthermore, we obtained 23.58–24.13 Mb clean reads in 18 cDNA libraries, with the clean data rate of 99.94–99.97%. Then, these clean reads were respectively mapped to the reference transcriptome database. Total mapped reads percentage and unique match percentage in the 18 libraries were similar. Total mapped reads percentage ranged from 86.39 to 88.44%, and unique match percentage ranged from 53.7 to 56.29% (Table S6). The correlation coefficients among biological replications are greater than 0.94 in all samples (Fig. 2, Table S7).
Screening and analysis of DEGs
We performed a detailed comparative analysis of the DEGs (P value ≤0.01 and |log2Ratio| ≥ 1) in three comparisons including L2 vs H2, L4 vs H4, and L6 vs H6. In the L2 vs H2 comparison, there were 7365 and 15,864 unigenes were down-regulated and up-regulated, respectively. Meanwhile, 8393 and 22,755 unigenes were down-regulated and up-regulated in the L4 vs H4 combination, respectively. In the L6 vs H6 comparison, 9551 unigenes were down-regulated and 22,011 unigenes were up-regulated, respectively (Fig. 3a). In total, there were more up-regulated unigenes than down-regulated unigenes during floral opening under relatively low temperature. Additionally, 7496 up-regulated DEGs were shared by three comparisons, while 4108, 6132, and 6913 up-regulated DEGs were only up-regulated in the L2 vs H2, L4 vs H4, and L6 vs H6 comparison, respectively (Fig. 3b). What’s more, 2968 down-regulated DEGs were shared by three comparisons, while 2954, 2189, and 3806 down-regulated unigenes were only down-expressed in the L2 vs H2, L4 vs H4, and L6 vs H6 comparison, respectively (Fig. 3c).
GO enrichment and KEGG pathway analysis of DEGs
In this work, GO analysis were used to classify the functions of the annotated DEGs under different temperatures treated for different number of days. In all comparisons, the most significantly enriched GO terms in biological process were “metabolic process”, “cellular process”, and “single-organism process”. In cellular component, five GO terms including “cell”, “cell part”, “membrane”, “membrane part”, and “organelle” were highly enriched. “Binding” and “catalytic activity” which were in the molecular function were significantly enriched in all comparison’s groups (Table 1). All DEGs were mapped to KEGG pathways in order to investigate the major pathways of them (Table 2). They were enriched in 132, 133, and 134 KEGG metabolic pathways in the L2 vs H2, L4 vs H4, and L6 vs H6 comparisons. In all the comparison group, “metabolic pathways”, “biosynthesis of secondary metabolites”, “plant-pathogen interaction”, “starch and sucrose metabolism”, and “plant hormone signal transduction” were the top five pathways with the largest number of DEGs (Table 2).
Identification of DEGs involved in the cell wall metabolism
The expansion of petal cells depends on the cell wall loosening and cellulose biosynthesis, soluble carbohydrate allocation, ion and water transport, and cytoskeleton rearrangement . Therefore, we screened the DEGs which may be involved in the cell wall metabolism under different temperatures. EXP and XTH, are two kinds of essential proteins involved in cell expansion, loosening and rearranging the cell wall fibers in growing tissues . In cucumber, EXPs consist of four sub families: α-EXP, β-EXP, EXP-like A and EXP-like B and were first discovered to loosen cell walls in pH dependent manner . In our study, four α-EXP and one β-EXP were significantly different expressed under 23 °C and 19 °C comparisons (Fig. 4, Table S8). Among them, four α-EXP were up-regulated under relatively low temperature (19 °C) and one β-EXP was down-regulated under the treatment of 19 °C.
Eleven unigenes which encoded as XTH protein were significantly expressed, seven of them (CL256.Contig3, CL3252.Contig2, CL10571.Contig1, Unigene1021, Unigene1631, Unigene24045, and Unigene25288) were up-regulated under the treatment of 19 °C and four (CL1547.Contig1, CL7741.Contig1, CL8777.Contig1, and Unigene26777) of them were down-regulated under 19 °C treatment (Fig. 4, Table S8).
The expression of unigenes involved in cell wall synthesis, modification or hydrolysis which probably participate in the process of flower opening were also assessed (Fig. 4, Table S8). There were three cellulose synthase (CES) genes, six xylosidase (XYL) genes, two pectin esterase (PE) genes, five polygalacturonase (PG) genes, and two pectate lyase (PL) genes with dramatic fold change in all three comparisons (Fig. 4, Table S8). The expression of all CESs and PLs genes increased when treated with relatively low temperature of 19 °C for different number of days. Most of the XYLs were up-regulated under relatively low temperatures except Unigene17352. Three PGs’ expression increased and two of them decreased when exposed to relatively low temperature. One PE gene was greatly up-regulated while one PE gene was down-regulated when the plants were exposed to 19 °C. Several genes associated with cell wall synthesis, modification or hydrolysis are up-regulated, suggesting that these genes may play vital roles in regulating petal cell expansion.
Aquaporins can facilitate the passage of water and/or small neutral solute fluxes across membranes and can be classified into four subclasses such as plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), nodulin-26-like intrinsic membrane proteins (NIPs) and small basic intrinsic proteins (SIPs) based on the sequence homology and cellular localization . Three kinds of aquaporins including two PIPs, two TIPs and three NIPs, were found in our transcriptomes and most of them were significantly up-regulated under relatively low temperature (Fig. 4, Table S8). However, only the expression level of two NIPs (Unigene22675, and Unigene40148) decreased under relatively low temperature. Dehydrins are hydrophilic, thermostable stress proteins and some of them were identified as osmotic stress-responsive protein which could facilitate the inflow and outflow of water . In our research, all the dehydrin genes were down-regulated when exposed to relatively low temperature (Fig. 4, Table S8).
Identification of DEGs involved in the phytohormone signal transduction pathway
There are some reports on phytohormones that can affect the floral opening , and therefore, we explored the DEGs involved in eight phytohormone signal transduction pathways. The levels of endogenous AUX were much higher in the basal petals and application of exogenous indole-3-butyric acid (IBA) can enlarge the cell length in chrysanthemum . Many unigenes involved in AUX signal transduction pathway were increasingly expressed with the flower opening of rose, suggesting that AUX signaling participated in petal development . In this study, nineteen unigenes were significantly differently expressed under different temperatures in AUX signal transduction pathway which included the auxin transporter protein (AUX1), auxin-induced protein/auxin-responsive protein (AUX/IAA), auxin response factor (ARF), auxin responsive GH3 family protein (GH3) and small auxin-up RNA family protein (SAUR) (Fig. 5). The majority of them had lower expression under relatively low temperature, especially three ARF members (CL661.Contig9, CL274.Contig1, and CL12908.Contig6), which decreased sharply as the treated days lengthened. Only one AUX1, two AUX/IAA, one ARF, and one GH3 were up-regulated when treated with relatively low temperature (Fig. 5, Table S9).
The phytohormone JA and its derivative methyl jasmonate (MeJA) can affect the flower opening of E. grandiflorum , rice , A. thaliana . The gene BIGPETALp (BPEp), which can be regulated by JA, may also play a role in JA-mediated petal growth . In JA signal transduction pathway, most of genes such as Coronatine-insensitive protein 1 (COI1), jasmonate ZIM domain protein (JAZ), and MYC2 were down-regulated in this study. Only two JAZ (Unigene43759 and CL8754.Contig3) genes and four MYC2 (CL37.Contig1, CL6240.Contig1, CL3114.Contig2, and CL2658.Contig2) genes were up-regulated under 19 °C treatment (Fig. 5, Table S9), suggesting that these unigenes may function as positive regulators of flower opening in O. fragrans (Fig. 5, Table S9). Only six TGACG-sequence-specific DNA-binding protein (TGA) genes in salicylic acid (SA) signal transduction pathway were expressed differently under different temperatures and two of them (CL2000.Contig1 and CL225.Contig1) were increased after 19 °C treatment (Fig. 5, Table S9).
Nineteen unigenes in the gibberellin acid (GA) signaling pathway were also differentially expressed under different temperatures (Fig. 5). Five of six GID1 were down-regulated under relatively low temperature, with the exception of Unigene22649, which was nearly increased eight thousand times after 2 d treatment of 19 °C. Three DELLA and five PIF3 genes were decreased while two DELLA and three PIF3 genes were increased after relatively low temperature treatment (Fig. 5, Table S9). Only four of twenty unigenes including two PYR/PYL and two ABF (ABA responsive element binding factor) genes in ABA signal transduction pathway were up-regulated after 19 °C treatment. However, the expression of one ABF (CL12556.Contig1) was decreased by two hundred times under the relatively low temperature at 19 °C compared with that under 23 °C (Fig. 5, Table S9).
Most of the unigenes in cytokinin (CTK) signal transduction pathway were down-regulated under relatively low temperature including all Arabidopsis histidine phosphotransfer Protein (AHP) genes, and all B-ARR members, especially the AHP (CL11134.Contig2), which was decreased severely after 6 d treatment of 19 °C (Fig. 5). However, most of the A-ARR members (CL7866.Contig2, Unigene18524, and CL8358.Contig1) were up-regulated when treated by relatively low temperature (Fig. 5, Table S9).
The effect of ETH on flower opening may be dose-dependent, species-dependent or cultivar-dependent . In ETH signal transduction pathway, the expression pattern of some ethylene receptor (ETR, Unigene14562), serine/threonine-protein kinase (CTR1, CL10648.Contig2), and ethylene insensitive 3 protein (EIN3, CL10633.Contig1) genes, were more complicated. For example, the expression level of Unigene14562 and CL10633.Contig1 were decreased when treated by 2 d of relatively low temperature, while they were increased when the treatment time was lengthened. Most of the genes in this pathway were down-regulated after the relatively low temperature’ treatment (Fig. 5, Table S9).
Cell elongation is also regulated by BR, which is a plant-specific steroid hormone and longitudinal cell expansion of BR mutants is greatly reduced in A. thaliana [33, 34]. Eight of fifteen unigenes in the BR signal transduction pathways were increased under 19 °C treatment, including all brassinosteroid receptor (BRI1), one BRI1 Associated receptor Kinase 1 (BSK), one TCH4, and one cyclin-D3 (CYCD3) (Fig. 5, Table S9). The expression level of BRI1, especially Unigene37674 and Unigene35684 were dramatically increased when treated with 19 °C. However, seven of fifteen unigenes in this signal pathway were decreased under 19 °C treatment. The expression level of BZR1/2 was eight-fold under 19 °C compared to 23 °C, and the expression level were continuously increased as the treatment time prolonged, suggesting the same role of BZR1/2 in petal development of O. fragrans and G. hybrida. Similarly, the expression level of TCH4 (CL8777.Contig2) was enhanced by relatively low temperature.
Identification of DEGs of transcription factors (TFs)
Many TFs changed dramatically under different temperatures, so we improved the screening criteria. The TFs with a |log2Ratio| ≥ 2 in any comparison were analyzed further. 79 TFs which were divided into eight gene family (Fig. 6, Table S10) were selected for further analysis because they may be involved in petal cell expansion [5, 30]. In the MYB TF family, only five out of twenty genes (Unigene36971, CL242.Contig2, CL10695.Contig2, Unigene11061, and Unigene30676) were up-expressed under relatively low temperature. Overexpression of the MYB62 gene results in a GA-deficient phenotype and can be partially alleviated with the application of exogenous GA . In our study, we did not find the homologous unigenes to MYB62 in O. fragrans, but some unigenes homologous to other gene family members, such as MYB4 (CL9450.Contig1) and MYB86 (Unigene19015), were expressed differently under different temperatures (Fig. 6, Table S10), indicating they may be involved in the regulation of relatively low temperature on flowering opening in O. fragrans.
In the bHLH TF family, three genes (Unigene30441, CL796.Contig2, and CL6865.Contig1) were up expressed after treated by 19 °C. Among them, especially the bHLH79, the expression level of bHLH79 isoform X2 (CL6865.Contig1) increased five hundred times when treated by 6 d of relatively low temperature. However, the expression level of bHLH 117 decreased sharply when exposed to 19 °C. None of the TF genes was up-regulated under relatively low temperature in the AP2 and NAC TF family. Plant-specific NAC family genes can regulate plant development, cell division, senescence, and responses to abiotic stress . In recent research, the N-terminal binding domain of RhNAC2 could be bound to the promoter of RhEXPA4 thus regulating of dehydration tolerance during the expansion of rose petals . The expression of all the members of NAC family decreased under the relatively low temperature (Fig. 6, Table S10), indicating that the suppression of NAC genes may facilitate the expression of EXPs to promote the elongation of flower petals in O. fragrans.
In the TCP TF family, two genes (CL9936.Contig2, and Unigene10400) were up-expressed under relatively low temperature treatment. miRNAs, such as miRNA319a, have also been found to be involved in petal elongation by regulating TCP4 which regulates AUX effects . Therefore, TCP4-like gene (CL9936.Contig2, and Unigene10400) increased with relatively low temperature treatment, implying they may regulate the petal development in O. fragrans (Fig. 6, Table S10).
Four out of nine genes were up-regulated in the MADS TF family, especially the AGAMOUS gene, which is the C-function genes, and its expression increased significantly with the 19 °C treatment. In the WRKY TF family, only two genes (CL11608.Contig1, and Unigene860) were up-regulated under relatively low temperature. In Zinc finger TF family, four out of eleven genes’ expression increased when exposed to 19 °C. Generally speaking, only a small proportion of TFs were up-regulated under relatively low temperature.
Ten differentially expressed unigenes comprising four genes involving in cell wall metabolism, four genes related to phytohormone signal transduction and two genes of TF were analyzed by qRT-PCR in order to confirm the results of RNA-Seq. All ten genes tested by qRT-PCR confirmed significant differential expression (P < 0.05) between 19 °C and 23 °C treatments (Fig. 7). Although fold changes of the ten genes between 19 °C and 23 °C treatments were not always similar in RNA-Seq and qRT-PCR results, the overall trend was consistent. Generally, the expression levels of CL10934.Contig2 (EXP), CL1473.Contig2 (XYL) and CL6542.Contig1 (PL) were significantly lower in 23 °C treatment than those in 19 °C treatment with qRT-PCR analysis, while the expression levels of other seven genes were higher in 23 °C treatment. Similarly, CL10934.Contig2, CL1473.Contig2 and CL6542.Contig1 were down-regulated, while other seven genes were up-regulated in 23 °C treatment, which verified the accuracy of RNA-Seq results.
The flowering time of once-flowering cultivars of O. fragrans is greatly affected by the relatively low temperature and they need a certain amount of low temperature accumulation at autumn before blooming . The sensitivity difference to the requirement of relatively low temperature among different cultivars leads to the variation of flowering time . In this study, 19 °C treatment in experimental condition for 6 d is enough for flower opening in the cultivar ‘Yanhong Gui’. The early flowering cultivar such as ‘Zao Yingui’ would blossom when the daily minimum temperature lower than 21 °C last for after 5–8 days, while the late blossoming cultivar, such as ‘Wan Yingui’, needs 5–8 days of the daily minimum temperature lower than 19 °C .
Flower opening is accompaning with cell division and cell expansion in petal cells . are up-regulated during floral open process in sandersonia , Mirabilis jalapa , petunia , carnation , and wintersweet . Similarly, during flower opening of O. fragrans under relatively low temperature, four α-EXP were up-regulated in this study (Fig. 4, Table S8). It was also reported that β-EXP is more efficient in disrupting cell wall polymers than α-EXP , however, one β-EXP was down-regulated under 19 °C treatment. The result is coincidence with the results in rose , suggesting that β-EXP may be involved in petal development in a special manner different from α-EXP. As to XTH, four rose XTH genes (RbXTH3, RbXTH5, RbXTH6, and RbXTH12) were up-regulated during flower opening and associated with petal movement . Rice OsXTH12 gene can cause rice glume-unclosing after anthesis under high temperature conditions . In our experiment, we cannot found the homologous genes of XTH3, XTH5, XTH6 and XTH12, but other family members were discovered. Most of XTH genes were up-regulated under relatively low temperature, indicating that they may be involved in the flower opening process of O. fragrans and different family members may have different functions.
These cell wall metabolism processes were also influenced by multiple kinds of phytohormones via transcription factors . The expression patterns of phytohormone pathway genes are different under different temperatures treated by different days in our research. In A. thaliana, both the AUX and JA signal transduction function synergistically during petal development . And ARF8 can interacts with BIGPETALp (BPEp) to regulate cell expansion in petals  and can trigger the expression of MYB21 and MYB24 to induce the production of JA to promote petal growth . However, most DEGs involved in AUX and JA signal pathways were down-regulated (Fig. 5, Table S9), suggesting that the regulation of relatively low temperature on flower opening in O. fragrans doesn’t depend on AUX and JA signaling.
GA can regulate seed germination, stem elongation, petal growth and flowering , and ABA has an antagonistically role such as in floral transition and fruit development . The same circumstance occurred in the G. hybrida petal growth . Compared to unigenes involved in ABA signal transduction pathway, there were much more unigenes in GA signal transduction pathway which increased sharply under the treatment of 19 °C. This result was consistent with the result in rose . However, whether this kind of antagonistic role exists in petal growth of O. fragrans still needs further demonstration.
ETH can promote flower opening in carnation, Phalaenopsis orchids, and petunia, but it inhibits flower opening in several rose cultivars . ETH can induce the expression of DELLA genes via EIN3–3 when it functions as inhibitors during flower opening in rose . In our study, EIN3 was down-regulated under relatively low temperature (Fig. 5, Table S9), suggesting that ETH may function negatively on flower opening in response to relatively low temperature in O. fragrans.
In G. hybrida, BL shows more efficient effect on lengthening the cells in the middle and basal regions of petals than GA . But the expression level of BES1/BZR1 was decreased after only 0.5 h BR treatment, indicating its vital role in petal cell expansion . In our study, expression level of the homologous unigene of BZR1 (Unigene45694), decreased dramatically when exposed to relatively low temperature (Fig. 5, Table S9), suggesting the same role in petal development of O. fragrans as that in G. hybrida.
The comprehensive transcriptomic dataset relating to flower opening of O. fragrans under relatively low temperature (19 °C) treatment was performed using RNA-seq technology, and a group of DEGs, which may regulate petal growth through pairwise comparison analysis of DEGs between different temperature treatments, had been identified, including unigenes involved in cell wall metabolism, phytohormone signal transduction pathway, and TFs. Overall, the resources generated in this study would lay foundation on fully revealing the molecular mechanism of relatively-low temperature-regulated flower opening of O. fragrans and provide theoretical reference for the prediction and regulation of flowering time and genetic modification in O. fragrans. Our future research will be focused on understanding the biological function of these candidate important DEGs and revealing how these DEGs response to relatively low temperature in O. fragrans.
Plant materials and treatments
Plants of O. fragrans ‘Yanhong Gui’ were potted and grown in the resource nursery of Zhejiang Agriculture and Forestry University. All plant materials were owned by Zhejiang Agriculture and Forestry University. The development of flower opening process are classified into S1-S6 as described in this paper . When the floral development stage was at S1, potted plants were moved into the growth cabinet with the temperature of 19 °C (relatively low temperature) or 23 °C (relative high temperature), respectively. The relative humidity was about 50–60% and the photoperiod was 12 h light/12 h dark regime with an illumination of 80 μmol·m− 2·s− 1. Petals were harvested every 2 days after treatment until the developmental stage of plants under 19 °C condition reached S2. The samples were named as H2, H4, H6 (treated for 2 d, 4 d, and 6 d under 23 °C treatment, respectively) and L2, L4, L6 (treated by 2 d, 4 d, and 6 d under 19 °C treatment, respectively). All samples were collected at 10:00 am in order to avoid the influence of diurnal rhythm. Three biological replicate samples were collected to generate 18 samples in total. In addition, petals at different developmental stages (S1-S6) were respectively sampled for SEM and TEM analysis.
SEM and TEM analysis
According to previous studies [25, 48], SEM and TEM were carried out and were observed and photographed in Hitachi Model TM-1000 SEM and Model H-7650 respectively in Zhejiang University. Cell size of adaxial and abaxial petal epidermal cells, and vacuole size were manually measured using Image J software (http://rsb.info.nih.gov/ij/, NIH, MD, USA).
RNA extraction, cDNA library preparation and sequencing
RNA isolation and RNA purification of 18 samples were carried out as described previously . The cDNA library preparation and sequencing of 18 samples of O. fragrans were respectively performed as described previously  for RNA-seq analysis. What’s more, a RNA pool mixed from 18 RNA samples was used to construct a cDNA library for reference transcriptome sequencing. The Illumina sequencing was performed at the Beijing Genomics Institute (BGI) (Shenzhen, China) (http://www.genomics.cn/index.php) according to the manufacturer’s instructions (Illumina Inc., San Diego, CA, USA). A total of 18 sets of raw reads were obtained for RNA-seq analysis, corresponding to H2–1, H2–2, H2–3, H4–1, H4–2, H4–3, H6–1, H6–2, H6–3, L2–1, L2–2, L2–3, L4–1, L4–2, L4–3, L6–1, L6–2, and L6–3, while raw reads were obtained for reference transcriptome corresponding to ALL.
Transcriptome data processing and analysis
Raw reads were cleaned by removing the sequence of adaptor, high content of unknown bases and low-quality reads before downstream analysis to decrease data noise. As O. fragrans did not have an appropriate reference genome sequence, the Trinity method  was used to de novo assemble all clean reads and the Tgicl method  was used to cluster the transcripts to obtain non-redundant unigenes. The assembled unigene sequences were aligned by Blastn  to nucleotide database NT, aligned by Blastx  or Diamond  to protein databases NR, Swiss-Prot, KEGG, and KOG, and aligned by InterProScan5  to protein database InterPro to get the annotations. With the NR annotation, GO annotations of unigenes were obtained using Blast2GO software . FPKM (reads mapped per 1000 bp per million sequenced reads) method was used in calculated expression level . Comparisons of FPKM between samples treated under different temperatures for the same number of days (L2 vs H2, L4 vs H4, and L6 vs H6 respectively) were performed for each unigene. To identify the DEGs in the two samples, the Audic and Claverie method was employed . Unigene with a P value ≤0.01 and a |log2Ratio| ≥ 1 were considered significant DEGs. All DEGs were mapped to each term of the GO and KEGG databases. For the analysis of transcription factors, |log2Ratio| ≥ 2 were marked to be significantly different between the two samples.
Quantitative real-time PCR (qRT-PCR) analysis
The extraction of total RNA, first-strand cDNA synthesis, and qRT-PCR were performed as previously described . Ten DEGs were randomly chosen and the primer sequences were listed in Table S11. Gene relative expression levels of those ten genes were normalized relatively to the expression level of an internal control OfACT  using the 2-△Ct method.
The statistical analysis of the size of adaxial petal epidermal cells, abaxial petal epidermal cells, and vacuole, as well as relative expression level and FPKM analysis was performed by one-way analysis of variance (ANOVA) using SPSS software version 18.0 (SPSS Inc., Chicago, IL, USA). Duncan’s multiple range test was employed and differences of P < 0.05 were considered significant.
Availability of data and materials
The RNA-Seq datasets used in the current study are available on the NCBI Short Read Archive Project - PRJNA642321 (http://www.ncbi.nlm.nih.gov/bioproject/PRJNA642321). All RNA-seq reads generated by this study are publicly available at the NCBI SRA under accession numbers SRR12109992-SRR12110009.
Auxin response factor
Arabidopsis histidine phosphotransfer protein
Auxin-induced protein/auxin-responsive protein
Auxin transporter protein
Coronatine-insensitive protein 1
Differentially expressed genes
Ethylene insentitive 3 protein
Reads mapped per 1000 bp per million sequenced reads
Auxin responsive GH3 family protein
Jasmonate ZIM domain protein
Kyoto Encyclopedia of Genes and Genomes
Eukaryotic Orthologous Groups of proteins
Nodulin-26-like intrinsic membrane protein
Non-Redundant Protein Sequence Database
Nucleotide Sequence Database
Plasma membrane intrinsic protein
Quantitative real-time PCR
Small auxin-up RNA family protein
Scanning electron microscopy
Small basic intrinsic protein
TGACG-sequence-specific DNA-binding protein
Transmission electron microscopy
Tonoplast intrinsic protein
Yang K. Chinese Osmanthus. Beijing: China Forestry Publication House; 2011.
Xiang Q, Liu Y. An illustrated monograph of the sweet osmanthus cultivars in China. Hangzhou: Zhejiang Science and Technology Press; 2008.
Wang Y, Zhang C, Fu J, Zhao H. Progresses on flower bud differentiation and flower opening in Osmanthus fragrans. J Zhejiang A & F Univ. 2016;33(2):340–7.
Dong L, Wang X, Ding Y. Study on the Osmanthus fragrans blooming season phenology. J Nanjing Forestry Univ (Nat Sci Ed). 2014;38(S1):51–6.
Wang J, Wang H, Ding L, Song A, Shen F, Jiang J, Chen S, Chen F. Transcriptomic and hormone analyses reveal mechanisms underlying petal elongation in Chrysanthemum morifolium 'Jinba'. Plant Mol Biol. 2017;93(6):593–606.
Varaud E, Brioudes F, Szecsi J, Leroux J, Brown S, Perrot-Rechenmann C, Bendahmane M. AUXIN RESPONSE FACTOR8 regulates Arabidopsis petal growth by interacting with the bHLH transcription FACTOR BIGPETALp. Plant Cell. 2011;23(3):973–83.
Christiaens A, De KE, Pauwels E, De RJ, Gobin B, Van Labeke MC. Suboptimal light conditions influence source-sink metabolism during flowering. Front Plant Sci. 2016;7:249.
Cosgrove DJ. Wall structure and wall loosening. A look backwards and forwards. Plant Physiol. 2001;125(1):131–4.
Zenoni S, Reale L, Tornielli GB, Lanfaloni L, Porceddu A, Ferrarini A, Moretti C, Zamboni A, Speghini A, Ferranti F. Downregulation of the Petunia hybrida alpha-expansin gene PhEXP1 reduces the amount of crystalline cellulose in cell walls and leads to phenotypic changes in petal limbs. Plant Cell. 2004;16(2):295–308.
Ma J, Li Z, Wang B, Sui S, Li M. Cloning of an expansin gene from Chimonanthus praecox flowers and its expression in flowers treated with ethephon or 1-methylcyclopropene. Hortscience. 2012;47(10):1472–7.
Panavas T, Reid PD, Rubinstein B. Programmed cell death of daylily petals: activities of wall-based enzymes and effects of heat shock. Plant Physiol Biochem. 1998;36(5):379–88.
Tarpley L, Sassenrath GF. Carbohydrate profiles during cotton floral bud (square) development. J Agron Crop Sci. 2006;192(5):363–72.
Ichimura KK, Koketsu M, Shimamura M, Ito A. Identification of pinitol as a main sugar constituent and changes in its content during flower bud development in carnation (Dianthus caryophyllus L.). J Plant Physiol. 1998;152(5):363–7.
Yamada K, Ichimura K. Determination of subcellular concentrations of soluble carbohydrates in rose petals during opening by nonaqueous fractionation method combined with infiltration–centrifugation method. Planta. 2009;230(6):1115–27.
Ichimura K, Kohata K, Yamaguchi Y, Douzono M, Ikeda H, Koketsu M. Identification of L-inositol and scyllitol and their distribution in various organs in chrysanthemum. Biosci Biotechnol Biochem. 2000;64(4):865–8.
Norikoshi R, Imanishi H, Ichimura K. Changes in cell number, osmotic potential and concentrations of carbohydrates and inorganic ions in Tweedia caerulea during flower opening. J Jpn Soc Hortic Sci. 2013;82(1):51–6.
Norikoshi R, Kohata K, Shimizu-Yumoto H, Goto R, Ichimura K. Identification of soluble carbohydrates and their subcellular concentrations in petals during flower opening in Eustoma grandiflorum. Horticult J. 2016;85(3):238–47.
Krizek BA, Anderson JT. Control of flower size. J Exp Bot. 2013;64(6):1427–37.
van Doorn WG, Kamdee C. Flower opening and closure: an update. J Exp Bot. 2014;65(20):5749–57.
Brioudes F, Joly C, Szécsi J, Varaud E, Leroux J, Bellvert F, Bertrand C, Bendahmane M. Jasmonate controls late development stages of petal growth in Arabidopsis thaliana. Plant J. 2009;60(6):1070–80.
Huang G, Han M, Yao W, Wang Y. Transcriptome analysis reveals the regulation of brassinosteroids on petal growth in Gerbera hybrida. Peerj. 2017;5(5):e3382.
Li L, Zhang W, Zhang L, Li N, Peng J, Wang Y, Zhong C, Yang Y, Sun S, Liang S. Transcriptomic insights into antagonistic effects of gibberellin and abscisic acid on petal growth in Gerbera hybrida. Front Plant Sci. 2015;6:168.
Koning RE. The roles of plant hormones in the growth of the corolla of Gaillardia grandiflora (Asteraceae) ray flowers. Am J Bot. 1984;71(1):1–8.
Kenis JD, Silvente ST, Trippi VS. Nitrogen metabolism and senescence-associated changes during growth of carnation flowers. Physiol Plantarum. 1985;65(4):455–9.
Yamada K, Norikoshi R, Suzuki K, Nishijima T, Imanishi H, Ichimura K. Cell division and expansion growth during rose petal development. J Jpn Soc Hortic Sci. 2009;78(3):356–62.
Norikoshi R, Shibata T, Ichimura K. Cell division and expansion in petals during flower development and opening in Eustoma grandiflorum. Horticult J. 2016;85(2):154–60.
Marowa P, Ding A, Kong Y. Expansins: roles in plant growth and potential applications in crop improvement. Plant Cell Rep. 2016;35(5):949–65.
Maurel C, Verdoucq L, Luu DT, Santoni V. Plant aquaporins: membrane channels with multiple integrated functions. Annu Rev Plant Biol. 2008;59:595–624.
Kumar M, Lee SC, Kim JY, Kim SJ, Aye SS, Kim SR. Over-expression of dehydrin gene, OsDhn1, improves drought and salt stress tolerance through scavenging of reactive oxygen species in rice (Oryza sativa L.). J Plant Biol. 2014;57(6):383–93.
Yu H, Wan H, Cheng T, Jia W, Yang W, Pan H, Zhang Q. Comparative RNA-seq analysis of transcriptome dynamics during petal development in Rosa chinensis. Sci Rep. 2017;7:43382.
Ochiai M, Matsumoto S, Yamada K. Methyl jasmonate treatment promotes flower opening of cut Eustoma by inducing cell wall loosening proteins in petals. Postharvest Biol Technol. 2013;82(21):1–5.
He YM, Lin YJ, Zeng XC. Dynamic changes of jasmonic acid biosynthesis in rice florets during natural anthesis. Acta Agron Sin. 2012;38(10):1891–9.
Gallego-Bartolomé J, Blázquez MA. Molecular mechanism for the interaction between gibberellin and brassinosteroid signaling pathways in Arabidopsis. Proc Natl Acad Sci U S A. 2012;109(33):13446–51.
Li QF, Wang C, Jiang L, Li S, Sun SS, He JX. An interaction between BZR1 and DELLAs mediates direct signaling crosstalk between brassinosteroids and gibberellins in Arabidopsis. Sci Signal. 2012;5(244):ra72.
Devaiah BN, Madhuvanthi R, Karthikeyan AS, Raghothama KG. Phosphate starvation responses and gibberellic acid biosynthesis are regulated by the MYB62 transcription factor in Arabidopsis. Mol Plant. 2009;2(1):43–58.
Olsen AN, Ernst HA, Leggio LL, Skriver K. NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005;10(2):79–87.
Dai F, Zhang C, Jiang X, Kang M, Yin X, Lü P, Zhang X, Zheng Y, Gao J. RhNAC2 and RhEXPA4 are involved in the regulation of dehydration tolerance during the expansion of rose petals. Plant Physiol. 2012;160(4):2064–82.
Nag A, King S, Jack T. miR319a targeting of TCP4 is critical for petal growth and development in Arabidopsis. Proc Natl Acad Sci U S A. 2009;106(52):22534–9.
Wang Y, Hu Y. Study of temperature influence on blooming of Osmanthus fragrans in Shanghai. In: Advances in Ornamental Horticulture of China: 2011: China Forestry Press; 2011. p. 768–71.
O'Donoghue EM, Somerfield SD, Heyes JA. Organization of cell walls in Sandersonia aurantiaca floral tissue. J Exp Bot. 2002;53(368):513–23.
Gookin TE, Hunter DA, Reid MS. Temporal analysis of alpha and beta-expansin expression during floral opening and senescence. Plant Sci. 2003;164(5):769–81.
Harada T, Torii Y, Morita S, Onodera R, Hara Y, Yokoyama R, Nishitani K, Satoh S. Cloning, characterization, and expression of xyloglucan endotransglucosylase/hydrolase and expansin genes associated with petal growth and development during carnation flower opening. J Exp Bot. 2011;62(2):815–23.
Singh AP, Dubey S, Lakhwani D, Pandey SP, Khan K, Dwivedi UN, Nath P, Sane AP. Differential expression of several xyloglucan endotransglucosylase/hydrolase genes regulates flower opening and petal abscission in roses. Aob Plants 2013;5(1):plt030.
Fu C, Wang F, Liu W, Liu D, Li J, Zhu M, Liao Y, Liu Z, Huang H, Zeng X. Transcriptomic analysis reveals new insights into high-temperature-dependent glume-unclosing in an elite rice male sterile line. Front Plant Sci. 2017;8:112.
Reeves PH, Ellis CM, Ploense SE, Wu MF, Yadav V, Tholl D, Chételat A, Haupt I, Kennerley BJ, Hodgens C. A regulatory network for coordinated flower maturation. PLoS Genet. 2012;8(2):e1002506.
Hedden P, Thomas SG. Gibberellin biosynthesis and its regulation. Biochem J. 2012;444(1):11–25.
Razem FA, Baron K, Hill RD. Turning on gibberellin and abscisic acid signaling. Curr Opin Plant Biol. 2006;9(5):454–9.
Zhong K, Li X, Le X, Kong X, Zhang Z. MoDnm1 dynamin mediating peroxisomal and mitochondrial fission in complex with MoFis1 and MoMdv1 is important for development of functional appressorium in Magnaporthe oryzae. PLoS Pathog. 2016;12(8):e1005823.
Zhang C, Wang Y, Fu J, Bao Z, Zhao H. Transcriptomic analysis and carotenogenic gene expression related to petal coloration in Osmanthus fragrans ‘Yanhong Gui’. Trees - Structure Function. 2016;30(4):1207–23.
Zhang C, Wang Y, Fu J, Dong L, Gao S, Du D. Transcriptomic analysis of cut tree peony with glucose supply using the RNA-Seq technique. Plant Cell Rep. 2014;33(1):111–29.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
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.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33:W116–20.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Mortazavi A, W BA, M K, S L, W B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.
Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.
Zhang C, Fu J, Wang Y, Bao Z, Zhao H. Identification of suitable reference genes for gene expression normalization in the quantitative real-time PCR analysis of sweet osmanthus (Osmanthus fragrans Lour.). PLoS ONE. 2015;10(8):e0136355.
Consent to publish
This research was supported by Zhejiang Provincial Natural Science Foundation of China (Grant Nos. LY19C160006, LY19C160002, and LQ16C160003), Zhejiang Provincial Major Program of New Cultivar Breeding (Grant No. 2016C02056–12), and China Scholarship Council (201808330145 and 201808330146). The funding organizations provided the financial support to the research projects, but were not involved in the design of the study, data collection, analysis of the data, or the writing of the manuscript.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A species-based distribution of BLASTX matches for unigenes from the reference transcriptome of Osmanthus fragrans. We used all the plant proteins in the NCBI NR database to perform the homology search and for each sequence we selected the closest match for analysis.
GO classification of unigenes from the reference transcriptome of Osmanthus fragrans. Results are summarized under three main GO categories: biological process, cellular component and molecular function. The right x-axis indicates the number of genes in the same category.
KOG classification of Osmanthus fragrans unigenes from the reference transcriptome. From a total of 96,920 de novo assembled unigenes, 43,496 transcripts with significant homologies in the KOG database were classified into 25 KOG categories.
The size of adaxial petal epidermal cells, abaxial petal epidermal cells, and vacuole of adaxial petal epidermal cells in sweet osmanthus flowers at different developmental stages
Summary of Illumina transcriptomic sequencing.
Summary of assembly quality of transcripts and unigenes.
Annotation results of the assembled unigenes from the reference transcriptome against public databases.
KEGG pathway annotation of Osmanthus fragrans unigenes obtained from the reference transcriptome. There were 40,764 unigenes mapped into 134 KEGG pathways.
Number of total clean reads of the 18 samples mapped to reference sequences.
Pearson correlation coefficients of six samples each with three biological replicates.
DEGs involved in the cell wall metabolism under different temperatures treated by different days in Osmanthus fragrans.
DEGs involved in the phytohormone signal transduction pathway under different temperatures treated by different days in Osmanthus fragrans.
DEGs involved in TFs under different temperatures treated by different days in Osmanthus fragrans.
Primer sequences and amplicon characteristics.
About this article
Cite this article
Fu, J., Zhang, C., Liu, Y. et al. Transcriptomic analysis of flower opening response to relatively low temperatures in Osmanthus fragrans. BMC Plant Biol 20, 337 (2020). https://doi.org/10.1186/s12870-020-02549-3
- Flower development
- Petal epidermal cells
- Sweet osmanthus
- RNA-seq analysis