Skip to main content

Integrated metabolic profiling and transcriptome analysis of pigment accumulation in Lonicera japonica flower petals during colour-transition



Plants have remarkable diversity in petal colour through the biosynthesis and accumulation of various pigments. To better understand the mechanisms regulating petal pigmentation in Lonicera japonica, we used multiple approaches to investigate the changes in carotenoids, anthocyanins, endogenous hormones and gene expression dynamics during petal colour transitions, i.e., green bud petals (GB_Pe), white flower petals (WF_Pe) and yellow flower petals (YF_Pe).


Metabolome analysis showed that YF_Pe contained a much higher content of carotenoids than GB_Pe and WF_Pe, with α-carotene, zeaxanthin, violaxanthin and γ-carotene identified as the major carotenoid compounds in YF_Pe. Comparative transcriptome analysis revealed that the key differentially expressed genes (DEGs) involved in carotenoid biosynthesis, such as phytoene synthase, phytoene desaturase and ζ-carotene desaturase, were significantly upregulated in YF_Pe. The results indicated that upregulated carotenoid concentrations and carotenoid biosynthesis-related genes predominantly promote colour transition. Meanwhile, two anthocyanins (pelargonidin and cyanidin) were significantly increased in YF_Pe, and the expression level of an anthocyanidin synthase gene was significantly upregulated, suggesting that anthocyanins may contribute to vivid yellow colour in YF_Pe. Furthermore, analyses of changes in indoleacetic acid, zeatin riboside, gibberellic acid, brassinosteroid (BR), methyl jasmonate and abscisic acid (ABA) levels indicated that colour transitions are regulated by endogenous hormones. The DEGs involved in the auxin, cytokinin, gibberellin, BR, jasmonic acid and ABA signalling pathways were enriched and associated with petal colour transitions.


Our results provide global insight into the pigment accumulation and the regulatory mechanisms underlying petal colour transitions during the flower development process in L. japonica.


Different plant lineages have adopted various mechanisms of flower colour determination to fulfil the requirements for pollinator attraction or non-pollinator-related traits [1]. Due to the importance of colour formation in angiosperms, especially in ornamentals, the biosynthesis pathways of pigments in colour formation have been extensively reported [2,3,4,5,6]. Three major classes of pigments, flavonoids/anthocyanins, carotenoids and chlorophylls, are distributed ubiquitously in plants. Among them, pigmented flavonoids, mainly anthocyanins, are major pigmentation compounds in flowering plants [5]. Carotenoids also widely participate in the yellow-to-red colouration of flowers [7]. Chlorophylls, which are a class of essential photosynthetic components, exist in almost all plants and are mainly involved in the formation of green colour in flowers [8]. Moreover, anthocyanins and carotenoids often coexist simultaneously, and their combination causes diversity in flower colour [9].

As one of the main subgroups of flavonoids, anthocyanins primarily cause the colour formation of red, orange, blue and purple flower colours [5]. The specific colouration of anthocyanin is also greatly dependent on various moieties, co-pigments and pH [9]. Anthocyanins in plants are synthesized via the flavonoid pathway, and their biosynthetic pathway and genes have been well characterized [5, 9, 10]. It was reported that anthocyanin synthesis shares the same upstream pathway as the formation of anthocyanidins (cyanidin, delphinidin, and pelargonidin), followed by specific downstream branch for anthocyanin modification. In this comprehensive synthesis process, the formation of anthocyanidins by the catalysis of anthocyanidin synthase (ANS; also known as leucoanthocyanidin dioxygenase, LDOX) is an important node that leads to flavonoid flux into the anthocyanin branch.

Carotenoids, a kind of C40 isoprenoids, are distributed in some flowers and provide distinct colours ranging from yellow/orange to red [3, 4]. In the initial steps of the carotenoid biosynthetic pathway, the key enzymes have been well characterized, including phytoene synthase (PSY), phytoene desaturase (PDS) and ζ-carotene desaturase (ZDS) [11,12,13]. Although the biosynthetic pathways of anthocyanins and carotenoids are well established, the balance of the expression dynamics of anthocyanins and carotenoids-related genes during petal colour transition in a single flower remains poorly understood.

Changes in flower colour are comprehensively regulated by physiological changes and transcriptional level fluctuations in related genes. To date, high-resolution mass spectrometry (MS)-based metabolomics represents an effective technique to detect the accumulation and dynamic changes of metabolites [14,15,16]. Furthermore, transcriptome analysis has developed into a powerful approach that provides abundant sequence resources to study the mechanisms regulating flower colour formation [17, 18]. To identify pigment accumulation, endogenous hormone changes and related gene fluctuations in petal colour transitions, global analysis of the metabolome combined with the transcriptional levels of pigment biosynthesis genes is required.

Lonicera japonica Thunb. is known as “gold and silver flower” in China and is widely cultivated in East Asian countries [19, 20]. It has excellent ornamental properties due to the dynamic petal colours of every single flower and provides plant materials to uncover the molecular mechanisms of petal colour transition [21,22,23,24]. In our study, we reported the metabolic profiling of carotenoids, anthocyanins, endogenous hormones and gene expression dynamics in L. japonica petals at various stages (i.e., green flower bud, white flower and yellow flower) using integrated analyses of the metabolome, physiology, and transcriptome. With this extensive analysis of multiple data points in L. japonica petals, we reveal changes in the key pigments, hormones, and related biosynthesis genes that are associated with petal colour transitions.


Plant materials sampling and color detection

According to a previous study [21], flower buds and opening flowers of Lonicera japonica Thunb. were collected from the Beijing Botanical Garden, Beijing, China, with permission. The colour of buds/flowers gradually changed from green to white and then to yellow during floral development. Three developmental stages of petal, i.e., green bud petals (GB_Pe), white flower petals (WF_Pe) and yellow flower petals (YF_Pe), were selected to perform pigment metabolome, transcriptome and plant hormone analyses. For each biological replication, petals were dissected from 3 to 6 uniform flowers, weighed, sampled, frozen immediately in liquid nitrogen and stored at − 80 °C until use. Meanwhile, the colour index of petals was measured using a CR-400 chroma meter (Konica Minolta Sensing Inc., Osaka, Japan). Hunter parameters of L*, a* and b* were mainly used according to the CIELAB colour model.

Carotenoid extraction and quantification

Petal samples from GB_Pe, WF_Pe and YF_Pe were used for carotenoid extraction. Three biological replicates were performed for each developmental stage. Each replicate (~ 1 g of fresh weight) was freeze-dried and ground into fine powder, and then about 100 mg of powder was dissolved in 1 mL of a solution of n-hexane: acetone: ethanol (2:1:1 (V:V:V)). The solution sample was vortexed (30 s), followed by ultrasound-assisted extraction at room temperature for 20 min, and centrifuged at 12,000 g for 5 min to collect the supernatant. The extraction steps above were repeated, and the supernatants from the two centrifugations were combined. Subsequently, the combined supernatant was evaporated to dryness using a nitrogen gas stream and then reconstituted in 200 μL of (acetonitrile:methanol = 3:1 (V:V)): methyl tert-butyl ether = 85:15 (V:V). Finally, the solution was centrifuged at 12,000 g for 2 min to collect the supernatant for LC-MS/MS analysis. To monitor the stability of the LC-MS/MS analytical conditions, a quality control (QC) sample was used with equal mixing of all measured samples and was run at intervals during the assay. Stock solutions of standards purchased from Sigma-Aldrich (St. Louis, MO, USA) or Olchemim Ltd. (Olomouc, Czech Republic) were prepared at a concentration of 10 mg/mL in HPLC-grade acetonitrile (ACN) and then diluted with ACN to working solutions. For each carotenoid, five successive concentration gradients were used to plot the standard curve.

The prepared sample solutions were analysed by an LC-APCI-MS/MS system (UPLC, Shim-pack UFLC SHIMADZU CBM30A system,; MS, Applied Biosystems 6500 Triple Quadrupole, equipped with an APCI Turbo Ion-Spray interface and operated in positive ion mode and controlled by Analyst 1.6.3 software (Applied Biosystems Company, Framingham, MA, USA) at Wuhan Metware Biotechnology Co., Ltd. (Wuhan, China). The APCI source operation parameter settings and multiple reaction monitoring (MRM) transitions were performed according to a previous study [25]. Finally, the carotenoid levels were calculated according to the following formula: carotenoid content (μg/g) = A*(B/1000)/C (A: the carotenoid concentration calculated by the standard curve using chromatographic peak area integrals, μg/mL; B: the volume of solution for redissolution, μL; C: the dry weight of plant tissue powder, g). SPSS 16.0 was used for the analysis of variance (ANOVA). The differences in means were compared using Duncan’s multiple range test.

Anthocyanin extraction and analyses

Nine petal samples (each sample was parallel to the above sample in the carotenoid analyses) were used for anthocyanin analyses. The freeze-dried petals were crushed into powder by a mixer mill (MM 400, Retsch). Each powder sample weighing ~ 100 mg was then extracted overnight at 4 °C with 70% (V/V) aqueous methanol and centrifuged at 10,000 g for 10 min. The supernatant extract was absorbed using a CNWBOND Carbon-GCB SPE Cartridge (250 mg, 3 mL; ANPEL, Shanghai, China) and filtered through a 0.22 μm SCAA-104 membrane (ANPEL, Shanghai, China). To examine the precision and repeatability of the instrumental assay system and analysis process, the QC sample was prepared by blending all of the samples equally and inserting test samples at intervals. The prepared samples were analysed by an LC-ESI-MS/MS system at Wuhan Metware Biotechnology Co., Ltd. The HPLC conditions and ESI source operation parameters were set according to previous studies [26, 27].

Qualitative and quantitative analyses of anthocyanins were basically consistent with the analyses of carotenoids above. Specifically, anthocyanin identification was performed based on the MWDB database (Wuhan MetWare Biotechnology Co., Ltd.) and publicly available metabolite databases following standard procedures if the standards were unavailable. Anthocyanin with a statistically significant difference in content was determined with thresholds of fold change ≥1.6 or ≤ 0.625, P-value < 0.05, and variable importance in projection (VIP) ≥ 0.8.

Determination of various hormones during flower development

The contents of indoleacetic acid (IAA), zeatin riboside (ZR), gibberellic acid (GA3), brassinosteroid (BR), methyl jasmonate (MeJA) and abscisic acid (ABA) from GB_Pe, WF_Pe and YF_Pe were measured by an indirect ELISA technique. The extraction, purification and determination of each hormone were performed according to previous study [28, 29] and the instructions of the corresponding kit produced by Shanghai Enzyme-linked Biotechnology Co., Ltd. (Shanghai, China). Each stage was prepared with three biological replicates. All data were analysed by ANOVA using SPSS 16.0, and the differences in means were compared by Duncan’s multiple range test.

RNA isolation, cDNA library construction and sequencing

Petal samples of GB_Pe, WF_Pe and YF_Pe were collected and three independent biological replicates were used. Total RNA was extracted from plant tissue using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol and digested with DNase I (Takara, Dalian, China). The quality and purity of total RNA were evaluated by stringent RNA quality control. cDNA library construction and sequencing were performed by Annoroad Gene Technology (Beijing, China). Each constructed cDNA library (10 ng) was subjected to paired-end 150 bp sequencing on an Illumina HiSeq™ 4000 system (San Diego, CA, USA) according to the manufacturer’s instructions.

Data assembly and annotation

The raw reads were filtered to remove adapter-polluted reads, low-quality reads, and reads with more than 5% ambiguous nucleotides. These clean reads with high-quality were subjected to the following analyses. Trinity software [30] was used to perform the de novo transcriptome assembly with default parameter values.

The assembled unigenes were annotated by homology search to publicly accessible databases using local BLAST programmes (version 2.2.28) with a significance threshold of E < 1e-5. Meanwhile, all unigenes were analysed with Blast2GO (version 3.0.8) to obtain the gene ontology (GO) annotations [31], which included BP, CC, and MF, with an E-value cut-off = 1e-5. Web Gene Ontology annotation software was adopted to perform GO functional classifications [32]. Furthermore, the sequences were searched against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database using KEGG Automatic Annotation Server (KAAS) with an E-value threshold <1e-10.

Differentially expressed gene (DEG) identification and analysis

The expression level of each unigene was calculated by reads per kilobase millon mapped reads (RPKM) to assess the length and depth of sequencing [33]. Then, the differences in the expression abundance of each gene between each pair of compared samples were calculated by DESeq 2 software (version 1.4.5) [34]. Each resulting p-value was adjusted to a q-value, following the Benjamini-Hochberg procedure for controlling the false discovery rate [35]. The DEGs were identified with q ≤ 0.05 and |log2(fold-change)| ≥ 1 as thresholds. Then, the GO and KEGG analyses were considered to be significantly enriched with q ≤ 0.05 [36]. A GO functional enrichment analysis was performed using the BiNGO plugin of Cytoscape [37].

Reverse transcription quantitative PCR (RT-qPCR) analysis

Two micrograms of total RNA from L. japonica petals was reverse transcribed using the M-MuLV Reverse Transcriptase Kit (Takara, Japan) and oligo (dT) primers according to the manufacturer’s instructions. The diluted cDNA reaction mixture was used as a template in a 20 μl PCR for transcript measurements. qPCR was carried out in a Bio-Rad CFX96™ Real-Time System. The qPCR programme was initiated with a preliminary step of 94 °C for 5 min, followed by 45 cycles of 94 °C for 10 s, 56 °C for 10 s, and 72 °C for 10 s. A melting curve was generated for each sample at the end of each run to ensure the purity of the amplified products with a temperature change of 0.5 °C/s from 65 °C to 95 °C and 5 s for each step afterwards for the melt curve. For each sample, three biological replicates were used. The gene-specific primers of qPCR were designed according to the selected sequences derived from RNA-seq (Table S1). The expression level of Actin was applied to normalize the mRNA levels for each sample [21]. All data were analysed by ANOVA, and the mean differences were compared utilizing Duncan’s multiple range test.


Morphology analysis of petal colour transitions

The petal colour of every single flower was transformed continuously from green to white and then to yellow during flower development in L. japonica. Early in the development of floral buds, primary buds with green petals grew to approximately 3.5 cm in length (Fig. 1a). At the early stage of anthesis, the petals turned from green to white (Fig. 1b). Then, the petals gradually transformed to yellow from white before the withering stage (Fig. 1c). During petal colour transitions, petals at the green bud, white flower and yellow flower stages were selected. The changes in the colour index of GB_Pe, WF_Pe and YF_Pe were significantly different (Table S2). The values of redness (a*) in GB_Pe, WF_Pe and YF_Pe were − 12.36, − 0.58 and 1.25, respectively. The parameter lightness (L*) in WF_Pe was 82.40, which was higher than that in GB_Pe and YF_Pe. The index of yellowness (b*) in YF_Pe was the highest (42.55).

Fig. 1

Morphological observation of L. japonica flowers. a Flower buds with green petals (GB_Pe). b Flowers with white petals (WF_Pe). c Flowers with yellow petals (YF_Pe)

Carotenoid and anthocyanin accumulation in L. japonica petals at various stages

To obtain an accurate understanding of carotenoid accumulation, carotenoid profiling was analysed in L. japonica petals using LC-MS/MS during petal colour transitions. A total of 13 carotenoids were detected from GB_Pe, WF_Pe and YF_Pe (Table 1). The major carotenoids of GB_Pe were lutein, violaxanthin, neoxanthin and zeaxanthin. The lutein content significantly decreased from 39.30 μg/g in GB_Pe to 1.46 μg/g in WF_Pe and slightly decreased to 1.09 μg/g in YF_Pe. The violaxanthin content first significantly decreased from 18.87 μg/g in GB_Pe to 0.60 μg/g in WF_Pe and then drastically increased to 43.81 μg/g in YF_Pe. The trends of neoxanthin and zeaxanthin were similar to those of violaxanthin. Compared with GB_Pe and WF_Pe, most of carotenoids, including α-carotene, antheraxanthin, lycopene, zeaxanthin, violaxanthin, γ-carotene, neoxanthin, β-carotene, β-cryptoxanthin and apocarotenal, were significantly upregulated in YF_Pe. Among them, violaxanthin (43.81 μg/g), zeaxanthin (27.45 μg/g), α-carotene (20.84 μg/g) and γ-carotene (19.97 μg/g) were the major carotenoid compounds in YF_Pe. Except for zeaxanthin, the contents of the remaining 12 carotenoids were all detected at lower levels in WF_Pe.

Table 1 Contents (μg/g DW) of carotenoids in GB_Pe, WF_Pe and YF_Pe

To better understand the content changes of anthocyanins, quantitative analysis of anthocyanins was further performed by LC-MS/MS technology. A total of 10 anthocyanins were identified from GB_Pe, WF_Pe and YF_Pe (Table 2). In GB_Pe, delphinidin was the predominant component of anthocyanins. Specifically, delphinidin was reduced by 3.55- and 2.85-fold in WF_Pe and YF_Pe, respectively, compared with GB_Pe (Table S3). In contrast, pelargonidin was not detected in GB_Pe but at higher levels in WF_Pe and YF_Pe. In the YF_Pe vs GB_Pe comparison, pelargonidin and cyanidin O-syringic acid significantly increased, while delphinidin, cyanidin O-malonyl-malonylhexoside and delphin chloride significantly decreased. Compared with WF_Pe, the contents of pelargonidin and cyanidin were increased by 2.11- and 2.36-fold in YF_Pe, respectively. However, cyanidin O-malonyl-malonylhexoside and delphin chloride were not detected in YF_Pe.

Table 2 Anthocyanins in GB_Pe, WF_Pe and YF_Pe

Effects of endogenous hormones during petal color transitions

To obtain the changes in endogenous hormones, the concentrations of IAA, ZR, GA3, BR, MeJA and ABA were analysed. During petal colour transitions, the concentrations of IAA, ZR, GA3, BR and MeJA decreased, but the content of ABA increased (Fig. 2). The IAA concentration decreased significantly from 717.3 ng·g− 1 FW (GB_Pe) to 191.0 ng·g− 1 FW (WF_Pe) and then to 118.8 ng·g− 1 FW (YF_Pe). The ZR and GA3 concentrations both first decreased significantly from GB_Pe to WF_Pe and then remained stable from WF_Pe to YF_Pe. The BR concentration was highest in GB_Pe. From GB_Pe to YF_Pe, the BR concentration decreased significantly from 9.2 ng·g− 1 FW in GB_Pe to 7.3 ng·g− 1 FW in WF_Pe and then increased slightly to 8.3 ng·g− 1 FW in YF_Pe. The level of MeJA first decreased significantly from GB_Pe to WF_Pe, reaching the lowest level, and then slightly increased from WF_Pe to YF_Pe. However, the ABA concentration increased significantly from 98.0 ng·g− 1 FW to 205.2 ng·g− 1 FW from GB_Pe to YF_Pe (Fig. 2).

Fig. 2

Concentrations of endogenous hormones in petal color transitions in L. japonica. a IAA concentration. b ZR concentration. c GA3 concentration. d BR concentration. e MeJA concentration. f ABA concentration. Significant differences are indicated by different letters at P < 0.05

Sequencing, de novo assembly and annotation

To identify key candidate genes for petal colour transitions, RNA sequencing was carried out from GB_Pe, WF_Pe and YF_Pe. Nine cDNA libraries were sequenced and 448,565,884 raw reads were generated. After data filtering, 408,576,816 (91.1%) clean reads were produced, and the Q30 values were greater than 96.7%. For each sample, clean reads were obtained from 6.6 to 7.1 Gb (Table S4). A total of 69,946 unigenes were generated with an average length of 871 bp and an N50 of 1636 bp (Table S5). Most unigenes (96.6%) were generated from 200 to 3200 bp in length, and 2383 (3.4%) unigenes were more than 3200 bp in length (Fig. S1).

A total of 34,068 assembled unigenes were annotated (Table S6). Based on sequence similarity, 22,662 (32.4%) unigenes were enriched into three groups (biological process, cellular component and molecular function) based on GO term analysis (Fig. S2). The biological processes were mainly focused on ‘cellular process’ and ‘metabolic process’. The cellular components were mainly involved in ‘cell part’. The molecular functions were mainly classified into ‘binding’ and ‘catalytic activity’. KEGG term analysis was used to identify the functional classifications of the unigenes. A total of 9309 (13.31%) unigenes were enriched in 32 KEGG pathway groups, of which ‘signal transduction’ represented the largest group, followed by ‘carbohydrate metabolism’, ‘translation’ and ‘folding, sorting and degradation’ (Fig. S3).

Identification and analysis of DEGs

To detect alterations in gene expression, transcriptomic analyses of WF_Pe vs GB_Pe, YF_Pe vs WF_Pe and YF_Pe vs GB_Pe were carried out to identify the key DEGs during petal colour transition in L. japonica (Fig. S4). A total of 29,679 DEGs were identified based on a 2-fold change at P < 0.05 (Fig. S4a). For each comparison, the numbers of total DEGs, upregulated DEGs and downregulated DEGs were counted, as shown in Fig. S4b.

All 29,679 identified DEGs were further classified into 8 clusters on the basis of expression alterations during petal colour transition (Fig. 3a). A total of 3470 DEGs were classified into two profiles based on expression changes across the three developmental stages: expression stable and then increased (profile 4) and expression stable and then decreased (profile 3). The opposite change patterns of gene expression during the petal colour transition from white to yellow suggest a tight linkage of these genes with petal colour transition in L. japonica.

Fig. 3

Analysis of DEGs. a Cluster analysis of all DEGs from the expression profiles. Black lines indicated the average expression level of unigenes grouped into the same profile. b Enrichment of selected GO terms for DEGs (> 1 RPKM) selected from profiles 3 and 4. The biological process with false discovery rate (adjusted P < 0.05, Student’s t-test) is shown

GO enrichment analysis was further performed to investigate the biological functions of these 1897 DEGs (RPKM > 1 in at least one sample from the 3470 DEGs) that showed higher or lower expression in YF_Pe. The hormone-mediated signalling pathway was significantly enriched in the biological process subcategory (Fig. 3b). DEGs involved in hormone-mediated signalling pathways, such as small auxin-up RNA (SAUR) and PYRABACTIN RESISTANCE1-like (PYL), were significantly differentially expressed between yellow petals and non-yellow petals and seemed relevant to the goal of this study.

Analysis of DEGs involved in hormone-mediated signalling pathways

GO enrichment analysis showed that DEGs were mainly enriched in hormone-mediated signalling pathways. To better investigate hormonal regulation in the colour transitions, we analysed the 67 DEGs (> 1 RPKM) that were enriched in the signalling pathways of auxin, cytokinin, gibberellin, BR, jasmonic acid, ABA and ethylene in YF_Pe vs GB_Pe and YF_Pe vs WF_Pe (Fig. S5 and Table S7).

In the auxin signalling pathway, 15 DEGs were identified, of which the AUX1, TIR1, ARF and SAUR genes were significantly downregulated from GB_Pe to YF_Pe, while three IAAs were upregulated at WF_Pe (Fig. S5a). A total of 18 DEGs were enriched in the cytokinine signalling pathway, including HKs, HPs, type-B RRs and type-A RRs. All of these DEGs were downregulated from GB_Pe to WF_Pe and YF_Pe (Fig. S5b). Meanwhile, in the gibberellin signalling pathway, GID1, GID2 and DELLA genes were identified and significantly downregulated in YF_Pe (Fig. S5c). In the BR signalling pathway, 13 DEGs were identified, most of which were first downregulated and then upregulated in the transition. Specifically, the expression of BRI1, BSK, BZR1_2, CYCD3 and TCH4 was significantly higher in YF_Pe than in WF_Pe (Fig. S5d). Four DEGs were enriched in the jasmonic acid signalling pathway, and their expression levels were higher in GB_Pe than in WF_Pe and YF_Pe (Fig. S5e). Furthermore, JAR1, COI-1 and MYC2 were expressed at higher levels in YF_Pe than in WF_Pe, while JAZ was expressed at lower levels in YF_Pe than in WF_Pe. However, seven DEGs were identified in the ABA signalling pathway, including PYL, PP2C, SNRK2 and ABF, of which PYLs and SNRK2 were significantly upregulated in YF_Pe (Fig. S5f). In the ethylene signalling pathway, five DEGs were identified, of which EIN3, ERS and ERF1 were significantly upregulated in YF_Pe (Fig. S5g).

Analysis of pigment-related DEGs during petal colour transitions

To investigate the pathways of pigment synthesis/degradation during the transitions, the expression levels of carotenoid, anthocyanin and chlorophyll metabolism-related genes were analysed. A total of 49 DEGs (> 1 RPKM) regulating carotenoid, anthocyanin and chlorophyll metabolism were identified and significantly differentially expressed between yellow petals and non-yellow petals (GB_Pe or WF_Pe) (Fig. 4 and Table S7).

Fig. 4

Expression levels of the pigments synthesis/catabolism in GB_Pe, WF_Pe and YF_Pe. a DEGs of carotenoid metabolism-related genes. b DEGs of porphyrin and chlorophyll metabolism-related genes. c DEGs of flavonoid/anthocyanin metabolism-related genes. High expression levels are represented in orchid. Low expression levels are represented in blue

In the carotenoid biosynthesis pathway, PSY1, PDS1, ZDS1 and ZDS2 were significantly upregulated in YF_Pe. However, three carotenoid degradation-related genes, carotenoid cleavage dioxygenase 4 (CCD4), CCD7 and abscisic-aldehyde oxidase 3 (AAO3), were significantly downregulated in YF_Pe (Fig. 4a). Meanwhile, the expression levels of chlorophyll metabolism-related genes showed significant differences. Among these genes, biosynthesis-related genes, including glutamyl-tRNA synthetase (GltX), protoporphyrinogen IX oxidase (PPO) and chlorophyll synthase (CHLG), were significantly upregulated in GB_Pe. However, pheophytinase (PPH), pheophorbide a oxygenase (PAO) and red chlorophyll catabolite reductase (RCCR) were significantly downregulated in GB_Pe (Fig. 4b).

In the basic upstream pathway of flavonoid/anthocyanin biosynthesis, some DEGs were identified in WF_Pe vs GB_Pe, YF_Pe vs WF_Pe and/or YF_Pe vs GB_Pe comparisons, such as phenylalanine ammonia-lyase (PAL), trans-cinnamate 4-monooxygenase (C4H), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3β-hydroxylase (F3H), flavonoid 3′-monooxygenase (F3’H), dihydroflavonol 4-reductase (DFR) and ANS. In the specific anthocyanin downstream branch, one anthocyanidin 3-O-glucosyltransferase (BZ1) gene and one anthocyanidin 3-O-glucoside 5-O-glucosyltransferase (UGT75C1) gene were identified as DEGs. The expression levels of ANS and BZ1 were significantly upregulated in YF_Pe. In contrast, some DEGs, including CHS2, DFRs and UGT75C1, were upregulated in GB_Pe. Although the expression levels of four DFRs and one UGT75C1 first declined in white petals, they then slightly rose from WF_Pe to YF_Pe (Fig. 4c).

Validation of the expression analysis of key genes

A total of fourteen unigenes in pigment metabolism pathways or hormone-mediated signalling pathways were randomly selected and identified by RT-qPCR. The expression patterns of these DEGs corresponded well with the RPKM values obtained by RNAseq (Fig. 5). Pearson correlation analysis showed high correlation coefficients between the RNA-seq and RT-qPCR data, suggesting that the sequencing data are reliable.

Fig. 5

Validation of the expression of pigment-related genes in L. japonica by RT-qPCR. Error bars indicate the standard deviation of three independent biological repeats. Significant differences are indicated by different letters at P < 0.05


Pigment accumulation and metabolism-related genes

During petal colour transitions, a previous study reported that the concentrations of total chlorophyll significantly decreased from green buds to white flowers and then remained stable to yellow flowers (22). In this study, the expression levels of chlorophyll biosynthesis-related genes, including GltX, PPO and CHLG, were significantly downregulated in WF_Pe and YF_Pe compared with GB_Pe (Fig. 4b). As known to us, CHLG is involved in the final step in chlorophyll synthesis [38], the level of CHLG expression showed a significant correlation with chlorophyll content in the petals of L. japonica. In addition to chlorophyll biosynthesis-related genes described above, the expression levels of chlorophyllide a oxygenase (CAO) and hydroxymethyl chlorophyll a reductase (HCAR), which regulate the interconversion of chlorophyll a and chlorophyll b and play a crucial role in greening processes [39,40,41], also drastically decreased in WF_Pe and YF_Pe. In a previous study, the expression level of HCAR showed a positive correlation with chlorophyll content in pale-green petals of carnation [42]. In contrast, chlorophyll degradation-related genes, such as PPH, PAO and RCCR, were significantly upregulated in YF_Pe. These results revealed that the decrease in chlorophyll accumulation mainly accounted for the loss of visual green colour in WF_Pe and YF_Pe.

The total contents of chlorophyll and carotenoid were both found to decrease to their lowest levels in WF_Pe (22). Accordingly, the expression of chlorophyll and carotenoid biosynthesis genes was significantly downregulated in WF_Pe. Interestingly, some anthocyanins were detected in WF_Pe, such as delphinidin and petunidin 3-O-glucoside, suggesting the presence of unobservable pigments in white petals.

In flowering plants, carotenoids mainly participate in petal colours ranging from yellow to red [3]. Several plant lineages have yellow flowers and contain pigments derived from carotenoids [5]. A previous study reported that the content of total carotenoids dramatically increased from the WF to YF stages [22, 24]. Similarly, most of the detected carotenoids were significantly upregulated in YF_Pe compared with WF_Pe in this study (Table 1). Among the 10 significantly higher carotenoids recorded in YF_Pe, the top four carotenoids were α-carotene, zeaxanthin, violaxanthin and γ-carotene. These four carotenoids may be the major contributors to yellow color of YF_Pe.

In accordance with the changes in carotenoid content, carotenoid biosynthesis-related genes, such as PSY1, PDS1, ZDS1, ZDS2, LCYB and VDE, were significantly upregulated in YF_Pe compared with WF_Pe (Fig. 4a), and similar results were reported in a previous study [24]. Because of their key roles in regulating carotenoid biosynthesis, PSY, PDS and ZDS have been subjected to intensive investigation. Previously, in the steps of carotenoid biosynthesis, PSY was shown to be involved in the condensation of two geranylgeranyl diphosphate molecules into phytoene, and the upregulation of PSY enhances carotenoid accumulation [43,44,45,46,47,48]. Then, the phytoene is subjected to a series of desaturation reactions catalysed by carotene desaturases, such as PDS and ZDS [4, 44]. Similar observations of petal colours have been reported in monocots, such as Lilium and Oncidium [7, 49]. In different cultivars of Asiatic hybrid lily, the petal colours are correlated well with the transcription levels of biosynthetic genes, including PSY, PDS and ZDS [49].

In contrast, the expression levels of the carotenoid catabolism-related genes CCD4 and CCD7 were significantly downregulated in YF_Pe. It has been found that the expression of CCD4 affects carotenoid levels in various plants. As reported for pigments of Brassica and Dendranthema, the increase in carotenoid content is related to disruption of a CCD4 gene involving the petal colour transition from white to yellow [50, 51]. In rose cultivars of yellow petals, carotenoid degradation has a high correlation with the expression of RhCCD4 [52].

Furthermore, we found that two anthocyanins (pelargonidin and cyanidin) were significantly upregulated in YF_Pe. In particular, the pelargonidin and cyanidin concentrations both tripled from WF_Pe to YF_Pe. In Gentiana lutea L. var. aurantiaca, the orange petal colour is predominantly caused by newly synthesized pelargonidin glycosides, which confer a reddish hue to the yellow background colour derived from carotenoids [53]. Pelargonidin and cyanidin accumulated at increasing concentrations from WF_Pe to YF_Pe, in sync with the variation in carotenoid content, which might assist the petal colour transition from a white to golden appearance.

Herein, the expression levels of some key anthocyanin biosynthesis-related genes, such as DFRs, ANS, BZ1 and UGT75C1, were upregulated in YF_Pe compared with WF_Pe (Fig. 4c). As an important step in the anthocyanin biosynthetic pathway, DFR is associated with anthocyanin colouration in flowers, such as Freesia hybrida [54], Andean genus Iochroma (Solanaceae) [55] and Cape Erica species [56]. In the present study, four DFR genes were identified and increased from WF_Pe to YF_Pe. To our knowledge, ANS is considered the first key enzyme that could lead to flavonoid flux into the anthocyanin branch. The expression level of ANS was upregulated in YF_Pe, which accounted for the increased concentrations of pelargonidin and cyanidin.

Changes in hormones concentrations associated with pigment accumulation

Plant hormones are involved in all stages of flower development. In our study, changes in the concentrations of IAA, ZR, GA3, BR, MeJA and ABA revealed that the transitions are regulated by endogenous hormones. Accordingly, the key DEGs involved in the hormone signal transduction pathways were significantly enriched in the auxin, cytokinin, gibberellin, BR, jasmonic acid, ABA and ethylene signalling pathways.

In this study, the concentration of IAA was significantly decreased from GB_Pe to YF_Pe (Fig. 2a). Accordingly, the genes from the auxin signalling pathway were mostly downregulated in YF_Pe (Fig. S5a). For example, SAUR genes, a family of auxin-responsive genes in auxin signalling pathways, were downregulated in the transitions. However, compared with WF_Pe, the expression levels of three IAAs were significantly downregulated in YF_Pe. It was previously reported that IAA represses the transcription of carotenoid biosynthesis-related genes, e.g., PSY, ZISO, PDS and CRTISO [57]. This indicated that downregulation of the IAA genes play important roles in carotenoid accumulation.

The concentrations of BR and MeJA were both slightly higher in YF_Pe than in WF_Pe (Fig. 2d, e). The expression levels of most genes from the BR and jasmonic acid signalling pathways, including BRI1, BSK, BZR1_2, TCH4, CYCD3, JAR1, COI-1 and MYC2, were upregulated in YF_Pe (Fig. S5d, e). In Solanum lycopersicum, the application of BR was previously reported to increase carotenoid accumulation [58]. Furthermore, ectopic expression of BZR1-1D results in an increase in carotenoid accumulation [59]. Meanwhile, the exogenous application of MeJA in the Never ripe mutant of Solanum lycopersicum significantly enhanced lycopene accumulation, as well as the expression levels of PSY1 and PDS [60]. In the present study, the expression levels of PSY1 and PDS1 were significantly higher in YF_Pe than in WF_Pe (Fig. 4a).

Genes from the ethylene and ABA signalling pathways were significantly upregulated during petal colour transitions (Fig. S5f, g). Previous studies indicate that upregulation of ethylene-related genes accelerates chlorophyll degradation [61] and carotenoid accumulation [57]. In our study, the upregulation of EIN3, ERS and ERF1 might lead to chlorophyll degradation and carotenoid accumulation in YF_Pe. In the ABA signalling pathway, the expression level of ABF reached the highest level in WF_Pe; PYL and SNRK2 were upregulated in WF_Pe and YF_Pe. In Arabidopsis, overexpression of ABF is associated with triggering chlorophyll degradation [62,63,64]. ABA signalling pathway related genes upregulated in WF_Pe and YF_Pe along with the increased concentration of ABA during flower development play roles in petal color transitions.


In this study, the most comprehensive metabolome, hormone, and transcriptome analyses investigated petal colour transitions in L. japonica. Analyses of key candidate genes, metabolites and hormones highlighted the effects of carotenoids, anthocyanins and endogenous hormones; this enabled us to clarify the regulatory mechanisms underlying the transitions. Based on our results and previously published studies, we provide a conceptual model for the regulatory network of the transitions in L. japonica (Fig. 6). In this model, the existing chlorophyll/carotenoid balance is disturbed during flower development. Overall, the content of total chlorophylls decreased in WF_Pe and YF_Pe, along with low expression of chlorophyll synthesis genes and high expression of chlorophyll degradation genes, resulting in green colour loss with flower development. In contrast, although the content of total carotenoids first dropped in WF_Pe, it then rapidly accumulated from WF_Pe to YF_Pe, in sync with the increase in the expression of genes related to carotenoid biosynthesis. The ten highest concentrations carotenoids were detected in YF_Pe, and the top four carotenoids were α-carotene, zeaxanthin, violaxanthin and γ-carotene, which are major contributors to yellow color of YF_Pe. It was also found that a few anthocyanins differentially accumulated during flower development, indicating that they may assist with to vivid colour of GB_Pe and YF_Pe. Meanwhile, this developmental process is regulated by endogenous hormones. These variations in key pigment and hormone-mediated signalling pathway-related genes, pigments and hormones promote petal colour transitions from green to white and then to yellow in L. japonica.

Fig. 6

Schematic of changes in the regulatory genes and metabolites in petal color- transition in L. japonica

Availability of data and materials

These sequence data have been submitted to the SRA database under accession number PRJNA574570. The datasets supporting the conclusions of this article are included within the article and its additional files.



Green bud petal


White flower petal


Yellow flower petal


Differentially expressed genes

GA3 :

Gibberellic acid




Abscisic acid


Indoleacetic acid


Zeatin riboside


Methyl jasmonate


Mass spectrometry


Quality control




Analysis of variance


Gene ontology


Kyoto Encyclopedia of Genes and Genomes


Reads per kilobase millon mapped reads


Reverse transcription quantitative PCR


Anthocyanidin synthase


Phytoene synthase


Phytoene desaturase


ζ-Carotene desaturase


Small auxin-up RNA




Carotenoid cleavage dioxygenase


Abscisic-aldehyde oxidase 3


Glutamyl-tRNA synthetase


Protoporphyrinogen IX oxidase


Chlorophyll synthase




Pheophorbide a oxygenase


Red chlorophyll catabolite reductase


Chlorophyllide a oxygenase


Hydroxymethyl chlorophyll a reductase


Phenylalanine ammonia-lyase


Trans-cinnamate 4-monooxygenase


4-Coumarate-CoA ligase


Chalcone synthase


Chalcone isomerase


Flavanone 3β-hydroxylase


Flavonoid 3′-monooxygenase


Dihydroflavonol 4-reductase


Anthocyanidin 3-O-glucosyltransferase


Anthocyanidin 3-O-glucoside 5-O-glucosyltransferase


  1. 1.

    Strauss SY, Whittall JB. Non-pollinator agents of selection on floral traits. In: Harder LD, Barrett SCH, editors. Ecology and Evolution of Flowers. Oxford: Oxford University Press; 2006. p. 120–38.

    Google Scholar 

  2. 2.

    Meng Y, Wang Z, Wang Y, Wang C, Zhu B, Liu H, et al. The MYB activator WHITE PETAL1 associates with MtTT8 and MtWD40-1 to regulate carotenoid-derived flower pigmentation in Medicago truncatula. Plant Cell. 2019;31(11):2751–67.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Nisar N, Li L, Lu S, Khin NC, Pogson BJ. Carotenoid metabolism in plants. Mol Plant. 2015;8(1):68–82.

    CAS  Article  Google Scholar 

  4. 4.

    Sun T, Yuan H, Cao H, Yazdani M, Tadmor Y, Li L. Carotenoid metabolism in plants: the role of plastids. Mol Plant. 2018;11(1):58–74.

    PubMed  Article  CAS  Google Scholar 

  5. 5.

    Grotewold E. The genetics and biochemistry of floral pigments. Annu Rev Plant Biol. 2006;57(1):761–80.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Zhao D, Tao J. Recent advances on the development and regulation of flower color in ornamental plants. Front Plant Sci. 2015;6:261.

    PubMed  PubMed Central  Google Scholar 

  7. 7.

    Chiou C-Y, Pan H-A, Chuang Y-N, Yeh K-W. Differential expression of carotenoid-related genes determines diversified carotenoid coloration in floral tissues of Oncidium cultivars. Planta. 2010;232(4):937–48.

    CAS  Article  Google Scholar 

  8. 8.

    Wang Z, Shen Y, Yang X, Pan Q, Ma G, Bao M, et al. Overexpression of particular MADS-box transcription factors in heat-stressed plants induces chloroplast biogenesis in petals. Plant Cell Environ. 2019;42(5):1545–60.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54(4):733–49.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Ferreyra MLF, Rius SP, Casati P. Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front Plant Sci. 2012;3:222.

    Google Scholar 

  11. 11.

    Hirschberg J. Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol. 2001;4(3):210–8.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Han Y, Wang X, Chen W, Dong M, Yuan W, Liu X, et al. Differential expression of carotenoid-related genes determines diversified carotenoid coloration in flower petal of Osmanthus fragrans. Tree Genet Genomes. 2014;10(2):329–38.

    Article  Google Scholar 

  13. 13.

    Moehs CP, Tian L, Osteryoung KW, Dellapenna D. Analysis of carotenoid biosynthetic gene expression during marigold petal development. Plant Mol Biol. 2001;45(3):281–93.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Dettmer K, Aronov PA, Hammock BD. Mass spectrometry-based metabolomics. Mass Spectrom Rev. 2007;26(1):51–78.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Bino RJ, Hall RD, Fiehn O, Kopka J, Saito K, Draper J, et al. Potential of metabolomics as a functional genomics tool. Trends Plant Sci. 2004;9(9):418–25.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Deng C, Li S, Feng C, Hong Y, Huang H, Wang J, et al. Metabolite and gene expression analysis reveal the molecular mechanism for petal colour variation in six Centaurea cyanus cultivars. Plant Physiol Biochem. 2019;142:22–33.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Lou Q, Liu Y, Qi Y, Jiao S, Tian F, Jiang L, et al. Transcriptome sequencing and metabolite analysis reveals the role of delphinidin metabolism in flower colour in grape hyacinth. J Exp Bot. 2014;65(12):3157–64.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Duan H, Wang L, Cui G, Zhou X, Duan X, Yang H. Identification of the regulatory networks and hub genes controlling alfalfa floral pigmentation variation using RNA-sequencing analysis. BMC Plant Biol. 2020;20(1):110.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Yan K, Cui M, Zhao S, Chen X, Tang X. Salinity stress is beneficial to the accumulation of chlorogenic acids in honeysuckle (Lonicera japonica Thunb.). Front Plant Sci. 2016;7:1563.

  20. 20.

    He SQ, Hu QF, Yang GY. Research of honeysuckle. Yunnan Chem Technol. 2010;37(3):72–5 (In Chinese).

    CAS  Google Scholar 

  21. 21.

    Wu J, Wang X, Liu Y, Du H, Shu Q, Su S, et al. Flavone synthases from Lonicera japonica and L. macranthoides reveal differential flavone accumulation. Sci Rep. 2016;6:19245.

  22. 22.

    Fu L, Li H, Li L, Yu H, Wang L. Reason of flower color change in Lonicera japonica. Scientia Silvae Sinicae. 2013;49(10):155–61 (In Chinese).

    CAS  Google Scholar 

  23. 23.

    Li J, Lian X, Ye C, Wang L. Analysis of flower color variations at different developmental stages in two honeysuckle (Lonicera Japonica Thunb.) cultivars. HortScience. 2019;54(5):779–82.

    CAS  Article  Google Scholar 

  24. 24.

    Pu X, Li Z, Tian Y, Gao R, Hao L, Hu Y, et al. The honeysuckle genome provides insight into the molecular mechanism of carotenoid metabolism underlying dynamic flower coloration. New Phytol. 2020;227(3):930–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Liu Y, Lv J, Liu Z, Wang J, Yang B, Chen W, et al. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in pepper fruit (Capsicum annuum L.). Food Chem. 2020;306:125629.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Wang D, Zhang L, Huang X, Wang X, Yang R, Mao J, et al. Identification of nutritional components in black sesame determined by widely targeted metabolomics and traditional Chinese medicines. Molecules. 2018;23(5):1180.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  27. 27.

    Huang D, Yuan Y, Tang Z, Huang Y, Kang C, Deng X, et al. Retrotransposon promoter of Ruby1 controls both light-and cold-induced accumulation of anthocyanins in blood orange. Plant Cell Environ. 2019;42(11):3092–104.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Yang Y-M, Xu C-N, Wang B-M, Jia J-Z. Effects of plant growth regulators on secondary wall thickening of cotton fibres. Plant Growth Regul. 2001;35(3):233–7.

    CAS  Article  Google Scholar 

  29. 29.

    Zeng Y-H, Zahng Y-P, Xiang J, Wu H, Chen H-Z, Zhang Y-K, et al. Effects of chilling tolerance induced by spermidine pretreatment on antioxidative activity,endogenous hormones and ultrastructure of indica-japonica hybrid rice seedlings. J Integr Agric. 2016;15(2):295–308.

    CAS  Article  Google Scholar 

  30. 30.

    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(7):644–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    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.

    CAS  Article  Google Scholar 

  32. 32.

    Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34(suppl_2):W293–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  35. 35.

    Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci. 2003;100(16):9440–5.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  37. 37.

    Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9.

    CAS  Article  Google Scholar 

  38. 38.

    Shalygo N, Czarnecki O, Peter E, Grimm B. Expression of chlorophyll synthase is also involved in feedback-control of chlorophyll biosynthesis. Plant Mol Biol. 2009;71(4):425–36.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Tanaka A, Ito H, Tanaka R, Tanaka NK, Yoshida K, Okada K. Chlorophyll a oxygenase (CAO) is involved in chlorophyll b formation from chlorophyll a. Proc Natl Acad Sci. 1998;95(21):12719–23.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Kunugi M, Takabayashi A, Tanaka A. Evolutionary changes in chlorophyllide a oxygenase (CAO) structure contribute to the acquisition of a new light-harvesting complex in micromonas. J Biol Chem. 2013;288(27):19330–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Meguro M, Ito H, Takabayashi A, Tanaka R, Tanaka A. Identification of the 7-hydroxymethyl chlorophyll a reductase of the chlorophyll cycle in Arabidopsis. Plant Cell. 2011;23(9):3442–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Ohmiya A, Hirashima M, Yagi M, Tanase K, Yamamizo C. Identification of genes associated with chlorophyll accumulation in flower petals. PLoS One. 2014;9(12):e113738.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. 43.

    Fraser PD, Truesdale MR, Bird CR, Schuch W, Bramley PM. Carotenoid biosynthesis during tomato fruit development (evidence for tissue-specific gene expression). Plant Physiol. 1994;105(1):405–13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Moise AR, Al-Babili S, Wurtzel ET. Mechanistic aspects of carotenoid biosynthesis. Chem Rev. 2014;114(1):164–93.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Shewmaker CK, Sheehy JA, Daley M, Colburn S, Ke DY. Seed-specific overexpression of phytoene synthase: increase in carotenoids and other metabolic effects. Plant J. 1999;20(4):401–12.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Ducreux LJM, Morris WL, Hedley PE, Shepherd T, Davies HV, Millam S, et al. Metabolic engineering of high carotenoid potato tubers containing enhanced levels of β-carotene and lutein. J Exp Bot. 2004;56(409):81–9.

    PubMed  PubMed Central  Google Scholar 

  47. 47.

    Ruiz-Sola MA, Rodriguez-Concepcion M. Carotenoid biosynthesis in Arabidopsis: a colorful pathway. Arabidopsis Book. 2012;10(10):e0158.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Yuan H, Zhang J, Nageswaran D, Li L. Carotenoid metabolism and regulation in horticultural crops. Horticult Res. 2015;2(1):15036.

    Article  CAS  Google Scholar 

  49. 49.

    Yamagishi M, Kishimoto S, Nakayama M. Carotenoid composition and changes in expression of carotenoid biosynthetic genes in tepals of Asiatic hybrid lily. Plant Breed. 2010;129(1):100–7.

    CAS  Article  Google Scholar 

  50. 50.

    Ohmiya A, Kishimoto S, Aida R, Yoshioka S, Sumitomo K. Carotenoid cleavage Dioxygenase (CmCCD4a) contributes to white color formation in chrysanthemum petals. Plant Physiol. 2006;142(3):1193–201.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Zhang B, Liu C, Wang Y, Yao X, Liu K. Disruption of a CAROTENOID CLEAVAGE DIOXYGENASE 4 gene converts flower colour from white to yellow in Brassica species. New Phytol. 2015;206(4):1513–26.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Glick A. Synthesis and degradation of carotenoids in cut rose petals during vase life, and characterization of the effect of Methyljasmonate treatment on the processes. Jerusalem: Hebrew University of Jerusalem; 2009.

    Google Scholar 

  53. 53.

    Diretto G, Jin X, Capell T, Zhu C, Gomez-Gomez L, Xu C. Differential accumulation of pelargonidin glycosides in petals at three different developmental stages of the orange-flowered gentian (Gentiana lutea L. var. aurantiaca). PLoS One. 2019;14(2):e0212062.

  54. 54.

    Yueqing L, Xingxue L, Xinquan C, Xiaotong S, Ruifang G, Song Y, et al. Dihydroflavonol 4-Reductase genes from Freesia hybrida play important and partially overlapping roles in the biosynthesis of flavonoids. Front Plant Sci. 2017;8:428.

    Google Scholar 

  55. 55.

    Smith SD, Shunqi W, Rausher MD. Functional evolution of an anthocyanin pathway enzyme during a flower color transition. Mol Biol Evol. 30(3):602–12.

  56. 56.

    Maitre NCL, Pirie MD, Bellstedt DU. Floral color, anthocyanin synthesis gene expression and control in cape Erica species. Front Plant Sci. 2019;10:1565.

    PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Su L, Diretto G, Purgatto E, Danoun S, Zouine M, Li Z, et al. Carotenoid accumulation during tomato fruit ripening is modulated by the auxin-ethylene balance. BMC Plant Biol. 2015;15(1):114.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  58. 58.

    Vidya Vardhini B, Rao SSR. Acceleration of ripening of tomato pericarp discs by brassinosteroids. Phytochemistry. 2002;61(7):843–7.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Liu L, Jia C, Zhang M, Chen D, Chen S, Guo R, et al. Ectopic expression of a BZR1-1D transcription factor in brassinosteroid signalling enhances carotenoid accumulation and fruit quality attributes in tomato. Plant Biotechnol J. 2014;12(1):105–15.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Liu L, Wei J, Zhang M, Zhang L, Li C, Wang Q. Ethylene independent induction of lycopene biosynthesis in tomato fruits by jasmonates. J Exp Bot. 2012;63(16):5751–61.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Zhang M, Yuan B, Leng P. The role of ABA in triggering ethylene biosynthesis and ripening of tomato fruit. J Exp Bot. 2009;60(6):1579–88.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Gao S, Gao J, Zhu X, Song Y, Li Z, Ren G, et al. ABF2, ABF3, and ABF4 promote ABA-mediated chlorophyll degradation and leaf senescence by transcriptional activation of chlorophyll catabolic genes and senescence-associated genes in Arabidopsis. Mol Plant. 2016;9(9):1272–85.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Yoshida T, Fujita Y, Sayama H, Kidokoro S, Maruyama K, Mizoi J, et al. AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. Plant J. 2010;61(4):672–85.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Fujita Y, Yoshida T, Yamaguchishinozaki K. Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants. Physiol Plant. 2013;147(1):15–27.

    CAS  PubMed  Article  Google Scholar 

Download references


Not Applicable.


This work was partially supported by National Natural Science Foundation of China (31800600), Fundamental Research Funds for the Central Universities (XDJK2019AA001 and XDJK2020B058), Innovation Research Group Funds for Chongqing Universities (CXQT19005). The funding bodies were not involved in the design of the study; collection, analysis, or interpretation of data; or manuscript writing.

Author information




GL and YX conceived the project. GL and QG supervised the project. YX, WC, and WX designed the experiments and analyzed the data. BX measured the color index of petals. DW and BX performed carotenoids extraction and analyses. XL and LX performed flavonoids extraction and analyses. WC and WX performed endogenous hormones extraction and determination. DW and SW performed RNA isolation and RT-qPCR analyses. YX, WC, and QG wrote the manuscript. All authors approved the final draft of the manuscript.

Corresponding authors

Correspondence to Qigao Guo or Guolu Liang.

Ethics declarations

Ethics approval and consent to participate

All the plant materials used in this study were provided by Beijing Botanical Garden, Chinese Academy of Science. The field experiments were conducted under local legislation and permissions.

Consent for publication

All authors agreed to publish.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Figure S1.

Length distribution of L. japonica transcriptome.

Additional file 2: Figure S2.

GO functional analysis and classification of L. japonica transcriptome.

Additional file 3: Figure S3.

KEGG classification of L. japonica transcriptome.

Additional file 4: Figure S4.

The numbers of DEGs among three L. japonica petals. a The relationships among the three DEGs datasets. b The number of DEGs in each comparison.

Additional file 5: Figure S5.

Expression profiles of DEGs that regulate plant hormone signal transduction in GB_Pe, WF_Pe, and YF_Pe. Orchid, high expression levels; blue, low expression levels. Genes encoding key enzymes of auxin (a), cytokinine (b), gibberellin (c), brassinosteroid (d), jasmonic acid (e), abscisic acid (f) and ethylene (g) signaling pathways were exhibited.

Additional file 6: Table S1.

Primers used for reverse transcription quantitative PCR (RT-qPCR). Table S2. Color parameters in GB_Pe, WF_Pe, and YF_Pe. Table S3. Differentially-accumulated anthocyanins in WF_Pe vs GB_Pe, YF_Pe vs GB_Pe and/or YF_Pe vs WF_Pe. Table S4. Throughput and quality of L. japonica transcriptome data. Table S5. Summary of L. japonica de novo transcriptome assembly. Table S6. Mapping results of L. japonica unigenes to various databases.

Additional file 7: Table S7.

Genes involving in plant hormone signal transduction and pigments metabolism.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xia, Y., Chen, W., Xiang, W. et al. Integrated metabolic profiling and transcriptome analysis of pigment accumulation in Lonicera japonica flower petals during colour-transition. BMC Plant Biol 21, 98 (2021).

Download citation


  • Lonicera japonica
  • Petal colour
  • Pigment
  • Gene expression
  • Endogenous hormones