Dihydroflavonol 4-reductase genes encode enzymes with contrasting substrate specificity and show divergent gene expression profiles in Fragaria species

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

FvDFR2 rather than FvDFR1 play key roles for anthocyanin synthesis in strawberry petioles

The accumulation of anthocyanins can be found in both the fruit and petioles of strawberries, but the fruit appears red while the petioles appear purple-red. Additionally, in the white-fruited diploid strawberries, the petioles can accumulate anthocyanins normally, suggesting a different synthesis pattern between the petioles and fruits. We screened the EMS mutagenized population of a red-fruited diploid strawberry 'Ruegen' and discovered a mutant which showed no anthocyanin accumulation in the petioles but normal accumulation in the fruit. Through BSA sequencing and allelic test, it was found that a mutation in FvDFR2 was responsible for this phenotype. Furthermore, the complex formed by the interaction between the petiole-specific FvMYB10L and FvTT8 only binds the promoter of FvDFR2 but not FvDFR1, resulting in the expression of only FvDFR2 in the petiole. FvDFR2 can catalyze the conversion of DHQ and eventually the formation of cyanidin and peonidin, giving the petiole a purplish-red color. In the fruit, however, both FvDFR1 and FvDFR2 can be expressed, which can mediate the synthesis of cyanidin and pelargonidin. Our study clearly reveals different regulation of FvDFR1 and FvDFR2 in mediating anthocyanin synthesis in petioles and fruits.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Petunia dihydroflavonol 4-reductase is only a few amino acids away from producing orange pelargonidin-based anthocyanins

Anthocyanins are responsible for the color spectrum of both ornamental and natural flowers. However, not all plant species produce all colors. For example, roses are not blue because they do not naturally possess a hydroxylase that opens the pathway for delphinidin and its derivatives. It is more intriguing why some plants do not carry orange or scarlet red flowers with anthocyanins based on pelargonidin, because the precursor for these anthocyanins should be available if anthocyanins are made at all. The key to this is the substrate specificity of dihydroflavonol 4-reductase (DFR), an enzyme located at the branch point between flavonols and anthocyanins. The most common example is petunia, which does not bear orange flowers unless the enzyme is complemented by biotechnology. We changed a few amino acids in the active site of the enzyme and showed that the mutated petunia DFR started to favor dihydrokaempferol, the precursor to orange pelargonidin, in vitro. When transferred to petunia, it produced an orange hue and dramatically more pelargonidin-based anthocyanins in the flowers.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

Biochemical Characterization of Black and Green Mutant Elderberry during Fruit Ripening

The content of sugars, organic acids, phenolic compounds and selected enzyme activities in the anthocyanin pathway were analyzed in NIGRA (Sambucus nigra var. nigra-black fruits) and VIRIDIS (S. nigra var. viridis-green fruits) fruits over four stages of ripening. The share of glucose and fructose in green fruits was higher than in colored fruits, and the sugar content increased significantly until the third developmental stage. Ripe NIGRA berries had 47% flavonol glycosides, 34% anthocyanins, 3% hydroxycinnamic acids and 14% flavanols, whereas the major phenolic group in the VIRIDIS fruits, making up 88% of the total analyzed polyphenols, was flavonols. NIGRA fruits were rich in anthocyanins (6020 µg g-1 FW), showing strong activation of the late anthocyanin pathway (dihydroflavonol 4-reductase, anthocyanidin synthase). In both color types, phenylalanine ammonia lyase and chalcone synthase/chalcone isomerase activities were highest in the first stage and decreased during ripening. In VIRIDIS fruit, no anthocyanins and only one flavanol (procyanidin dimer) were found. This was most likely caused by a lack of induction of the late anthocyanin pathway in the last period of fruit ripening. The VIRIDIS genotype may be useful in studying the regulatory structures of anthocyanin biosynthesis and the contribution of distinct flavonoid classes to the health benefits of elderberries.

Insights into the catalytic and regulatory mechanisms of dihydroflavonol 4-reductase, a key enzyme of anthocyanin synthesis in Zanthoxylum bungeanum

Accumulation of anthocyanins largely determines the fruit color, and dihydroflavonol 4-reductase (DFR) is a key enzyme involved in the formation of anthocyanins. However, the catalytic and regulatory mechanisms of DFR are unclear. In this study, the gene encoding DFR from Zanthoxylum bungeanum Maxim. was cloned and ZbDFR was analyzed in detail. The ZbDFR accepted dihydrokaempferol, dihydroquercetin and dihydromyricetin as substrates. Flavonols such as myricetin, quercetin and kaempferol significantly inhibited the activity of ZbDFR, while quercitrin and isoquercitrin slightly increased the activity. Quercetin was a competitive inhibitor at low concentrations, and it had a combined effect of competitive and noncompetitive inhibition at high concentrations, which was consistent with ZbDFR having two inhibitor binding sites. In addition, the content of different types of flavonoids in Z. bungeanum peel at green, semi-red and red stage was analyzed, and the in vivo results could be explained by the regulation of ZbDFR activity in vitro. Site-directed mutagenesis combined with enzyme activity experiments showed that Ser128, Tyr163, Phe164 and Lys167 are the key catalytic amino acid residues. The Ser128, Tyr163 and Lys167 were crucial for the hydrogen transfer reaction, and mutation of these amino acids resulted in the loss of all or most of the activity. Phe164 was found to be important for the regulation of ZbDFR by flavonols. Accordingly, ZbDFR is a node at which flavonoids regulate the synthesis of anthocyanins and proanthocyanins.

Functional Analysis of Genes GlaDFR1 and GlaDFR2 Encoding Dihydroflavonol 4-Reductase (DFR) in Gentiana lutea L. Var. Aurantiaca (M. Laínz) M. Laínz

Anthocyanins are important pigments for flower color, determining the ornamental and economic values of horticultural plants. As a key enzyme in the biosynthesis of anthocyanidins, dihydroflavonol 4-reductase (DFR) catalyzes the reduction of dihydroflavonols to generate the precursors for anthocyanidins (i.e., leucoanthocyanidins) and anthocyanins. To investigate the functions of DFRs in plants, we cloned the GlaDFR1 and GlaDFR2 genes from the petals of Gentiana lutea var. aurantiaca and transformed both genes into Nicotiana tabacum by Agrobacterium-mediated leaf disc method. We further investigated the molecular and phenotypic characteristics of T1 generation transgenic tobacco plants selected based on the hygromycin resistance and verified by both PCR and semiquantitative real-time PCR analyses. The phenotypic segregation was observed in the flower color of the transgenic tobacco plants, showing petals darker than those in the wild-type (WT) plants. Results of high-performance liquid chromatography (HPLC) analysis showed that the contents of gentiocyanin derivatives were decreased in the petals of transgenic plants in comparison to those of WT plants. Ours results revealed the molecular functions of GlaDFR1 and GlaDFR2 in the formation of coloration, providing solid theoretical foundation and candidate genes for further genetic improvement in flower color of plants.

Exploring Dihydroflavonol-4-Reductase Reactivity and Selectivity by QM/MM-MD Simulations

Flavonoids are secondary metabolites ubiquitously found in plants. Their antioxidant properties make them highly interesting natural compounds for use in pharmacology. Therefore, unravelling the mechanisms of flavonoid biosynthesis is an important challenge. Among all the enzymes involved in this biosynthetic pathway, dihydroflavonol-4-reductase (DFR) plays a key role in the production of anthocyanins and proanthocyanidins. Here, we provide new information on the mechanism of action of this enzyme by using QM/MM-MD simulations applied to both dihydroquercetin (DHQ) and dihydrokaempferol (DHK) substrates. The consideration of these very similar compounds shed light on the major role played by the enzyme on the stabilization of the transition state but also on the activation of the substrate before the reaction through near-attack conformer effects.

Transcriptome sequencing and metabolomics analyses provide insights into the flavonoid biosynthesis in Torreya grandis kernels

The kernel of Torreya grandis (T. grandis) is a rare nut with a variety of bioactive compounds. Flavonoids are a very important class of bioactive compounds with high antioxidant activity in T. grandis kernels. However, the flavonoid compositions which mainly contribute to antioxidant capacity and the molecular basis of flavonoid biosynthesis in T. grandis remain unclear. Here, transcriptome sequencing and metabolomics analysis for kernels were performed. In total, 124 flavonoids were identified. Among them, 9 flavonoids were highly correlated with antioxidant activity. Furthermore, unigenes encoding CHS, DFR and ANS showed significant correlation with the 9 flavonoids. Transient overexpression of TgDFR1 in tobacco leaves resulted in increased antioxidant activity. Moreover, several transcription factors from MYB, bHLH and bZIP families were identified by co-expression assay, suggesting that they may regulate flavonoid biosynthesis. Our findings provide a molecular basis and new insights into the flavonoid biosynthesis in T. grandis kernels.

Color variation in young and senescent leaves of Formosan sweet gum (Liquidambar formosana) by the gene regulation of anthocyanidin biosynthesis

In certain plants, leaf coloration occurs in young and senescent leaves; however, it is unclear whether these two developmental stages are controlled by the same regulatory mechanisms. Formosan sweet gum (Liquidambar formosana Hance) is a subtropical deciduous tree species that possesses attractive autumnal leaf coloration. The color of young leaves is closer to purplish red, while senescent leaves are more orange-red to dark red. It was confirmed that delphinidin and cyanidin are the two anthocyanidins that contribute to the color of Formosan sweet gum leaves, and the content of different anthocyanins influences the appearance of color. To elucidate the regulation of anthocyanidin biosynthesis, recombinant DIHYDROFLAVONOL-4-REDUCTASEs (LfDFR1 and LfDFR2) (EC 1.1.1.234) were produced, and their substrate acceptability was investigated both in vitro and in planta. The functions of flavanones and dihydroflavonols modification by FLAVONOID 3' HYDROXYLASE (LfF3'H1) (EC 1.14.14.82) and FLAVONOID 3'5' HYDROXYLASE (LfF3'5'H) (EC 1.14.14.81) were verified using a transient overexpression experiment in Nicotiana benthamiana. The results showed that LfMYB5 induced LfF3'5'H and LfMYB123 induced both LfF3'H1 and LfDFR1 in spring when the leaves were expanding, whereas LfMYB113 induced LfF3'H1, LfDFR1, and LfDFR2 in late autumn to winter when the leaves were undergoing leaf senescence. In conclusion, the color variation of Formosan sweet gum in young and senescent leaves was attributed to the composition of anthocyanidins through the transcriptional regulation of LfF3'H1 and LfF3'5'H by LfMYB5, LfMYB113, and LfMYB123.

Transcriptome profiling provides insights into the fruit color development of wild Lycium ruthenicum Murr. from Qinghai-Tibet Plateau

Lycium ruthenicum Murr. is an important ecological and economic species in the Qaidam Basin of Qinghai-Tibet Plateau. Its black fruits (BF) are rich in anthocyanins, which have health-promoting properties for humans and thus provide nutritional benefits for this plant. Although the fruit quality of natural white fruit (WF) is affected by the disappearance of pigmentation in phenotypes, this phenomenon provides an opportunity to unravel the complex color metabolic networks. In this study, anthocyanin profiling confirmed that WF was formed due to anthocyanin loss. Transcriptome analysis of BF and WF revealed 101,466 unigenes, 261 of which were identified as the putative homologs of color-related genes in other species. Genes encoding the enzymes involved in flavonoid biosynthesis were also identified systematically. The structural gene expression levels of chalcone synthase (CHS), chalcone isomerase (CHI), flavonoid 3'5'-hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), and anthocyanidin 3-O-glucosyltransferase (UFGT) were highly similar and significantly positively correlated with anthocyanin accumulation rate in BF. In particular, F3'5'H, UFGT, ANS, and DFR expression levels in BF were 2391, 119, 96, and 85 times higher than those in WF at S3 (35 days after anthesis), respectively. This result strongly suggests that the low expression of these genes in WF is responsible for the anthocyanin loss. Meanwhile, the expression patterns of the anthocyanin regulatory genes were also investigated by qRT-PCR. Mass sequencing data were obtained and annotated by deep sequencing and provided a platform for future function and molecular biological research on L. ruthenicum Murr.

Anthocyanin Genes Involved in the Flower Coloration Mechanisms of Cymbidium kanran

The Orchidaceae, otherwise known as orchids, is one of the largest plant families and is renowned for its spectacular flowers and ecological adaptations. Various polymorphisms of orchid flower colour can attract pollinators and be recognised as valuable horticultural ornamentals. As one of the longest historic cultured orchids, Cymbidium kanran has been domesticated for more than 2,500 years and is an ideal species to study coloration mechanisms because of plentiful variations in floral coloration and abundant traditional varieties. In this study, we used two distinct colour-type flowers of C. kanran as experimental materials to elucidate the mechanism of flower coloration. High-performance liquid chromatography (HPLC) analysis revealed that anthocyanins in purple-red-type flowers include three types of anthocyanidin aglycones, peonidin, malvidin, and cyanidin, whereas anthocyanins are lacking in white-type flowers. Through comparative transcriptome sequencing, 102 candidate genes were identified as putative homologues of colour-related genes. Based on comprehensive correlation analysis between colour-related compounds and gene expression profiles, four candidates from 102 captured genes showed a positive correlation with anthocyanidin biosynthesis. Furthermore, transient expression of CkCHS-1, CkDFR, and CkANS by particle bombardment confirmed that recovery of their expression completed the anthocyanin pathway and produced anthocyanin compounds in white-type flowers. Collectively, this study provided a comprehensive transcriptomic dataset for Cymbidium, which significantly facilitate our understanding of the molecular mechanisms of regulating floral pigment accumulation in orchids.

Allelic Variation of MYB10 Is the Major Force Controlling Natural Variation in Skin and Flesh Color in Strawberry (Fragaria spp.) Fruit

The fruits of diploid and octoploid strawberry (Fragaria spp) show substantial natural variation in color due to distinct anthocyanin accumulation and distribution patterns. Anthocyanin biosynthesis is controlled by a clade of R2R3 MYB transcription factors, among which MYB10 is the main activator in strawberry fruit. Here, we show that mutations in MYB10 cause most of the variation in anthocyanin accumulation and distribution observed in diploid woodland strawberry (F. vesca) and octoploid cultivated strawberry (F ×ananassa). Using a mapping-by-sequencing approach, we identified a gypsy-transposon in MYB10 that truncates the protein and knocks out anthocyanin biosynthesis in a white-fruited F. vesca ecotype. Two additional loss-of-function mutations in MYB10 were identified among geographically diverse white-fruited F. vesca ecotypes. Genetic and transcriptomic analyses of octoploid Fragaria spp revealed that FaMYB10-2, one of three MYB10 homoeologs identified, regulates anthocyanin biosynthesis in developing fruit. Furthermore, independent mutations in MYB10-2 are the underlying cause of natural variation in fruit skin and flesh color in octoploid strawberry. We identified a CACTA-like transposon (FaEnSpm-2) insertion in the MYB10-2 promoter of red-fleshed accessions that was associated with enhanced expression. Our findings suggest that cis-regulatory elements in FaEnSpm-2 are responsible for enhanced MYB10-2 expression and anthocyanin biosynthesis in strawberry fruit flesh.

Genetic analysis of phenylpropanoids and antioxidant capacity in strawberry fruit reveals mQTL hotspots and candidate genes

Phenylpropanoids are a large class of plant secondary metabolites, which play essential roles in human health mainly associated with their antioxidant activity. Strawberry (Fragaria × ananassa) is a rich source of phytonutrients, including phenylpropanoids, which have been shown to have beneficial effects on human health. In this study, using the F. × ananassa '232' × '1392' F₁ segregating population, we analyzed the genetic control of individual phenylpropanoid metabolites, total polyphenol content (TPC) and antioxidant capacity (TEAC) in strawberry fruit over two seasons. We have identified a total of 7, 9, and 309 quantitative trait loci (QTL) for TPC, TEAC and for 77 polar secondary metabolites, respectively. Hotspots of stable QTL for health-related antioxidant compounds were detected on linkage groups LG IV-3, LG V-2 and V-4, and LG VI-1 and VI-2, where associated markers represent useful targets for marker-assisted selection of new varieties with increased levels of antioxidant secondary compounds. Moreover, differential expression of candidate genes for major and stable mQTLs was studied in fruits of contrasting lines in important flavonoids. Our results indicate that higher expression of FaF3'H, which encodes the flavonoid 3'-hydroxylase, is associated with increased content of these important flavonoids.

Genome-wide interologous interactome map (TeaGPIN) of Camellia sinensis

Tea, prepared from the young leaves of Camellia sinensis, is a non-alcoholic beverage globally consumed due to its antioxidant properties, strong taste and aroma. Although, the genomic data of this medicinally and commercially important plant is available, studies related to its sub-cellular interactomic maps are less explored. In this work, we propose a genome-wide interologous protein-protein interaction (PPI) network of tea, termed as TeaGPIN, consisting of 12,033 nodes and 216,107 interactions, developed using draft genome of tea and known PPIs exhaustively collected from 49 template plants. TeaGPIN interactions are prioritized using domain-domain interactions along with the interolog information. A high-confidence TeaGPIN consisting of 5983 nodes and 58,867 edges is reported and its interactions are further evaluated using protein co-localization similarities. Based on three network centralities (degree, betweenness and eigenvector), 1302 key proteins are reported in tea to have p-value <0.01 by comparing the TeaGPIN with 10,000 realizations of Erdős-Rényi and Barabási-Albert based corresponding random network models. Functional content of TeaGPIN is assessed using KEGG and GO annotations and its modular architecture is explored. Network based characterization is carried-out on the transcription factors, and proteins involved flavonoid biosynthesis and photosynthesis pathways to find novel candidates involved in various regulatory processes. We believe the proposed TeaGPIN will impart useful insights in understanding various mechanisms related to growth and development as well as defence against biotic and abiotic perturbations.

Automatic Identification of Players in the Flavonoid Biosynthesis with Application on the Biomedicinal Plant Croton tiglium

The flavonoid biosynthesis is a well-characterised model system for specialised metabolism and transcriptional regulation in plants. Flavonoids have numerous biological functions such as UV protection and pollinator attraction, but also biotechnological potential. Here, we present Knowledge-based Identification of Pathway Enzymes (KIPEs) as an automatic approach for the identification of players in the flavonoid biosynthesis. KIPEs combines comprehensive sequence similarity analyses with the inspection of functionally relevant amino acid residues and domains in subjected peptide sequences. Comprehensive sequence sets of flavonoid biosynthesis enzymes and knowledge about functionally relevant amino acids were collected. As a proof of concept, KIPEs was applied to investigate the flavonoid biosynthesis of the medicinal plant Croton tiglium on the basis of a transcriptome assembly. Enzyme candidates for all steps in the biosynthesis network were identified and matched to previous reports of corresponding metabolites in Croton species.

Arg-type dihydroflavonol 4-reductase genes from the fern Dryopteris erythrosora play important roles in the biosynthesis of anthocyanins

Dihydroflavonol 4-reductase (DFR), a key enzyme involved in the biosynthesis of anthocyanins, has been cloned from various species. However, little research has been conducted on this enzyme in ferns, which occupy a unique evolutionary position. In this study, we isolated two novel DFR genes from the fern Dryopteris erythrosora. In vitro enzymatic analysis revealed that DeDFR1 and DeDFR2 enzymes can catalyze dihydrokaempferol and dihydroquercetin but cannot catalyze dihydromyricetin. Amino acid sequence analysis showed that DeDFR1 and DeDFR2 have an arginine at the same substrate-specificity-determining site as that in the ferns Salvinia cucullata and Azolla filiculoides. Thus, we speculate that the Arg-type DFR is a new DFR functional type. To further verify the substrate preferences of the Arg-type DFR, an amino acid substitution assay was conducted. When N133 was mutated to R133, Arabidopsis DFR protein completely lost its catalytic activity for dihydromyricetin, as observed for DeDFR1 and DeDFR2. Additionally, heterologous expression of DeDFR2 in the Arabidopsis tt3-1 mutant resulted in increasing anthocyanin accumulation. In summary, DeDFR1 and DeDFR2 are considered to be a new type of DFR with unique structures and functions. The discovery of the Arg-type DFR provides new insights into the anthocyanin biosynthesis pathway in ferns.

Comprehensive analysis of putative dihydroflavonol 4-reductase gene family in tea plant

One identified dihydroflavonol 4-reductases (DFR) encoding gene (named as CsDFRa herein) and five putative DFRs (named as CsDFRb1, CsDFRb2, CsDFRb3, CsDFRc and CsDFRd) in tea (Camellia sinensis) have been widely discussed in recent papers concerning multi-omics data. However, except for CsDFRa, their function and biochemical characteristics are not clear. This study aims to compare all putative CsDFRs and preliminarily evaluate their function. We investigated the sequences of genes (coding and promoter regions) and predicted structures of proteins encoded, and determined the activities of heterologously expressed CsDFRs under various conditions. The results showed that the sequences of five putative CsDFRs were quite different from CsDFRa, and had lower expression levels as well. The five putative CsDFRs could not catalyze three dihydroflavonol substrates. The functional CsDFRa had the strongest affinity with dihydroquercetin, and performed best at pH around 7 and 35°C but was not stable at lower pHs or higher temperatures. Single amino acid mutation at position 141 modified the preference of CsDFRa for dihydroquercetin and dihydromyricetin, and also weakened its stability. These data suggest that only CsDFRa works in the pathway for generating anthocyanidins and catechins. This study provides new insights into the function of CsDFRs and may assist to develop new strategies to manipulate the composition of tea flavonoids in the future.

Cloning and Functional Characterization of Dihydroflavonol 4-Reductase Gene Involved in Anthocyanidin Biosynthesis of Grape Hyacinth

Grape hyacinth (Muscari spp.) is a popular ornamental plant with bulbous flowers noted for their rich blue color. Muscari species have been thought to accumulate delphinidin and cyanidin rather than pelargonidin-type anthocyanins because their dihydroflavonol 4-reductase (DFR) does not efficiently reduce dihydrokaempferol. In our study, we clone a novel DFR gene from blue flowers of Muscari. aucheri. Quantitative real-time PCR (qRT-PCR) and anthocyanin analysis showed that the expression pattern of MaDFR had strong correlations with the accumulation of delphinidin, relatively weak correlations with cyanidin, and no correations with pelargonidin. However, in vitro enzymatic analysis revealed that the MaDFR enzyme can reduce all the three types of dihydroflavonols (dihydrokaempferol, dihydroquercetin, and dihydromyricetin), although it most preferred dihydromyricetin as a substrate to produce leucodelphinidin, the precursor of blue-hued delphinidin. This indicated that there may be other functional genes responsible for the loss of red pelargonidin-based pigments in Muscari. To further verify the substrate-specific selection domains of MaDFR, an assay of amino acid substitutions was conducted. The activity of MaDFR was not affected whenever the N135 or E146 site was mutated. However, when both of them were mutated, the catalytic activity of MaDFR was lost completely. The results suggest that both the N135 and E146 sites are essential for the activity of MaDFR. Additionally, the heterologous expression of MaDFR in tobacco (Nicotiana tabacum) resulted in increasing anthocyanin accumulation, leading to a darker flower color, which suggested that MaDFR was involved in color development in flowers. In summary, MaDFR has a high preference for dihydromyricetin, and it could be a powerful candidate gene for genetic engineering for blue flower colour modification. Our results also make a valuable contribution to understanding the basis of color variation in the genus Muscari.

Transcriptional regulation of abscisic acid biosynthesis and signal transduction, and anthocyanin biosynthesis in 'Bluecrop' highbush blueberry fruit during ripening

Highbush blueberry (Vaccinium corymbosum) fruit accumulate high levels of anthocyanins during ripening, which might be controlled by abscisic acid (ABA), a signal molecule in non-climacteric fruits. For an integrated view of the ripening process from ABA to anthocyanin biosynthesis, we analyzed the transcriptomes of ‘Bluecrop’ highbush blueberry fruit using RNA-Seq at three ripening stages, categorized based on fruit skin coloration: pale green at ca. 30 days after full bloom (DAFB), reddish purple at ca. 40 DAFB, and dark purple at ca. 50 DAFB. Mapping the trimmed reads against the reference sequences yielded 25,766 transcripts. Of these, 143 transcripts were annotated to encode five ABA biosynthesis enzymes, four ABA signal transduction regulators, four ABA-responsive transcription factors, and 12 anthocyanin biosynthesis enzymes. The analysis of differentially expressed genes between the ripening stages revealed that 11 transcripts, including those encoding nine-cis-epoxycarotenoid dioxygenase, SQUAMOSA-class MADS box transcription factor, and flavonoid 3′,5′-hydroxylase, were significantly up-regulated throughout the entire ripening stages. In fruit treated with 1 g L⁻¹ ABA, at least nine transcripts of these 11 transcripts as well as one transcript encoding flavonoid 3′-hydroxylase were up-regulated, presumably promoting anthocyanin accumulation and fruit skin coloration. These results will provide fundamental information demonstrating that ABA biosynthesis and signal transduction, and anthocyanin biosynthesis are closely associated with anthocyanin accumulation and skin coloration in highbush blueberry fruit during ripening.

Biochemical and Functional Characterization of Anthocyanidin Reductase (ANR) from Mangifera indica L

Mango (Mangifera indica L.) is abundant in proanthocyanidins (PAs) that are important for human health and plant response to abiotic stresses. However, the molecular mechanisms involved in PA biosynthesis still need to be elucidated. Anthocyanidin reductase (ANR) catalyzes a key step in PA biosynthesis. In this study, three ANR cDNAs (MiANR1-1,1-2,1-3) were isolated from mango, and expressed in Escherichia coli. In vitro enzyme assay showed MiANR proteins convert cyanidin to their corresponding flavan-3-ols, such as (−)-catechin and (−)-epicatechin. Despite high amino acid similarity, the recombinant ANR proteins exhibited differences in enzyme kinetics and cosubstrate preference. MiANR1-2 and MiANR1-3 have the same optimum pH of 4.0 in citrate buffer, while the optimum pH for MiANR1-1 is pH 3.0 in phosphate buffer. MiANR1-1 does not use either NADPH or NADH as co-substrate while MiANR1-2/1-3 use only NADPH as co-substrate. MiANR1-2 has the highest Km and Vmax for cyanidin, followed by MiANR1-3 and MiANR1-1. The overexpression of MiANRs in ban mutant reconstructed the biosynthetic pathway of PAs in the seed coat. These data demonstrate MiANRs can form the ANR pathway, leading to the formation of two types of isomeric flavan-3-ols and PAs in mango.

The rare orange-red colored Euphorbia pulcherrima cultivar 'Harvest Orange' shows a nonsense mutation in a flavonoid 3'-hydroxylase allele expressed in the bracts

BACKGROUND

Commercially available poinsettia (Euphorbia pulcherrima) varieties prevalently accumulate cyanidin derivatives and show intense red coloration. Orange-red bract color is less common. We investigated four cultivars displaying four different red hues with respect to selected enzymes and genes of the anthocyanin pathway, putatively determining the color hue.

RESULTS

Red hues correlated with anthocyanin composition and concentration and showed common dark red coloration in cultivars 'Christmas Beauty' and 'Christmas Feeling' where cyanidin derivatives were prevalent. In contrast, orange-red bract color is based on the prevalent presence of pelargonidin derivatives that comprised 85% of the total anthocyanin content in cv. 'Premium Red' and 96% in cv. 'Harvest Orange' (synonym: 'Orange Spice'). cDNA clones of flavonoid 3'-hydroxylase (F3'H) and dihydroflavonol 4-reductase (DFR) were isolated from the four varieties, and functional activity and substrate specificity of the corresponding recombinant enzymes were studied. Kinetic studies demonstrated that poinsettia DFRs prefer dihydromyricetin and dihydroquercetin over dihydrokaempferol, and thus, favor the formation of cyanidin over pelargonidin. Whereas the F3'H cDNA clones of cultivars 'Christmas Beauty', 'Christmas Feeling', and 'Premium Red' encoded functionally active enzymes, the F3'H cDNA clone of cv. 'Harvest Orange' contained an insertion of 28 bases, which is partly a duplication of 20 bases found close to the insertion site. This causes a frameshift mutation with a premature stop codon after nucleotide 132 and, therefore, a non-functional enzyme. Heterozygosity of the F3'H was demonstrated in this cultivar, but only the mutated allele was expressed in the bracts. No correlation between F3'H-expression and the color hue could be observed in the four species.

CONCLUSIONS

Rare orange-red poinsettia hues caused by pelargonidin based anthocyanins can be achieved by different mechanisms. F3'H is a critical step in the establishment of orange red poinsettia color. Although poinsettia DFR shows a low substrate specificity for dihydrokaempferol, sufficient precursor for pelargonidin formation is available in planta, in the absence of F3'H activity.

Transcriptomic Profiling of Fruit Development in Black Raspberry Rubus coreanus

The wild Rubus species R. coreanus, which is widely distributed in southwest China, shows great promise as a genetic resource for breeding. One of its outstanding properties is adaptation to high temperature and humidity. To facilitate its use in selection and breeding programs, we assembled de novo 179,738,287 R. coreanus reads (125 bp in length) generated by RNA sequencing from fruits at three representative developmental stages. We also used the recently released draft genome of R. occidentalis to perform reference-guided assembly. We inferred a final 95,845-transcript reference for R. coreanus. Of these genetic resources, 66,597 (69.5%) were annotated. Based on these results, we carried out a comprehensive analysis of differentially expressed genes. Flavonoid biosynthesis, phenylpropanoid biosynthesis, plant hormone signal transduction, and cutin, suberin, and wax biosynthesis pathways were significantly enriched throughout the ripening process. We identified 23 transcripts involved in the flavonoid biosynthesis pathway whose expression perfectly paralleled changes in the metabolites. Additionally, we identified 119 nucleotide-binding site leucine-rich repeat (NBS-LRR) protein-coding genes, involved in pathogen resistance, of which 74 were in the completely conserved domain. These results provide, for the first time, genome-wide genetic information for understanding developmental regulation of R. coreanus fruits. They have the potential for use in breeding through functional genetic approaches in the near future.

Great Cause-Small Effect: Undeclared Genetically Engineered Orange Petunias Harbor an Inefficient Dihydroflavonol 4-Reductase

A recall campaign for commercial, orange flowering petunia varieties in spring 2017 caused economic losses worldwide. The orange varieties were identified as undeclared genetically engineered (GE)-plants, harboring a maize dihydroflavonol 4-reductase (DFR, A₁), which was used in former scientific transgenic breeding attempts to enable formation of orange pelargonidin derivatives from the precursor dihydrokaempferol (DHK) in petunia. How and when the A₁ cDNA entered the commercial breeding process is unclear. We provide an in-depth analysis of three orange petunia varieties, released by breeders from three countries, with respect to their transgenic construct, transcriptomes, anthocyanin composition, and flavonoid metabolism at the level of selected enzymes and genes. The two possible sources of the A₁ cDNA in the undeclared GE-petunia can be discriminated by PCR. A special version of the A₁ gene, the A₁ type 2 allele, is present, which includes, at the 3'-end, an additional 144 bp segment from the non-viral transposable Cin4-1 sequence, which does not add any functional advantage with respect to DFR activity. This unequivocally points at the first scientific GE-petunia from the 1980s as the A₁ source, which is further underpinned e.g., by the presence of specific restriction sites, parts of the untranslated sequences, and the same arrangement of the building blocks of the transformation plasmid used. Surprisingly, however, the GE-petunia cannot be distinguished from native red and blue varieties by their ability to convert DHK in common in vitro enzyme assays, as DHK is an inadequate substrate for both the petunia and maize DFR. Recombinant maize DFR underpins the low DHK acceptance, and, thus, the strikingly limited suitability of the A₁ protein for a transgenic approach for breeding pelargonidin-based flower color. The effect of single amino acid mutations on the substrate specificity of DFRs is demonstrated. Expression of the A₁ gene is generally lower than the petunia DFR expression despite being under the control of the strong, constitutive p35S promoter. We show that a rare constellation in flavonoid metabolism-absence or strongly reduced activity of both flavonol synthase and B-ring hydroxylating enzymes-allows pelargonidin formation in the presence of DFRs with poor DHK acceptance.

A wild 'albino' bilberry (Vaccinium myrtillus L.) from Slovenia shows three bottlenecks in the anthocyanin pathway and significant differences in the expression of several regulatory genes compared to the common blue berry type

Relative expressions of structural genes and a number of transcription factors of the anthocyanin pathway relevant in Vaccinium species, and related key enzyme activities were compared with the composition and content of metabolites in skins of ripe fruits of wild albino and blue bilberry (Vaccinium myrtillus) found in Slovenia. Compared to the common blue type, the albino variant had a 151-fold lower total anthocyanin and a 7-fold lower total phenolic content in their berry skin, which correlated with lower gene expression of flavonoid 3-O-glycosyltransferase (FGT; 33-fold), flavanone 3-hydroxylase (FHT; 18-fold), anthocyanidin synthase (ANS; 11-fold), chalcone synthase (CHS, 7.6-fold) and MYBPA1 transcription factor (22-fold). The expression of chalcone isomerase (CHI), dihydroflavonol 4-reductase (DFR), leucoanthocyanidin reductase (LAR), anthocyanidin reductase (ANR) and MYBC2 transcription factor was reduced only by a factor of 1.5-2 in the albino berry skins, while MYBR3 and flavonoid 3',5'-hydroxylase (F3'5'H) were increased to a similar extent. Expression of the SQUAMOSA class transcription factor TDR4, in contrast, was independent of the color type and does therefore not seem to be correlated with anthocyanin formation in this variant. At the level of enzymes, significantly lower FHT and DFR activities, but not of phenylalanine ammonia-lyase (PAL) and CHS/CHI, were observed in the fruit skins of albino bilberries. A strong increase in relative hydroxycinnamic acid derivative concentrations indicates the presence of an additional bottleneck in the general phenylpropanoid pathway at a so far unknown step between PAL and CHS.

A new buckwheat dihydroflavonol 4-reductase (DFR), with a unique substrate binding structure, has altered substrate specificity

BACKGROUND

Dihydroflavonol 4-reductase (DFR) is the key enzyme committed to anthocyanin and proanthocyanidin biosynthesis in the flavonoid biosynthetic pathway. DFR proteins can catalyse mainly the three substrates (dihydrokaempferol, dihydroquercetin, and dihydromyricetin), and show different substrate preferences. Although relationships between the substrate preference and amino acids in the region responsible for substrate specificity have been investigated in several plant species, the molecular basis of the substrate preference of DFR is not yet fully understood.

RESULTS

By using degenerate primers in a PCR, we isolated two cDNA clones that encoded DFR in buckwheat (Fagopyrum esculentum). Based on sequence similarity, one cDNA clone (FeDFR1a) was identical to the FeDFR in DNA databases (DDBJ/Gen Bank/EMBL). The other cDNA clone, FeDFR2, had a similar sequence to FeDFR1a, but a different exon-intron structure. Linkage analysis in an F₂ segregating population showed that the two loci were linked. Unlike common DFR proteins in other plant species, FeDFR2 contained a valine instead of the typical asparagine at the third position and an extra glycine between sites 6 and 7 in the region that determines substrate specificity, and showed less activity against dihydrokaempferol than did FeDFR1a with an asparagine at the third position. Our 3D model suggested that the third residue and its neighbouring residues contribute to substrate specificity. FeDFR1a was expressed in all organs that we investigated, whereas FeDFR2 was preferentially expressed in roots and seeds.

CONCLUSIONS

We isolated two buckwheat cDNA clones of DFR genes. FeDFR2 has unique structural and functional features that differ from those of previously reported DFRs in other plants. The 3D model suggested that not only the amino acid at the third position but also its neighbouring residues that are involved in the formation of the substrate-binding pocket play important roles in determining substrate preferences. The unique characteristics of FeDFR2 would provide a useful tool for future studies on the substrate specificity and organ-specific expression of DFRs.

The draft genome of Ruellia speciosa (Beautiful Wild Petunia: Acanthaceae)

The genus Ruellia (Wild Petunias; Acanthaceae) is characterized by an enormous diversity of floral shapes and colours manifested among closely related species. Using Illumina platform, we reconstructed the draft genome of Ruellia speciosa, with a scaffold size of 1,021 Mb (or ∼1.02 Gb) and an N50 size of 17,908 bp, spanning ∼93% of the estimated genome (∼1.1 Gb). The draft assembly predicted 40,124 gene models and phylogenetic analyses of four key enzymes involved in anthocyanin colour production [flavanone 3-hydroxylase (F3H), flavonoid 3'-hydroxylase (F3'H), flavonoid 3',5'-hydroxylase (F3'5'H), and dihydroflavonol 4-reductase (DFR)] found that most angiosperms here sampled harboured at least one copy of F3H, F3'H, and DFR. In contrast, fewer than one-half (but including R. speciosa) harboured a copy of F3'5'H, supporting observations that blue flowers and/or fruits, which this enzyme is required for, are less common among flowering plants. Ka/Ks analyses of duplicated copies of F3'H and DFR in R. speciosa suggested purifying selection in the former but detected evidence of positive selection in the latter. The genome sequence and annotation of R. speciosa represents only one of only four families sequenced in the large and important Asterid clade of flowering plants and, as such, will facilitate extensive future research on this diverse group, particularly with respect to floral evolution.

Early metabolic and transcriptional variations in fruit of natural white-fruited Fragaria vesca genotypes

Strawberry fruits (Fragaria vesca) are valued for their sweet fruity flavor, juicy texture, and characteristic red color caused by anthocyanin pigments. To gain a deeper insight into the regulation of anthocyanin biosynthesis, we performed comparative metabolite profiling and transcriptome analyses of one red-fruited and two natural white-fruited strawberry varieties in two tissues and three ripening stages. Developing fruit of the three genotypes showed a distinctive pattern of polyphenol accumulation already in green receptacle and achenes. Global analysis of the transcriptomes revealed that the ripening process in the white-fruited varieties is already affected at an early developmental stage. Key polyphenol genes showed considerably lower transcript levels in the receptacle and achenes of both white genotypes, compared to the red genotype. The expression of the anthocyanidin glucosyltransferase gene and a glutathione S-transferase, putatively involved in the vacuolar transport of the anthocyanins, seemed to be critical for anthocyanin formation. A bHLH transcription factor is among the differentially expressed genes as well. Furthermore, genes associated with flavor formation and fruit softening appear to be coordinately regulated and seem to interact with the polyphenol biosynthesis pathway. This study provides new information about polyphenol biosynthesis regulators in strawberry, and reveals genes unknown to affect anthocyanin formation.

Transcriptomics and Metabolite Analysis Reveals the Molecular Mechanism of Anthocyanin Biosynthesis Branch Pathway in Different Senecio cruentus Cultivars

The cyanidin (Cy), pelargonidin (Pg), and delphinidin (Dp) pathways are the three major branching anthocyanin biosynthesis pathways that regulate flavonoid metabolic flux and are responsible for red, orange, and blue flower colors, respectively. Different species have evolved to develop multiple regulation mechanisms that form the branched pathways. In the current study, five Senecio cruentus cultivars with different colors were investigated. We found that the white and yellow cultivars do not accumulate anthocyanin and that the blue, pink, and carmine cultivars mainly accumulate Dp, Pg, and Cy in differing densities. Subsequent transcriptome analysis determined that there were 43 unigenes encoding anthocyanin biosynthesis genes in the blue cultivar. We also combined chemical and transcriptomic analyses to investigate the major metabolic pathways that are related to the observed differences in flower pigmentation in the series of S. cruentus. The results showed that mutations of the ScbHLH17 and ScCHI1/2 coding regions abolish anthocyanin formation in the white and the yellow cultivars; the competition of the ScF3'H1, ScF3'5'H, and ScDFR1/2 genes for naringenin determines the differences in branching metabolic flux of the Cy, Dp, and Pg pathways. Our findings provide new insights into the regulation of anthocyanin branching and also supplement gene resources (including ScF3'5 'H, ScF3'H, and ScDFRs) for flower color modification of ornamentals.

Red Anthocyanins and Yellow Carotenoids Form the Color of Orange-Flower Gentian (Gentiana lutea L. var. aurantiaca)

Flower color is an important characteristic that determines the commercial value of ornamental plants. Gentian flowers occur in a limited range of colors because this species is not widely cultivated as a cut flower. Gentiana lutea L. var. aurantiaca (abbr, aurantiaca) is characterized by its orange flowers, but the specific pigments responsible for this coloration are unknown. We therefore investigated the carotenoid and flavonoid composition of petals during flower development in the orange-flowered gentian variety of aurantiaca and the yellow-flowered variety of G. lutea L. var. lutea (abbr, lutea). We observed minor varietal differences in the concentration of carotenoids at the early and final stages, but only aurantiaca petals accumulated pelargonidin glycosides, whereas these compounds were not found in lutea petals. We cloned and sequenced the anthocyanin biosynthetic gene fragments from petals, and analyzed the expression of these genes in the petals of both varieties to determine the molecular mechanisms responsible for the differences in petal color. Comparisons of deduced amino acid sequences encoded by the isolated anthocyanin cDNA fragments indicated that chalcone synthase (CHS), chalcone isomerase (CHI), anthocyanidin synthase 1 (ANS1) and ANS2 are identical in both aurantiaca and lutea varieties whereas minor amino acid differences of the deduced flavonone 3-hydroxylase (F3H) and dihydroflavonol 4-reductase (DFR) between both varieties were observed. The aurantiaca petals expressed substantially higher levels of transcripts representing CHS, F3H, DFR, ANS and UDP-glucose:flavonoid-3-O-glucosyltransferase genes, compared to lutea petals. Pelargonidin glycoside synthesis in aurantiaca petals therefore appears to reflect the higher steady-state levels of pelargonidin synthesis transcripts. Moreover, possible changes in the substrate specificity of DFR enzymes may represent additional mechanisms for producing red pelargonidin glycosides in petals of aurantiaca. Our report describing the exclusive accumulation of pelargonidin glycosides in aurantiaca petals may facilitate the modification of gentian flower color by the production of red anthocyanins.

Functional Characterization of a Dihydroflavanol 4-Reductase from the Fiber of Upland Cotton (Gossypium hirsutum)

Dihydroflavanol 4-reductase (DFR) is a key later enzyme involved in two polyphenols' (anthocyanins and proanthocyanidins (PAs)) biosynthesis, however it is not characterized in cotton yet. In present reports, a DFR cDNA homolog (designated as GhDFR1) was cloned from developing fibers of upland cotton. Silencing GhDFR1 in cotton by virus-induced gene silencing led to significant decrease in accumulation of anthocyanins and PAs. More interestingly, based on LC-MS analysis, two PA monomers, (-)-epicatachin and (-)-epigallocatachin, remarkably decreased in content in fibers of GhDFR1-silenced plants, but two new monomers, (-)-catachin and (-)-gallocatachin were present compared to the control plants infected with empty vector. The ectopic expression of GhDFR1 in an Arabidopsis TT3 mutant allowed for reconstruction of PAs biosynthesis pathway and led to accumulation of PAs in seed coat. Taken together, these data demonstrate that GhDFR1 contributes to the biosynthesis of anthocyanins and PAs in cotton.

Transcript Quantification by RNA-Seq Reveals Differentially Expressed Genes in the Red and Yellow Fruits of Fragaria vesca

Fragaria vesca (2n = 2x = 14), the woodland strawberry, is a perennial herbaceous plant with a small sequenced genome (240 Mb). It is commonly used as a genetic model plant for the Fragaria genus and the Rosaceae family. Fruit skin color is one of the most important traits for both the commercial and esthetic value of strawberry. Anthocyanins are the most prominent pigments in strawberry that bring red, pink, white, and yellow hues to the fruits in which they accumulate. In this study, we conducted a de novo assembly of the fruit transcriptome of woodland strawberry and compared the gene expression profiles with yellow (Yellow Wonder, YW) and red (Ruegen, RG) fruits. De novo assembly yielded 75,426 unigenes, 21.3% of which were longer than 1,000 bp. Among the high-quality unique sequences, 45,387 (60.2%) had at least one significant match to an existing gene model. A total of 595 genes, representing 0.79% of total unigenes, were differentially expressed in YW and RG. Among them, 224 genes were up-regulated and 371 genes were down-regulated in the fruit of YW. Particularly, some flavonoid biosynthetic pathway genes, including C4H, CHS, CHI, F3H, DFR and ANS, as well as some transcription factors (TFs), including MYB (putative MYB86 and MYB39), WDR and MADS, were down-regulated in YW fruit, concurrent with a reduction in anthocyanin accumulation in the yellow pigment phenotype, whereas a putative transcription repressor MYB1R was up-regulated in YW fruit. The altered expression levels of the genes encoding flavonoid biosynthetic enzymes and TFs were confirmed by quantitative RT-PCR. Our study provides important insights into the molecular mechanisms underlying the yellow pigment phenotype in F. vesca.