Time-Series Transcriptome Sheds Light on Anthocyanin Metabolism and Species Diversification Between two Lonicera Japonica Thunb Cultivars

Time-Series Transcriptome Sheds Light on Anthocyanin Metabolism and Species Diversification Between two Lonicera Japonica Thunb Cultivars | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Time-Series Transcriptome Sheds Light on Anthocyanin Metabolism and Species Diversification Between two Lonicera Japonica Thunb Cultivars Jianjun Li, Chenglin Ye, Jingxiao Ma, Ting Cheng, Yan Lv, Cuifang Chang This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-104094/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background：‘Yujin 2’ is new variety of Lonicera japonica Thunb and its flower color can change from red to yellow; hence, it is a good model for investigating flower color development mechanisms. Results：High throughput transcriptome sequencing of seven flower development stages of Yujin No.2 was carried out, and 133,487 unigenes were annotated, among which 73,088 were differentially expressed. Then the real-time PCR analysis was carried out. Further, the number of up-regulated DEGs was higher than those that were down-regulated. Of these annotated DEGs, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis were active throughout the flowering process during each stage, whereas carotenoid biosynthesis was inactive in the S1-6 stages. Furthermore, phenylalanine synthesis was enhanced in the S1 phase; however, anthocyanin synthesis was weakened in the S5 and S6 phases, which may be consistent with the changes in petal color of ‘Yujin 2’ from red (S1) to white (S5) and gold (S6). The results showed that 114 unigenes were associated with anthocyanin metabolism, and 72 were significantly upregulated or downregulated. According to the analysis of TFs in anthocyanin metabolism, we obtained 47 transcription factors, which belonged to 18 families. The LjDFR, LjABCB1, LjMYC6, LjDDB2, and LjANS genes rapidly increased during the first three stages. However, only LjF3'5'H expression was significantly down-regulated at S5, which was consistent with anthocyanin accumulation. Conclusions：This study developed a transcriptome profile of flower color generation for L. japonica as well as annotated unigene sets of seven anthesis phases, thereby providing possibilities for improving the germplasm of L. japonica with genetic engineering technologies and cultivating new varieties of different colors. Plant Physiology and Morphology Plant Molecular Biology and Genetics Lonicera Japonica Thunb Anthocyanin Yujin 2 Comparative transcriptomics RNA-seq Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Background Lonicera japonica Thunb ( L. japonica ) is a perennial semi-evergreen twining stoloniferous shrub. Traditionally, the L. japonica flower bud, which has been listed in the Chinese Pharmacopoeia as Lonicera japonica (Jinyinhua in Chinese) is considered to be the primary medicinal component [ 1 ]. Recent studies have found that L. japonica has many medicinal properties due to its antioxidant [ 2 ], hypoglycemic and hypolipidemic [ 3 ], anti-allergic [ 4 ], anti-inflammatory [ 5 ], and antibacterial effects [ 6 ]. Further, the demand for L. japonica in pharmaceuticals [ 1 ], cosmetics [ 7 ], and health foods [ 8 ] has increased. Although L. japonica is widely distributed in China, its yields and quality are the highest in Fengqiu of Henan and Shandong Provinces. Therefore, our research group selected and bred a new variety (Yujin 2). It has the characteristics of red buds, strong resistance, and a high content of chlorogenic acid (CGA) and luteoloside. The red buds endow ‘Yujin 2’ with enhanced beauty, which is more valuable in utility [ 9 ]. As one of the most important phenotypic traits of ornamental plants, the formation of flower color is the result of the combined effect of environmental influences and genetic characteristics [ 10 ]. The modification of petal color regulation genes is indispensable in modern flower breeding [ 11 ]. Through various studies it has been established that changes in flower color are primarily due to alterations in plant pigment metabolism [ 12 ]. The key pigments to adjust the color are carotenoids, flavonoids, and alkaloids [11; 13; 14]. Anthocyanins, which are produced through flavonoid biosynthesis pathways during plant metabolism, are important water-soluble natural food pigment, which is widely present in plant epidermal cell vacuoles [ 15 ]. In the first stage of anthocyanin biosynthesis, phenylalanine transforms into 4-coumaryl:CoA, referred to as the phenylpropanoid metabolic pathway. The next stage, 4-coumaryl:CoA transforms into multiple flavonoid compounds [ 16 ]. The biosynthesis of anthocyanins is related to a variety of structures and regulatory genes [ 17 ]. The expression of chalcone synthase (CHS), Flavanone 3-hydroxylase (F3H), dihydroflavonol-4-reductase (DFR), anthocyanin synthase (ANS) and UDP glucoflavonoid 3-o-glucosyltransferase (UFGT) directly affect the accumulation of anthocyanins in plants [ 18 ]. These structural genes have been cloned in different plants, and the expression pathway is common to different types of plants. Key regulatory genes serve as regulators that combine individually, or in the form of dimers and trimers, which play a critical role in the regulation of anthocyanin synthesis. So far, MYB, bHLH and WD40 transcription factors are widely regarded as the three main regulatory factors [ 19 ]. The MYB family transcription factor is one of the largest transcription factor families in plants, which has a MYB domain that is composed of 1 ~ 4 MYB repeating units in its structure [ 20 ]. In conjunction with MYB, the bHLH transcription factors play active and critical roles in anthocyanin biosynthesis. The high anthocyanin content of many crops results from the increased expression of bHLH transcription factors [ 10 ]. WD40 is an ancient and highly conserved family of proteins that is involved in myriad plant cellular processes, such as transcriptional regulation, cell division, vesicle formation and transport, processing, signal transduction and so on [ 21 ]. There were great differences in the corolla color of ‘Yujin 2’ at the various developmental stages. The upper half of the flower bud during the early stage was purple-red, whereas the whole flower bud from the three green periods to the great white period was red, but gradually became lighter; whereas during the silver period, the exterior was red and the interior was white. During the golden period, the exterior of ‘Yujin 2’ was red and the interior was white, but for ‘Damaohua’ it was yellow [ 9 ]. Color is considered to be an important indicator in the assessment of ornamental value and offers a basis for plant classification [ 11 ]. However, the molecular mechanisms of anthocyanin metabolism in L. japonica flowers have not yet been discovered. Therefore, this paper employed ‘Yujin 2 and ‘Damaohua’ flowers at different stages of development to study the molecular mechanisms of anthocyanin formation. This, by means of morphology, transcriptome, and bioinformatics, etc., to provide possibilities for enhancing the L. japonica germplasm using genetic engineering technologies, and cultivating new varieties with different colors. Results Characteristics of ‘Damaohua’ and ‘Yujin 2’ transcriptome A total of 42 samples (‘Damaohua’ and ‘Yujin 2’, respectively, have seven stages, each with three independent replicates) were sequenced using the Illumina HiSeq X Ten sequencer (Illumina Inc., USA), from which 311 Gb of raw data was harvested. Using a pair-end sequencing strategy, 2G raw reads with mean lengths of 150 nt were obtained from the 42 samples with 47.62 M reads for each sample. After the filtering of raw reads, 1.95 G reads were counted, with an average of 46.43 M reads per sample (Additional file 1: Table S1). Following the removal of redundant sequences and quantity control, 133,487 unigenes were assembled. The average length of the unigenes was 1,047 bp with 45.75% GC and 1642nt N50. The longest transcript was 16,864 bp, and the shortest was 301 bp. The unigenes were annotated to obtain functional information using seven separate databases, including NCBI non-redundant protein (NR), Clusters of Orthologous Groups (COG) for eukaryotic complete genomes, Gene Ontology (GO), Swiss-Prot, evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases by Diamond Software, and mapped to Pfam databases by HMMER. Overall, 66.79% (133,487 unigenes) were successfully annotated, among which most unigenes (66.46%) were assigned to annotation terms in the NR database (Table 1 ). Further, 61.43% of the unigenes were annotated on eggNOG, 46.62% on Swiss-prot, 43.46% on GO, 35.76% on KOG, and 26.02% on the KEGG pathway (Table 1 ). Table 1 Percentage of annotated unigenes. Database Count of annotated unigenes% NR 66.46% KOG 35.76% GO 43.46% Swiss-prot 46.62% eggNOG 61.43% KEGG 26.02% Pfam 0.13% Total 66.79% According to the NR annotation results, the current study compared L. japonica with taxonomic species stored in the NR database. A total of 97247 sequences were matched with L. japonica unigenes, which covered 1512 taxonomy categories (Additional file 2: Table S2). The highest ranked species was Alternaria alternate , which took up 12.82%, whereas the second highest-ranked was Vitis vinifera , which contributed 8.80% (Fig. 1 ). Differentially expressed genes ( DEGs ) When comparing seven flowering stages of expressed genes between ‘Damaohua’ and ‘Yujin 2’, 73,088 differentially expressed genes exhibited remarkable changes. A total of 5,355 genes were differentially expressed, with 2,864 genes up-regulated and 2491 genes down-regulated at the S1 stage. At the S2 stage, 21,159 genes were differentially expressed, with 12,734 genes up-regulated and 8,425 genes down-regulated. At the S3 stage, 19,126 genes were differentially expressed, with 10,343 genes up-regulated and 8,783 genes down-regulated. At the S4 stage, 11,153 genes were differentially expressed, with 6,029 genes up-regulated and 5,124 genes down-regulated. At the S5 stage, 27,718 genes were differentially expressed, with 12,102 genes up-regulated and 15,616 genes down-regulated. At the S6 stage, 18,652 genes were differentially expressed, with 8,486 genes up-regulated and 10,166 genes down-regulated. Finally, at the S7 stage, 38,162 genes were differentially expressed, with 21,153 genes up-regulated and 17,009 genes down-regulated (Fig. 2 ). As shown in the figures, the number of up-regulated DEGs was higher than those that were down-regulated. KEGG-pathway analysis of differentially expressed genes Unigenes were assigned to KEGG pathways based on annotation data. DEGs were enriched into 181, 202, 202, 191, 207, 204, and 210 pathways (S1, S2, S3, S4, S5, S6, and S7) (Additional file 3: Table S3). At the S1 stage, 53 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), and phenylpropanoid biosynthesis (ko00940) during the entire process, as well as up-regulation and down-regulation processes, played active roles, and the expression was significant. Amino sugar and nucleotide sugar metabolism (ko00520) were also active during the entire process and up-regulation, while they were not active in the down-regulation process. Starch and sucrose metabolism (ko00500) were active during the entire process and down-regulation process, while they were not active in the up-regulation process (Figs. 1 A). At the S2 stage, 72 pathways were enriched with a statistically significant enrichment p-value, which was more than the other six paired comparison phases. Endocytosis (ko04144), plant hormone signal transduction (ko04075), and starch and sucrose metabolism (ko00500) during the entire process, as well as up- and down-regulation processes played an active role, and the expression was significant. Phenylpropanoid biosynthesis (ko00940) was active during the entire process and the down-regulation process, while it decreased in the up-regulation process (Figs. 1 B). At the S3 stage, 66 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500) during the entire process, and the up- and down-regulation processes played an active role, and the expression is significant. Glycolysis/Gluconeogenesis (ko00010) was active during the entire process and the down-regulation process, while it was not active in the up-regulation process (Figs. 1 C). At the S4 stage, 71 pathways were enriched with a statistically significant enrichment p-value, and plant hormone signal transduction (ko04075) during the entire process, the up- and down-regulation processes played an active role, and the expression was significant. Phenylpropanoid biosynthesis (ko00940) and cell cycle (ko04110) were active during the entire process and the up-regulation process, while they were not active in the down-regulation process. Starch and sucrose metabolism (ko00500) were active during the entire process and the down-regulation process, while they were not active in the up-regulation process. Cell cycle-yeast (ko04111) was triggered in the up-regulation process, while carbon metabolism (ko01200) was triggered in the down-regulation process (Figs. 1 D). At the S5 stage, 63 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075) and phenylpropanoid biosynthesis (ko00940) played active role during the entire process as well as up- and down-regulation processes, and the expression was significant. Glycerophospholipid metabolism (ko00564), starch and sucrose metabolism (ko00500) were active in the total process and the up-regulation process, which were not active in the down-regulation process. Glutathione metabolism (ko00480), drug metabolism - cytochrome P450 (ko00982), and metabolism of xenobiotics by cytochrome P450 (ko00980) were triggered in the down -regulation process (Figs. 1 E). At the S6 stage, 31 pathways were enriched a with statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500), glycolysis/gluconeogenesis (ko00010) played an active role during the entire process, the up- and down-regulation process, and the expression was significant. Valine, leucine and isoleucine degradation (ko00280) was triggered in the up-regulation process (Figs. 1 F). At the S7 stage, 39 pathways were enriched with a statistically significant enrichment p-value. Fatty acid metabolism (ko01212) was active during the entire process, but was not in the up- and down-regulation process. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500) were active during the entire process and up-regulation process, while they were not in the down-regulation process. Tyrosine metabolism (ko00350) was active during the entire process and the down-regulation process, while it was not in the up-regulation process. Spliceosome (ko003040) was triggered in the up-regulation process, whereas MAPK signaling pathway-yeast (ko04011), as well as glycine, serine and threonine metabolism (ko00260) were triggered in the down -regulation process (Figs. 1 G). Summarizing the above KEGG pathways enriched during the flowering period, it was found that phenylpropanoid biosynthesis and flavonoid biosynthesis were active through the entire flowering process and nearly up-regulated at each stage, Peroxisome(ko04146), Carbon metabolism༈ko01200༉ and Amino sugar and nucleotide sugar metabolism༈ko00520༉can be found in each period(Fig. 3 A) ,and carotenoid biosynthesis was down-regulated from S1-6. The anthocyanin precursor, phenylalanine, was enhanced at the S1 phase; however, the synthesis of anthocyanin was weakened at the S5 and S6 phases (Table 2 ). Table 2 Enriched KEGG pathways by up-regulated and down-regulated genes. KEGG Pathways S1 S2 S3 S4 S5 S6 S7 Phenylpropanoid biosynthesis ↑↓ ↑↓ ↑↓ ↑ ↑↓ ↑↓ ↑ Carotenoid biosynthesis ↓ ↑↓ ↑↓ ↓ ↑↓ ↑ Flavonoid biosynthesis ↑ ↑↓ ↑ ↑ ↑ ↑ ↑ Phenylalanine metabolism ↑ ↓ Phenylalanine, tyrosine and tryptophan biosynthesis ↑ Anthocyanin biosynthesis ↓ ↓ *↑represents KEGG pathways enriched by up-regulated genes;↓represents KEGG pathways enriched by down-regulated genes. Gene Ontology (GO) analysis of differentially expressed genes Differentially expressed genes were annotated by Gene Ontology terms (Additional file 4: Table S4). Genes were classified into three GO categories, that is, Biological Process, Cellular Component, and Molecular Function. From the Biological Process, auxin-activated signaling pathway (GO:0009734), response to abscisic acid (GO:0009737), response to wounding (GO:0009611) were remarkable at the S1 stage. The negative regulation of the hydrogen peroxide metabolic process (GO:0010727), positive regulation of the receptor biosynthetic process (GO:0010870), cellular response to topologically incorrect protein (GO:0035967) were remarkable at the S2 stage. The regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836), nucleolar chromatin organization (GO:1990700), and regulation of superoxide anion generation (GO:0032928) were remarkable at the S3 stage. The priming of the cellular response to stress (GO:0080136), cellular response to sucrose starvation (GO:0043617), and the cellulose metabolic process (GO:0030243) were remarkable at the S4 stage. The negative regulation of the gibberellin biosynthetic process (GO:0010373), guanosine tetraphosphate metabolic process (GO:0015969), and the regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836) were remarkable at the S5 stage. The regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836), nucleolar chromatin organization (GO:1990700), and the establishment of protein localization to endoplasmic reticulum membrane (GO:0097051) were remarkable at the S6 stage. The negative regulation of the extrinsic apoptotic signaling pathway in absence of ligand (GO:2001240), the phosphatidylinositol-3- phosphate biosynthetic process (GO:0036092), and the regulation of the embryo sac egg cell differentiation (GO:0045694) were remarkable at the S7 stage (Figs. 2 ). From the Cellular Component process, apoplast (GO:0048046), plant-type cell wall (GO:0009505), and cell wall (GO:0005618) were remarkable at the S1 stage. Inclusion body (GO:0016234), nitrite reductase complex [NAD(P)H] (GO:0009344), and plasma membrane (GO:0005886) were remarkable at the S2 stage. Nucleolar chromatin (GO:0030874) vacuole (GO:0005773), and chloroplast (GO:0009507) were remarkable at the S3 stage. Anchored component of plasma membrane (GO:0046658), apoplast (GO:0048046), and plasma membrane (GO:0005886) were remarkable at the S4 stage. Nucleolar chromatin (GO:0030874), P4 peroxisome (GO:0019822), chloroplast (GO:0009507) were remarkable at the S5 stage. Nucleolar chromatin (GO:0030874), chloroplast (GO:0009507), vacuole (GO:0005773) were remarkable at the S6 stage. The DNA topoisomerase complex (ATP-hydrolyzing) (GO:0009330), nuclear ubiquitin ligase complex (GO:0000152), and TRAPPII protein complex (GO:1990071) were remarkable at the S7 stage (Figs. 2 ). From the Molecular Function process, DNA binding transcription factor activity (GO:0003700), aspartic-type endopeptidase activity (GO:0004190), hydrolase activity, and acting on ester bonds (GO:0016788) were remarkable at the S1 stage. Tau protein binding (GO:0048156), euchromatin binding (GO:1990188), and IAA-Phe conjugate hydrolase activity (GO:0010210) were remarkable at the S2 stage. IAA-Phe conjugate hydrolase activity (GO:0010210), IAA-Leu conjugate hydrolase activity (GO:0010211), and dehydroquinate synthase activity (GO:0102042) were remarkable at the S3 stage. Glucuronoxylan 4-O-methyltransferase activity (GO:0030775), DNA binding transcription factor activity (GO:0003700), sequence-specific DNA binding (GO:0043565) were remarkable at the S4 stage. Euchromatin binding (GO:1990188), malonyl-CoA decarboxylase activity (GO:0050080), and glutamine-phenylpyruvate transaminase activity (GO:0047316) were remarkable in the S5 stage; euchromatin binding (GO:1990188), IAA-Phe conjugate hydrolase activity (GO:0010210), IAA-Leu conjugate hydrolase activity (GO:0010211) were remarkable at the S6 stage. Arylacetonitrilase activity (GO:0047428), euchromatin binding (GO:1990188), and calcium-independent phospholipase A2 activity (GO:0047499) were remarkable at the S7 stage (Figs. 2 ). From seven periods,defense response to bacterium (GO:0042742)༌transcription, DNA-templated (GO:0006351)༌regulation of transcription, DNA-templated (GO:0006355)༌vacuole (GO:0005773)༌vacuolar membrane (GO:0005774)༌extracellular region (GO:0005576)༌chloroplast (GO:0009507)༌nucleus (GO:0005634)༌integral component of membrane (GO:0016021)༌plasma membrane (GO:0005886)༌cytoplasm (GO:0005737)༌sequence-specific DNA binding (GO:0043565)༌DNA binding transcription factor activity (GO:0003700)༌metal ion binding (GO:0046872)༌ATP binding (GO:0005524) were found at the every period(Fig. 3 B) Anthocyanin metabolism analysis of differentially expressed genes Sequencing results revealed that 114 unigenes were associated with anthocyanin metabolism, and 72 were significantly upregulated or downregulated across the seven flowering stages. The number of up-regulated unigenes was highest at the S7 stage, while the number of down-regulated unigenes was highest at the S5 stage (Fig. 4 A). According to the analysis of TFs in anthocyanin metabolism, we obtained 47 genes relating to transcription factors in anthocyanin metabolism that belonged to 18 families, of which bHLH had the largest number, followed by GeBP (Additional file 5: Table S5, Fig. 4 B). According to the anthocyanin heatmap analysis, the expression of genes at S7 was different from that at other stages, with the highest up-regulated expression. The expression of genes at S6 was similar to that at S5, while those at S4, S3, S2, and S1 were similar, among which the expression of genes in S1 and S2 were the most similar (Fig. 4 C). qRT-PCR verification To validate the accuracy and repeatability of our RNA-Seq data, the current study selected six genes for quantitative real-time polymerase chain reaction (qRT-PCR) analysis, with gene-specific primers designed using Primer software (version 5.0), as shown in Additional file 6: Table S6. As shown in Fig. 5 , the results of qRT-PCR indicated that most of these genes had expression patterns that correlated with the RNA-Seq data, which confirmed the reliability of our data. The genes LjDFR (TRINITY_DN37612_c0_g1_i2_29), LjABCB1 (TRINITY_DN23592_c0_g2_i2_8), LjMYC6 (TRINITY_DN17985_c0_g1_i2_8), LjDDB2 (TRINITY_DN25019_c0_g2_i6_8), and LjANS (TRINITY_DN27734_c0_g1_i1_12) rapidly increased during the first three stages. However, the expression of only LjF3'5'H (TRINITY_DN27086_c0_g1_i1_29) was significantly down-regulated at S5, which was consistent with anthocyanin accumulation. Discussion The biological regulation of flower color had garnered much attention from researchers, particularly in terms of ornamental plant research [ 11 ]. This present study focused on a new variety of L. japonica (Yujin no. 2). There are a variety of designs and colors for ‘Yujin 2’ that has a high ornamental value, which may play a significant role in landscaping [ 9 ]. It is very important to regulate anthocyanin synthesis in L. japonica breeding [ 22 ]. With respect to molecular breeding, the identification of the key genes involved in anthocyanin synthesis will enable the overexpression or knockout of these genes to facilitate color modification [ 23 ]. Therefore, this study began with the aim of identifying the essential genes that regulated the alteration of flower color and subsequently exploring the potential molecular mechanisms behind the development of flower color. However, little evidence has been reported as relates to L. japonica petal color. The current study results provide a foundation for research into petal color mechanisms at the transcriptome level. Transcriptome analysis based on an EST sequencing dataset illustrated the characteristics of the carnation transcriptome and revealed that carotenoid, chlorophyll, and anthocyanin biosynthesis played roles in the regulation of carnation flower color [ 24 ]. Correlation analysis showed that the changes in Paeonia flower color, from coral to pink to pale yellow were due to a significant decrease in anthocyanin content [ 25 ]. It was deduced that the disequilibrium of expression levels in flavonol synthases and dihydroflavonol-4-reductases resulted in different levels of anthocyanin accumulation and coloration in white and pink tea flowers via metabolome and transcriptome analysis [ 26 ]. Besides, the transcriptome dataset assisted with the establishment of linkages between gene expression and other aspects, including the phenotypical, biochemical and metabolic features of Meconopsis petal development [ 10 ]. Petal color change is a significant feature of gene expression. As a whole, the number of differentially expressed genes at each flowering development stage was variable, which revealed that transcriptome expression changes were more active during color alteration. For example, when the petal color changed from purple-red (S1) to red (S2), 15804 genes were differentially expressed; while only 2033 genes changed from S2 to S3, when the petal color at both were red. Further, when the petal colors were red, but gradually became lighter from the three-green period (S2) to the great white period (S4), the number of differentially expressed genes were eventually lower, which were 21159, 19126, and 11153, respectively. Besides, gene expression changes occurred when the petal color transitioned from red (S4), to the outer part being red and the inner part being white (S5), where 6073 genes were up-regulated from the bud stage to cracking stage (Fig. 2 ). Petal color phenotypes result from plant pigment metabolism [ 13 ]. Anthocyanin, flavonoid, and alkaloid are well-known plant pigments that play active roles in the regulation of flower color [ 11 ]. Flavonoids are involved in the production of red, pink, purple and blue in plants and are widely distributed in petals [ 27 ]. Additionally, and most importantly, flavonoid accumulation in a spontaneous cotton mutant resulted in red coloration and enhanced disease resistance [ 28 ]. Anthocyanins are a group of water-soluble pigments that confer the blue, purple and red color to many fruits [ 29 ]. Carotenoids primarily affect yellow, orange and red coloration [ 30 ]. To identify the potential mechanisms that regulate the modification of flowering colors, the current study annotated differentially expressed genes via KEGG pathway analysis and Gene Ontology. In the pathway analysis, several pigment metabolism pathways were highlighted, including plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), anthocyanin biosynthesis (ko00942), flavonoid biosynthesis (ko00941), and carotenoid biosynthesis (ko00906). Furthermore, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis were active throughout the flowering process and nearly up-regulated at each stage. These results corroborated the findings of several previous investigations into the linkages between petal color and plant pigment. Previous studies found that the flowers of ‘Yujin 2’ were red due to their higher carotenoid and flavonoid content in contrast to ‘Damaohua’ at all developmental stages [ 9 ]. However, the results of our study revealed that carotenoid biosynthesis was down-regulated throughout S1-6, which was different from those of previous studies to some extent, due to the formation of mutants that varied with sampling times. Phenylalanine is an important anthocyanin biosynthesis pathway [ 13 ]. In this study, phenylalanine synthesis was enhanced at the S1 phase, and anthocyanin synthesis was weakened at S5 and S6. This result may be consistent with the changes of ‘Yujin 2’ petal color from red (S1) to white (S5) and gold (S6). The results obtained in the current study revealed that the phenylpropanoid biosynthetic process (GO:0009699), carotenoid biosynthetic process (GO:0016117), and anthocyanin 5-O-glucoside-4'''-O-malonyltransferase activity (GO:0102801) were triggered. Anthocyanin is the stable configuration of anthocyanidins that is always combined with a sugar moiety [ 31 ], which may explain the active glucose-related metabolism pathways in our research. GO annotation also supported the KEGG evidence. Furthermore, GO annotation showed that plasma membrane (GO:0005886) was active throughout the flowering process. Chloroplast (GO:0009507) and vacuole (GO:0005773) were active through most of the flowering processes, and could be explained by the results of previous experiments. The colors of plants are primarily codetermined by the different distributions of chlorophyll, carotenoid, flavonoid, and betaine, which exist in the cytoplasms or vacuoles of plants [ 9 ]. Moreover, GO analysis revealed that environment-related responses (such as to cold) were triggered (GO:0009409). Generally, the binding of upstream transcription factors, initiation of gene transcription and expression, final synthesis of anthocyanins, and coloration of fruits and peels are closely related to the expression of structural genes in the anthocyanin biosynthetic pathway [ 19 ]. At present, MYB, BHLH and WD40 are the three major transcription factors currently known. Anthocyanin biological regulation includes single-structure gene regulation, two-gene interaction regulation, and some MBW (MYB-bHLH-WD40) complexes regulate anthocyanin accumulation [ 16 ]. In the first three stages of flower development, LjDFR, LjABCB1, LjMYC6, LjDDB2, and LjANS genes increased rapidly. However, only the expression of LjF3'5'H was significantly down-regulated at S5, which was consistent with the accumulation of anthocyanin. The expression of other transcription factors did not increase during anthocyanin accumulation, which may not be directly related to anthocyanin synthesis. To verify the RNA-Seq results, we randomly selected six genes and tested their expression patterns via qPCR, which revealed the same expression trend with the sequencing results. Conclusions In summary, this study developed a transcriptome profile of flower color generation for L. japonica as well as annotated unigene sets of seven anthesis phases. Bioinformatics analysis revealed a number of distinct clues in regard to color development mechanisms. Further, candidate genes and potential pathways were obtained, which provided a foundation for further research at the transcriptome level. However, the current results still did not fully elucidate the overall regulatory mechanisms of plant pigments and the formation of the red color for L. japonica . Therefore, future studies relating to the current topic, such as transcript regulation, protein expression, and biochemical metabolism, are recommended. Candidate key genes involved in anthocyanin metabolism should also be further investigated. Methods Plant material preparation Five year old ‘Damaohua’ and ‘Yujin 2’ were planted in the resource garden of the College of Life Sciences, at Henan Normal University (N35°18 '13.71 ", E113°55' 15.05). The land owner allows us to sample on their private land. Professor Jingyuan Li identifies ‘Damaohua’ and ‘Yujin 2’ as honeysuckle [ 9 ] (Fig. 6 ). Our research group collected fresh flower buds of of ‘Damaohua’ and ‘Yujin 2’ in seven periods: young bud stage (S1), three green stage (S2), two white stage (S3), great white stage (S4), silver stage (S5), golden stage (S6), fade stage (S7). During sample collection, the flowers were combined and regarded as one biological replicate representing each stage, and three independent replicates were prepared. Partial flower materials were immediately frozen in liquid nitrogen following collection and stored at -80 °C. RNA isolation and library construction The total RNA was extracted by mirVana miRNA Isolation Kit (Ambion) [ 32 ]. RNA integrity was evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), and samples with the RNA Integrity Number (RIN) ≥ 7 were subjected to subsequent analysis. Separate mRNA from total RNA with oligonucleotide beads (dT) and add fragment buffer to cut it into short fragments, and then use random hexamer primers as templates to synthesize first strand cDNA. The library was constructed according to the manufacturer's instructions for TruSeq chain mRNA LTSample Prep Kit (Illumina, San Diego, California, USA). Sequencing, de novo assembly and annotation The above libraries were sequenced using an Illumina HiSeq X Ten sequencer (Illumina Inc., USA) and 150 bp paired-end reads were generated. The preparation of the cDNA library and sequencing were performed at Shanghai OE Biotech. Co., Ltd., Shanghai, China, and raw data (raw reads) were processed using Trimmomatic [ 33 ]. The reads containing ploy-N and low quality reads were removed to obtain clean reads. Following the removal of adaptor and low quality sequences, the clean reads were assembled into expressed sequence tag clusters (contigs) and assembled de novo into transcripts using Trinity [ 34 ] (version: trinityrnaseq_r20131110) through a paired-end method. The longest transcript was selected as a unigene based on the similarity and length of a sequence for subsequent analysis. The function of the unigenes was annotated to the NCBI non-redundant protein (NR), Clusters of Orthologous Groups (COG) for eukaryotic complete genomes, Gene Ontology (GO), Swiss-Prot, and evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases by diamond software, and mapped to Pfam databases by HMMER. The search was conducted using Blastx with a threshold E-value cut-off of 1e-5. Unigene quantification, analysis of differentially expressed unigenes (DEGs), cluster analysis, GO and KEGG enrichment FPKM [ 35 ], and read count values of each unigene were calculated using bowtie2 [ 36 ]and eXpress[ 37 ]. We use the DESeq function to estimate the size factor and the nbinom test to determine DEG [ 38 ]. The threshold value of P 2 or multiple change < 0.5 was set as the significant difference threshold.To explore the transcriptional expression patterns of DEGs, cluster analysis was performed. We performed GO analysis and KEGG pathway enrichment analysis for DEGs based on hypergeometric distribution of R. Expression analysis of transcription factors (TFs) PlantTFDB ( http://planttfdb.cbi.pku.edu.cn/index.php ) contains sequences of 58 plant transcription factor families from 165 plants, and is a plant transcription factor database [ 39 ]. Blastx was used to compare the sequence of unigenes with transcription factor database, and the best gene with E value less than 1e-5 was selected as annotation information of unigene. Candidates that contained DNA binding domains were recognized by GO annotation for the final TF identification. Differentially expressed TFs (DETFs) between samples were identified using the value of Fragments Per Kilobase of transcript Per Million fragments mapped (FPKM) with |log2(fold change)| > 1, p value ≤ 0.05 and q value ≤ 0.05 (Bolger et al. 2014). Multiexperiment Viewer software (v4.9) was employed to exhibit the expression profiles of DETFs by clustering [ 34 ] and the DETFs were clustered via the K-means method. Using the total TFs in L. japonica as reference, GO enrichment analysis of DETFs was performed using Agrigo ( http:// bioinfo.cau.edu.cn/agriGO/) with the hypergeometric test and FDR cut-off of 0.05 [ 40 ]. qRT-PCR The same RNA samples used for the RNA-seq experiments were employed for qRT-PCR. We use NanoDrop 2000 spectrophotometer to measure RNA yield (Thermo Scientific, USA), whereas the integrity was evaluated using agarose gel electrophoresis stained with ethidium bromide. Quantification was performed using a two-step reaction process: reverse transcription (RT) and PCR. Each sample was run in triplicate for analysis. At the end of the PCR cycles, melting curve analysis was performed to validate the specific generation of the expected PCR product. Quantification was performed using the 2 − ΔΔCT method, and data were normalized to the ACT2/7 transcript [ 41 ]. The sequences of the primers used are listed in Additional file 1: Table S1. Statistical analysis Significant differences were calculated using a one-way ANOVA analysis with a Turkey test and a significance level at p ≤ 0.05 and p ≤ 0.01 in SPSS 19.0 software [ 9 ]. All expression analyses were repeated three times. The reported value represents the arithmetic mean of three replicates. Data was expressed as a mean plus or minus a standard deviation (mean ± SD). Abbreviations L. japonica Lonicera japonica Thunb DEGs differentially expressed genes CGA chlorogenic acid CHS chalcone synthase, F3H flavanone 3-hydroxylase, DFR dihydroflavonol-4-reductase, ANS anthocyanidin synthase, UFGT UDP -glucose flavonoid 3-O-glucosyltransferase Declarations Ethics approval and consent to participate Not applicable Consent for publication Not applicable Availability of data and material All data generated or analyzed for this study are included in this article and its supplementary information files. Competing interests The authors declare that they have no competing interests Funding The work was carried out under the auspices of the Public Health Services of Traditional Chinese Medicine in 2017 [(2017)66] and Henan Province Enterprise Technology Innovation Guidance Special Project (172107000031). Authors’ contributions LJ.J. designed the research methodology, C C.F. and L Y. performed data analysis and wrote the article, and Y C.L,MJ.X. and C T prepared the plant samples. All authors have read and approved the manuscript. Acknowledgment At the point of finishing this paper, I’d like to express my sincere thanks to all those who have lent me hands in the course of my writing this paper. References Shang X, Pan H, Li M, Miao X, Ding H. Lonicera japonica Thunb.: ethnopharmacology, phytochemistry and pharmacology of an important traditional Chinese medicine. J Ethnopharmacol. 2011;138(1):1–21. doi: 10.1016/j.jep.2011.08.016 . Hsu HF, Hsiao PC, Kuo TC, Chiang ST, Chen SL, Chiou SJ, Ling XH, Liang MT, Cheng WY, Houng JY. Antioxidant and anti-inflammatory activities of Lonicera japonica Thunb. var. sempervillosa Hayata flower bud extracts prepared by water, ethanol and supercritical fluid extraction techniques. Ind Crop Prod. 2016;89:543–9. doi: 10.1016/j.indcrop.2016.05.010 . .3 Wang DY, Zhao XM, Liu YL. (2017) Hypoglycemic and hypolipidemic effects of a polysaccharide from flower buds of Lonicera japonica in streptozotocin-induced diabetic rats. International journal of biological macromolecules 102:396–404. doi: 10.1016/j.ijbiomac.2017.04.056 . Tian J, Che HL, Ha D, Wei YP, Zheng SY. Characterization and anti-allergic effect of a polysaccharide from the flower buds of Lonicera japonica. Carbohyd Polym. 2012;90(4):1642–7. doi: 10.1016/j.carbpol.2012.07.044 . Ryu KH, Rhee HI, Kim JH, Yoo H, Lee BY, Um KA, Kim K, Noh JY, Lim KM, Chung JH. Anti-Inflammatory and Analgesic Activities of SKLJI, a Highly Purified and Injectable Herbal Extract of Lonicera japonica. Biosci Biotech Bioch. 2010;74(10):2022–8. doi: 10.1271/bbb.100279 . Han J, Lv QY, Jin SY, Zhang TT, Jin SX, Li XY, Yuan HL. Comparison of anti-bacterial activity of three types of di-O-caffeoylquinic acids in Lonicera japonica flowers based on microcalorimetry. Chin J Nat Medicines. 2014;12(2):108–13. doi: 10.1016/S1875-5364(14)60017-0 . Ning W, Peng X, Ma LY, Cui L, Lu XP, Wang J, Tian JK, Li XM, Wang W, Zhang L. Enhanced secondary metabolites production and antioxidant activity in postharvest Lonicera japonica Thunb. in response to UV radiation. Innov Food Sci Emerg. 2012;13:231–43. doi: 10.1016/j.ifset.2011.10.005 . Wang Z, Clifford MN, Sharp P. Analysis of chlorogenic acids in beverages prepared from Chinese health foods and investigation, in vitro, of effects on glucose absorption in cultured Caco-2 cells. Food Chem. 2008;108(1):369–73. doi: 10.1016/j.foodchem.2007.10.083 . Li JJ, Lian XY, Ye CL, Wang L. Analysis of Flower Color Variations at Different Developmental Stages in Two Honeysuckle ( Lonicera Japonica Thunb.) Cultivars. Hortscience. 2019;54(5):779–82. doi: 10.21273/Hortsci13819-18 . Qu Y, Ou Z, Yang FS, Wang S, Peng J. The study of transcriptome sequencing for flower coloration in different anthesis stages of alpine ornamental herb (Meconopsis 'Lingholm'). Gene. 2019;689:220–6. doi: 10.1016/j.gene.2018.12.017 . Zhao D, Tao J. Recent advances on the development and regulation of flower color in ornamental plants. Front Plant Sci. 2015;6:261. doi: 10.3389/fpls.2015.00261 . Yoshida K, Mori M, Kondo T. Blue flower color development by anthocyanins: from chemical structure to cell physiology. Nat Prod Rep. 2009;26(7):884–915. doi: 10.1039/b800165k . Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54(4):733–49. doi: 10.1111/j.1365-313X.2008.03447.x . Tanaka Y, Brugliera F, Kalc G, Senior M, Dyson B, Nakamura N, Katsumoto Y, Chandler S. Flower Color Modification by Engineering of the Flavonoid Biosynthetic Pathway: Practical Perspectives. Biosci Biotech Bioch. 2010;74(9):1760–9. doi: 10.1271/bbb.100358 . Sadilova E, Stintzing FC, Carle R. Anthocyanins, colour and antioxidant properties of eggplant (Solanum melongena L.) and violet pepper (Capsicum annuum L.) peel extracts. Z Naturforsch C. 2006;61(7–8):527–35. Yang Y, Cui BH, Tan ZW, Song BX, Cao HN, Zong CW. (2018) RNA sequencing and anthocyanin synthesis-related genes expression analyses in white-fruited Vaccinium uliginosum. Bmc Genomics 19. doi:ARTN 93010.1186/s12864-018-5351-0. Gao LX, Yang HX, Liu HF, Yang J, Hu YH. (2016) Extensive Transcriptome Changes Underlying the Flower Color Intensity Variation in Paeonia ostii. Frontiers in Plant Science 6. doi:ARTN 120510.3389/fpls.2015.01205. Tripathi AM, Niranjan A, Roya S. Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna. Plant Physiol Bioch. 2018;127:1–10. doi: 10.1016/j.plaphy.2018.03.008 . Zhao L, Gao LP, Wang HX, Chen XT, Wang YS, Yang H, Wei CL, Wan XC, Xia T. The R2R3-MYB, bHLH, WD40, and related transcription factors in flavonoid biosynthesis. Funct Integr Genomic. 2013;13(1):75–98. doi: 10.1007/s10142-012-0301-4 . Bi H, Luang S, Li Y, Bazanova N, Morran S, Song Z, Perera MA, Hrmova M, Borisjuk N, Lopato S. Identification and characterization of wheat drought-responsive MYB transcription factors involved in the regulation of cuticle biosynthesis. J Exp Bot. 2016;67(18):5363–80. doi: 10.1093/jxb/erw298 . Zheng J, Wu H, Zhu HB, Huang CY, Liu C, Chang YS, Kong ZC, Zhou ZH, Wang GW, Lin YJ, Chen H. Determining factors, regulation system, and domestication of anthocyanin biosynthesis in rice leaves. New Phytol. 2019;223(2):705–21. doi: 10.1111/nph.15807 . Wiczkowski W, Szawara-Nowak D, Topolska J. Changes in the content and composition of anthocyanins in red cabbage and its antioxidant capacity during fermentation, storage and stewing. Food Chem. 2015;167:115–23. doi: 10.1016/j.foodchem.2014.06.087 . Tasaki K, Higuchi A, Watanabe A, Sasaki N, Nishihara M. Effects of knocking out three anthocyanin modification genes on the blue pigmentation of gentian flowers. Sci Rep. 2019;9(1):15831. doi: 10.1038/s41598-019-51808-3 . Tanase K, Nishitani C, Hirakawa H, Isobe S, Tabata S, Ohmiya A, Onozaki T. Transcriptome analysis of carnation (Dianthus caryophyllus L.) based on next-generation sequencing technology. BMC Genom. 2012;13:292. doi: 10.1186/1471-2164-13-292 . Guo L, Wang Y, da Silva JAT, Fan Y, Yu X. Transcriptome and chemical analysis reveal putative genes involved in flower color change in Paeonia 'Coral Sunset'. Plant Physiol Biochem. 2019;138:130–9. doi: 10.1016/j.plaphy.2019.02.025 . Zhou C, Mei X, Rothenberg DO, Yang Z, Zhang W, Wan S, Yang H, Zhang L. (2020) Metabolome and Transcriptome Analysis Reveals Putative Genes Involved in Anthocyanin Accumulation and Coloration in White and Pink Tea (Camellia sinensis) Flower. Molecules 25 (1). doi: 10.3390/molecules25010190 . Williams CA, Grayer RJ. Anthocyanins and other flavonoids. Nat Prod Rep. 2004;21(4):539–73. doi: 10.1039/b311404j . Long L, Liu J, Gao Y, Xu FC, Zhao JR, Li B, Gao W. Flavonoid accumulation in spontaneous cotton mutant results in red coloration and enhanced disease resistance. Plant Physiol Biochem. 2019;143:40–9. doi: 10.1016/j.plaphy.2019.08.021 . Fang J. Classification of fruits based on anthocyanin types and relevance to their health effects. Nutrition. 2015;31(11–12):1301–6. doi: 10.1016/j.nut.2015.04.015 . Rosas-Saavedra C, Stange C. Biosynthesis of Carotenoids in Plants: Enzymes and Color. Subcell Biochem. 2016;79:35–69. doi: 10.1007/978-3-319-39126-7_2 . Veitch NC, Grayer RJ. Flavonoids and their glycosides, including anthocyanins. Nat Prod Rep. 2011;28(10):1626–95. doi: 10.1039/c1np00044f . Yang B, Guan Q, Tian J, Komatsu S. Transcriptomic and proteomic analyses of leaves from Clematis terniflora DC. under high level of ultraviolet-B irradiation followed by dark treatment. Journal of proteomics. 2017;150:323–40. doi: 10.1016/j.jprot.2016.10.001 . Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. doi: 10.1093/bioinformatics/btu170 . Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, Chen ZH, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644-U130. doi: 10.1038/nbt.1883 . Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5. doi: 10.1038/nbt.1621 . Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357-U354. doi: 10.1038/Nmeth.1923 . Roberts A, Pachter L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods. 2013;10(1):71–3. doi: 10.1038/nmeth.2251 . Gupta R, Dewan I, Bharti R, Bhattacharya A. (2012) Differential Expression Analysis for RNA-Seq Data. ISRN Bioinform 2012:817508. doi: 10.5402/2012/817508 . Jin J, Tian F, Yang DC, Meng YQ, Kong L, Luo J, Gao G. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017;45(D1):D1040–5. doi: 10.1093/nar/gkw982 . 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. doi: 10.1038/nmeth.1226 . Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. doi: 10.1006/meth.2001.1262 . Supplementary Files Additionalfile.xlsx Additional file 1: Table S1. The preprocessing results of sequencing data quality. Additional file 2: Table S2. Specie taxonomy based on Nucleotide (NR) database. Additional file 3: Table S3. Enriched KEGG pathways of the L.japonica flower transcriptome. Additional file 4: Table S4. Enriched GO of the L.japonica flower transcriptome. Additional file 5: Table S5. Differentially expressed TFs involved in anthocyanin metabolism. Additional file 6: Table S6. The primers used for analysis of gene expression by qRT-PCR. Additional file 7: Fig. 1 Summary of top 20 KEGG pathway assignments for the L. japonica flower transcriptome A-G: The top 20 KEGG pathways in S1, S2, S3, S4, S5, S6, and S7. Additional file 8:Fig. 2 Summary of top 10 GO terms for the L. japonica flower transcriptome. A-G: The top 10 GO terms in S1, S2, S3, S4, S5, S6, and S7. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-104094","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research","associatedPublications":[],"authors":[{"id":4459891,"identity":"5262be41-d5da-4652-87d5-582fb097db04","order_by":0,"name":"Jianjun Li","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAoklEQVRIiWNgGAWjYBACNobjBw58MLBJIF4LP+OZxIczCtISeIjWItl8wNiY58NhErQYHDuQJjnD4HyevUQC4+cCorScOXhM4oPB7WIeiQRm6RlEabkBtuV2Yo9EAhszUa4zuP/ATJrH4BwJWiQbQN43OECCFn4GUCAbJCf2nHnYLE2UFkhU/rFLbG9PPviZ+JCGAMYGEjWMglEwCkbBKMAJAGA2NfrM4qEWAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-6034-2411","institution":"Henan Normal University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Jianjun","middleName":"","lastName":"Li","suffix":""},{"id":4459892,"identity":"7413b6ea-29b5-4b79-b408-03c497291266","order_by":1,"name":"Chenglin Ye","email":"","orcid":"","institution":"Henan Normal University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Chenglin","middleName":"","lastName":"Ye","suffix":""},{"id":4459893,"identity":"6a08d145-8d9f-4451-a724-1d6d21351089","order_by":2,"name":"Jingxiao Ma","email":"","orcid":"","institution":"Henan Normal University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jingxiao","middleName":"","lastName":"Ma","suffix":""},{"id":4459894,"identity":"26ee0aee-140a-4b27-99b0-da6aded148db","order_by":3,"name":"Ting Cheng","email":"","orcid":"","institution":"Henan Normal University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Ting","middleName":"","lastName":"Cheng","suffix":""},{"id":4459895,"identity":"bb3ea9d2-cdc9-434a-8ad2-d3d6b8952de4","order_by":4,"name":"Yan Lv","email":"","orcid":"","institution":"Henan Normal University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yan","middleName":"","lastName":"Lv","suffix":""},{"id":4459896,"identity":"8bcc3d97-8436-49c5-acc0-9884310f87e6","order_by":5,"name":"Cuifang Chang","email":"","orcid":"","institution":"Henan Normal University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Cuifang","middleName":"","lastName":"Chang","suffix":""}],"badges":[],"createdAt":"2020-11-06 14:59:04","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-104094/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-104094/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":3541005,"identity":"f532dd45-32ba-4ed6-88c8-34c61516396c","added_by":"auto","created_at":"2020-11-12 17:04:15","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":87722,"visible":true,"origin":"","legend":"Species taxonomy based on the Nucleotide (NR) database.","description":"","filename":"1.PNG","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/a8551fc15fa61fb218615461.PNG"},{"id":3541007,"identity":"4362c4f8-bfff-4315-aa00-5fbc823b971d","added_by":"auto","created_at":"2020-11-12 17:04:15","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":18643,"visible":true,"origin":"","legend":"Number of up-regulated and down-regulated genes across seven flowering stages.","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/93e0f3a6c08938a865cc6603.png"},{"id":3541008,"identity":"a266639e-bae3-4dbb-aa9e-e436c5fd717c","added_by":"auto","created_at":"2020-11-12 17:04:15","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":110906,"visible":true,"origin":"","legend":"Enriched KEGG pathway that existed in seven periods；B: Enriched GO pathway that existed in seven periods.","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/11bcd7da9687c9a540b63433.png"},{"id":3541009,"identity":"9c2fda93-8d95-46c4-8899-89192ada927c","added_by":"auto","created_at":"2020-11-12 17:04:16","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":5641641,"visible":true,"origin":"","legend":"Identification of potential unigenes associated with anthocyanin metabolism, and their expression levels at different stages of L. japonica. A: The numbers of upregulated or downregulated across the seven flowering stages; B: DEGs related to TFs in anthocyanin development; C: Transcript expression analysis for unigenes associated with anthocyanin development. Changes in expression level are represented by color; green indicates a lower expression level, whereas red indicates a higher expression level.","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/8d9ff28a4a5fc3e168dfe67a.png"},{"id":3541010,"identity":"2af57db6-92cd-4ec6-a1f0-15b5240f673e","added_by":"auto","created_at":"2020-11-12 17:04:17","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":87299,"visible":true,"origin":"","legend":"Real-time qPCR validation of six unigenes involved in anthocyanin metabolism. The orange line indicates qPCR results, blue line indicates RNA-Seq results.","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/67c1d4e730cfd69729d09ef8.png"},{"id":3541011,"identity":"37452af6-0ab2-491b-8d2a-addc645dfe24","added_by":"auto","created_at":"2020-11-12 17:04:17","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":393955,"visible":true,"origin":"","legend":"The morphology of flowers at different developmental stages in two Lonicera Japonica Thunb.\n\n\n","description":"","filename":"6.PNG","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/46e7eeeb0034b79b2f464ab1.PNG"},{"id":3541012,"identity":"123f0f1d-b13a-4473-8d89-13868dd943bc","added_by":"auto","created_at":"2020-11-12 17:04:17","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":5641641,"visible":true,"origin":"","legend":"","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/34864a86469dfdddeed0dd13.png"},{"id":13613478,"identity":"22b45721-03f0-47e1-9531-b5e7400b61c6","added_by":"auto","created_at":"2021-09-17 06:37:13","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2633023,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/8c33f48e-4c91-4852-b51b-36ef90c3f9ea.pdf"},{"id":3541006,"identity":"fe764ce3-299f-4938-aa4f-dceff71a514c","added_by":"auto","created_at":"2020-11-12 17:04:15","extension":"xlsx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":4636836,"visible":true,"origin":"","legend":"Additional file 1: Table S1. The preprocessing results of sequencing data quality.\nAdditional file 2: Table S2. Specie taxonomy based on Nucleotide (NR) database.\nAdditional file 3: Table S3. Enriched KEGG pathways of the L.japonica flower transcriptome.\nAdditional file 4: Table S4. Enriched GO of the L.japonica flower transcriptome.\nAdditional file 5: Table S5. Differentially expressed TFs involved in anthocyanin metabolism.\nAdditional file 6: Table S6. The primers used for analysis of gene expression by qRT-PCR.\nAdditional file 7: Fig. 1 Summary of top 20 KEGG pathway assignments for the L. japonica flower transcriptome \nA-G: The top 20 KEGG pathways in S1, S2, S3, S4, S5, S6, and S7. \nAdditional file 8:Fig. 2 Summary of top 10 GO terms for the L. japonica flower transcriptome. \nA-G: The top 10 GO terms in S1, S2, S3, S4, S5, S6, and S7.\n","description":"","filename":"Additionalfile.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-104094/v1/8d27b107d5345345ceb8c484.xlsx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eTime-Series Transcriptome Sheds Light on Anthocyanin Metabolism and Species Diversification Between two Lonicera Japonica Thunb Cultivars\u003c/p\u003e","fulltext":[{"header":"Background","content":" \u003cp\u003e \u003cem\u003eLonicera japonica\u003c/em\u003e Thunb (\u003cem\u003eL. japonica\u003c/em\u003e) is a perennial semi-evergreen twining stoloniferous shrub. Traditionally, the \u003cem\u003eL. japonica\u003c/em\u003e flower bud, which has been listed in the Chinese Pharmacopoeia as Lonicera japonica (Jinyinhua in Chinese) is considered to be the primary medicinal component [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. Recent studies have found that \u003cem\u003eL. japonica\u003c/em\u003e has many medicinal properties due to its antioxidant [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e], hypoglycemic and hypolipidemic [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e], anti-allergic [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e], anti-inflammatory [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e], and antibacterial effects [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Further, the demand for \u003cem\u003eL. japonica\u003c/em\u003e in pharmaceuticals [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e], cosmetics [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e], and health foods [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e] has increased. Although \u003cem\u003eL. japonica\u003c/em\u003e is widely distributed in China, its yields and quality are the highest in Fengqiu of Henan and Shandong Provinces. Therefore, our research group selected and bred a new variety (Yujin 2). It has the characteristics of red buds, strong resistance, and a high content of chlorogenic acid (CGA) and luteoloside. The red buds endow \u0026lsquo;Yujin 2\u0026rsquo; with enhanced beauty, which is more valuable in utility [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eAs one of the most important phenotypic traits of ornamental plants, the formation of flower color is the result of the combined effect of environmental influences and genetic characteristics [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. The modification of petal color regulation genes is indispensable in modern flower breeding [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Through various studies it has been established that changes in flower color are primarily due to alterations in plant pigment metabolism [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. The key pigments to adjust the color are carotenoids, flavonoids, and alkaloids [11; 13; 14].\u003c/p\u003e \u003cp\u003eAnthocyanins, which are produced through flavonoid biosynthesis pathways during plant metabolism, are important water-soluble natural food pigment, which is widely present in plant epidermal cell vacuoles [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. In the first stage of anthocyanin biosynthesis, phenylalanine transforms into 4-coumaryl:CoA, referred to as the phenylpropanoid metabolic pathway. The next stage, 4-coumaryl:CoA transforms into multiple flavonoid compounds [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe biosynthesis of anthocyanins is related to a variety of structures and regulatory genes [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. The expression of chalcone synthase (CHS), Flavanone 3-hydroxylase (F3H), dihydroflavonol-4-reductase (DFR), anthocyanin synthase (ANS) and UDP glucoflavonoid 3-o-glucosyltransferase (UFGT) directly affect the accumulation of anthocyanins in plants [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. These structural genes have been cloned in different plants, and the expression pathway is common to different types of plants.\u003c/p\u003e \u003cp\u003eKey regulatory genes serve as regulators that combine individually, or in the form of dimers and trimers, which play a critical role in the regulation of anthocyanin synthesis. So far, MYB, bHLH and WD40 transcription factors are widely regarded as the three main regulatory factors [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. The MYB family transcription factor is one of the largest transcription factor families in plants, which has a MYB domain that is composed of 1\u0026thinsp;~\u0026thinsp;4 MYB repeating units in its structure [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. In conjunction with MYB, the bHLH transcription factors play active and critical roles in anthocyanin biosynthesis. The high anthocyanin content of many crops results from the increased expression of bHLH transcription factors [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. WD40 is an ancient and highly conserved family of proteins that is involved in myriad plant cellular processes, such as transcriptional regulation, cell division, vesicle formation and transport, processing, signal transduction and so on [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThere were great differences in the corolla color of \u0026lsquo;Yujin 2\u0026rsquo; at the various developmental stages. The upper half of the flower bud during the early stage was purple-red, whereas the whole flower bud from the three green periods to the great white period was red, but gradually became lighter; whereas during the silver period, the exterior was red and the interior was white. During the golden period, the exterior of \u0026lsquo;Yujin 2\u0026rsquo; was red and the interior was white, but for \u0026lsquo;Damaohua\u0026rsquo; it was yellow [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eColor is considered to be an important indicator in the assessment of ornamental value and offers a basis for plant classification [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. However, the molecular mechanisms of anthocyanin metabolism in \u003cem\u003eL. japonica\u003c/em\u003e flowers have not yet been discovered. Therefore, this paper employed \u0026lsquo;Yujin 2 and \u0026lsquo;Damaohua\u0026rsquo; flowers at different stages of development to study the molecular mechanisms of anthocyanin formation. This, by means of morphology, transcriptome, and bioinformatics, etc., to provide possibilities for enhancing the \u003cem\u003eL. japonica\u003c/em\u003e germplasm using genetic engineering technologies, and cultivating new varieties with different colors.\u003c/p\u003e "},{"header":"Results","content":" \u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eCharacteristics of \u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo; transcriptome\u003c/h2\u003e \u003cp\u003eA total of 42 samples (\u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo;, respectively, have seven stages, each with three independent replicates) were sequenced using the Illumina HiSeq X Ten sequencer (Illumina Inc., USA), from which 311\u0026nbsp;Gb of raw data was harvested. Using a pair-end sequencing strategy, 2G raw reads with mean lengths of 150 nt were obtained from the 42 samples with 47.62\u0026nbsp;M reads for each sample. After the filtering of raw reads, 1.95 G reads were counted, with an average of 46.43\u0026nbsp;M reads per sample (Additional file 1: Table S1). Following the removal of redundant sequences and quantity control, 133,487 unigenes were assembled. The average length of the unigenes was 1,047\u0026nbsp;bp with 45.75% GC and 1642nt N50. The longest transcript was 16,864\u0026nbsp;bp, and the shortest was 301\u0026nbsp;bp.\u003c/p\u003e \u003cp\u003eThe unigenes were annotated to obtain functional information using seven separate databases, including NCBI non-redundant protein (NR), Clusters of Orthologous Groups (COG) for eukaryotic complete genomes, Gene Ontology (GO), Swiss-Prot, evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases by Diamond Software, and mapped to Pfam databases by HMMER. Overall, 66.79% (133,487 unigenes) were successfully annotated, among which most unigenes (66.46%) were assigned to annotation terms in the NR database (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Further, 61.43% of the unigenes were annotated on eggNOG, 46.62% on Swiss-prot, 43.46% on GO, 35.76% on KOG, and 26.02% on the KEGG pathway (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003ePercentage of annotated unigenes.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"2\"\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eDatabase\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCount of annotated unigenes%\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e66.46%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKOG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e35.76%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGO\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e43.46%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSwiss-prot\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e46.62%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eeggNOG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e61.43%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKEGG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e26.02%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePfam\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.13%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e66.79%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eAccording to the NR annotation results, the current study compared \u003cem\u003eL. japonica\u003c/em\u003e with taxonomic species stored in the NR database. A total of 97247 sequences were matched with \u003cem\u003eL. japonica\u003c/em\u003e unigenes, which covered 1512 taxonomy categories (Additional file 2: Table S2). The highest ranked species was \u003cem\u003eAlternaria alternate\u003c/em\u003e, which took up 12.82%, whereas the second highest-ranked was \u003cem\u003eVitis vinifera\u003c/em\u003e, which contributed 8.80% (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eDifferentially expressed genes (\u003c/b\u003eDEGs\u003cb\u003e)\u003c/b\u003e\u003c/p\u003e \u003cp\u003eWhen comparing seven flowering stages of expressed genes between \u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo;, 73,088 differentially expressed genes exhibited remarkable changes. A total of 5,355 genes were differentially expressed, with 2,864 genes up-regulated and 2491 genes down-regulated at the S1 stage. At the S2 stage, 21,159 genes were differentially expressed, with 12,734 genes up-regulated and 8,425 genes down-regulated. At the S3 stage, 19,126 genes were differentially expressed, with 10,343 genes up-regulated and 8,783 genes down-regulated. At the S4 stage, 11,153 genes were differentially expressed, with 6,029 genes up-regulated and 5,124 genes down-regulated. At the S5 stage, 27,718 genes were differentially expressed, with 12,102 genes up-regulated and 15,616 genes down-regulated. At the S6 stage, 18,652 genes were differentially expressed, with 8,486 genes up-regulated and 10,166 genes down-regulated. Finally, at the S7 stage, 38,162 genes were differentially expressed, with 21,153 genes up-regulated and 17,009 genes down-regulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). As shown in the figures, the number of up-regulated DEGs was higher than those that were down-regulated.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eKEGG-pathway analysis of differentially expressed genes\u003c/h2\u003e \u003cp\u003eUnigenes were assigned to KEGG pathways based on annotation data. DEGs were enriched into 181, 202, 202, 191, 207, 204, and 210 pathways (S1, S2, S3, S4, S5, S6, and S7) (Additional file 3: Table S3). At the S1 stage, 53 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), and phenylpropanoid biosynthesis (ko00940) during the entire process, as well as up-regulation and down-regulation processes, played active roles, and the expression was significant. Amino sugar and nucleotide sugar metabolism (ko00520) were also active during the entire process and up-regulation, while they were not active in the down-regulation process. Starch and sucrose metabolism (ko00500) were active during the entire process and down-regulation process, while they were not active in the up-regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA).\u003c/p\u003e \u003cp\u003eAt the S2 stage, 72 pathways were enriched with a statistically significant enrichment p-value, which was more than the other six paired comparison phases. Endocytosis (ko04144), plant hormone signal transduction (ko04075), and starch and sucrose metabolism (ko00500) during the entire process, as well as up- and down-regulation processes played an active role, and the expression was significant. Phenylpropanoid biosynthesis (ko00940) was active during the entire process and the down-regulation process, while it decreased in the up-regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB).\u003c/p\u003e \u003cp\u003eAt the S3 stage, 66 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500) during the entire process, and the up- and down-regulation processes played an active role, and the expression is significant. Glycolysis/Gluconeogenesis (ko00010) was active during the entire process and the down-regulation process, while it was not active in the up-regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003eAt the S4 stage, 71 pathways were enriched with a statistically significant enrichment p-value, and plant hormone signal transduction (ko04075) during the entire process, the up- and down-regulation processes played an active role, and the expression was significant. Phenylpropanoid biosynthesis (ko00940) and cell cycle (ko04110) were active during the entire process and the up-regulation process, while they were not active in the down-regulation process. Starch and sucrose metabolism (ko00500) were active during the entire process and the down-regulation process, while they were not active in the up-regulation process. Cell cycle-yeast (ko04111) was triggered in the up-regulation process, while carbon metabolism (ko01200) was triggered in the down-regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD).\u003c/p\u003e \u003cp\u003eAt the S5 stage, 63 pathways were enriched with a statistically significant enrichment p-value. Plant hormone signal transduction (ko04075) and phenylpropanoid biosynthesis (ko00940) played active role during the entire process as well as up- and down-regulation processes, and the expression was significant. Glycerophospholipid metabolism (ko00564), starch and sucrose metabolism (ko00500) were active in the total process and the up-regulation process, which were not active in the down-regulation process. Glutathione metabolism (ko00480), drug metabolism - cytochrome P450 (ko00982), and metabolism of xenobiotics by cytochrome P450 (ko00980) were triggered in the down -regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eE).\u003c/p\u003e \u003cp\u003eAt the S6 stage, 31 pathways were enriched a with statistically significant enrichment p-value. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500), glycolysis/gluconeogenesis (ko00010) played an active role during the entire process, the up- and down-regulation process, and the expression was significant. Valine, leucine and isoleucine degradation (ko00280) was triggered in the up-regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eF).\u003c/p\u003e \u003cp\u003eAt the S7 stage, 39 pathways were enriched with a statistically significant enrichment p-value. Fatty acid metabolism (ko01212) was active during the entire process, but was not in the up- and down-regulation process. Plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500) were active during the entire process and up-regulation process, while they were not in the down-regulation process. Tyrosine metabolism (ko00350) was active during the entire process and the down-regulation process, while it was not in the up-regulation process. Spliceosome (ko003040) was triggered in the up-regulation process, whereas MAPK signaling pathway-yeast (ko04011), as well as glycine, serine and threonine metabolism (ko00260) were triggered in the down -regulation process (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eG).\u003c/p\u003e \u003cp\u003eSummarizing the above KEGG pathways enriched during the flowering period, it was found that phenylpropanoid biosynthesis and flavonoid biosynthesis were active through the entire flowering process and nearly up-regulated at each stage, Peroxisome(ko04146), Carbon metabolism༈ko01200༉ and Amino sugar and nucleotide sugar metabolism༈ko00520༉can be found in each period(Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA) ,and carotenoid biosynthesis was down-regulated from S1-6. The anthocyanin precursor, phenylalanine, was enhanced at the S1 phase; however, the synthesis of anthocyanin was weakened at the S5 and S6 phases (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eEnriched KEGG pathways by up-regulated and down-regulated genes.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"8\"\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKEGG Pathways\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eS1\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eS2\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS3\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS4\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eS5\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eS6\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003eS7\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePhenylpropanoid biosynthesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eCarotenoid biosynthesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFlavonoid biosynthesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u0026uarr;\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePhenylalanine metabolism\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePhenylalanine, tyrosine and tryptophan biosynthesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026uarr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAnthocyanin biosynthesis\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u0026darr;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"8\"\u003e*\u0026uarr;represents KEGG pathways enriched by up-regulated genes;\u0026darr;represents KEGG pathways enriched by down-regulated genes.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"8\"\u003e\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eGene Ontology (GO) analysis of differentially expressed genes\u003c/h2\u003e \u003cp\u003eDifferentially expressed genes were annotated by Gene Ontology terms (Additional file 4: Table S4). Genes were classified into three GO categories, that is, Biological Process, Cellular Component, and Molecular Function. From the Biological Process, auxin-activated signaling pathway (GO:0009734), response to abscisic acid (GO:0009737), response to wounding (GO:0009611) were remarkable at the S1 stage. The negative regulation of the hydrogen peroxide metabolic process (GO:0010727), positive regulation of the receptor biosynthetic process (GO:0010870), cellular response to topologically incorrect protein (GO:0035967) were remarkable at the S2 stage. The regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836), nucleolar chromatin organization (GO:1990700), and regulation of superoxide anion generation (GO:0032928) were remarkable at the S3 stage.\u003c/p\u003e \u003cp\u003eThe priming of the cellular response to stress (GO:0080136), cellular response to sucrose starvation (GO:0043617), and the cellulose metabolic process (GO:0030243) were remarkable at the S4 stage. The negative regulation of the gibberellin biosynthetic process (GO:0010373), guanosine tetraphosphate metabolic process (GO:0015969), and the regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836) were remarkable at the S5 stage. The regulation of the transcription of nuclear large rRNA transcript from RNA polymerase I promoter (GO:1901836), nucleolar chromatin organization (GO:1990700), and the establishment of protein localization to endoplasmic reticulum membrane (GO:0097051) were remarkable at the S6 stage. The negative regulation of the extrinsic apoptotic signaling pathway in absence of ligand (GO:2001240), the phosphatidylinositol-3- phosphate biosynthetic process (GO:0036092), and the regulation of the embryo sac egg cell differentiation (GO:0045694) were remarkable at the S7 stage (Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFrom the Cellular Component process, apoplast (GO:0048046), plant-type cell wall (GO:0009505), and cell wall (GO:0005618) were remarkable at the S1 stage. Inclusion body (GO:0016234), nitrite reductase complex [NAD(P)H] (GO:0009344), and plasma membrane (GO:0005886) were remarkable at the S2 stage. Nucleolar chromatin (GO:0030874) vacuole (GO:0005773), and chloroplast (GO:0009507) were remarkable at the S3 stage. Anchored component of plasma membrane (GO:0046658), apoplast (GO:0048046), and plasma membrane (GO:0005886) were remarkable at the S4 stage. Nucleolar chromatin (GO:0030874), P4 peroxisome (GO:0019822), chloroplast (GO:0009507) were remarkable at the S5 stage. Nucleolar chromatin (GO:0030874), chloroplast (GO:0009507), vacuole (GO:0005773) were remarkable at the S6 stage. The DNA topoisomerase complex (ATP-hydrolyzing) (GO:0009330), nuclear ubiquitin ligase complex (GO:0000152), and TRAPPII protein complex (GO:1990071) were remarkable at the S7 stage (Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFrom the Molecular Function process, DNA binding transcription factor activity (GO:0003700), aspartic-type endopeptidase activity (GO:0004190), hydrolase activity, and acting on ester bonds (GO:0016788) were remarkable at the S1 stage. Tau protein binding (GO:0048156), euchromatin binding (GO:1990188), and IAA-Phe conjugate hydrolase activity (GO:0010210) were remarkable at the S2 stage. IAA-Phe conjugate hydrolase activity (GO:0010210), IAA-Leu conjugate hydrolase activity (GO:0010211), and dehydroquinate synthase activity (GO:0102042) were remarkable at the S3 stage. Glucuronoxylan 4-O-methyltransferase activity (GO:0030775), DNA binding transcription factor activity (GO:0003700), sequence-specific DNA binding (GO:0043565) were remarkable at the S4 stage. Euchromatin binding (GO:1990188), malonyl-CoA decarboxylase activity (GO:0050080), and glutamine-phenylpyruvate transaminase activity (GO:0047316) were remarkable in the S5 stage; euchromatin binding (GO:1990188), IAA-Phe conjugate hydrolase activity (GO:0010210), IAA-Leu conjugate hydrolase activity (GO:0010211) were remarkable at the S6 stage. Arylacetonitrilase activity (GO:0047428), euchromatin binding (GO:1990188), and calcium-independent phospholipase A2 activity (GO:0047499) were remarkable at the S7 stage (Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFrom seven periods,defense response to bacterium (GO:0042742)༌transcription, DNA-templated (GO:0006351)༌regulation of transcription, DNA-templated (GO:0006355)༌vacuole (GO:0005773)༌vacuolar membrane (GO:0005774)༌extracellular region (GO:0005576)༌chloroplast (GO:0009507)༌nucleus (GO:0005634)༌integral component of membrane (GO:0016021)༌plasma membrane (GO:0005886)༌cytoplasm (GO:0005737)༌sequence-specific DNA binding (GO:0043565)༌DNA binding transcription factor activity (GO:0003700)༌metal ion binding (GO:0046872)༌ATP binding (GO:0005524) were found at the every period(Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB)\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eAnthocyanin metabolism analysis of differentially expressed genes\u003c/h2\u003e \u003cp\u003eSequencing results revealed that 114 unigenes were associated with anthocyanin metabolism, and 72 were significantly upregulated or downregulated across the seven flowering stages. The number of up-regulated unigenes was highest at the S7 stage, while the number of down-regulated unigenes was highest at the S5 stage (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). According to the analysis of TFs in anthocyanin metabolism, we obtained 47 genes relating to transcription factors in anthocyanin metabolism that belonged to 18 families, of which bHLH had the largest number, followed by GeBP (Additional file 5: Table S5, Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). According to the anthocyanin heatmap analysis, the expression of genes at S7 was different from that at other stages, with the highest up-regulated expression. The expression of genes at S6 was similar to that at S5, while those at S4, S3, S2, and S1 were similar, among which the expression of genes in S1 and S2 were the most similar (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eqRT-PCR verification\u003c/h2\u003e \u003cp\u003eTo validate the accuracy and repeatability of our RNA-Seq data, the current study selected six genes for quantitative real-time polymerase chain reaction (qRT-PCR) analysis, with gene-specific primers designed using Primer software (version 5.0), as shown in Additional file 6: Table S6. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, the results of qRT-PCR indicated that most of these genes had expression patterns that correlated with the RNA-Seq data, which confirmed the reliability of our data. The genes LjDFR (TRINITY_DN37612_c0_g1_i2_29), LjABCB1 (TRINITY_DN23592_c0_g2_i2_8), LjMYC6 (TRINITY_DN17985_c0_g1_i2_8), LjDDB2 (TRINITY_DN25019_c0_g2_i6_8), and LjANS (TRINITY_DN27734_c0_g1_i1_12) rapidly increased during the first three stages. However, the expression of only LjF3'5'H (TRINITY_DN27086_c0_g1_i1_29) was significantly down-regulated at S5, which was consistent with anthocyanin accumulation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e "},{"header":"Discussion","content":" \u003cp\u003eThe biological regulation of flower color had garnered much attention from researchers, particularly in terms of ornamental plant research [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. This present study focused on a new variety of \u003cem\u003eL. japonica\u003c/em\u003e (Yujin no. 2). There are a variety of designs and colors for \u0026lsquo;Yujin 2\u0026rsquo; that has a high ornamental value, which may play a significant role in landscaping [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. It is very important to regulate anthocyanin synthesis in \u003cem\u003eL. japonica\u003c/em\u003e breeding [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. With respect to molecular breeding, the identification of the key genes involved in anthocyanin synthesis will enable the overexpression or knockout of these genes to facilitate color modification [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. Therefore, this study began with the aim of identifying the essential genes that regulated the alteration of flower color and subsequently exploring the potential molecular mechanisms behind the development of flower color.\u003c/p\u003e \u003cp\u003eHowever, little evidence has been reported as relates to \u003cem\u003eL. japonica\u003c/em\u003e petal color. The current study results provide a foundation for research into petal color mechanisms at the transcriptome level. Transcriptome analysis based on an EST sequencing dataset illustrated the characteristics of the carnation transcriptome and revealed that carotenoid, chlorophyll, and anthocyanin biosynthesis played roles in the regulation of carnation flower color [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eCorrelation analysis showed that the changes in Paeonia flower color, from coral to pink to pale yellow were due to a significant decrease in anthocyanin content [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. It was deduced that the disequilibrium of expression levels in flavonol synthases and dihydroflavonol-4-reductases resulted in different levels of anthocyanin accumulation and coloration in white and pink tea flowers via metabolome and transcriptome analysis [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Besides, the transcriptome dataset assisted with the establishment of linkages between gene expression and other aspects, including the phenotypical, biochemical and metabolic features of \u003cem\u003eMeconopsis\u003c/em\u003e petal development [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e].\u003c/p\u003e \u003cp\u003ePetal color change is a significant feature of gene expression. As a whole, the number of differentially expressed genes at each flowering development stage was variable, which revealed that transcriptome expression changes were more active during color alteration. For example, when the petal color changed from purple-red (S1) to red (S2), 15804 genes were differentially expressed; while only 2033 genes changed from S2 to S3, when the petal color at both were red. Further, when the petal colors were red, but gradually became lighter from the three-green period (S2) to the great white period (S4), the number of differentially expressed genes were eventually lower, which were 21159, 19126, and 11153, respectively. Besides, gene expression changes occurred when the petal color transitioned from red (S4), to the outer part being red and the inner part being white (S5), where 6073 genes were up-regulated from the bud stage to cracking stage (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003ePetal color phenotypes result from plant pigment metabolism [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Anthocyanin, flavonoid, and alkaloid are well-known plant pigments that play active roles in the regulation of flower color [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Flavonoids are involved in the production of red, pink, purple and blue in plants and are widely distributed in petals [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. Additionally, and most importantly, flavonoid accumulation in a spontaneous cotton mutant resulted in red coloration and enhanced disease resistance [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. Anthocyanins are a group of water-soluble pigments that confer the blue, purple and red color to many fruits [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e]. Carotenoids primarily affect yellow, orange and red coloration [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. To identify the potential mechanisms that regulate the modification of flowering colors, the current study annotated differentially expressed genes via KEGG pathway analysis and Gene Ontology.\u003c/p\u003e \u003cp\u003eIn the pathway analysis, several pigment metabolism pathways were highlighted, including plant hormone signal transduction (ko04075), phenylpropanoid biosynthesis (ko00940), anthocyanin biosynthesis (ko00942), flavonoid biosynthesis (ko00941), and carotenoid biosynthesis (ko00906). Furthermore, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis were active throughout the flowering process and nearly up-regulated at each stage. These results corroborated the findings of several previous investigations into the linkages between petal color and plant pigment. Previous studies found that the flowers of \u0026lsquo;Yujin 2\u0026rsquo; were red due to their higher carotenoid and flavonoid content in contrast to \u0026lsquo;Damaohua\u0026rsquo; at all developmental stages [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eHowever, the results of our study revealed that carotenoid biosynthesis was down-regulated throughout S1-6, which was different from those of previous studies to some extent, due to the formation of mutants that varied with sampling times. Phenylalanine is an important anthocyanin biosynthesis pathway [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. In this study, phenylalanine synthesis was enhanced at the S1 phase, and anthocyanin synthesis was weakened at S5 and S6. This result may be consistent with the changes of \u0026lsquo;Yujin 2\u0026rsquo; petal color from red (S1) to white (S5) and gold (S6).\u003c/p\u003e \u003cp\u003eThe results obtained in the current study revealed that the phenylpropanoid biosynthetic process (GO:0009699), carotenoid biosynthetic process (GO:0016117), and anthocyanin 5-O-glucoside-4'''-O-malonyltransferase activity (GO:0102801) were triggered. Anthocyanin is the stable configuration of anthocyanidins that is always combined with a sugar moiety [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e], which may explain the active glucose-related metabolism pathways in our research.\u003c/p\u003e \u003cp\u003eGO annotation also supported the KEGG evidence. Furthermore, GO annotation showed that plasma membrane (GO:0005886) was active throughout the flowering process. Chloroplast (GO:0009507) and vacuole (GO:0005773) were active through most of the flowering processes, and could be explained by the results of previous experiments. The colors of plants are primarily codetermined by the different distributions of chlorophyll, carotenoid, flavonoid, and betaine, which exist in the cytoplasms or vacuoles of plants [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. Moreover, GO analysis revealed that environment-related responses (such as to cold) were triggered (GO:0009409).\u003c/p\u003e \u003cp\u003eGenerally, the binding of upstream transcription factors, initiation of gene transcription and expression, final synthesis of anthocyanins, and coloration of fruits and peels are closely related to the expression of structural genes in the anthocyanin biosynthetic pathway [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. At present, MYB, BHLH and WD40 are the three major transcription factors currently known. Anthocyanin biological regulation includes single-structure gene regulation, two-gene interaction regulation, and some MBW (MYB-bHLH-WD40) complexes regulate anthocyanin accumulation [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn the first three stages of flower development, LjDFR, LjABCB1, LjMYC6, LjDDB2, and LjANS genes increased rapidly. However, only the expression of LjF3'5'H was significantly down-regulated at S5, which was consistent with the accumulation of anthocyanin. The expression of other transcription factors did not increase during anthocyanin accumulation, which may not be directly related to anthocyanin synthesis. To verify the RNA-Seq results, we randomly selected six genes and tested their expression patterns via qPCR, which revealed the same expression trend with the sequencing results.\u003c/p\u003e "},{"header":"Conclusions","content":" \u003cp\u003eIn summary, this study developed a transcriptome profile of flower color generation for \u003cem\u003eL. japonica\u003c/em\u003e as well as annotated unigene sets of seven anthesis phases. Bioinformatics analysis revealed a number of distinct clues in regard to color development mechanisms. Further, candidate genes and potential pathways were obtained, which provided a foundation for further research at the transcriptome level. However, the current results still did not fully elucidate the overall regulatory mechanisms of plant pigments and the formation of the red color for \u003cem\u003eL. japonica\u003c/em\u003e. Therefore, future studies relating to the current topic, such as transcript regulation, protein expression, and biochemical metabolism, are recommended. Candidate key genes involved in anthocyanin metabolism should also be further investigated.\u003c/p\u003e "},{"header":"Methods","content":" \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003ePlant material preparation\u003c/h2\u003e \u003cp\u003eFive year old \u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo; were planted in the resource garden of the College of Life Sciences, at Henan Normal University (N35\u0026deg;18 '13.71 \", E113\u0026deg;55' 15.05). The land owner allows us to sample on their private land. Professor Jingyuan Li identifies \u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo; as honeysuckle [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e6\u003c/span\u003e). Our research group collected fresh flower buds of of \u0026lsquo;Damaohua\u0026rsquo; and \u0026lsquo;Yujin 2\u0026rsquo; in seven periods: young bud stage (S1), three green stage (S2), two white stage (S3), great white stage (S4), silver stage (S5), golden stage (S6), fade stage (S7). During sample collection, the flowers were combined and regarded as one biological replicate representing each stage, and three independent replicates were prepared. Partial flower materials were immediately frozen in liquid nitrogen following collection and stored at -80\u0026nbsp;\u0026deg;C.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eRNA isolation and library construction\u003c/h2\u003e \u003cp\u003eThe total RNA was extracted by mirVana miRNA Isolation Kit (Ambion) [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. RNA integrity was evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), and samples with the RNA Integrity Number (RIN)\u0026thinsp;\u0026ge;\u0026thinsp;7 were subjected to subsequent analysis. Separate mRNA from total RNA with oligonucleotide beads (dT) and add fragment buffer to cut it into short fragments, and then use random hexamer primers as templates to synthesize first strand cDNA. The library was constructed according to the manufacturer's instructions for TruSeq chain mRNA LTSample Prep Kit (Illumina, San Diego, California, USA).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eSequencing, de novo assembly and annotation\u003c/h2\u003e \u003cp\u003eThe above libraries were sequenced using an Illumina HiSeq X Ten sequencer (Illumina Inc., USA) and 150\u0026nbsp;bp paired-end reads were generated. The preparation of the cDNA library and sequencing were performed at Shanghai OE Biotech. Co., Ltd., Shanghai, China, and raw data (raw reads) were processed using Trimmomatic [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. The reads containing ploy-N and low quality reads were removed to obtain clean reads. Following the removal of adaptor and low quality sequences, the clean reads were assembled into expressed sequence tag clusters (contigs) and assembled de novo into transcripts using Trinity [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e] (version: trinityrnaseq_r20131110) through a paired-end method. The longest transcript was selected as a unigene based on the similarity and length of a sequence for subsequent analysis.\u003c/p\u003e \u003cp\u003eThe function of the unigenes was annotated to the NCBI non-redundant protein (NR), Clusters of Orthologous Groups (COG) for eukaryotic complete genomes, Gene Ontology (GO), Swiss-Prot, and evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases by diamond software, and mapped to Pfam databases by HMMER. The search was conducted using Blastx with a threshold E-value cut-off of 1e-5. Unigene quantification, analysis of differentially expressed unigenes (DEGs), cluster analysis, GO and KEGG enrichment FPKM [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e], and read count values of each unigene were calculated using bowtie2 [\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]and eXpress[\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eWe use the DESeq function to estimate the size factor and the nbinom test to determine DEG [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e]. The threshold value of P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and multiple change\u0026thinsp;\u0026gt;\u0026thinsp;2 or multiple change\u0026thinsp;\u0026lt;\u0026thinsp;0.5 was set as the significant difference threshold.To explore the transcriptional expression patterns of DEGs, cluster analysis was performed. We performed GO analysis and KEGG pathway enrichment analysis for DEGs based on hypergeometric distribution of R.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eExpression analysis of transcription factors (TFs)\u003c/h2\u003e \u003cp\u003ePlantTFDB (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://planttfdb.cbi.pku.edu.cn/index.php\u003c/span\u003e\u003c/span\u003e) contains sequences of 58 plant transcription factor families from 165 plants, and is a plant transcription factor database [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. Blastx was used to compare the sequence of unigenes with transcription factor database, and the best gene with E value less than 1e-5 was selected as annotation information of unigene. Candidates that contained DNA binding domains were recognized by GO annotation for the final TF identification.\u003c/p\u003e \u003cp\u003eDifferentially expressed TFs (DETFs) between samples were identified using the value of Fragments Per Kilobase of transcript Per Million fragments mapped (FPKM) with |log2(fold change)| \u0026gt; 1, p value\u0026thinsp;\u0026le;\u0026thinsp;0.05 and q value\u0026thinsp;\u0026le;\u0026thinsp;0.05 (Bolger et al. 2014). Multiexperiment Viewer software (v4.9) was employed to exhibit the expression profiles of DETFs by clustering [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e] and the DETFs were clustered via the K-means method. Using the total TFs in \u003cem\u003eL. japonica\u003c/em\u003e as reference, GO enrichment analysis of DETFs was performed using Agrigo (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://\u003c/span\u003e\u003c/span\u003e bioinfo.cau.edu.cn/agriGO/) with the hypergeometric test and FDR cut-off of 0.05 [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e].\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eqRT-PCR\u003c/h2\u003e \u003cp\u003eThe same RNA samples used for the RNA-seq experiments were employed for qRT-PCR. We use NanoDrop 2000 spectrophotometer to measure RNA yield (Thermo Scientific, USA), whereas the integrity was evaluated using agarose gel electrophoresis stained with ethidium bromide. Quantification was performed using a two-step reaction process: reverse transcription (RT) and PCR. Each sample was run in triplicate for analysis. At the end of the PCR cycles, melting curve analysis was performed to validate the specific generation of the expected PCR product. Quantification was performed using the 2\u0026thinsp;\u0026minus;\u0026thinsp;ΔΔCT method, and data were normalized to the ACT2/7 transcript [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. The sequences of the primers used are listed in Additional file 1: Table S1.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eStatistical analysis\u003c/h2\u003e \u003cp\u003eSignificant differences were calculated using a one-way ANOVA analysis with a Turkey test and a significance level at p \u0026le; 0.05 and p \u0026le; 0.01 in SPSS 19.0 software [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. All expression analyses were repeated three times. The reported value represents the arithmetic mean of three replicates. Data was expressed as a mean plus or minus a standard deviation (mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD).\u003c/p\u003e "},{"header":"Abbreviations","content":"\u003cp\u003e\u003cem\u003eL. japonica Lonicera japonica \u003c/em\u003eThunb\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eDEGs differentially expressed genes\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eCGA \u003cem\u003e chlorogenic acid \u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eCHS chalcone synthase, \u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eF3H flavanone 3-hydroxylase, \u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eDFR dihydroflavonol-4-reductase, \u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eANS anthocyanidin synthase, \u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eUFGT UDP -glucose flavonoid 3-O-glucosyltransferase \u003c/em\u003e\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and material\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data generated or analyzed for this study are included in this article and its supplementary information files.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe work was carried out under the auspices of the Public Health Services of Traditional Chinese Medicine in 2017 [(2017)66] and Henan Province Enterprise Technology Innovation Guidance Special Project (172107000031).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLJ.J. designed the research methodology, C C.F. and L Y. performed data analysis and wrote the article, and Y C.L,MJ.X. and C T prepared the plant samples. All authors have read and approved the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAt the point of finishing this paper, I\u0026rsquo;d like to express my sincere thanks to all those who have lent me hands in the course of my writing this paper.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eShang X, Pan H, Li M, Miao X, Ding H. \u003cem\u003eLonicera japonica\u003c/em\u003e Thunb.: ethnopharmacology, phytochemistry and pharmacology of an important traditional Chinese medicine. J Ethnopharmacol. 2011;138(1):1\u0026ndash;21. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.jep.2011.08.016\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHsu HF, Hsiao PC, Kuo TC, Chiang ST, Chen SL, Chiou SJ, Ling XH, Liang MT, Cheng WY, Houng JY. Antioxidant and anti-inflammatory activities of Lonicera japonica Thunb. var. sempervillosa Hayata flower bud extracts prepared by water, ethanol and supercritical fluid extraction techniques. Ind Crop Prod. 2016;89:543\u0026ndash;9. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.indcrop.2016.05.010\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e.3 Wang DY, Zhao XM, Liu YL. (2017) Hypoglycemic and hypolipidemic effects of a polysaccharide from flower buds of Lonicera japonica in streptozotocin-induced diabetic rats. International journal of biological macromolecules 102:396\u0026ndash;404. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.ijbiomac.2017.04.056\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTian J, Che HL, Ha D, Wei YP, Zheng SY. Characterization and anti-allergic effect of a polysaccharide from the flower buds of Lonicera japonica. Carbohyd Polym. 2012;90(4):1642\u0026ndash;7. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.carbpol.2012.07.044\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRyu KH, Rhee HI, Kim JH, Yoo H, Lee BY, Um KA, Kim K, Noh JY, Lim KM, Chung JH. Anti-Inflammatory and Analgesic Activities of SKLJI, a Highly Purified and Injectable Herbal Extract of Lonicera japonica. Biosci Biotech Bioch. 2010;74(10):2022\u0026ndash;8. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1271/bbb.100279\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHan J, Lv QY, Jin SY, Zhang TT, Jin SX, Li XY, Yuan HL. Comparison of anti-bacterial activity of three types of di-O-caffeoylquinic acids in Lonicera japonica flowers based on microcalorimetry. Chin J Nat Medicines. 2014;12(2):108\u0026ndash;13. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/S1875-5364(14)60017-0\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNing W, Peng X, Ma LY, Cui L, Lu XP, Wang J, Tian JK, Li XM, Wang W, Zhang L. Enhanced secondary metabolites production and antioxidant activity in postharvest \u003cem\u003eLonicera japonica\u003c/em\u003e Thunb. in response to UV radiation. Innov Food Sci Emerg. 2012;13:231\u0026ndash;43. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.ifset.2011.10.005\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang Z, Clifford MN, Sharp P. Analysis of chlorogenic acids in beverages prepared from Chinese health foods and investigation, in vitro, of effects on glucose absorption in cultured Caco-2 cells. Food Chem. 2008;108(1):369\u0026ndash;73. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.foodchem.2007.10.083\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi JJ, Lian XY, Ye CL, Wang L. Analysis of Flower Color Variations at Different Developmental Stages in Two Honeysuckle (\u003cem\u003eLonicera Japonica\u003c/em\u003e Thunb.) Cultivars. Hortscience. 2019;54(5):779\u0026ndash;82. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.21273/Hortsci13819-18\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eQu Y, Ou Z, Yang FS, Wang S, Peng J. The study of transcriptome sequencing for flower coloration in different anthesis stages of alpine ornamental herb (Meconopsis 'Lingholm'). Gene. 2019;689:220\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.gene.2018.12.017\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao D, Tao J. Recent advances on the development and regulation of flower color in ornamental plants. Front Plant Sci. 2015;6:261. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fpls.2015.00261\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYoshida K, Mori M, Kondo T. Blue flower color development by anthocyanins: from chemical structure to cell physiology. Nat Prod Rep. 2009;26(7):884\u0026ndash;915. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1039/b800165k\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54(4):733\u0026ndash;49. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1111/j.1365-313X.2008.03447.x\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTanaka Y, Brugliera F, Kalc G, Senior M, Dyson B, Nakamura N, Katsumoto Y, Chandler S. Flower Color Modification by Engineering of the Flavonoid Biosynthetic Pathway: Practical Perspectives. Biosci Biotech Bioch. 2010;74(9):1760\u0026ndash;9. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1271/bbb.100358\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSadilova E, Stintzing FC, Carle R. Anthocyanins, colour and antioxidant properties of eggplant (Solanum melongena L.) and violet pepper (Capsicum annuum L.) peel extracts. Z Naturforsch C. 2006;61(7\u0026ndash;8):527\u0026ndash;35.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Y, Cui BH, Tan ZW, Song BX, Cao HN, Zong CW. (2018) RNA sequencing and anthocyanin synthesis-related genes expression analyses in white-fruited Vaccinium uliginosum. Bmc Genomics 19. doi:ARTN 93010.1186/s12864-018-5351-0.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGao LX, Yang HX, Liu HF, Yang J, Hu YH. (2016) Extensive Transcriptome Changes Underlying the Flower Color Intensity Variation in Paeonia ostii. Frontiers in Plant Science 6. doi:ARTN 120510.3389/fpls.2015.01205.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTripathi AM, Niranjan A, Roya S. Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna. Plant Physiol Bioch. 2018;127:1\u0026ndash;10. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.plaphy.2018.03.008\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao L, Gao LP, Wang HX, Chen XT, Wang YS, Yang H, Wei CL, Wan XC, Xia T. The R2R3-MYB, bHLH, WD40, and related transcription factors in flavonoid biosynthesis. Funct Integr Genomic. 2013;13(1):75\u0026ndash;98. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s10142-012-0301-4\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBi H, Luang S, Li Y, Bazanova N, Morran S, Song Z, Perera MA, Hrmova M, Borisjuk N, Lopato S. Identification and characterization of wheat drought-responsive MYB transcription factors involved in the regulation of cuticle biosynthesis. J Exp Bot. 2016;67(18):5363\u0026ndash;80. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/jxb/erw298\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZheng J, Wu H, Zhu HB, Huang CY, Liu C, Chang YS, Kong ZC, Zhou ZH, Wang GW, Lin YJ, Chen H. Determining factors, regulation system, and domestication of anthocyanin biosynthesis in rice leaves. New Phytol. 2019;223(2):705\u0026ndash;21. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1111/nph.15807\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWiczkowski W, Szawara-Nowak D, Topolska J. Changes in the content and composition of anthocyanins in red cabbage and its antioxidant capacity during fermentation, storage and stewing. Food Chem. 2015;167:115\u0026ndash;23. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.foodchem.2014.06.087\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTasaki K, Higuchi A, Watanabe A, Sasaki N, Nishihara M. Effects of knocking out three anthocyanin modification genes on the blue pigmentation of gentian flowers. Sci Rep. 2019;9(1):15831. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41598-019-51808-3\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTanase K, Nishitani C, Hirakawa H, Isobe S, Tabata S, Ohmiya A, Onozaki T. Transcriptome analysis of carnation (Dianthus caryophyllus L.) based on next-generation sequencing technology. BMC Genom. 2012;13:292. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1186/1471-2164-13-292\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuo L, Wang Y, da Silva JAT, Fan Y, Yu X. Transcriptome and chemical analysis reveal putative genes involved in flower color change in Paeonia 'Coral Sunset'. Plant Physiol Biochem. 2019;138:130\u0026ndash;9. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.plaphy.2019.02.025\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou C, Mei X, Rothenberg DO, Yang Z, Zhang W, Wan S, Yang H, Zhang L. (2020) Metabolome and Transcriptome Analysis Reveals Putative Genes Involved in Anthocyanin Accumulation and Coloration in White and Pink Tea (Camellia sinensis) Flower. Molecules 25 (1). doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3390/molecules25010190\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilliams CA, Grayer RJ. Anthocyanins and other flavonoids. Nat Prod Rep. 2004;21(4):539\u0026ndash;73. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1039/b311404j\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLong L, Liu J, Gao Y, Xu FC, Zhao JR, Li B, Gao W. Flavonoid accumulation in spontaneous cotton mutant results in red coloration and enhanced disease resistance. Plant Physiol Biochem. 2019;143:40\u0026ndash;9. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.plaphy.2019.08.021\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFang J. Classification of fruits based on anthocyanin types and relevance to their health effects. Nutrition. 2015;31(11\u0026ndash;12):1301\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.nut.2015.04.015\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRosas-Saavedra C, Stange C. Biosynthesis of Carotenoids in Plants: Enzymes and Color. Subcell Biochem. 2016;79:35\u0026ndash;69. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/978-3-319-39126-7_2\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVeitch NC, Grayer RJ. Flavonoids and their glycosides, including anthocyanins. Nat Prod Rep. 2011;28(10):1626\u0026ndash;95. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1039/c1np00044f\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang B, Guan Q, Tian J, Komatsu S. Transcriptomic and proteomic analyses of leaves from Clematis terniflora DC. under high level of ultraviolet-B irradiation followed by dark treatment. Journal of proteomics. 2017;150:323\u0026ndash;40. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.jprot.2016.10.001\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114\u0026ndash;20. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/bioinformatics/btu170\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGrabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, Chen ZH, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644-U130. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/nbt.1883\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTrapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511\u0026ndash;5. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/nbt.1621\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLangmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357-U354. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/Nmeth.1923\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRoberts A, Pachter L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods. 2013;10(1):71\u0026ndash;3. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/nmeth.2251\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGupta R, Dewan I, Bharti R, Bhattacharya A. (2012) Differential Expression Analysis for RNA-Seq Data. ISRN Bioinform 2012:817508. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.5402/2012/817508\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJin J, Tian F, Yang DC, Meng YQ, Kong L, Luo J, Gao G. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017;45(D1):D1040\u0026ndash;5. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/nar/gkw982\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-SEq. Nat Methods. 2008;5(7):621\u0026ndash;8. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/nmeth.1226\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLivak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402\u0026ndash;8. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1006/meth.2001.1262\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Lonicera Japonica Thunb, Anthocyanin, Yujin 2, Comparative transcriptomics, RNA-seq","lastPublishedDoi":"10.21203/rs.3.rs-104094/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-104094/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eBackground：‘Yujin 2’ is new variety of \u003cem\u003eLonicera japonica \u003c/em\u003eThunb and its flower color can change from red to yellow; hence, it is a good model for investigating flower color development mechanisms. \u003c/p\u003e\u003cp\u003eResults：High throughput transcriptome sequencing of seven flower development stages of Yujin No.2 was carried out, and 133,487 unigenes were annotated, among which 73,088 were differentially expressed. Then the real-time PCR analysis was carried out. Further, the number of up-regulated DEGs was higher than those that were down-regulated. Of these annotated DEGs, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis were active throughout the flowering process during each stage, whereas carotenoid biosynthesis was inactive in the S1-6 stages. Furthermore, phenylalanine synthesis was enhanced in the S1 phase; however, anthocyanin synthesis was weakened in the S5 and S6 phases, which may be consistent with the changes in petal color of ‘Yujin 2’ from red (S1) to white (S5) and gold (S6). The results showed that 114 unigenes were associated with anthocyanin metabolism, and 72 were significantly upregulated or downregulated. According to the analysis of TFs in anthocyanin metabolism, we obtained 47 transcription factors, which belonged to 18 families. The LjDFR, LjABCB1, LjMYC6, LjDDB2, and LjANS genes rapidly increased during the first three stages. However, only LjF3'5'H expression was significantly down-regulated at S5, which was consistent with anthocyanin accumulation. \u003c/p\u003e\u003cp\u003eConclusions：This study developed a transcriptome profile of flower color generation for L. japonica as well as annotated unigene sets of seven anthesis phases, thereby providing possibilities for improving the germplasm of \u003cem\u003eL. japonica\u003c/em\u003e with genetic engineering technologies and cultivating new varieties of different colors.\u003c/p\u003e","manuscriptTitle":"Time-Series Transcriptome Sheds Light on Anthocyanin Metabolism and Species Diversification Between two Lonicera Japonica Thunb Cultivars","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2020-11-12 17:04:13","doi":"10.21203/rs.3.rs-104094/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"a186d82f-67b4-48b5-9ca9-9045fa3b786b","owner":[],"postedDate":"November 12th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":1059484,"name":"Plant Physiology and Morphology"},{"id":1059485,"name":"Plant Molecular Biology and Genetics"}],"tags":[],"updatedAt":"2020-12-04T17:40:31+00:00","versionOfRecord":[],"versionCreatedAt":"2020-11-12 17:04:13","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-104094","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-104094","identity":"rs-104094","version":["v1"]},"buildId":"WrCJVZZCHTDjtuVLN7oU0","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}