Skip to main content

Study on the causes of changes in colour during Hibiscus syriacus flowering based on transcriptome and metabolome analyses

Abstract

Background

The flower colour of H. syriacus ‘Qiansiban’ transitions from fuchsia to pink–purple and finally to pale purple, thereby enhancing the ornamental value of the cultivars. However, the molecular mechanism underlying this change in flower colour in H. syriacus has not been elucidated. In this study, the transcriptomic data of H. syriacus ‘Qiansiban’ at five developmental stages were analysed to investigate the impact of flavonoid components on flower colour variation. Additionally, five cDNA libraries were constructed from H. syriacus ‘Qiansiban’ during critical blooming stages, and the transcriptomes were sequenced to investigate the molecular mechanisms underlying changes in flower colouration.

Results

High-performance liquid chromatography‒mass spectrometry detected five anthocyanins in H. syriacus ‘Qiansiban’, with malvaccin-3-O-glucoside being the predominant compound in the flowers of H. syriacus at different stages, followed by petunigenin-3-O-glucoside. The levels of these five anthocyanins exhibited gradual declines throughout the flowering process. In terms of the composition and profile of flavonoids and flavonols, a total of seven flavonoids were identified: quercetin-3-glucoside, luteolin-7-O-glucoside, Santianol-7-O-glucoside, kaempferol-O-hexosyl-C-hexarbonoside, apigenin-C-diglucoside, luteolin-3,7-diglucoside, and apigenin-7-O-rutinoside. A total of 2,702 DEGs were identified based on the selected reference genome. Based on the enrichment analysis of differentially expressed genes, we identified 9 structural genes (PAL, CHS, FLS, DRF, ANS, CHI, F3H, F3’5’H, and UFGT) and 7 transcription factors (3 MYB, 4 bHLH) associated with flavonoid biosynthesis. The qRT‒PCR results were in good agreement with the high-throughput sequencing data.

Conclusion

This study will establish a fundamental basis for elucidating the mechanisms underlying alterations in the flower pigmentation of H. syriacus.

Peer Review reports

Introduction

Hibiscus syriacus L., a perennial deciduous shrub or small tree, belongs to the Malvaceae family and the genus Hibiscus. Its cultivars are rich in colour and include red, pink, purple, blue, and white. The trees are fully flowered in full bloom and have high ornamental value. The flowering period of individual H. syriacus is short and lasts only one day, a phenomenon also known as blooming in the morning and falling in the evening. Although the flowers of these plants have a short lifespan, the flower colour changes in some cultivars, such as ‘Qiansiban’, during the flowering process from fuchsia to pink‒purple and then to pale purple. This characteristic confers high ornamental and research value to this cultivar. Many plants change flower colour during flowering. Paeonia ‘Coral Sunset’ and ‘Pink Hawaiian Coral’ flowers change from coral to pink and then to pale yellow during flowering [1]. The corolla of purple‒red Mirabilis jalapa changes from light green to purple‒red between budding and blooming [2].

Flower colour, an important ornamental trait of plants, can be affected by anthocyanin content, cell sap pH, and metal ions. Among these, the anthocyanin components and contents are the most important factors. There are three main types of pigments related to flower colour formation: flavonoids [3], carotenoids [4], and alkaloids [5]. Flavonoids include mainly anthocyanins, flavones, and flavonols [6]. At present, 27 categories of anthocyanins have been identified in nature, of which only 6 are the most common, namely, cyanidin, geranium, delphinidin, paeoniflorin, petunia, and mallow [5]. Geranium is usually orange to red, cyanidin is pink, delphinidin is blue to purple, and flavonoids and flavonols are colourless or light yellow, which generally affects the formation of flower colour as pigments [7]. P Z et al. [8] detected major anthocyanins, including cyanidin, geranium, delphinidin, paeoniflorin, and petunia, in H. syriacus. In addition, the flavonoid biosynthetic pathway was studied in detail. In the flavonoid metabolism pathway, structural genes such as CHS, CHI, and DRF regulate the synthesis of anthocyanins and affect the colouration of H. syriacus [9].

However, the mechanism underlying the colour change in H. syriacus during flowering has not been determined. To investigate the underlying mechanism, H. syriacus ‘Qiansiban’, which changes flower colour during flowering, was selected as the experimental material. In the present study, flavonoids in petals at five developmental stages were identified via high-performance liquid chromatography‒mass spectrometry (HPLC‒MS), and the relationship between flower colouration and pigment content was analysed. Furthermore, a comparative transcriptomic analysis of five critical flowering states in ‘Qiansiban’ was conducted to identify the differentially expressed genes (DEGs) related to flavonoid biosynthesis. This research will help elucidate the mechanism of flower colour change in H. syriacus.

Materials and methods

Plant materials

The petals of ‘Qiansiban’ were used as experimental materials. Samples were collected at five different flowering stages (S1, S2, S3, S4, S4, and S5) in July 2021 at the Hebei Normal University of Science and Technology, and 15 samples were collected at each stage. Materials for the tests were immediately placed in liquid nitrogen and stored at −80 °C.

Anthocyanin quantification

Flower petals were collected at five sequential developmental stages: the splitting stage (S1), the early stage (S2), the initial opening period (S3), the blooming period (S4), and the final flowering stage (S5). A total of 0.1 g of H. syriacus corolla from each of the five periods was ground into powder with liquid nitrogen, poured into a centrifuge tube, and soaked in 5 mL of a 1% hydrochloric acid–methanol solution in a dark refrigerator at 4 °C for 24 h (during several electric shocks), and the supernatant obtained after centrifugation at 4 °C for 15 min (13,000 r·min− 1) was submitted for determination by HPLC‒MS.

The samples were analysed with an HPLC‒MS system. The analysis methods used were described by C Y [10] and H C [11] and were slightly modified. The liquid-phase conditions were as follows: the column temperature was set at 40 °C, and 10.0 µL of sample was injected at 0.50 mL/min. An ultrapure aqueous solution (with 1% formic acid) was used as mobile phase A, and acetonitrile was used as mobile phase B. After 95% A + 5% B elution for 0 min, the volume fraction of B was increased to 40%, and the mixture was eluted for 40 min. Finally, the volume fraction of B was reduced to 5%, and the mixture was eluted for 40 min. The absorption spectra were scanned at 350 and 520 nm.

The MS conditions for the positive ion mode, electrospray ionization source, dryer, and ion source temperature were set at 400 °C and 120 °C, respectively. The capillary pyrolysis and cone hole voltages are 3.5 kV and 40 V, respectively. The scanning points were in the range of 100–1200 m/z.

With standard cyanidin and rutin as the controls, the concentration as the abscissa, and the peak area as the ordinate, a standard curve was generated to calculate the contents of anthocyanins and flavonoids.

The standard equation for cyanidin content was y = 528.57x−9.6762 (R2 = 0.9989); the standard equation for rutin content was y = 534.58x + 499.02 (R2 = 0.996).

RNA extraction, cDNA library construction, and transcriptome sequencing

Total RNA was extracted using an RNAprep Pure Plant Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. Transcriptome sequencing and library construction were completed by BMK Biotechnology Co., Ltd. (China). First, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. First-strand cDNA was synthesized, and second-strand cDNA was subsequently generated. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of the 3’ ends of the DNA fragments, an NEBNext Adaptor with a hairpin loop structure was ligated to prepare the fragments for hybridization. The library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, USA). Then, 3 µl of USER Enzyme (NEB, USA) was incubated with size-selected, adaptor-ligated cDNA at 37°C for 15 min followed by 5 min at 95°C before PCR. PCR was subsequently performed with Phusion High-Fidelity DNA polymerase, universal PCR primers, and Index (X) Primer. Finally, the PCR products were purified (AMPure XP system), and the library quality was assessed on an Agilent Bioanalyzer 2100 system. The transcriptome library of H. syriacus ‘Qiansiban’ was subsequently sequenced on the Illumina NovaSeq 6000 sequencing platform.

De novo assembly and annotation

The raw reads were filtered by removing adaptor and low-quality sequences to obtain high-quality reads. The gene functions were annotated based on the following databases: Nr (NCBI nonredundant protein sequences); Pfam (protein family); KOG/COG (Clusters of Orthologous Groups of proteins); Swiss-Prot (a manually annotated and reviewed protein sequence database); KO (KEGG Orthologue database); and GO (Gene Ontology).

Verification by qRT‒PCR

The expression levels of key genes in the petals of ‘Qiansiban’ were measured using qRT‒qPCR. GAPDH was used as an internal reference gene for real-time fluorescence-based quantitative reactions. The 10 µl reaction mixture was formulated based on the instructions provided for the TB Green TaKaRa Premix Ex Taq™ II (Tli RNaseH Plus; TaKaRa Bio, Inc., Japan). The reaction procedure involved predenaturation at 95 °C for 30 s. The cycle stage was as follows: denaturation at 95 °C for 5 s and annealing at 60 °C for 30 s for 40 cycles. Dissolution curves were recorded from 60 °C to 95 °C with a 0.5 °C increase every 1 min. Each reaction was repeated three times. The relative expression levels of the target genes were calculated by the 2 − ΔΔCt method [12], and the relevant data were analysed with SPSS version 26.

Results

Qualitative and quantitative analyses of anthocyanins

The flower colour of ‘Qiansiban’ is purplish red at the S1 stage, and as the flowers bloom, the colour gradually changes to pinkish purple and purple (Fig. 1A).

Fig. 1
figure 1

Anthocyanin contents in the ‘Qiansiban’ petals at different stages. (A) The colour of the flower samples varied from S1 to S5. (B) Total contents of individual anthocyanidins from S1 to S5. The bars represent the average data of the three HPLC‒MS/MS analyses

Anthocyanins and flavones were identified by comparison with the HPLC retention times, elution orders, UV–vis spectra, and MS fragmentation patterns in published data (Table 1). Through the analysis of molecular ions and fragment ions via mass spectrometry and in combination with previous reports [13,14,15], five anthocyanins were detected in the petals of ‘Qiansiban’ in the five stages, namely, delphinidin-3-O-glucoside (De3G), cyanidin-3-O-glucoside (Cy3G), petunidin-3-O-glucoside (Pe3G), pelargonidin-3-O-glucoside (Pg3G), and malvacin-3-O-glucoside (Ma3G) (Table 1). The remaining 7 compounds were speculated to be quercetin-3-glucoside, luteolin-7-O-glucoside, santianol-7-O-glucoside, kaempferol-O-hexose-C-hexoside, apigenin-C-glucoside, luteolin-3,7-diglucoside, and apigenin-7-O-rutoside (Table 1) [10, 11, 16].

H. syriacus petals undergo major changes in flavonoid and anthocyanin accumulation from S1 to S5. Anthocyanidin is the main contributor to the colour of H. syriacus; thus, its type, content, and biosynthesis process in H. syriacus petals were analysed thoroughly.

The flavonoid accumulation trend was S1 > S2 > S3 > S4 > S5, which gradually decreased with increasing flower opening. The analysis based on HPLC‒MS/MS revealed that all anthocyanidins were highest in content in S1 and lowest in content in S5. The contents of these five anthocyanins gradually decreased over the course of the flowering process, decreasing by 46%, 52%, 40%, 64%, and 39%, from S1 to S5, respectively. The anthocyanidin with the highest content in S1 was Ma3G (Fig. 1B), and the flavone with the highest content was apigenin-C-glucoside (Supplementary Table S1).

In summary, the accumulation of flavonoids and anthocyanidin reached the highest and lowest levels, respectively, at S1 and S5. Ma3G (the anthocyanidin with the highest content), De3G, Pe3G, Cy3G, and Pg3G were the main chromogenic anthocyanidins. With decreasing anthocyanidin concentration, the colour of the petals changed from purplish red to purple (Fig. 1).

Table 1 Chromatographic and spectral data of flavonoids from ‘Qiansiban’

Transcriptome profiling of H. syriacus petals in different states

Overview of the transcriptome sequence

Transcriptome sequencing was performed on 15 ‘Qiansiban’ samples at 5 stages. A total of 89 GB of clean data were obtained, and the amount of clean data obtained for each sample reached 6.616 GB. The GC content ranged from 44.34 to 45.62%, and the percentage of Q30 bases ranged from 93.62 to 95.03% (Supplementary Table S2).

Fig. 2
figure 2

Overview of the H. syriacus petal transcriptome. (A) Statistics of DEGs between petals harvested at different developmental stages. (B) Venn diagram representing the common and specific DEGs identified in the three comparisons. G0: S1 vs. S2; G1: S1 vs. S3; G2: S1 vs. S4; G3: S1 vs. F5; S1–S5 represent flower development stages

Functional annotation of new genes

The DIAMOND [17] code was used to match the newly discovered sequences with those in the NR [18], Swiss professional [19], COG [20], KOG [21], and KEGG [22] databases, and the KEGG origins and alternative results of the new sequences were obtained. InterProScan [23] used the information component of InterPro to investigate the GO Orthology [24] results of the latest genes and foreseen aminoalkanoic acid sequences of the latest genes and subsequently compared them with the Pfam [25] database victimization HMMER [26] code to obtain annotation data from the latest genes.

Based on the chosen order, the mapped reads were spliced with the NeighborTie code and were compared with the first-order annotation data. A total of 32,956 new genes were discovered, of which 23,322 were annotated. The numbers of genes annotated in different databases ranged from 4319 to 23,235, with the TrEMBL annotation having the highest number of genes and the COG annotation having the lowest number of genes (Supplementary Table S7).

Analysis of DEGs identified in the four libraries

With the thresholds FDR < 0.01 and FC > 2, 3456 DEGs (3,053 upregulated and 403 downregulated) were identified between the FS1 and FS2 libraries. Moreover, the number of DEGs between the FS1 and FS3 libraries was the highest, with 27,736 DEGs (12,655 upregulated and 15,081 downregulated). A total of 25,968 DEGs (12,124 upregulated and 13,844 downregulated) were identified between the FS1 and FS4 libraries, and 20,937 DEGs (10,013 upregulated and 10,924 downregulated) were identified between the FS1 and FS5 libraries (Fig. 2A). There were 2,702 common DEGs among the 5 comparisons (Fig. 2B).

GO enrichment analysis was also conducted on the screened differentially expressed genes, and Fig. 3A shows the top 20 pathways associated with the genes with the highest GO enrichment. The DEGs were distributed in three primary categories: biological process, cell composition, and molecular operation. The genes were enriched in cell process, biological process, and single-organism process in the biological process category; catalytic activity, binding, and transporter activity in the cell composition category; and cell, cell part, and organelle in the molecular operation category.

A differential factor KEGG metabolic pathway map was integrated into a complete metabolic pathway by classifying the DEGs [27,28,29]. The KEGG metabolic pathway enrichment of the ‘Qiansiban’ differential factors primarily included essential amino acid synthesis, flavonoid synthesis, cyanide metabolism, flavonoid and flavonol synthesis, essential amino acid metabolism, and anthocyanin synthesis. The metabolic pathways associated with flower colour formation included essential amino acid synthesis and metabolism, flavonoid and flavonol synthesis, and anthocyanin synthesis (Fig. 3B).

By analysing the transcriptome sequences of H. syriacus, combined with the KEGG annotation, 64 differentially expressed genes were obtained from the 5 stages of ‘Qiansiban’. The gene with the greatest variety was UFGT, with fourteen, and the gene with the smallest variety was F3’5’H, with one (Fig. 3C).

Fig. 3
figure 3

Differential gene analysis (A) Differential expression gene GO annotation classification statistical chart. (B) Metabolic pathway map of KEGG differential genes in H. syriacus ‘Qiansiban’. (C) Number of differentially expressed genes

Identification of structural genes involved in the flavonoid synthesis pathway

Transcriptome analysis revealed that sixty-four genes were considerably enriched in the flavonoid biogenesis pathway. The accumulation patterns of these anthocyanidins at different stages were consistent. The relative expression level of PAL sequentially decreased after its highest expression was reached at S2, and its expression then increased again at S5. The relative expression levels of the CHS1 and ANS1 genes also reached their highest levels in S2, with the expression of the CHS1 gene gradually decreasing throughout the S2–S5 stages and the expression of the ANS1 gene increasing in the S4 stage and then decreasing again in the S5 stage. The relative expression level of the CHI1 gene gradually decreased in the S1–S4 stages and then increased in the S5 stage. The relative expression level of the FLS1 gene increased incrementally through the S1–S4 stages and decreased in the S5 stage. The relative expression of the DFR1 gene decreased incrementally through the S1–S5 stages. The relative expression of the UFGT1 gene was the lowest in the S2 stage and increased in the S2–S5 stages. The relative expression of F3H, which was associated with the F3’5’H gene, showed an increasing decreasing trend, reaching its highest values in the S4 and S2 stages (Figs. 4 and 5, Supplementary Tables S3–S6).

According to the GO and KEGG enrichment maps constructed for the DEGs, nine structural genes were found to play key roles in the change in flower colour in ‘Qiansiban’, namely, UFGT1 (gene-LOC120143098), DFR1 (gene-LOC120216676), CHI1 (gene-LOC120193176), FLS1 (gene-LOC120188023), ANS1 (gene-LOC120214326), CHS1 (Hibiscus_syriacus_newGene_21723), PAL1 (Hibiscus_syriacus_newGene_32513), F3’5’H (Hibiscus_syriacus_newGene_29100) and F3H (gene-LOC120154124).

Changes in transcription factor expression

Transcription factors such as MYB, bHLH, and WD40 (MBW) affect the colour of H. syriacus flowers by regulating the synthesis of anthocyanins. According to transcriptome sequencing and classification statistics, the expression levels of five transcription factors, namely, AP2-ERF-EFR, MYB, bHLH, C2H2, and NAC (Supplementary Figure S1), were relatively high in ‘Qiansiban’. Among them, there were 279 MYB transcription factors and 203 bHLH transcription factors (Supplementary Table S4).

Based on the GO and KEGG enrichment results, the transcription factors (TFs) encoded by the DEGs and nine structural genes were subjected to protein‒protein interaction analysis to screen for crucial TFs that control anthocyanin synthesis. The results were as follows (Supplementary Table S4): bHLH2 (gene-LOC120149108), bHLH4 (gene-LOC120171442), bHLH6 (gene-LOC120209233), bHLH7 (gene-LOC120215242), MYB1 (Hibiscus_syriacus_newGene_28382), MYB3 (gene-LOC120154943), and MYB5 (gene-LOC120214348) may participate in the regulation of the seven structural genes.

Fig. 4
figure 4

Heatmaps (log 2(fold change)) of the anthocyanin biosynthesis genes. S1 (splitting stage), S2 (early stage), S3 (initial opening period), S4 (blooming period), and S5 (final flowering stage) represent the flower development stages

qRT‒PCR analysis and validation of anthocyanin synthesis structural genes

To verify the accuracy of the RNA-Seq results, we selected seven structural genes and seven regulatory genes involved in flavonoid metabolism and analysed their relative expression levels in the five flowering stages through qRT‒PCR (Fig. 5). The qRT‒PCR results were consistent with the trends in the relative expression levels of the genes identified in the sequencing data (Supplementary Tables S3, S4), indicating that the RNA-Seq results were reliable.

Fig. 5
figure 5

qRT‒PCR expression of anthocyanin-related genes. (A) The relative expression levels of the structural genes; (B) The relative expression levels of the transcription factors. The lowercase letters on the histogram represent differences at the 0.5 level

Fig. 6
figure 6

Differential regulation of anthocyanin biosynthesis in H. syriacus petals. The genes and metabolites that were differentially expressed and accumulated, respectively, are notated in red text. The dotted arrows/lines indicate that the genes were not found or are not known. The transcriptome and metabolome data were analysed according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis, and a flow chart was drawn. The Kyoto Encyclopedia of Genes and Genomes has authorized the use of these images. Each of the 9 colour bars includes 4 blocks, and the 4 blocks from left to right represent the trends of change of the relative expression levels of related genes in the S2, S3, S4, and S5 stages compared with the S1 stage. The colours of the blocks from light to deep represent downregulated, normal, and upregulated relative expression, respectively. The gene name abbreviations correspond to the full names given in Supplementary Table S6

Combined analysis of the transcriptome and metabolome

As summarized in Supplementary Table S5, bHLH2, bHLH4, and bHLH6 were significantly negatively correlated with UFGT. bHLH6 and bHLH7 were significantly positively correlated with PAL, DFR, and CHI. These findings indicate that increases in the expression levels of these two bHLH genes will also promote the synthesis of anthocyanins. bHLH2 and bHLH4 were significantly positively correlated with ANS, F3’5’H, and DFR, respectively, indicating that the expression levels of some structural genes also increased with increasing flower height, thereby promoting the accumulation of anthocyanins. In terms of correlations with total anthocyanins, four bHLH and MYB3 genes were significantly positively correlated with total anthocyanins, indicating that they can also promote the generation of total anthocyanins, thereby affecting petal colour.

Figure 5 shows that the changes in the expression levels of DFR and bHLH4 were the same as the changes in the anthocyanin contents during flowering (Fig. 1), indicating that these structural genes may be induced by transcription factors to regulate the colour of ‘Qiansiban’ petals.

Discussion

Petal colour transition is attributed to changes in the most highly expressed anthocyanins

By using a combination of transcriptomic and metabolomic methods, the reasons for the changes in colour of ‘Qiansiban’ during flowering were investigated. The results revealed that the changes in colour of the ‘Qiansiban’ fruits were closely related to the pigment composition and contents in the petals. The composition and contents of anthocyanins change continuously during development, which in turn affects the external colour of the petals. The anthocyanin contents decreased from S1 to S5, and the colour of the petals changed from purplish red to purple (Fig. 1). Y Q et al. [30] reported that the petal colours of two peony varieties changed from red to yellow during the flowering process, mainly due to the drastic decrease in anthocyanin content. Similarly, purple jasmine (Mirabilis jalapa L.) flowers change from light green in the bud stage to purple‒red in the blooming stage, which is mainly related to a decrease in chlorophyll content and an increase in anthocyanin content [2]. These studies further support our conclusions. In this study, a total of five anthocyanins (Ma3G, De3G, Pe3G, Cy3G, and Pg3G) were detected in ‘Qiansiban’ (Fig. 1; Table 1), which are the same as the types of anthocyanins previously identified in H. syriacus petals [8]. The proportion of total anthocyanins in these plants changes constantly, and because the anthocyanins can exhibit different colours, this change in total anthocyanin content also affects petal colour.

The petals of ‘Qiansiban’ contain flavone and flavonol substances throughout the flowering process. However, their colours are not visible in the flower. Flavones and flavonols are light yellow and are speculated to be chiefly used as auxiliary pigments in H. syriacus petals. This phenomenon is also common in other plants, as luteolin, apigenin-7-O-glucoside, kaempferol-3-O-glucoside, and different substances are expressed at high levels in Camellia nitidissima C. W. Chi [31] and florist’s chrysanthemum [32].

The regulatory effects of structural genes on the petal colour transition of H. Syriacus

The early biogenic genes of anthocyanins include mainly genes located upstream of the anthocyanin synthesis pathway, chiefly PAL, CHS, and CHI, which are responsible for encoding flavonoids and flavonols and different flavonoid synthetases. DFR, ANS, UFGT, and other downstream genes are mainly responsible for encoding anthocyanin synthase [33]. In this study, we conducted a transcriptomic analysis of five stages (S1, S2, S3, S4, and S5) of ‘Qiansiban’ and identified 9 genes that were significantly differentially expressed during these processes. These functional genes were significantly overexpressed in the S2 stage and contributed to the accumulation of anthocyanins. However, after the S3 stage, the expression levels of these genes decreased, leading to significant decreases in the anthocyanin contents in the later stages of development (Fig. 6). These findings are correspond with the overall changes in the expression levels of related genes in Paeonia lactiflora Pall [3] and Paeonia × suffruticosa Andrews [10]. The expression levels of genes in P. lactiflora and peony reach their maximum in the early stages of flowering and gradually decrease as their anthocyanin contents decrease throughout the entire flowering process.

FLS is the key factor involved in flavonol synthesis. Additionally, to aid in anthocyanin colouration, flavonols appear as a very light yellow once they are highly abundant. However, whether the expression of FLS causes an increase in flavonol content has yet to be verified. The DFR gene chiefly catalyses dihydroflavonol to synthesize colourless anthocyanins, whereas the ANS gene chiefly catalyses colourless anthocyanins to make coloured anthocyanins [34]. However, the high expression levels of most genes in S2, S3, and different early flowering stages promoted the synthesis of anthocyanins in the later stages. Reductions in organic phenomena within the later 2 stages resulted in a reduction in anthocyanin content, leading to changes in flower colour [1]. The analysis of Chrysanthemum × morifolium Ramat [35] and different decorative plants verified these findings.

Relationships between transcription factors and regulatory genes

The transcription factors MYB and bHLH can affect flower colour by regulating the expression of structural genes in the anthocyanin metabolism pathway. This study identified a total of 3 MYBs and 4 bHLHs. Correlation analysis revealed that the MYB3 and 4 bHLH genes were closely related to functional genes. These genes interact with structural genes to increase or decrease their expression, thereby affecting the flower colour of ‘Qiansiban’. In recent years, studies on Lilium davidii var. willmottiae (E. H. Wilson), Raffill [36], Prunus persica (L.) Batsch [37], Morus alba L. [38], and different plants have shown that MYBs and bHLHs are key transcription factors involved in plant anthocyanin synthesis. These findings also confirm the role of these genes in ‘Qiansiban’.

Although the MYB and bHLH transcription factors are involved in anthocyanin regulation, the regulatory patterns of different transcription factors are very different. Some MYBs directly bind to the promoters and promote the expression of structural genes and subsequently activate the entire anthocyanin expression network [39]. This type of interaction is very common between MYB and genes regulating the early stages of anthocyanin biosynthesis, such as CHS, CHI, and F3’H [40]. However, some MYB transcription factors require the interaction of other transcription factors, such as WD40 and bHLH, to induce the expression of structural genes, thereby affecting anthocyanin accumulation [41, 42]. In particular, late biosynthesis genes such as DFR and UFGT are generally regulated by the MBW-TF ternary protein complex of MYB-bHLH-WD40 [43]. Therefore, simply detecting the correlations between MYB and structural genes may not yield the desired results. This distinction may also explain why MYB5 was not significantly correlated with any of the structural genes. Niu et al. reported that the activation of the DFR promoter and the promotion of anthocyanin synthesis by MrMYB1 in Myrica rubra (Lour.) S. et al. are dependent on the coexpression of bHLH [44]. Similarly, in the absence of the FhbHLH cofactor, FhMYB5 had no significant effect on activating FhDFR3, FhCHI2, FhCHS1, FhF3’5’H, or FhF3H in Freesia hybrida Klatt [45].

Conclusion

Through metabolome and transcriptome sequencing and qRT‒PCR analysis, the colour changes of H. syriacus ‘Qiansiban’ petals from S1 to S5 were studied. The colour of its petals changes from red to pink‒purple and eventually to purple, which is attributed to the reduction in anthocyanin content. The anthocyanins expressed include delphinidin-3-O-glucoside, cyanidin-3-O-glucoside, petunidin-3-O-glucoside, pelargonidin-3-O-glucoside, and malvacin-3-O-glucoside. According to the transcriptome sequencing results, the interactions of structural genes (DFR, CHI, and FLS) and regulatory genes (MYB3 and bHLH2, bHLH4, bHLH6, and bHLH7) regulates the petal colour changing process of H. syringae.

Data availability

All transcriptomic sequencing data associated with this study have been submitted to the NCBI SRA under the accession number PRJNA995230.

References

  1. Guo LP, Wang YJ, et al. Transcriptome and chemical analysis reveal putative genes involved in flower color change in Paeonia ‘Coral Sunset’. Plant Physiol Biochem. 2019;138130–9. https://doi.org/10.1016/j.plaphy.2019.02.025.

  2. Zhou B, Yan XH, Lei YH, et al. Variation of flower color and pigment contents of Mirabilis jalapa during flowering. Bull Bot Res. 2022;42(3):475–82. https://doi.org/10.7525/j.issn.1673-5102.2022.03.017.

    Article  CAS  Google Scholar 

  3. Zhao D, Tao J, Han C, et al. Flower color diversity was revealed by differential expression of flavonoid biosynthesis-related genes and flavonoid accumulation in herbaceous peony (Paeonia lactiflora Pall). Mol Biol Rep. 2012;39(12):11263–75. https://doi.org/10.1007/s11033-012-2036-7.

    Article  CAS  PubMed  Google Scholar 

  4. Zheng XD, Pan XX, Meng XM, et al. Research progress in the determination of carotenoids in Lycium chinense by HPLC. China Fruit Vegetable. 2013;1027–30. https://doi.org/10.3969/j.issn.1008-1038.2013.10.013.

  5. Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54(4):733–49. https://doi.org/10.1111/j.1365-313X.2008.03447.x.

    Article  CAS  PubMed  Google Scholar 

  6. Deng J, Chen S, Yin XJ, et al. Systematic qualitative and quantitative assessment of anthocyanins, flavones, and flavonols in the petals of 108 lotus (Nelumbo nucifera) cultivars-ScienceDirect. Food Chem. 2013;139(1–4):307–12. https://doi.org/10.1016/j.foodchem.2013.02.010.

    Article  CAS  PubMed  Google Scholar 

  7. Li X, Duan JJ, Luo XN, et al. Formation mechanism of different tree peony flower colors by anatomy and biochemistry. J Northeast Forestry Univ. 2019;43(3):38–43. https://doi.org/10.3969/j.issn.1000-5382.2019.03.008.

    Article  Google Scholar 

  8. Zhang P, Li Y, Chong SL et al. Identification and quantitative analysis of anthocyanin composition and stability from different strains of Hibiscus syriacus L. flowers. Industrial Crops Prod. 2022; (177):114457. http://doi.org/10.1016/j.indcrop.2021.114457

  9. Chen JL, Hou RM, Zhu JJ, et al. RNA sequence analysis of petals of two cultivars of Hibiscus syriacus with different flower colors. Mol Plant Breed. 2022;20(8):2507–16. https://doi.org/10.13271/j.mpb.020.002507.

    Article  Google Scholar 

  10. Yang C. Analysis of petal pigment composition and related gene expression during flower color change of ‘Fengdan’ Peony. Yangling: Northwest University of Agriculture and Forestry Science and Technology; 2020. http://doi.org/10.27409/d.cnki.gxbnu.2020.001010

    Google Scholar 

  11. Chen H. Analysis of anthocyanin components and antioxidant capacity of Hibiscus syriacus with different flower color. Qinhuangdao: Hebei Normal University of Science and Technology; 2021. http://doi.org/10.27741/d.cnki.ghbkj.2021.000025

    Google Scholar 

  12. Wu C, Dai MY, Zhang C et al. Study on the Mechanism of the FLS Gene Regulating Flower Color Formation in Yulania denudata and Y. liliiflora. Journal of Nuclear Agricultural Sciences. 2023,37(10):1947–1956. https://doi.org/10.11869/j.issn.1000-8551.2023.10.1947.

  13. He YF, Zhang CY, Liu QL. Characterization of flower coloration in petals of Rhododendron Simsii Planch. J Agricultural Sci Technol. 2021;23(2):50–6. https://doi.org/10.13304/j.nykjdb.2019.0693.

    Article  CAS  Google Scholar 

  14. Gao S, Li TL. Relationship between the composition of anthocyanins and flower color variation in lisianthus(Eustoma grandiflorum). J Shenyang Agricultural Univ. 2020;51(4):410–6. https://doi.org/10.3969/j.issn.1000-1700.2020.04.004.

    Article  CAS  Google Scholar 

  15. Cui HL, He X, Zhang Q. Dynamic changes of anthocyanin and flavonoid composition in petals of different peony varieties during account opening. Scientia Agricultura Sinica. 2021;54(13):2858–69. https://doi.org/10.3864/j.issn.0578-1752.2021.13.014.

    Article  Google Scholar 

  16. Wang XT, Li XG, Li DL, et al. Identification and analysis of nonanthocyanin phenolics in various raspberry (Rubus spp.) cultivars. North Hortic. 2021;1136–43. https://doi.org/10.11937/bfyy.20203753.

  17. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;1259–60. https://doi.org/10.1038/NMETH.3176.

  18. Deng YY, Li JQ, Wu SF, Zhu YP, et al. Integrated Nr Database in protein annotation system and its localization. Comput Eng. 2006;32(5):71–4. https://doi.org/10.1109/INFOCOM.2006.241.

    Article  Google Scholar 

  19. Apweiler R, Bairoch A, Wu CH, et al. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32:D115–9. https://doi.org/10.1093/nar/gkw1099.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tatusov RL, Galperin MY, Natale DA. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28(1):33–6. https://doi.org/10.1093/nar/28.1.33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Koonin EV, Fedorova ND, Jackson JD, et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 2004;5(2):R7. https://doi.org/10.1186/gb-2004-5-2-r7.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Kanehisa M, Goto S, Kawashima S, Okuno Y, et al. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:D277–80. https://doi.org/10.1093/NAR/GKH063.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Jones P, Binns D, Chang HY, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40. https://doi.org/10.1093/bioinformatics/btu031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ashburner M, Ball CA, Blake JA et al. Gene ontology: a tool for the unification of biology. Nature genetics. 2000; 25(1): 25–29. https://doi.org/10.1038/75556.

  25. Finn RD, Bateman A, Clements J et al. Pfam: the protein families database. Nucleic Acids Res. 2014;(42):D222-D230. http://doi.org/10.1093/nar/gkt1223.

  26. Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14(9):755–63. https://doi.org/10.1198/016214502388618870.

    Article  CAS  PubMed  Google Scholar 

  27. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30. https://doi.org/10.1093/nar/27.1.29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51. https://doi.org/10.1002/pro.3715.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–592. http://doi.org/10.1093/nar/gkac963.

    Article  CAS  PubMed  Google Scholar 

  30. Yang Q, Yuan T, Sun XB. Preliminary studies on the changes of flower color during the flowering period in two tree peony cultivars. Acta Horticulturae Sinica. 2015;42(5):930–8. https://doi.org/10.16420/j.issn.0513-353x.2014-0996.

    Article  CAS  Google Scholar 

  31. Li XL, Wang JT, Sun ZY, et al. The relationship between flavonoid components and flower color of Camellia japonica L, Camellia Alba, and their three hybrid varieties. J Hortic. 2019;46(6):1145–54. https://doi.org/10.1139/G82-039.

    Article  Google Scholar 

  32. Xie ZF, Zhang QQ, Zhu LJ, Han G. Research progress on chemical components and pharmacological activities of Chrysanthemum × morifolium Ramat. J Henan Univ (Medical Edition). 2015;34(4):290–300. http://doi.org/10.15991/j.cnki.41-1361/r.2015.04.018

    CAS  Google Scholar 

  33. Luo TJ. Research progress on genes related to plant Anthocyanidin metabolism pathway. South Agric. 2022;16(16):202–6. https://doi.org/10.19415/j.cnki.1673-890x.2022.16.064.

    Article  Google Scholar 

  34. Dai SL, Hong Y. Molecular breeding for flower color improvement of ornamental plants based on the mechanism of Anthocyanidin glycoside synthesis and color rendering. China Agricultural Sci. 2016;49(3):529–42. https://doi.org/10.3864/j.issn.0578-1752.2016.03.011.

    Article  CAS  Google Scholar 

  35. Han KT. Effects of the introduction of key structural genes for anthocyanidin glycoside synthesis on Chrysanthemum × morifolium Ramat flower color. Beijing: Beijing Forestry University; 2010.

    Google Scholar 

  36. Fan WG, Li BY, Tian H et al. Screening and analysis of MYB and bHLH transcription factors related to Anthocyanidin synthesis in lily bulbs of Lanzhou [J/OL]. Food and fermentation industry: 1–12 [2023-08-03]. https://doi.org/10.13995/j.cnki.11-1802/ts.033270https://kns.cnki.net/kns8/defaultresult/index.

  37. Hua Y, Lv XY, Zhang Y, et al. Research progress of peach MYB transcription factor regulating Anthocyanidin synthesis. Mod Agricultural Sci Technol. 2022;2175–9. https://doi.org/10.3969/j.issn.1007-5739.2022.21.019.

  38. Fan YT. Research on the regulation and mechanism of Mulberry MYB61 Gene on Anthocyanidin and Lignin Synthesis. Tai’an: Shandong Agricultural University; 2022. http://doi.org/ 10.27277/d.cnki.gsdnu.2022.000245.

    Google Scholar 

  39. Li W, Ding Z, Ruan M, Yu X, Peng M, Liu Y, Kiwifruit. R2R3-MYB transcription factors and contribution of the novel AcMYB75 to red kiwifruit anthocyanin biosynthesis. Sci Rep. 2017;7(1):16861. https://doi.org/10.1038/s41598-017-16905-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Pelletier MK, Burbulis IE, Winkel-Shirley B. Disruption of specific flavonoid genes enhances the accumulation of flavonoid enzymes and end-products in Arabidopsis seedlings. Plant Mol Biol. 1999;40(1):45–54. https://doi.org/10.1023/a:1026414301100.

    Article  CAS  PubMed  Google Scholar 

  41. Albert NW, Lewis DH, Zhang H, Schwinn KE, Jameson PE, Davies KM. Members of an R2R3-MYB transcription factor family in Petunia are developmentally and environmentally regulated to control complex floral and vegetative pigmentation patterning. Plant J. 2011;65(5):771–84. https://doi.org/10.1111/j.1365-313X.2010.04465.x.

    Article  CAS  PubMed  Google Scholar 

  42. Zhou H, Liao L, Xu S, Ren F, Zhao J, Ogutu C, Wang L, Jiang Q, Han Y. Two amino acid changes in the R3 repeat cause functional divergence of two clustered MYB10 genes in peach. Plant Mol Biol. 2018;98(1–2):169–83. https://doi.org/10.1007/s11103-018-0773-2.

    Article  CAS  PubMed  Google Scholar 

  43. Ilk N, Ding J, Ihnatowicz A, Koornneef M, Reymond M. Natural variation for anthocyanin accumulation under high-light and low-temperature stress is attributable to the ENHANCER OF AG-4 2 (HUA2) locus in combination with PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2. New Phytol. 2015;206(1):422–35. https://doi.org/10.1111/nph.13177.

    Article  CAS  PubMed  Google Scholar 

  44. Niu SS, Xu CJ, Zhang WS, Zhang B, Li X, Lin-Wang K, Ferguson IB, Allan AC, Chen KS. Coordinated regulation of anthocyanin biosynthesis in Chinese bayberry (Myrica rubra) fruit by an R2R3 MYB transcription factor. Planta. 2010;231(4):887–99. https://doi.org/10.1007/s00425-009-1095-z.

    Article  CAS  PubMed  Google Scholar 

  45. Li Y, Shan X, Zhou L, Gao R, Yang S, Wang S, Wang L, Gao X. The R2R3-MYB factor FhMYB5 from Freesia Hybrida contributes to the regulation of anthocyanin and Proanthocyanidin Biosynthesis. Front Plant Sci. 2019;9(1935). https://doi.org/10.3389/fpls.2018.01935.

Download references

Funding

The work was supported financially by Hebei Agriculture Research System(HBCT2024200202)the Survey of Herbaceous Germplasm Resources in Shandong Province (Grant No. Lu Financial [2021]1); Hebei Provincial Department of Science and Technology Project (21326301D); and the Germplasm Resources Bank of the National Forestry and Grassland Administration (2005DKA210003).

Author information

Authors and Affiliations

Authors

Contributions

Zhezhe LI prepared samples for analysis, determined and analyzed the data, interpreted the results, and wrote the original draft paper. Dan LIU determined and analyzed the data and interpreted the results. Meng SUN, Dongsheng WANG, and Yu WU measured supplementary experimental data. Yidan ZHANG collected samples and assisted. Dongsheng WANG and Beibei CHENG designed the study, analyzed the data, interpreted the results, and revised the manuscript. Guojun ZHANG made revisions to the first draft of the article. All authors reviewed the manuscript.

Corresponding author

Correspondence to Beibei Cheng.

Ethics declarations

Ethics approval and consent to participate

Our research did not involve any human or animal subjects, material, or data. We declare that the plant material used in the experiment was collected and studied by relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Conflict of interest

The authors hereby declare no conflicts of interest.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Z., Liu, D., Wang, D. et al. Study on the causes of changes in colour during Hibiscus syriacus flowering based on transcriptome and metabolome analyses. BMC Plant Biol 24, 431 (2024). https://doi.org/10.1186/s12870-024-05142-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05142-0

Keywords