Next Article in Journal
Comparative Analyses of Chloroplast Genomes for Parasitic Species of Santalales in the Light of Two Newly Sequenced Species, Taxillus nigrans and Scurrula parasitica
Next Article in Special Issue
Floral Development Stage-Specific Transcriptomic Analysis Reveals the Formation Mechanism of Different Shapes of Ray Florets in Chrysanthemum
Previous Article in Journal
Effect of Dietetic Obesity on Testicular Transcriptome in Cynomolgus Monkeys
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium ‘Siberia’

1
Fujian Academy of Agricultural Sciences Biotechnology Institute, Fuzhou 350003, China
2
College of Horticulture, Fujian Agriculture and Forestry University, Fuzhou 350002, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2023, 14(3), 559; https://doi.org/10.3390/genes14030559
Submission received: 24 November 2022 / Revised: 21 February 2023 / Accepted: 22 February 2023 / Published: 23 February 2023
(This article belongs to the Special Issue Genetic Studies of Ornamental Horticulture and Floriculture)

Abstract

:
The bulb formation of Lilium is affected by many physiological and biochemical phenomena, including flower bud differentiation, starch and sucrose accumulation, photoperiod, carbon fixation, plant hormone transduction, etc. The transcriptome analysis of flower buds of Lilium hybrid ‘Siberia’ at different maturity stages showed that floral bud formation is associated with the accumulation of anthocyanins. The results of HPLC-MS showed that cyanidin is the major anthocyanin found in Lilium ‘Siberia’. Transcriptome KEGG enrichment analysis and qRT-PCR validation showed that two genes related to flavonoid biosynthesis (LhANS-rr1 and LhDFR) were significantly up-regulated. The functional analysis of differential genes revealed that LhMYB114 was directly related to anthocyanin accumulation among 19 MYB transcription factors. Furthermore, the qRT-PCR results suggested that their expression patterns were very similar at different developmental stages of the lily bulbs. Virus-induced gene silencing (VIGS) revealed that down-regulation of LhANS-rr1, LhDFR, and LhMYB114 could directly lead to a decrease in anthocyanin accumulation, turning the purple phenotype into a white color. Moreover, this is the first report to reveal that LhMYB114 can regulate anthocyanin accumulation at the mature stage of lily bulbs. The accumulation of anthocyanins is an important sign of lily maturity. Therefore, these findings have laid a solid theoretical foundation for further discussion on lily bulb development in the future.

1. Introduction

Lilies (Lilium, spp.), a group of monocotyledonous ornamental plants, are widely grown for commercial purposes. They significantly contribute to the global ornamental industry and are utilized for commercial bulb and flower production, including outdoor and indoor fresh-cut flowers and potted plants, and for landscaping in private gardens [1,2].
Anthocyanins, a type of flavonoid, are secondary metabolites, exist in many horticultural plants, and have multiple biological roles [3], such as growth inhibition of tumors [4], as well as anti-inflammation and -oxidation [5]. Furthermore, anthocyanins act as plant-coloring substances, play a vital role in defense mechanisms, and protect the plants from UV damage [5]. The rate of anthocyanin synthesis in the plant, membrane transport, and utilization or degradation are the main factors to influence the final anthocyanin concentration in ripe fruits [6,7]. Anthocyanin biosynthesis in fruits may be regulated by the activities of their metabolism-related enzymes, such as L-phenylalanine ammonia-lyase (PAL), cinnamate 4-hydrogenase (C4H), 4-coumarate: coenzyme A Ligase (4CL), chalcone synthase (CHS), UPD-3-O- glycosyltransferase (UFGT), and glutathione S-transferase (GST) [8,9,10,11]. In the anthocyanin biosynthesis pathway, the phenylalanine is converted into cinnamic acid by PAL [12,13]. The 4CL and CHS enzymes are derived from polyketone synthases (PKSs) [14] and are involved in the synthesis of naringin chalcone [15]. Likewise, the main function of the UFGT enzyme is to stabilize anthocyanin biosynthesis by attaching the sugar portion to anthocyanin glycogens, and it is a rate-limiting gene involved in anthocyanin biosynthesis [16]. The GST enzyme is very important for anthocyanin transportation in plants, and the results show that it is a key factor in the expression of anthocyanin-based red coloration in bracts [17].
With the cloning of anthocyanin-related transcription factors, anthocyanin biosynthesis in plants has been formally studied [18]. In plants, anthocyanin synthesis is mainly related to three elements, including two R2R3-MYB type TF, two basic helix-loop-helix type TF (bHLH), and one WD40 repeat TF [19], and their interaction [20]. Anthocyanins have been widely studied in Arabidopsis, pear, rice, maize, apple, petunia, and lily [21]. The overexpression of R2R3-MYB type TFs in Arabidopsis, including MYB114, MYB113, MYB90, and MYB75, enhanced anthocyanin accumulation [22]. Furthermore, in an experiment on red-skinned pear, it was revealed that PyMYB114 could inhibit anthocyanin biosynthesis, which was similar to the function in tobacco and strawberry [23]. Further research suggests that ethylene response factors interacted with PyMYB114 to regulate anthocyanin biosynthesis in red-skinned pears [24]. The transcriptome data on Hydrangea macrophylla revealed that MYB114 acted as a negative regulation transcription factor in anthocyanin synthesis [25]. A WD40 repeat gene (OsTTG1) regulated anthocyanin biosynthesis in rice [26]. ZmC1 (MYB-type C1) and ZmR (bHLH R) are two TFs; their up-regulation controlled anthocyanin biosynthesis in maize [27]. In apples, most MYB TFs positively regulated anthocyanin biosynthesis, except for MdMYB16 [28]. MdMYB114’s transcription level was increased with the deepening of apple color, and overexpression promoted the production and accumulation of anthocyanin in callus [29].
The buildup of anthocyanin in lily varieties causes the difference in floral color [30]. Several genes and TFs are involved in the anthocyanin synthesis of lily. In the Lilium cultivar ‘Montreux’, LhMYB6 and LhMYB12 were isolated, which interacted with LhbHLH2 protein and positively regulated anthocyanin biosynthesis [31]. During flower development, LhMYB12 could directly activate the promoters of chalcone synthase and dihydroflavonol 4-reductase in Asiatic hybrid lily [32]. The new allele of the LhMYB12, named LhMYB12-Lat, determined the presence or absence of splatters on the tepal of the Asiatic hybrid lily [33]. However, R3-MYB TFs usually play a native role in repressing anthocyanin biosynthesis [34]. For example, the transient expression of the LvMYB1 gene suggests that it inhibits anthocyanin biosynthesis in lily flowers [35].
In this study, a LhMYB114 transcription factor was isolated from Oriental hybrid Lilium ‘Siberia’ based on the transcriptome data. Our study suggested that it can effectively regulate the genetic expression of ANS and DFR via the VIGS system, thereby affecting the accumulation of anthocyanins in the floral bud of the Lilium bulb.

2. Materials and Methods

2.1. Experimental Material

The floral buds of Lilium hybrid ‘Siberia’ were sampled at six development stages according to the size of bulb circumference (Figure 1), where anthocyanin accumulation gradually increased from stage C to stage F. We finally focused on the last two stages (Figure 1E,F) for transcriptomics sequencing, named stage E and F.

2.2. Anthocyanidin Identification

Samples were more than 5 g. All samples were homogenized in a homogenizer and stored at −18 °C. The samples were mixed with the HCl (0.5% v/v) and 70% methanol buffer to a volume of 50 mL, and this process was repeated 3 times. Then, the samples were shook for 1 min and extracted for 30 min by ultrasonication. It was hydrolyzed in a boiling water bath for 1 h, and then cooled in constant volume. Standing and taking the supernatant, it was filtered with 0.45 aqueous phase filter membrane, and stored at 4 degrees for test. The injection volume of 20 μL was passed through the chromatographic column C18. The mobile phase A is an aqueous solution containing 1% formic acid, and the mobile phase B is a 100% acetonitrile solution at 0.25 mL/min flow rate. It was detected by a wavelength of 530 nm at LC-QTOF (USA, Agilent 6545B). Delphinidin, cyanidin, petunidin, pelargonin, and mallow pigment were used as an internal standard for quantitation.

2.3. Transcriptomic Sequencing and Differential Gene Expression Analysis

Total RNA was extracted from the samples of two developmental stages (Figure 1E,F) using RNAprep pure Plant Kit (Tiangen, Beijing, China; Code:DP360) according to the manufacturer’s protocol. HiScript® II 1st Strand cDNA Synthesis Kit (+gDNA wiper) (Vazyme, Nanjing, China) was used as reverse transcription kit. The quantity and quality of the cDNA were assessed by using a NanoDrop ND2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). cDNA library construction and transcriptome sequencing were performed by Novogene Technologies Co., Ltd. (Beijing, China). The quality of the sample libraries was determined by an Agilent 2100 Bioanalyzer and Qubit2.0. Finally, the well-constructed library was sequenced using an Illumina NovaSeq 6000 platform (Beijing Novogene Technology, Beijing, China). The R package DESeq2 (version 1.10.1) was used to identify differential expression analysis of two groups (E&F), and the gene expression values were calculated as reads aligned to the fold change of the normalized (RPKM) expression values by Ballgown (version 2.14.1). FDR was used to judge the threshold of p value. DEGs mainly depended on multiples hypotheses, such as the p-value threshold, false discovery rate (FDR) (≤0.05) and (log 2 FoldChange ≥ 1), which were considered as the significance of each gene expression difference.

2.4. Assay of Quantitative PCR

The real-time quantitative PCR (qRT-PCR) was performed on an ABI Q1 Real-Time PCR System with 10 µL mixture according to the following procedures: 95 °C for 1 min, 40 cycles of 95 °C for 15 s, and 60 °C for 30 s. The primers were designed using the software Primer Premier 6.0 and are listed in Table S1. Lilium × formolongi EF-1a was selected as a reference gene and the relative expression was analyzed with the 2−∆∆Ct method [36]. Three independent biological as well as technical replications were used.

2.5. Virus-Induced Gene Silencing in Lilium

Silencing of the targeted gene (LhMYB114, LhANS-rr1, or LhDFR) by VIGS was performed as described earlier [37,38]. According to the experiment protocol, a 200–500 base pair (bp) fragment, specific to the target gene, was designed and cloned into the pTRV2 vector (primers are listed in Table S1) using a One-step Cloning Kit (Vazyme, Nanjing, China), which was named 35S:LhMYB115, 35S:LhANS-rr1, and 35S:LhDFR. The blank vector of TRV1, TRV2, and the target genes inserted into TRV2 were transformed individually into the Agrobacterium GV3101 strain. Then, monoclonal bacteria were picked up and inoculated in a 100 mL LB medium containing antibiotics, then grown overnight at 28 °C in a shaker. Agrobacterium cells were collected and resuspended in infiltration buffer (10 mM MgCl2, 200 μM acetosyringone, and 10 mM 2-(N morpholino) ethane sulfonic acid (pH 5.6) to a final OD600 of 1.8). Equal volumes of TRV1 and TRV2 (the control), as well as the TRV1 and TRV2 with target genes, were mixed together and kept in the dark for 3–6 h at room temperature before infiltration. The buds were infiltrated using a needleless 1 mL syringe and grew at room temperature for seven days. Then, the flower buds continued vernalization in a refrigerator at 4 °C under dark for 15 days.

2.6. Electron Microscopy and Morphological Observation

To evaluate the features of apical meristem structure, the apical meristem from Lilium bulbs was cut for scanning electron microscopy (SEM) at 6 developmental stages. The cutting sections were soon fixed in 2.5% glutaraldehyde, and then dehydrated in a series of graded ethanol, as previously described [39]. The samples were observed and photographed with an SEM SU-8010 (Hitachi Ltd., Tokyo, Japan) (Figure S1).

2.7. Statistical Analysis

All data are presented as the mean ± standard deviation (SD) of at least three independent replicates. All statistical analyses in this paper were performed with SPSS (version 22.0, USA) software. The significant and extreme differences were confirmed if p-values were <0.05 and <0.01, respectively.

3. Results

3.1. The Chromatographic Analysis of Purple Anthocyanin in Flower Buds

To identify the major anthocyanin compound in Lilium flower buds, five types of anthocyanins were used as standard, including delphinidin, cyanidin, petunia, pelargonin, and mallow (Figure 2). Cyanidin was detected as the major anthocyanin compound in Lilium flower buds, as suggested by the results of HPLC-MS.

3.2. Identification of DEGs Related to Flavonoid Biosynthesis

Based on transcriptome data, putative DEGs from E and F groups identified about 8944 DEGs (Figure 3C), including 4557 transcripts showing up-regulation and 4387 transcripts showing up-regulation. The Arabidopsis genome acted as an internal reference. Interestingly, we detected that the sixteen expression genes related to flavonoid biosynthesis were dramatically changed. All the differential expression patterns and gene annotations are shown in the heatmap between E and F groups (Figure 3A). Lilium anthocyanins mainly comprised cyanidins; the metabolic pathway is presented in Figure 3B. The quantitative PCR verification showed that the expressions of anthocyanin-related genes were highly consistent with transcriptome analysis, and F3H, LhDFR, and LhANS-rr1 were up-regulated (Figure 3D).

3.3. MYB Gene and Cluster Analysis

Based on full-length transcriptomics, 28 kinds of TFs were enriched, such as C3H, MYB-ralated, C2H2, bHLH, GRAS, FAR1, and AP2/EFF-ERF, etc. Among them, MYB and bHLH were mainly involved in anthocyanin biosynthesis. The number of MYBs reached 33, and bHLH reached 54 (Figure 4A). We finally selected the MYB transcription factor as the research object, which was the largest one related to anthocyanin synthesis. The results showed that 19 major MYB TFs were related to the differential expression in DEGs, including 8 up-regulated genes and 11 down-regulated genes. Furthermore, only LhMYB114 was involved in the regulation of anthocyanin biosynthesis (Figure 4B). The result of the evolutionary tree with Arabidopsis showed that it was the homologous genes, with many anthocyanin-biosynthesis-related genes, such as Atmyb123, Atmyb114, Atmyb113, Atmyb90, and Atmyb75, in Arabidopsis (Figure 4C).

3.4. LhMYB114 Gene Cloning and Phylogenetic Analysis

We cloned the LhMYB114 gene from Lilium ‘Siberia’, which includes 214 amino acid sequences and TA insert to pGXT vector for sequencing. The results of structural alignment showed that it belonged to the members of MYB gene family (PLN03212), including R2-R3 repeat domain and the bHLH binding domain (Figure 5). The R2R3-MYB related to anthocyanin pigmentation patterns was highly homologous with LhMYB114 (Figure 5A).

3.5. The Expression Analysis of Three Key Genes Associated with Anthocyanin Metabolism in Flower Buds

We analyzed the expression patterns of three key genes in different development stages, i.e., LhMYB114, LhANS-rr1, and LhDFR. The results showed the expression patterns of the three genes were significantly similar (Figure 6).

3.6. Silencing Flavonoid-Biosynthesis-Related Genes Reduced Anthocyanin Accumulation

To determine whether three key genes were essential for anthocyanin accumulation in flower buds, we knocked down LhMYB114, LhANS-rr1, and LhDFR by VIGS. Agrobacterium strains of pTRV1 and pTRV2 were mixed in a ratio of 1:1, and buds were infected by syringe and grown for approximately 7 dpi at 25 °C in the dark. Then, we turned it to 7–8 °C ice storage to continue germination for the experiment. The leaf of gene-silenced buds was white around the injection site at 15 dpi, whereas that of control leaf continued to be purple in color (Figure 7A). The PCR amplification results suggested that the target band of RNA1 and RNA2 could be detected in new grown buds of the positive control (empty vector) and gene-silenced plants (pTRV2-LhMYB114, LhANS-rr1, and LhDFR vector), but there were no bands in the negative control (buffer) (Figure 6B). The phenotype and testing results showed that TRV had successfully invaded the buds of Lilium and replicated and transferred in vivo. The qRT-PCR results suggested that the expression of three key genes was evidently reduced in gene-silenced buds compared with control buds (empty vector) at 15 dpi (Figure 8A–F). Other genes in the cyanidin pathway were up-regulated, except for the three major genes, which was the opposite of positive regulation (Figure 8G). These results suggested that three key genes played a vital role in anthocyanin accumulation in Lilium buds.

4. Discussion

The color is an important attribute that determines the quality of flowers [40,41]. The purple color with the development and vernalization in lily buds, which is due to the accumulation of anthocyanins, has become an important target in lily breeding [42]. The anthocyanin accumulation of the bud is not only an indicator of bulb maturity, but also crucial for the development of the bulb and improvements in cold storage resistance [43]. Therefore, it is necessary to analyze the biosynthetic process involved in anthocyanin accumulation in buds, to improve the quality of Lilium buds. In the present study, the results reveal that three genes of LhANS-rr1, LhDFR, and LhMYB114 are closely related to anthocyanin biosynthesis in lily buds. Our results (according to VIGS) suggested that three key genes might be closely related to anthocyanin accumulation in floral buds of lilies.
During the process of bulb development, the flower bud of a lily becomes larger with the expansion of bulb size, and anthocyanin accumulation simultaneously occurs. The coordinated expression of genes, which encode the anthocyanin biosynthetic pathway enzymes, usually controls the anthocyanin accumulation through a ternary MYB-bHLH-WD40 (MBW) transcription complex [44]. Anthocyanin accumulation mainly happens in lily tepals and causes bicolor, bud blush, splatter spots, or raised spots [45,46]. The anthocyanin-accumulation-pathway-related genes in flower buds were significantly enriched in all DEGs. Previous studies showed that a single enzyme-coding structural gene might be involved in anthocyanin biosynthesis, including C4H, 4CL, CHS, CHI, UFGT, F3H, DFR, and ANS genes. Among the unigenes encoding six CCoAOMT, UDPG, HCT, C4H, 4CL, CHS, CHI, F3H, LhDFR, LhANS-rr1, and UFGT, the results suggested that the FPKM values of the three CCoAOMT, LhANS-rr1, and LhDFR unigenes were higher in the purple buds of bulbs than in the white buds, whereas the FPKM values of three CCoAOMT, UDPG, HCT, C4H, 4CL, CHS, CHI, F3H, and UFGT were not (Figure 3A). Meanwhile, the qRT PCR results confirmed that the expression of F3H, LhANS-rr1, and LhDFR was significantly up-regulated (Figure 3D). Furthermore, we clone the LhDFR and LhANS-rr1 genes from Lilium hybrid ‘Siberia’. The cDNAs in NCBI blast show that both have their alleles in Lilium speciosum (LsDFR-ws and LsANS-rr1), which proved that they were related to anthocyanin biosynthesis [47].
There are three main downstream metabolic pathways of anthocyanins, i.e., cyanidin, pelargonium, and delphinidin, in plants. Though ANS and DFR genes do not specifically regulate anthocyanin biosynthesis, they are the key genes in the cyanidin biosynthesis pathway [48]. Previous studies proved that the structural genes, including DFR, ANS, and 3GT, directly determine the flower color formation in lilies [49]. In our findings, the results of qRT-PCR showed that their expression levels were significantly different at different development stages of Lilium bulbs (Figure 6B,C). Further explanation suggested that they might be involved in anthocyanin accumulation in floral lily buds.
MYB TFs regulate the phenylpropanoid pathway, including lignin, flavonoids, and other metabolites, in plants [50]. The MBW complex related to anthocyanin biosynthesis is regulated by MYB TFs [51]. It is important that MYB is the most abundant TF related to anthocyanins. We analyzed all the DEGs and found that there were 19 major MYB transcription factors. One of them, named LhMYB114, was involved in anthocyanin synthesis (Figure 4). The LhMYB114 sequencing blast results suggest that it is the allele gene of Lilium regale LrMYB15 for transcription factor R2R3-MYB [52]. The Cluster analysis results also show that LhMYB114 is homologous with anthocyanin-related genes in Arabidopsis thaliana, including AtMYB123, AtMYB5, AtMYB114, AtMYB113, AtMYB90, and AtMYB75. Most R2R3-MYB TFs can promote anthocyanin biosynthesis by up-regulating the expression level in plants [53]. Our results show that the LhMYB114 expression level increases significantly in the mature stage of lily bulb development (Figure 6A).
The structural genes may be regulated by a single MYB transcription factor or the MBW complex [54]. VIGS mediated by Agrobacterium tumefaciens is a powerful tool to reveal gene functions in plants, which are difficult to transform, such as lily, potato, etc. [55]. Thus, to prove the regulatory response relationship among them, gene transient silencing is carried out by virus-induced gene silencing (VIGS). The pTRV1 and pTRV2 gene vector-infiltrated flower buds suggested an obvious white phenotype compared with empty vector (pTRV1 and pTRV2) at 15 dpi, including 35S:LhMYB114, 35S:LhANS-rr1, 35S:LhDFR, 35S:LhMYB114+35S:LhANS-rr1, 35S:LhMYB114+35S:LhDFR, and 35S:LhMYB114+35S:LhANS-rr1+LhDFR (Figure 8A). The fading phenotype of LhMYB114 and LhANS-rr1 was more typical than that of LhDFR after VIGS experiments (Figure 7A). At the same time, the qRT-PCR results from six groups showed that the gene expression levels were significantly down-regulated in the new growth buds (Figure 8A–F). This indicated that the three genes were positive regulators of anthocyanin synthesis. Their upstream genes were correspondingly up-regulated, such as C4H, 4CL, CHS, CHI, and F3H (Figure 8G). The quantitative data proved that three genes were the key factors in regulating anthocyanin accumulation through the cyanidin pathway in the floral buds of lilies. It also laid a foundation for us to analyze the relationship between anthocyanins and vernalization in the future.

5. Conclusions

In this study, it was revealed that the maturation of lily bulbs is closely associated with the amount of anthocyanin accumulation. Therefore, anthocyanin accumulation in flower buds is a very important marker in the development of lily bulbs, and the concentration of anthocyanins directly decides on the ripeness of bulbs. In addition, the maturity of lily bulbs will also directly decide on the time of bulb quality and lily vernalization. The transcriptomics and VIGS analysis firstly revealed that LhMYB114, LhANS-rr1, and LhDFR were the key genes for regulating the accumulation of cyanidin in floral lily buds. Among them, LhMYB114 is an important regulatory transcription factor during lily bulb maturation. In a word, the data in our article contribute to the understanding for the molecular mechanism of anthocyanin biosynthesis and the judging of bulb maturity in Lilium flower bud development.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14030559/s1, Figure S1: The microscopic view of epical meristem of Lilium hybrid ‘Siberia’ at E and F stage, Table S1: The sequences of primers used in the study, Technical terms: full name of genes.

Author Contributions

Conceptualization, Z.L.; methodology, S.F., M.L., X.Y. and S.W.; software, M.L.; validation, M.L., M.M.A. and Z.L.; formal analysis, S.F.; investigation, Z.L.; resources, Z.L.; data curation, Y.Z., M.L. and M.M.A.; writing—original draft preparation, S.F., M.L. and M.M.A.; writing—review and editing, M.L., X.Y. and Z.L.; supervision, Z.L.; project administration, S.W. and Z.L.; funding acquisition, F.C. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Fujian Forestry Science and Technology Project (ZMGG-0711 to Z.L. and F.C.) and Science and Technology Innovation Platform Project of Fujian Academy of Agricultural Sciences (CXPT202204 to Z.L.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bakhshaie, M.; Khosravi, S.; Azadi, P.; Bagheri, H.; van Tuyl, J.M. Biotechnological Advances in Lilium. Plant Cell Rep. 2016, 35, 1799–1826. [Google Scholar] [CrossRef] [PubMed]
  2. Arzate-Fernández, A.-M.; Nakazaki, T.; Okumoto, Y.; Tanisaka, T. Efficient Callus Induction and Plant Regeneration from Filaments with Anther in Lily (Lilium longiflorum Thunb.). Plant Cell Rep. 1997, 16, 836–840. [Google Scholar] [CrossRef]
  3. Ali, M.M.; Anwar, R.; Yousef, A.F.; Li, B.; Luvisi, A.; De Bellis, L.; Aprile, A.; Chen, F. Influence of Bagging on the Development and Quality of Fruits. Plants 2021, 10, 358. [Google Scholar] [CrossRef] [PubMed]
  4. Caltagirone, S.; Rossi, C.; Poggi, A.; Ranelletti, F.O.; Natali, P.G.; Brunetti, M.; Aiello, F.B.; Piantelli, M. Flavonoids Apigenin and Quercetin Inhibit Melanoma Growth and Metastatic Potential. Int. J. Cancer 2000, 87, 595–600. [Google Scholar] [CrossRef] [PubMed]
  5. Harborne, J.B.; Williams, C.A. Advances in Flavonoid Research since 1992. Phytochemistry 2000, 55, 481–504. [Google Scholar] [CrossRef] [PubMed]
  6. Cao, Y.; Han, Y.; Li, D.; Lin, Y.; Cai, Y. Systematic Analysis of the 4-Coumarate:Coenzyme A Ligase (4CL) Related Genes and Expression Profiling during Fruit Development in the Chinese Pear. Genes 2016, 7, 89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Haskill, J.S.; Häyry, P.; Radov, L.A. Systemic and Local Immunity in Allograft and Cancer Rejection. Contemp. Top. Immunobiol. 1978, 8, 107–170. [Google Scholar] [CrossRef]
  8. Fu, H.; Qiao, Y.; Wang, P.; Mu, X.; Zhang, J.; Fu, B.; Du, J. Changes of Bioactive Components and Antioxidant Potential during Fruit Development of Prunus humilis Bunge. PLoS ONE 2021, 16, e0251300. [Google Scholar] [CrossRef] [PubMed]
  9. Lister, C.E.; Lancaster, J.E.; Sutton, K.H.; Walker, J.R.L. Developmental Changes in the Concentration and Composition of Flavonoids in Skin of a Red and a Green Apple Cultivar. J. Sci. Food Agric. 1994, 64, 155–161. [Google Scholar] [CrossRef]
  10. Zoratti, L.; Karppinen, K.; Luengo Escobar, A.; HÃggman, H.; Jaakola, L. Light-Controlled Flavonoid Biosynthesis in Fruits. Front. Plant Sci. 2014, 5, 534. [Google Scholar] [CrossRef]
  11. Shi, M.; Ali, M.M.; He, Y.; Ma, S.; Rizwan, H.M.; Yang, Q.; Li, B.; Lin, Z.; Chen, F. Flavonoids Accumulation in Fruit Peel and Expression Profiling of Related Genes in Purple (Passiflora edulis f. Edulis) and Yellow (Passiflora edulis f. Flavicarpa) Passion Fruits. Plants 2021, 10, 2240. [Google Scholar] [CrossRef] [PubMed]
  12. Kong, J.-Q. Phenylalanine Ammonia-Lyase, a Key Component Used for Phenylpropanoids Production by Metabolic Engineering. RSC Adv. 2015, 5, 62587–62603. [Google Scholar] [CrossRef]
  13. Singh, R.; Rastogi, S.; Dwivedi, U.N. Phenylpropanoid Metabolism in Ripening Fruits. Compr. Rev. Food Sci. Food Saf. 2010, 9, 398–416. [Google Scholar] [CrossRef]
  14. Lu, Y.; Song, Y.; Zhu, J.; Xu, X.; Pang, B.; Jin, H.; Jiang, C.; Liu, Y.; Shi, J. Potential Application of CHS and 4CL Genes from Grape Endophytic Fungus in Production of Naringenin and Resveratrol and the Improvement of Polyphenol Profiles and Flavour of Wine. Food Chem. 2021, 347, 128972. [Google Scholar] [CrossRef]
  15. Flores-Sanchez, I.J.; Verpoorte, R. Plant Polyketide Synthases: A Fascinating Group of Enzymes. Plant Physiol. Biochem. 2009, 47, 167–174. [Google Scholar] [CrossRef]
  16. Li, X.-J.; Zhang, J.-Q.; Wu, Z.-C.; Lai, B.; Huang, X.-M.; Qin, Y.-H.; Wang, H.-C.; Hu, G.-B. Functional Characterization of a Glucosyltransferase Gene, LcUFGT1, Involved in the Formation of Cyanidin Glucoside in the Pericarp of Litchi Chinensis. Physiol. Plant. 2016, 156, 139–149. [Google Scholar] [CrossRef] [PubMed]
  17. Vilperte, V.; Boehm, R.; Debener, T. A Highly Mutable GST Is Essential for Bract Colouration in Euphorbia pulcherrima Willd. Ex Klotsch. BMC Genom. 2021, 22, 208. [Google Scholar] [CrossRef] [PubMed]
  18. LaFountain, A.M.; Yuan, Y. Repressors of Anthocyanin Biosynthesis. New Phytol. 2021, 231, 933–949. [Google Scholar] [CrossRef] [PubMed]
  19. Wen, C.; Tsao, N.; Wang, S.; Chu, F. Color Variation in Young and Senescent Leaves of Formosan Sweet Gum (Liquidambar formosana) by the Gene Regulation of Anthocyanidin Biosynthesis. Physiol. Plant. 2021, 172, 1750–1763. [Google Scholar] [CrossRef]
  20. Outchkourov, N.S.; Karlova, R.; Hölscher, M.; Schrama, X.; Blilou, I.; Jongedijk, E.; Simon, C.D.; van Dijk, A.D.J.; Bosch, D.; Hall, R.D.; et al. Transcription Factor-Mediated Control of Anthocyanin Biosynthesis in Vegetative Tissues. Plant Physiol. 2018, 176, 1862–1878. [Google Scholar] [CrossRef] [Green Version]
  21. Tanaka, Y.; Sasaki, N.; Ohmiya, A. Biosynthesis of Plant Pigments: Anthocyanins, Betalains and Carotenoids. Plant J. 2008, 54, 733–749. [Google Scholar] [CrossRef]
  22. Wada, T.; Kunihiro, A.; Tominaga-Wada, R. Arabidopsis CAPRICE (MYB) and GLABRA3 (BHLH) Control Tomato (Solanum lycopersicum) Anthocyanin Biosynthesis. PLoS ONE 2014, 9, e109093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yao, G.; Ming, M.; Allan, A.C.; Gu, C.; Li, L.; Wu, X.; Wang, R.; Chang, Y.; Qi, K.; Zhang, S.; et al. Map-Based Cloning of the Pear Gene MYB114 Identifies an Interaction with Other Transcription Factors to Coordinately Regulate Fruit Anthocyanin Biosynthesis. Plant J. 2017, 92, 437–451. [Google Scholar] [CrossRef] [Green Version]
  24. Ni, J.; Bai, S.; Zhao, Y.; Qian, M.; Tao, R.; Yin, L.; Gao, L.; Teng, Y. Ethylene Response Factors Pp4ERF24 and Pp12ERF96 Regulate Blue Light-Induced Anthocyanin Biosynthesis in ‘Red Zaosu’ Pear Fruits by Interacting with MYB114. Plant Mol. Biol. 2019, 99, 67–78. [Google Scholar] [CrossRef] [PubMed]
  25. Peng, J.; Dong, X.; Xue, C.; Liu, Z.; Cao, F. Exploring the Molecular Mechanism of Blue Flower Color Formation in Hydrangea Macrophylla Cv. “Forever Summer”. Front. Plant Sci. 2021, 12, 585665. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, X.; Wang, J.; Xia, X.; Zhang, Z.; He, J.; Nong, B.; Luo, T.; Feng, R.; Wu, Y.; Pan, Y.; et al. OsTTG1, a WD40 Repeat Gene, Regulates Anthocyanin Biosynthesis in Rice. Plant J. 2021, 107, 198–214. [Google Scholar] [CrossRef] [PubMed]
  27. Riaz, B.; Chen, H.; Wang, J.; Du, L.; Wang, K.; Ye, X. Overexpression of Maize ZmC1 and ZmR Transcription Factors in Wheat Regulates Anthocyanin Biosynthesis in a Tissue-Specific Manner. Int. J. Mol. Sci. 2019, 20, 5806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Gao, H.-N.; Jiang, H.; Cui, J.-Y.; You, C.-X.; Li, Y.-Y. Review: The Effects of Hormones and Environmental Factors on Anthocyanin Biosynthesis in Apple. Plant Sci. 2021, 312, 111024. [Google Scholar] [CrossRef] [PubMed]
  29. Jiang, S.; Sun, Q.; Zhang, T.; Liu, W.; Wang, N.; Chen, X. MdMYB114 Regulates Anthocyanin Biosynthesis and Functions Downstream of MdbZIP4-like in Apple Fruit. J. Plant Physiol. 2021, 257, 153353. [Google Scholar] [CrossRef] [PubMed]
  30. Tasaki, K.; Terada, H.; Masuta, C.; Yamagishi, M. Virus-Induced Gene Silencing (VIGS) in Lilium leichtlinii Using the Cucumber Mosaic Virus Vector. Plant Biotechnol. 2016, 33, 373–381. [Google Scholar] [CrossRef] [Green Version]
  31. Yamagishi, M.; Shimoyamada, Y.; Nakatsuka, T.; Masuda, K. Two R2R3-MYB Genes, Homologs of Petunia AN2, Regulate Anthocyanin Biosyntheses in Flower Tepals, Tepal Spots and Leaves of Asiatic Hybrid Lily. Plant Cell Physiol. 2010, 51, 463–474. [Google Scholar] [CrossRef] [Green Version]
  32. Lai, Y.-S.; Shimoyamada, Y.; Nakayama, M.; Yamagishi, M. Pigment Accumulation and Transcription of LhMYB12 and Anthocyanin Biosynthesis Genes during Flower Development in the Asiatic Hybrid Lily (Lilium spp.). Plant Sci. 2012, 193–194, 136–147. [Google Scholar] [CrossRef]
  33. Yamagishi, M.; Toda, S.; Tasaki, K. The Novel Allele of the Lh MYB 12 Gene Is Involved in Splatter-type Spot Formation on the Flower Tepals of Asiatic Hybrid Lilies (Lilium spp.). New Phytol. 2014, 201, 1009–1020. [Google Scholar] [CrossRef] [PubMed]
  34. Sakai, M.; Yamagishi, M.; Matsuyama, K. Repression of Anthocyanin Biosynthesis by R3-MYB Transcription Factors in Lily (Lilium spp.). Plant Cell Rep. 2019, 38, 609–622. [Google Scholar] [CrossRef]
  35. Yin, X.; Zhang, Y.; Zhang, L.; Wang, B.; Zhao, Y.; Irfan, M.; Chen, L.; Feng, Y. Regulation of MYB Transcription Factors of Anthocyanin Synthesis in Lily Flowers. Front. Plant Sci. 2021, 12, 761668. [Google Scholar] [CrossRef] [PubMed]
  36. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  37. Baum, S. The PAS Reaction for Staining Cell Walls. Cold Spring Harb. Protoc. 2008, 2008, pdb.prot4956. [Google Scholar] [CrossRef]
  38. Li, J.; Seng, S.; Li, D.; Zhang, F.; Liu, Y.; Yao, T.; Liang, J.; Yi, M.; Wu, J. Antagonism between Abscisic Acid and Gibberellin Regulates Starch Synthesis and Corm Development in Gladiolus hybridus. Hortic. Res. 2021, 8, 155. [Google Scholar] [CrossRef]
  39. Chehregani, A.; Kouhkan, F. Diesel Exhaust Particles and Allergenicity of Pollen Grains of Lilium martagon. Ecotoxicol. Environ. Saf. 2008, 69, 568–573. [Google Scholar] [CrossRef]
  40. Sadique, S.; Ali, M.M.; Usman, M.; Hasan, M.U.; Yousef, A.F.; Adnan, M.; Gull, S.; Nicola, S. Effect of Foliar Supplied PGRs on Flower Growth and Antioxidant Activity of African Marigold (Tagetes erecta L.). Horticulturae 2021, 7, 378. [Google Scholar] [CrossRef]
  41. Ali, M.M.; Yousef, A.F.; Li, B.; Chen, F. Effect of Environmental Factors on Growth and Development of Fruits. Trop. Plant Biol. 2021, 14, 226–238. [Google Scholar] [CrossRef]
  42. Salachna, P.; Byczyńska, A.; Zawadzińska, A.; Piechocki, R.; Mizielińska, M. Stimulatory Effect of Silver Nanoparticles on the Growth and Flowering of Potted Oriental Lilies. Agronomy 2019, 9, 610. [Google Scholar] [CrossRef] [Green Version]
  43. Dufoo-Hurtado, M.D.; Zavala-Gutiérrez, K.G.; Cao, C.-M.; Cisneros-Zevallos, L.; Guevara-González, R.G.; Torres-Pacheco, I.; Vázquez-Barrios, M.E.; Rivera-Pastrana, D.M.; Mercado-Silva, E.M. Low-Temperature Conditioning of “Seed” Cloves Enhances the Expression of Phenolic Metabolism Related Genes and Anthocyanin Content in ‘Coreano’ Garlic (Allium sativum) during Plant Development. J. Agric. Food Chem. 2013, 61, 10439–10446. [Google Scholar] [CrossRef]
  44. Allan, A.C.; Hellens, R.P.; Laing, W.A. MYB Transcription Factors That Colour Our Fruit. Trends Plant Sci. 2008, 13, 99–102. [Google Scholar] [CrossRef]
  45. Yamagishi, M. How Genes Paint Lily Flowers: Regulation of Colouration and Pigmentation Patterning. Sci. Hortic. 2013, 163, 27–36. [Google Scholar] [CrossRef]
  46. Yamagishi, M.; Uchiyama, H.; Handa, T. Floral Pigmentation Pattern in Oriental Hybrid Lily (Lilium spp.) Cultivar ‘Dizzy’ Is Caused by Transcriptional Regulation of Anthocyanin Biosynthesis Genes. J. Plant Physiol. 2018, 228, 85–91. [Google Scholar] [CrossRef]
  47. Suzuki, K.; Tasaki, K.; Yamagishi, M. Two Distinct Spontaneous Mutations Involved in White Flower Development in Lilium speciosum. Mol. Breed. 2015, 35, 193. [Google Scholar] [CrossRef]
  48. Jaakola, L.; Määttä, K.; Pirttilä, A.M.; Törrönen, R.; Kärenlampi, S.; Hohtola, A. Expression of Genes Involved in Anthocyanin Biosynthesis in Relation to Anthocyanin, Proanthocyanidin, and Flavonol Levels during Bilberry Fruit Development. Plant Physiol. 2002, 130, 729–739. [Google Scholar] [CrossRef] [Green Version]
  49. Yin, X.; Lin, X.; Liu, Y.; Irfan, M.; Chen, L.; Zhang, L. Integrated Metabolic Profiling and Transcriptome Analysis of Pigment Accumulation in Diverse Petal Tissues in the Lily Cultivar ‘Vivian’. BMC Plant Biol. 2020, 20, 446. [Google Scholar] [CrossRef]
  50. Ma, D.; Constabel, C.P. MYB Repressors as Regulators of Phenylpropanoid Metabolism in Plants. Trends Plant Sci. 2019, 24, 275–289. [Google Scholar] [CrossRef]
  51. Jin, W.; Wang, H.; Li, M.; Wang, J.; Yang, Y.; Zhang, X.; Yan, G.; Zhang, H.; Liu, J.; Zhang, K. The R2R3 MYB Transcription Factor PavMYB10.1 Involves in Anthocyanin Biosynthesis and Determines Fruit Skin Colour in Sweet Cherry (Prunus avium L.). Plant Biotechnol. J. 2016, 14, 2120–2133. [Google Scholar] [CrossRef] [Green Version]
  52. Yamagishi, M. A Novel R2R3-MYB Transcription Factor Regulates Light-Mediated Floral and Vegetative Anthocyanin Pigmentation Patterns in Lilium regale. Mol. Breed. 2016, 36, 3. [Google Scholar] [CrossRef]
  53. Li, C.; Yu, W.; Xu, J.; Lu, X.; Liu, Y. Anthocyanin Biosynthesis Induced by MYB Transcription Factors in Plants. Int. J. Mol. Sci. 2022, 23, 11701. [Google Scholar] [CrossRef] [PubMed]
  54. Lloyd, A.; Brockman, A.; Aguirre, L.; Campbell, A.; Bean, A.; Cantero, A.; Gonzalez, A. Advances in the MYB–BHLH–WD Repeat (MBW) Pigment Regulatory Model: Addition of a WRKY Factor and Co-Option of an Anthocyanin MYB for Betalain Regulation. Plant Cell Physiol. 2017, 58, 1431–1441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Gao, X.; Britt, R.C., Jr.; Shan, L.; He, P. Agrobacterium-Mediated Virus-Induced Gene Silencing Assay in Cotton. J. Vis. Exp. 2011, 54, 2938. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The appearance and morphology of Lilium bulb and bud at different developmental stages (AF). (A) The bulb diameter is 3 cm; (B) the bulb diameter is 6 cm; (C) the bulb diameter is 10 cm; (D) the bulb diameter is 14 cm; (E) the bulb diameter is 16 cm; (F) the bulb diameter is 22 cm. The data represent the average of three biological replicates.
Figure 1. The appearance and morphology of Lilium bulb and bud at different developmental stages (AF). (A) The bulb diameter is 3 cm; (B) the bulb diameter is 6 cm; (C) the bulb diameter is 10 cm; (D) the bulb diameter is 14 cm; (E) the bulb diameter is 16 cm; (F) the bulb diameter is 22 cm. The data represent the average of three biological replicates.
Genes 14 00559 g001
Figure 2. The concentrations of anthocyanin-related compounds in the flower buds of oriental lily ‘Siberia’. High-performance liquid chromatography-mass spectrometry (HPLC-MS) chromatograms suggested the presence of anthocyanin aglycones was completed at 530nm; the compounds were identified by comparing with the mass spectra of the standards. A—delphinidin, B/B1—cyanidin, C—petunidin, D—pelargonin, E—mallow pigment. Black color represents the standard, and red or green represents the text sample.
Figure 2. The concentrations of anthocyanin-related compounds in the flower buds of oriental lily ‘Siberia’. High-performance liquid chromatography-mass spectrometry (HPLC-MS) chromatograms suggested the presence of anthocyanin aglycones was completed at 530nm; the compounds were identified by comparing with the mass spectra of the standards. A—delphinidin, B/B1—cyanidin, C—petunidin, D—pelargonin, E—mallow pigment. Black color represents the standard, and red or green represents the text sample.
Genes 14 00559 g002
Figure 3. Transcriptome analysis and qRT-PCR validation of the candidate DEGs involved in anthocyanin biosynthesis in Lilium ′Siberia′. (A) A heatmap showing the expression patterns of the DEGs involved in anthocyanin biosynthesis. The heatmap was generated based on the normalized log2 (FoldChange) > 1 values for each DEG. Green and red scales represent relatively low or high expression, respectively. (B) Schematic diagram of metabolic pathway for cyanidin. (C) Volcano plot of the total differentially expressed genes in both E and F groups using the threshold of p < 0.05 and |log2FoldChange| > 1. The x-axis represents the log2 (fold change) values for gene expression, and the y-axis represents the −log10 (p value). (D) The qRT–PCR validation of the DEGs involved in anthocyanin biosynthesis. Expression data are the represent as the mean values of three biological replicates ± SD. Yellow color represents significant up-regulation of genes, while blue color represents insignificant changes.
Figure 3. Transcriptome analysis and qRT-PCR validation of the candidate DEGs involved in anthocyanin biosynthesis in Lilium ′Siberia′. (A) A heatmap showing the expression patterns of the DEGs involved in anthocyanin biosynthesis. The heatmap was generated based on the normalized log2 (FoldChange) > 1 values for each DEG. Green and red scales represent relatively low or high expression, respectively. (B) Schematic diagram of metabolic pathway for cyanidin. (C) Volcano plot of the total differentially expressed genes in both E and F groups using the threshold of p < 0.05 and |log2FoldChange| > 1. The x-axis represents the log2 (fold change) values for gene expression, and the y-axis represents the −log10 (p value). (D) The qRT–PCR validation of the DEGs involved in anthocyanin biosynthesis. Expression data are the represent as the mean values of three biological replicates ± SD. Yellow color represents significant up-regulation of genes, while blue color represents insignificant changes.
Genes 14 00559 g003
Figure 4. Genetic expressions of MYB genes and cluster analysis in transcriptome data. (A) The numbers and families of top 28 TFs enriched in flower bud. (B) The heatmap analysis of MYB DEGs through transcriptome data. The heatmap was generated based on the normalized log2 (FoldChange) > 1 values for each DEG. The color scale at the right represents gene expression values (the red corresponds to genes with high expression and the green corresponds to genes with low expression). (C) The cluster analysis map of R2R3-MYB proteins. The cluster analysis tree was completed using MEGA 5.0 software with the Arabidopsis MYB genes as reference. The aim gene of LhMYB114 is highlighted in red color.
Figure 4. Genetic expressions of MYB genes and cluster analysis in transcriptome data. (A) The numbers and families of top 28 TFs enriched in flower bud. (B) The heatmap analysis of MYB DEGs through transcriptome data. The heatmap was generated based on the normalized log2 (FoldChange) > 1 values for each DEG. The color scale at the right represents gene expression values (the red corresponds to genes with high expression and the green corresponds to genes with low expression). (C) The cluster analysis map of R2R3-MYB proteins. The cluster analysis tree was completed using MEGA 5.0 software with the Arabidopsis MYB genes as reference. The aim gene of LhMYB114 is highlighted in red color.
Genes 14 00559 g004
Figure 5. Amino acid sequence alignment of LhMYB114. (A) The phylogenetic analysis of LhMYB114 with homologous genes related to anthocyanin biosynthesis. (B) The important domain analysis of LhMYB114. Color represents amino acid sequence is high similarity, including pink, blue, cyan.
Figure 5. Amino acid sequence alignment of LhMYB114. (A) The phylogenetic analysis of LhMYB114 with homologous genes related to anthocyanin biosynthesis. (B) The important domain analysis of LhMYB114. Color represents amino acid sequence is high similarity, including pink, blue, cyan.
Genes 14 00559 g005
Figure 6. The qRT–PCR analysis of 3 key genes responsible for anthocyanin metabolism in flower buds of lily at different development stages. (A) The expression of LhMYB114 genes at different stages. (B) The expression of LhDFR genes at different stages. (C) The expression of LhANS-rr1 genes at different stages. The values are means ± SDs (n = 6). The data were normalized to a value of 1 for the A group. Significant differences were determined by Student′s t test (** p < 0.01; *** p < 0.001).
Figure 6. The qRT–PCR analysis of 3 key genes responsible for anthocyanin metabolism in flower buds of lily at different development stages. (A) The expression of LhMYB114 genes at different stages. (B) The expression of LhDFR genes at different stages. (C) The expression of LhANS-rr1 genes at different stages. The values are means ± SDs (n = 6). The data were normalized to a value of 1 for the A group. Significant differences were determined by Student′s t test (** p < 0.01; *** p < 0.001).
Genes 14 00559 g006
Figure 7. Transient silencing phenotypes among LhMYB114, LhANS-rr1, and LhDFR in Lilium flower buds by virus-induced gene silencing (VIGS). (A) Photographs of Lilium flower buds from TRV-infected plants. (B) RT-PCR detection of RNA1 and RNA2 of TRV1 and TRV2 in Lilium.
Figure 7. Transient silencing phenotypes among LhMYB114, LhANS-rr1, and LhDFR in Lilium flower buds by virus-induced gene silencing (VIGS). (A) Photographs of Lilium flower buds from TRV-infected plants. (B) RT-PCR detection of RNA1 and RNA2 of TRV1 and TRV2 in Lilium.
Genes 14 00559 g007
Figure 8. Relative expression levels of three key genes under different VIGS treatments. (A) Down-regulation of LhANS-rr1 and LhDFR expression in LhMYB114-silenced bud. (B) Down-regulation of LhMYB114 and LhDFR expression in LhANS-rr1-silenced bud. (C) Down-regulation of LhANS-rr1 and LhMYB114 expression in LhDFR-silenced bud. (D) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhANS-rr1-silenced bud. (E) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhDFR-silenced bud. (F) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhDFR+LhANS-rr1-silenced bud. (G) The six main genes (C4H, 4CL, CHS, CHI, F3H, and UFGT) expressed in TRV-LhMYB114, TRV-LhANS-rr1, TRV-LhDFR, and empty vector. The values are means ± SDs (n = 3). Lowercase letters (a,b) indicate statistically significant differences at p < 0.05.
Figure 8. Relative expression levels of three key genes under different VIGS treatments. (A) Down-regulation of LhANS-rr1 and LhDFR expression in LhMYB114-silenced bud. (B) Down-regulation of LhMYB114 and LhDFR expression in LhANS-rr1-silenced bud. (C) Down-regulation of LhANS-rr1 and LhMYB114 expression in LhDFR-silenced bud. (D) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhANS-rr1-silenced bud. (E) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhDFR-silenced bud. (F) Down-regulation of LhMYB114, LhANS-rr1 and LhDFR expression in LhMYB114+LhDFR+LhANS-rr1-silenced bud. (G) The six main genes (C4H, 4CL, CHS, CHI, F3H, and UFGT) expressed in TRV-LhMYB114, TRV-LhANS-rr1, TRV-LhDFR, and empty vector. The values are means ± SDs (n = 3). Lowercase letters (a,b) indicate statistically significant differences at p < 0.05.
Genes 14 00559 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fang, S.; Lin, M.; Ali, M.M.; Zheng, Y.; Yi, X.; Wang, S.; Chen, F.; Lin, Z. LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium ‘Siberia’. Genes 2023, 14, 559. https://doi.org/10.3390/genes14030559

AMA Style

Fang S, Lin M, Ali MM, Zheng Y, Yi X, Wang S, Chen F, Lin Z. LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium ‘Siberia’. Genes. 2023; 14(3):559. https://doi.org/10.3390/genes14030559

Chicago/Turabian Style

Fang, Shaozhong, Mi Lin, Muhammad Moaaz Ali, Yiping Zheng, Xiaoyan Yi, Shaojuan Wang, Faxing Chen, and Zhimin Lin. 2023. "LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium ‘Siberia’" Genes 14, no. 3: 559. https://doi.org/10.3390/genes14030559

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop