Next Article in Journal
Peanut Genotypes with Reduced Content of Immunogenic Proteins by Breeding, Biotechnology, and Management: Prospects and Challenges
Previous Article in Journal
Morpho-Physicochemical, Nutritional Composition and Phenolic Compound Profile of Two Avocado Landraces in Different Ripening Stages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Revealing the Molecular Regulatory Mechanism of Flavonoid Accumulation in Tender Leaves of Tea Plants by Transcriptomic and Metabolomic Analyses

by
Ruiyang Shan
1,†,
Yongheng Zhang
2,†,
Xiaomei You
1,
Xiangrui Kong
1,
Yazhen Zhang
1,
Xinlei Li
1,
Lu Wang
2,
Xinchao Wang
2,* and
Changsong Chen
1,*
1
Tea Research Institute, Fujian Academy of Agricultural Science, Fujian Branch of National Center for Tea Improvement, Fuzhou 350013, China
2
Key Laboratory of Biology, Genetics and Breeding of Special Economic Animals and Plants, Ministry of Agriculture and Rural Affairs, National Center for Tea Plant Improvement, Tea Research Institute, Chinese Academy of Agricultural Sciences, Hangzhou 310008, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2025, 14(4), 625; https://doi.org/10.3390/plants14040625
Submission received: 30 December 2024 / Revised: 3 February 2025 / Accepted: 13 February 2025 / Published: 19 February 2025
(This article belongs to the Special Issue Genetic Breeding and Quality Improvement of Tea)

Abstract

:
Flavonoids are secondary metabolites that are beneficial to life activities and are mainly concentrated in buds and leaves in the form of glycosides. Flavonoid glycosides have important effects on the properties and quality of tea plants. Research has shown that the abundance of flavonoid glycosides varies greatly among different cultivars, but research on the regulatory mechanisms that cause their differential accumulation among tea plant cultivars with different leaf colors is lacking. In this study, an integrated analysis of metabolomics and transcriptomics was conducted to determine the regulatory networks regulating astringency and color-related flavonoids in tea plant cultivars with diverse leaf colors. A total of five anthocyanidins, four catechins, and nine flavonol glycosides were found to partially contribute to the differences in taste and leaf color among tea plant cultivars with diverse leaf colors. Furthermore, 15 MYB genes and 5 Dof genes were identified as potential regulators controlling the expression of eight key structural genes, resulting in differences in the accumulation of specific compounds, including epicatechin (EC), catechin (C), cyanidin, cyanidin 3-O-glucoside, pelargonidin 3-O-glucoside, and quercetin 3-O-glucoside, in tea plant cultivars with diverse leaf colors. These findings provide insights into the development and utilization of resources from tea plants with diverse leaf colors.

1. Introduction

Tea has emerged as one of the most significant and widely consumed nonalcoholic beverages globally because of its distinctive flavor profile and inherent nutritional value [1,2]. Tea is roughly divided into six categories based on the processing process, namely white tea, green tea, black tea, yellow tea, dark tea, and oolong tea, and the unique tastes and aromas of these teas are attributed to the differing levels of amino acids, flavonoids, caffeine, and other active compounds present in them [3,4]. Flavonoids are primarily concentrated in buds and young leaves and constitute between 12% and 24% of the dry weight of tea [5]. The principal flavonoid categories include isoflavones, flavanones, flavanols, flavonols, flavones, and anthocyanins [6]. Extensive research has demonstrated that flavonoids possess free radical scavenging and antioxidant properties [5,6,7,8], exerting a profound influence on plant defense mechanisms, growth, and development [9]. Specifically, they influence color formation to attract pollinators, provide UV protection, confer pathogen resistance, mediate plant–microbial signaling interactions, impact pollen fertility, and regulate plant growth [10,11]. Within the human body, the intestinal microflora hydrolyzes glycosylated flavonoids into their respective aglycones, resulting in anti-inflammatory, immunomodulatory, and potent anticancer activities and assisting in the management of cardiovascular diseases [12,13,14].
Furthermore, flavonols and catechins are pivotal flavor constituents in tea liquor, exhibiting dual sensory attributes of bitterness and astringency [15]. Under the catalysis of polyphenol oxidase (PPO), catechins undergo oxidation to form theaflavins, thearubigins, and tea brown pigments, which not only significantly influence the color of tea liquor but also exhibit much lower astringency thresholds compared to catechins [16]. Anthocyanins, a notable subclass of flavonoids, have gained increasing attention due to their roles in cancer treatment and health care [17]. Eight types of anthocyanins have been isolated from purple-leaf tea plant varieties for investigation [18]. These compounds undergo intramolecular coloration via glycosylation and exhibit diverse hues, such as the red appearance of cyanidin-3-O-glucose and cyanidin-3-O-rutin and the dark blue of delphinidin 3,5,3′-O-triglucoside in gentian [19,20]. High concentrations of delphinidin 3-O-(6-O-p-coumaroyl) galactoside and cyanidin 3-O-(6-O-p-coumaroyl) galactoside were detected in ‘Zijuan’ purple-leaf tea in contrast to the minimal levels found in ‘Yunkang 10’ green-leaf tea [21].
Flavonoids accumulate in plant organs like flowers, fruits, and leaves as glycosides, consisting of an aglycone linked to sugar units (e.g., glucose, galactose, and rhamnose). In tea, flavonol glycosides are abundant, accounting for 3–4% of the dry weight, with over 20 aglycones identified [22,23]. These can be categorized as quercetin glycosides, myricetin glycosides, and kaempferol glycosides [23]. Purple-leaf plants are rich in anthocyanins, primarily pelargonidin, cyanidin, delphinidin, peonidin, petunidin, and malvidin [24]. Tea contains anthocyanin glycosides, such as delphinidin and cyanidin derivatives, with β-D-galactoside and β-D-(6-(E)-p-coumaroyl)-galactopyranoside moieties [25,26,27]. All flavonoids originate from the phenylpropanoid pathway, synthesized from L-phenylalanine under the catalysis of a series of enzymes, such as phenylalanine am-monia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), flavonoid 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3′H), and flavonoid 3′,5′-hydroxylase (F3′5′H), which regulate the early biosynthetic steps of the flavonoid pathway [28,29,30]. Glycosylation is crucial in plants, involving flavonoid UDP-glycosyltransferase (UGT) that transfers sugar moieties to phenolic compounds [22,31,32]. Currently, a series of UGTs have been found to play crucial roles in the glycosylation process of flavonoids in tea plants. For instance, CsUGT73A17 catalyzes the glycosylation of 17 flavonoids specifically at the 7-O position [33], while CsUGT72AM1 is responsible for promoting the glucosylation of kaempferol, quercetin, myricetin, and naringenin, among others, at the 3-O position [18,34]. Additionally, CsUGT73A20 can catalyze the production of various flavonoid glycosides, including flavonoid 3-O-glucosides, flavonoid 7-O-glycosides, flavonoid 3-O-rhamnosides, and flavonoid 3,7-di-O-rhamnosides [22,35]. These studies suggest that UGTs exhibit a complex catalytic role in the glycosylation of flavonoids in tea plants, possessing a diverse range of flavonoid substrates that contribute to the abundant glycosylated flavonoids found in these plants.
In addition to the aforementioned structural genes, transcription factors (TFs) exert crucial regulatory functions in the synthesis of flavonoids. A considerable amount of evidence supports the involvement of MYB TFs in modulating flavonoid synthesis in tea plants. Specifically, CsMYB196 and CsMYB184 exhibit substantial activation of CsANR and CsANS expression, thereby orchestrating the biosynthesis of catechins, anthocyanins, and flavonols [36]. Similarly, CsMYB34 has been reported to participate in the biosynthesis of galloylated catechins [37]. Notably, CsMYB1 interacts with CsGL3 and CsWD40 to form the MYB-bHLH-WD40 (MBW) transcriptional complex, which subsequently activates the expression of CsANR and CsSCPL1A, genes involved in galloylated cis-catechin biosynthesis [38]. Furthermore, other TF types, including CsbZIP1 [39], CsWRKY12 [40], CsWRKY48 [41], MYC2 [42], and CsMADSL1 [43], regulate flavonoid synthesis by modulating the transcription of specific flavonoid biosynthetic genes. Interestingly, Dof TFs have demonstrated their capacity to regulate flavonoid synthesis in other plants [44,45]. In tea plants, Dof members are implicated in the regulation of N remobilization [46], chlorophyll metabolism [47], and responses to abiotic stress [48]. However, the role of Dof TFs in regulating flavonoid synthesis within tea plants remains unknown and warrants further elucidation.
A series of tea germplasm resources exhibiting unique leaf colors, distinct from the conventional green color, were observed. These resources are characterized by varying accumulation levels of flavonoids, leading to diversity in tea quality. Recently, numerous studies have focused on the variation in flavonoid metabolism among tea plants with different leaf colors. However, the underlying regulatory mechanisms remain elusive. To address this, our study integrated transcriptome and metabolome analyses to elucidate the potential regulatory network governing specific flavonoid accumulation in tea plants with diverse leaf colors. Special emphasis was placed on the regulatory roles of MYB and DOF transcription factors in flavonoid accumulation in tea plants of different colors. These findings provide crucial foundations for the resource utilization and molecular breeding of tea plants.

2. Results

2.1. Metabolite Profiles of Leaf Samples from Three Tea Plant Cultivars

By observing the phenotypes of shoots from three tea plant cultivars at different developmental stages, we found that the buds and one bud with two leaves of ‘Mingguan’ (CO) were green and yellowish, whereas those of ‘Lvyafoshou’ (LF) were simply green. In contrast, the buds and one bud with two leaves of ‘Zijuan’ (ZJ) were purple (Figure 1A). A widely targeted metabolite analysis was applied to analyze the differences in their metabolites. A total of 863 components were detected, including 35 amino acids and derivatives, 56 phenylpropanoids and lignans, 21 benzene and derivatives, 89 phenols, 23 nucleotides and derivatives, 138 flavonoids, 28 coumarins, 14 anthraquinones, 129 alkaloids, 26 carbohydrates and organooxygen compounds, 93 terpenoids, 40 organic acids and derivatives, 17 steroids and derivatives, 61 lipids, and 93 others (Figure 1B and Table S2). A heatmap of all 863 components reveals that they had different accumulation patterns (Figure 1C), and the PCA results indicate that the metabolite profiles of the six samples were separated (Figure 1D). A heatmap was generated to assess the total levels of 14 metabolite categories in different leaf samples of CO, LF, and ZJ. The results reveal that the total levels of anthraquinones, organic acids and derivatives, and alkaloids were greater in CO1, whereas LF1 accumulated more phenols, terpenoids, and lipids. The most abundant phenylpropanoids and lignans, benzene and derivatives, flavonoids, steroids and derivatives, and carbohydrates and organooxygen compounds were found in ZJ2, whereas the most abundant amino acids and derivatives were found in LF2 (Figure 1E). These results collectively indicate that the samples of the buds and one bud with two leaves of the three tea plant cultivars have distinct metabolite profiles.

2.2. Differentially Abundant Metabolite Analysis

Six comparison groups, namely CO1_vs_LF1, CO1_vs_ZJ1, LF1_vs_ZJ1, CO2_vs_LF2, CO2_vs_ZJ2, and LF2_vs_ZJ2, were set to identify significantly differentially accumulated metabolites (DAMs), resulting in 47, 69, 75, 80, 73, and 76 DAMs, respectively (Figure 2A and Table S3). DAM1 (consisting of CO1_vs_LF1, CO1_vs_ZJ1, and LF1_vs_ZJ1) and DAM2 (consisting of CO2_vs_LF2, CO2_vs_ZJ2, and LF2_vs_ZJ2) were subsequently defined to reflect the total DAMs of the samples of buds and one bud with two leaves, respectively. A KEGG enrichment analysis revealed significant differences in the top 10 enriched pathways between DAM1 and DAM2 (Figure 2B), indicating distinct metabolic characteristics of the samples of buds and one bud with two leaves. Notably, ‘flavonoid biosynthesis’ was significantly enriched in both DAM1 and DAM2, with ‘flavone and flavonol biosynthesis’ and ‘anthocyanin biosynthesis’, the different continuations of ‘flavonoid biosynthesis’, enriched in DAM1 and DAM2, respectively (Figure 2B). These results reveal significant differences in the accumulation of flavonoid compounds in the samples of both the buds and one bud with two leaves.

2.3. Accumulation Patterns of Anthocyanidins, Flavan-3-Ols, and Flavonol Glycosides

A total of 47 differentially accumulated flavonoids were found in all of the compared groups and clustered into seven different accumulation patterns (Figure 3). The compounds clustered into groups 1, 2, 3, and 7 presented the greatest accumulation of CO2, LF1, LF2, and ZJ2, respectively. Cluster 4 compounds, with the exception of kaempferol-3,7-O-alpha-L-dirhamnoside, were more highly accumulated in CO1 and ZJ2. The cluster 5 compounds accumulated more in ZJ1 and ZJ2, whereas the cluster 6 compounds accumulated more in CO2 and ZJ2.
Considering that anthocyanidins are important compounds that result in differences in leaf color phenotypes, flavan-3-ols alongside flavonol glycosides are important flavonoids that contribute to the formation of characteristic flavors in tea. Therefore, these metabolites attracted our attention, and their biosynthetic pathways (Figure 4A) and accumulation patterns (Figure 4B) were further investigated. Five differentially accumulated anthocyanidins were detected, among which cyanidin was more highly accumulated in LF1 and LF2, whereas pelargonidin-3,5-O-diglucoside and pelargonidin-3-O-glucoside presented the greatest accumulation in CO2 and ZJ2. Notably, similar to procyanidin B2 (Figure 3), cyanidin 3-rutinoside and cyanidin 3-O-glucoside were significantly more abundant in ZJ1 and ZJ2. Four flavan 3-ols and 9 flavonol glycosides were differentially accumulated among these samples. Three flavan 3-ols (EC, C, and CG) and 3 flavonol glycosides (quercetin 3-rutinoside, quercetin 3-O-glucoside, and kaempferol 3-O-rutinoside) exhibited the greatest accumulation in ZJ2, whereas 1 flavan 3-ol (epiafzelechin) and 2 flavonol glycosides (myricetin 3-O-rhamnoside and quercetin 3-L-rhamnoside) exhibited the greatest accumulation in CO2. Additionally, quercetin 3-O-neohesperidoside and kaempferol 3-O-beta-sophoroside presented significantly greater accumulations in ZJ1 and ZJ2, whereas kaempferol-3,7-O-a-L-dirhamnoside and kaempferol 3-O-glucoside exhibited the greatest accumulations in CO1 and LF1, respectively. These results reveal complex accumulation patterns of flavan 3-ols and flavonol glycosides in the CO, LF, and ZJ cultivars. Notably, more flavan 3-ols and flavonol glycosides accumulated in the buds than in the one bud and two leaves stages of the CO, LF, and ZJ cultivars (Figure S1) and may be important factors affecting the taste of tea made from the buds or the one bud and two leaves stage of these cultivars.

2.4. Analysis of RNA Sequencing

To gain insight into the important genes underlying alterations in flavonoid accumulation across different leaf samples and cultivars, an RNA-seq analysis was conducted (Table S4). A Pearson analysis indicated that there was high repeatability among the 3 biological replicate samples (Figure S2). The PCA revealed that PC1 and PC2 of the gene expression data together accounted for 58.95% of the total variance in the expressed genes and clearly distinguished six samples (Figure 5A). Nine genes were randomly selected to verify the accuracy of the transcriptome data via a qRT-PCR analysis, and the results reveal that the expression patterns of these genes detected via qRT-PCR analysis were consistent with those detected via transcriptome sequencing (Figure S3), indicating that the transcriptome sequencing results could be reliable. Subsequently, 5349, 4655, 5670, 5452, 4831, and 5702 DEGs were identified in the CO1_vs_LF1, CO1_vs_ZJ1, LF1_vs_ZJ1, CO2_vs_LF2, CO2_vs_ZJ2, and LF2_vs_ZJ2 comparison groups, respectively (Figure 5B and Table S5), indicating significant variations in the gene expression profiles of the samples of the buds and one bud with two leaves among the CO, LF, and ZJ cultivars. These DEGs were subjected to a clustering analysis, which revealed six significant distinct clusters according to their expression levels (Figure 5C). The genes in cluster 3, cluster 2, and cluster 5 presented relatively high expression levels in the CO, LF, and ZJ cultivars, respectively, whereas those in cluster 6, cluster 4, and cluster 1 presented relatively low expression levels in the CO, LF, and ZJ cultivars, respectively. A KEGG enrichment analysis revealed significantly distinct characteristics of different cluster members; specifically, ‘anthocyanin biosynthesis’ was enriched in genes whose expression was relatively high in ZJ, whereas ‘sesquiterpenoid and triterpenoid biosynthesis’ and ‘flavonoid biosynthesis’ were enriched in genes whose expression was relatively low in CO and LF, respectively, which was consistent with the metabolite results (Figure 1E).

2.5. Correlation Analysis of Flavonoid Levels and Related Synthetic Gene Expression

There were 21 DEGs related to the biosynthesis of anthocyanidins, flavan 3-ols, and flavonol glycosides identified in this study, including 1 F3H, 1 F3′5′H, 6 DFR, 1 ANR, 3 ANS, 1 LAR, and 8 UGT, respectively (Table S6). The expression of these genes was cultivar specific, and the genes clustered into three patterns, with cluster 1, cluster 2, and cluster 3 members being highly expressed in CO, ZJ, and LF, respectively (Figure 6A). An association analysis revealed that some of these DEGs were highly correlated with specific flavonoids (r > 0.8 and p < 0.05) and thus may contribute to their accumulation (Figure 6B and Table S7). For example, flavonol 3-O-glucosyltransferase genes (CsUGT78A15-1/CSS0020068, CsUGT72AM1/CSS0004725, CsBZ1-1/CSS0004477, and CsBZ1-2/CSS0026138) were found to be associated with cyanidin 3-O-glucoside, and CsUGT72AM1 was also associated with pelargonidin 3-O-glucoside and quercetin 3-O-glucoside. CsUGT78A15-1 encodes a protein highly conserved with CsUGT78A15 (Figure S4), which is a broad substrate that functionally catalyzes the formation of a series of flavonol 3-O-glucosides in tea plants; CsUGT72AM1 has been demonstrated to catalyze the formation of quercetin 3-O-glucoside and cyanidin 3-O-glucoside; and CsBZ1-1 and CsZB1-2 are referred to as anthocyanidin 3-O-glucosyltransferase (EC: 2.4.1.115) members that catalyze the formation of anthocyanidin 3-O-glucoside. In addition, our results reveal that CsDFRb1 (CSS0016543) was associated with catechin (C) and epicatechin (EC), whereas CsDFRb3 (CSS0033342), CsANS (CSS0018498), and CsANSa (CSS0010687) were associated with cyanidin. Together, these findings suggest the complexity of the accumulation of specific flavonoids that involve multiple genes.

2.6. Identification of Transcription Factors Involved in Flavonoid Biosynthesis

Transcription factors (TFs) are important regulators that are involved in plant metabolism accumulation by binding to the promoters of structural genes. Therefore, to investigate the potential TFs involved in the accumulation of the specific flavonoids mentioned above, TF binding site scanning of their related structural genes was first performed. The results reveal that multiple TF binding sites were present in all of the gene promoters, whereas MYB and Dof TF binding sites were present in all of the investigated genes (Figure 7A and Table S8). A total of 62 MYBs and 12 Dofs were subsequently identified from the DEGs (Figure 7B and Table S9). Among them, 15 MYBs and 5 Dofs were highly associated with some investigated structural genes (r > 0.9 and p < 0.05) and were therefore considered to be involved in the accumulation of specific flavonoids (Figure 7C and Table S10). Specifically, eight MYBs (CSS0002706, CSS0008558, CSS0033018, CSS0048639, CSS0005971, CSS0023980, CSS0039199, and CSS0046853) and two Dofs (CSS0009351 and CSS0017752) were predicted to be involved in cyanidin accumulation by regulating CsDFRb3, CsANS, CsANSa, or CsF3Hb expression; two MYB (CSS0003042 and CSS0022533) and two Dofs (CSS0031872 and CSS0048616) were predicted to be involved in catechin and epicatechin accumulation by regulating CsDFRb1; and six MYBs (CSS0039604, CSS0016431, CSS0028980, CSS0042913, CSS0022533, and CSS0038191) and three Dofs (CSS0037953, CSS0031872, and CSS0048616) were predicted to be involved in cyanidin 3-O-glucoside by regulating CsUGT78A7815-1, CsBZ1-1, or CsBZ1-2 expression. CSS0028980, CSS0038191, CSS0031872, and CSS0048616 were also predicted to be involved in the accumulation of cyanidin 3-O-glucoside, pelargonidin 3-O-glucoside, and quercetin 3-O-glucoside through regulating CsUGT72AM1 expression. Notably, most specific structural genes were predicted to be regulated by multiple TFs, suggesting a complex regulatory network related to specific flavonoids that is mediated by multiple TFs and structural genes.

3. Discussion

An abundant array of tea germplasm resources, characterized by substantial diversity in the accumulation of biochemical metabolites, particularly flavonoids such as catechins (flavan-3-ols), anthocyanidins, and flavonol glycosides, exert a definitive influence on both tea processing procedures and the ultimate quality of the tea product. These flavonoids have attracted significant attention because of their health benefits and pivotal roles in imparting bitterness and astringency to tea. Additionally, anthocyanidins are crucial factors influencing the coloration of tea plant leaves. Consequently, uncovering the underlying molecular mechanisms responsible for the variations in the accumulation of these compounds among tea cultivars exhibiting distinct color phenotypes represents a promising avenue for the utilization of tea plant resources. To achieve this goal, the flavonoid levels in tea cultivars with these color characteristics were investigated, and a transcriptional regulatory network with flavonoid synthesis structural genes as its core was constructed.
Catechin, anthocyanidin, and flavonol glycoside compounds exhibit distinct variational patterns of accumulation among tea cultivars with different color phenotypes, which are attributed to alterations in the expression of related synthetase genes, including F3′H, DFR, ANS, ANR, LAR, FLS, and UGT [18,21,49,50]. For example, the green-leaf tea cultivar ‘Longjing 43′ presented lower expression levels of F3′H, FLS, DFR, and UGT than the purple-leaf mutant ‘Mooma 1’, resulting in decreased levels of anthocyanidin compounds, several glycosylated flavonols, and two catechin compounds [18], underscoring the complexity of the accumulation of these flavonoids in tea plants, which is orchestrated by a network of related biosynthetic genes.
In addition to the direct involvement of ANR and LAR in the biosynthesis of monomeric catechins, DFR and ANS also play synergistic roles [51,52,53]. For example, compared with green shoots, yellow mutant shoots of ‘Danzicha’ have decreased expression of DFR and ANS members, accompanied by lower catechin levels [54]. In our study, the levels of C and EC in ZJ were greater than those in CO and LF, which may be related to DFRb1. Moreover, other DFR, ANS, ANR, and LAR members presented relatively high expression levels in LF or CO, and these genes are primarily involved in the formation of other flavonoids not focused on in this study. UGTs play an important role in the diversity of metabolites by catalyzing the transfer of an activated sugar donor to acceptor molecules. In tea plants, several UGTs, such as CsUGT14, CsUGT15, CsUGT73A20, and CsUGT72AM1, were found to be responsible for the biosynthesis of flavonol 3-O-glucosides and anthocyanidin 3-O-glucosides [18,22,54,55]. Furthermore, BZ1 is a UGT that was suggested to catalyze the formation of anthocyanidin 3-glucoside in plants. In the present study, CsBZ1-1, CsBZ1-2, CsUGT78A15-1, and CsUGT72AM1 were found to be associated with cyanidin 3-O-glucoside. Moreover, CsUGT72AM1 was also associated with pelargonidin 3-O-glucoside and quercetin 3-O-glucoside, suggesting their synergistic role in the abundance of these compounds in ZJ. Several studies have indicated that cyanidin levels are notably greater in purple-leaf tea cultivars [21,50], whereas other reports indicate that cyanidin levels are not significantly different among yellow-leaf, green-leaf, and purple-leaf tea cultivars [56]. Interestingly, we found that the cyanidin levels in LF were greater than those in CO and ZJ because of the high expression of CsANS, CsANSa, and CsDFRb3 in LF, which encode enzymes directly responsible for its production and upstream reactions.
Transcription factors (TFs) in plants function as critical regulators involved in the flux of flavonoids. In tea plants, several MYB members function as critical regulators of the biosynthesis of catechins and anthocyanins, affecting the expression of DFR, ANS, ANR, and LAR either directly or indirectly [38,39]. Other TFs, such as bHLH, WRKY, and HSF, also have the ability to regulate the synthesis of specific flavonoids [40,49,57]. These findings indicate a complex network of flavonoid synthesis involving multiple structural genes and transcription factors in tea plants. Our study revealed that the promoters of all of the genes responsible for the biosynthesis of specific catechins, anthocyanidins, and flavonol 3-O-glycosides contain MYB-binding cis-elements. Furthermore, multiple MYB genes were highly correlated with the expression of some of these genes, supporting the critical role of MYB genes in regulating flavonoid biosynthesis. Notably, despite the lack of evidence of the involvement of Dof genes in flavonoid synthesis in tea plants, the present study revealed that Dof-binding cis-elements also exist in all investigated gene promoters and that the expression of five Dof genes was highly correlated with the expression of some of these genes. These findings suggest the potential involvement of Dof genes in regulating specific flavonoids, which is worthy of further study. In addition, although the function of UGT in tea plants has been well verified, the upstream regulatory factors of UGT are less well studied. Our study revealed that multiple MYB and Dof genes may synergistically regulate their expression, thereby leading to variations in anthocyanidin 3-glucoside and flavonol glycoside accumulation in different tea cultivars.

4. Materials and Methods

4.1. Plant Materials

Three tea plant cultivars, ‘Mingguan’ (CO), ‘Lvyafoshou’ (LF), and ‘Zijuan’ (ZJ), which were grown in the tea garden of the Tea Research Institute, Fujian Academy of Agricultural Sciences in Fu’an, China (119° 350′ E, 27° 100′ N), were used in this study. Samples of the buds and one bud with two leaves of each cultivar were collected from at least 60 individual healthy plants and were randomly divided into 3 groups representing three biological replicates. All of the samples were immediately frozen with liquid nitrogen after being removed from the plants and stored at −80 °C for further use.

4.2. Metabolite Extraction and UHPLC-MS Analysis

The freeze-dried samples were crushed with a mixer mill for 30 s at 60 Hz. A 100 mg aliquot of each individual sample was precisely weighed and extracted overnight at 4 °C with 500 μL of extraction solution (methanol/water = 3:1, containing an internal standard: 0.3 mg/mL 2-Chloro-L-phenylalanine) on a shaker. The supernatant was filtered through a 0.22 μm microporous membrane, and the resulting supernatants were subsequently diluted 10 times with extraction solution and vortexed for 30 s before UHPLC-MS analysis.
UHPLC-MS analysis was carried out using a 6500+ Triple Quad LC-MS/MS System equipped with an EXIONLC UHPLC unit (AB Sciex, Framingham, MA, USA). An ACQUITY UPLC HSS T3 column (1.8 μm, 2.1 mm × 100 mm) was used for the separation of metabolite compounds. The mobile phase consisted of water containing 1% formic acid (A) and 100% acetonitrile (B), with a flow rate of 0.4 mL/min. The typical ion source parameters were set as follows: ion spray voltage, +5500 and −4500 V; curtain gas, 35 psi; temperature, 400 °C; ion source gas, 1, 60 psi; ion source gas, 2, 60 psi; and DP, ±100 V.

4.3. Metabolite Data Preprocessing and Analysis

SCIEX Analyst Work Station Software (version 1.6.3) was employed for multiple reaction monitoring (MRM) data acquisition and processing. The qualitative analyses of the compounds were performed by comparing the obtained mass spectra information with a self-constructed database (Metware database: MWDB). During the analysis, isotope signals, as well as repeated signals containing K+, Na+, and NH4+ ions, and signals from fragments of other substances with larger molecular weights, were excluded. Metabolites were quantified by calculating the area of each individual peak. A principal component analysis (PCA) of the identified metabolites was performed using the R package gmodels (version 2.19.1). On the basis of variable importance in the projection (VIP) values of the orthogonal projection to latent structures discriminant analysis (OPLS-DA) model, metabolites with VIP ≥ 1.0, a p value < 0.05, and|Log2 (fold change)| > 1 were defined as significantly differentially accumulated metabolites (DAMs). The R package ‘pheatmap’ (version 1.0.12) was used to generate a heatmap for visualization of DAMs.

4.4. RNA Sequencing and Data Analysis

Total RNA was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA). The quality and integrity of the RNA from each sample were assessed by a NanoDrop spectrophotometer (Thermo Scientific, Waltham, MA, USA). The RNA library construction and sequencing were performed by Allwegene Company (Beijing, China). The paired-end 150 bp reads of each library were generated using the Illumina NovaSeq 6000 platform.
The raw sequences were transformed into clean reads by removing reads containing adapters, reads containing poly-N sequences, and low-quality reads. The clean reads were subsequently mapped to the ‘Shuchazao’ reference genome [58] by STAR. HTSeq (version 2.0) was used to count the number of reads mapped to each gene. Gene expression levels were estimated as fragments per kilobase of transcript per million fragments (FPKM). Genes with |Log2 (fold change)| >2 and an adjusted p value < 0.05 were considered to be differentially expressed genes (DEGs). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed via TBtools software (version 2.112). The R package ‘ClusterGVis’ was employed to integrate the expression clusters and KEGG enrichment results of the DEGs.

4.5. Real-Time Quantitative PCR (qRT-PCR) Analysis

Nine genes were randomly selected for qRT-PCR analysis to verify the accuracy of the transcriptome results. GAPDH and PTB were used as reference genes to normalize gene expression using the 2−ΔΔCt method. The specific primers used are listed in Table S1. The first-strand cDNA of each sample was prepared using the HiScript III 1st-Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). qRT-PCR was performed using the ChamQ universal SYBR qPCR Master Mix (Vazyme, China) on a LightCycler 480 II platform (Roche, Mannheim, Germany).

4.6. Correlation Analysis and Transcriptional Regulatory Network Construction

The correlation coefficients between the selected DAM levels and the corresponding synthetic gene expressions were calculated and visualized via the R package Hmisc (version 5.1.3) and the qcorrplot function in the R package Hy4m/linkET, respectively. The expression of the Pearson correlation coefficient (PCC) between transcription factor genes and selected structural genes was calculated via the R package Hmisc (version 5.1.3). Promoters of selected genes, 2000 bp sequences upstream of the start codon, were subjected to PlantRegMap [59] to scan the TF binding sites. The transcriptional regulatory networks were constructed by correlating the PCC  >  0.9 between each TF gene and selected structural gene pairs with the corresponding TF binding sites in the promoters of the selected structural genes.

5. Conclusions

In this study, we conducted an integrated analysis of metabolomics and transcriptomics to investigate the networks regulating astringency and color-related flavonoids in tea cultivars with diverse leaf colors. Through this analysis, we found that EC, C, cyanidin 3-O-glucoside, pelargonidin 3-O-glucoside, and quercetin 3-O-glucoside accumulated to higher levels in the purple-leaf tea cultivar ZJ compared to the yellow-leaf tea cultivar CO and the green-leaf tea cultivar LF. However, cyanidin accumulated more abundantly in the green-leaf tea cultivar LF. These different accumulation patterns of specific flavonoids, affecting the taste and leaf color of tea plants, may be contributed by the expression of eight key structural genes controlled by 15 MYB genes and 5 Dof genes as potential regulators. Despite the established role of MYB genes, there is scant evidence regarding the involvement of Dof genes in flavonoid synthesis within tea plants, which aroused our interest and will be the focus of our future research. Overall, our findings provide insights into the development and utilization of resources from tea plants with diverse leaf colors.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants14040625/s1, Figure S1. Comparison of flavan 3-ols and flavonol glycoside levels between the buds and one bud with two leaves of CO, LF, and ZJ; Figure S2. Correlation analysis of average expression levels of 18 samples (n = 3); Figure S3. qRT-PCR verification of the RNA-seq results for 9 differentially expressed genes (DEGs) in CO, LF, and ZJ samples. The data are shown as the mean ± standard deviations (SD) (n = 3); Figure S4. Sequence alignment of CsUGT78A15 and CsUGT78A15-1; Table S1. Primers used in this study; Table S2. All 863 metabolites were altered; Table S3. DAMs in different comparison groups; Table S4. Summary of RNA-seq results of 18 samples; Table S5. Information on DEGs in different comparison groups; Table S6. Information on flavonoid-related DEGs; Table S7. Pearson correlation coefficients between structural genes and specific flavonoids; Table S8. Information on TF binding sites of investigated structural genes; Table S9. Information on differentially expressed MYB and Dof genes; Table S10. Pearson correlation coefficient between each TF gene and selected structural gene pairs.

Author Contributions

The authors confirm their contributions to this paper as follows: conceptualization, X.W. and C.C.; methodology, R.S. and Y.Z. (Yongheng Zhang); software, R.S.; validation, R.S. and Y.Z. (Yongheng Zhang); formal analysis, Y.Z. (Yazhen Zhang) and X.Y.; investigation, X.K.; resources, X.L.; data curation, L.W. and Y.Z. (Yongheng Zhang); writing—original draft preparation, R.S.; writing—review and editing, Y.Z. (Yongheng Zhang), X.W. and C.C.; supervision, X.W. and C.C.; funding acquisition, X.W. and C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding-Tea Plant (2021C02067-5), Fujian public welfare projects (2022R1029003), the Ministry of Agriculture and Rural Affairs of P. R. China (CARS-19), and the ‘5511’ Collaborative Innovation Project for agricultural high-quality development and transcendence (XTCXGC2021004 and XTCXGC2021019-CYS01).

Data Availability Statement

The original contributions presented in this study are included in the article/Supplementary Materials. Further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fan, F.Y.; Huang, C.S.; Tong, Y.L.; Guo, H.W.; Zhou, S.J.; Ye, J.H.; Gong, S.Y. Widely targeted metabolomics analysis of white peony teas with different storage time and association with sensory attributes. Food Chem. 2021, 362, 130257. [Google Scholar] [CrossRef] [PubMed]
  2. Daglia, M.; Antiochia, R.; Sobolev, A.P.; Mannina, L. Untargeted and targeted methodologies in the study of tea (Camellia sinensis L.). Food Res. Int. 2014, 63, 275–289. [Google Scholar] [CrossRef]
  3. Kaczyński, P.; Iwaniuk, P.; Jankowska, M.; Orywal, K.; Socha, K.; Perkowski, M.; Farhan, J.A.; Łozowicka, B. Pesticide residues in common and herbal teas combined with risk assessment and transfer to the infusion. Chemosphere 2024, 367, 143550. [Google Scholar] [CrossRef] [PubMed]
  4. Feng, X.; Yang, S.; Pan, Y.; Zhou, S.; Ma, S.; Ou, C.; Fan, F.; Gong, S.; Chen, P.; Chu, Q. Yellow tea: More than turning green leaves to yellow. Crit. Rev. Food Sci. Nutr. 2024, 64, 7836–7853. [Google Scholar] [CrossRef]
  5. Wang, Y.; Ho, C.T. Polyphenolic chemistry of tea and coffee: A century of progress. J. Agric. Food Chem. 2009, 57, 8109–8114. [Google Scholar] [CrossRef]
  6. Mira, L.; Tereza Fernandez, M.; Santos, M.; Rocha, R.; Helena Florêncio, M.; Jennings, K.R. Interactions of flavonoids with iron and copper ions: A mechanism for their antioxidant activity. Free Radic. Res. 2002, 36, 1199–1208. [Google Scholar] [CrossRef]
  7. Juadjur, A.; Mohn, C.; Schantz, M.; Baum, M.; Winterhalter, P.; Richling, E. Fractionation of an anthocyanin-rich bilberry extract and in vitro antioxidative activity testing. Food Chem. 2015, 167, 418–424. [Google Scholar] [CrossRef]
  8. Matsuura, R.; Moriyama, H.; Takeda, N.; Yamamoto, K.; Morita, Y.; Shimamura, T.; Ukeda, H. Determination of antioxidant activity and characterization of antioxidant phenolics in the plum vinegar extract of cherry blossom (Prunus lannesiana). J. Agric. Food Chem. 2007, 56, 544–549. [Google Scholar] [CrossRef]
  9. Nabavi, S.M.; Šamec, D.; Tomczyk, M.; Milella, L.; Russo, D.; Habtemariam, S.; Suntar, I.; Rastrelli, L.; Daglia, M.; Xiao, J.; et al. Flavonoid biosynthetic pathways in plants: Versatile targets for metabolic engineering. Biotechnol. Adv. 2020, 38, 107316. [Google Scholar] [CrossRef]
  10. Jiang, X.L.; Shi, Y.F.; Fu, Z.P.; Li, W.W.; Lai, S.Y.; Wu, Y.H.; Wang, Y.S.; Liu, Y.J.; Gao, L.P.; Xia, T. Functional characterization of three flavonol synthase genes from camellia sinensis: Roles in flavonol accumulation. Plant Sci. 2020, 300, 110632. [Google Scholar] [CrossRef]
  11. Harborne, J.B.; Williams, C.A. Advances in flavonoid research since 1992. Phytochemistry 2000, 55, 481–504. [Google Scholar] [CrossRef] [PubMed]
  12. Abotaleb, M.; Samuel, S.; Varghese, E.; Varghese, S.; Kubatka, P.; Liskova, A.; Büsselberg, D. Flavonoids in cancer and apoptosis. Cancers 2018, 11, 28. [Google Scholar] [CrossRef] [PubMed]
  13. Yahfoufi, N.; Alsadi, N.; Jambi, M.; Matar, C. The immunomodulatory and anti-inflammatory role of polyphenols. Nutrients 2018, 10, 1618. [Google Scholar] [CrossRef]
  14. Choy, K.W.; Murugan, D.; Leong, X.F.; Abas, R.; Alias, A.; Mustafa, M.R. Flavonoids as natural anti-inflammatory agents targeting nuclear factor-kappa B (NFΚB) signaling in cardiovascular diseases: A mini review. Front. Pharmacol. 2019, 10, 1295. [Google Scholar] [CrossRef] [PubMed]
  15. Xu, Y.Q.; Zhang, Y.N.; Chen, J.X.; Wang, F.; Du, Q.Z.; Yin, J.F. Quantitative analyses of the bitterness and astringency of catechins from green tea. Food Chem. 2018, 258, 16–24. [Google Scholar] [CrossRef] [PubMed]
  16. Scharbert, S.; Hofmann, T. Molecular definition of black tea taste by means of quantitative studies, taste reconstitution, and omission experiments. J. Agric. Food Chem. 2005, 53, 5377–5384. [Google Scholar] [CrossRef]
  17. Hsu, C.P.; Shih, Y.T.; Lin, B.R.; Chiu, C.F.; Lin, C.C. Inhibitory effect and mechanisms of an anthocyanins- and anthocyanidins-rich extract from purple-shoot tea on colorectal carcinoma cell proliferation. J. Agric. Food Chem. 2012, 60, 3686–3692. [Google Scholar] [CrossRef]
  18. He, X.J.; Zhao, X.C.; Gao, L.P.; Shi, X.X.; Dai, X.L.; Liu, Y.J.; Xia, T.; Wang, Y.S. Isolation and characterization of key genes that promote flavonoid accumulation in purple-leaf tea (Camellia sinensis L.). Sci. Rep. 2018, 8, 130. [Google Scholar] [CrossRef]
  19. Morita, Y.; Ishiguro, K.; Tanaka, Y.; Iida, S.; Hoshino, A. Spontaneous mutations of the UDP-glucose: Flavonoid 3-O-glucosyltransferase gene confers pale-and dull-colored flowers in the Japanese and common morning glories. Planta 2015, 242, 575–587. [Google Scholar] [CrossRef]
  20. Kang, X.F.; Mikami, R.; Akita, Y. Characterization of 5-O-glucosyltransferase involved in anthocyanin biosynthesis in Cyclamen purpurascens. Plant Biotechnol. 2021, 38, 263–268. [Google Scholar] [CrossRef]
  21. Mei, Y.; Xie, H.; Liu, S.R.; Zhu, J.Y.; Zhao, S.Q.; Wei, C.L. Metabolites and transcriptional profiling analysis reveal the molecular mechanisms of the anthocyanin metabolism in the “Zijuan” tea Plant (Camellia sinensis var. assamica). J. Agric. Food Chem. 2021, 69, 414–427. [Google Scholar] [CrossRef] [PubMed]
  22. Zhao, X.C.; Wang, P.J.; Li, M.Z.; Wang, Y.R.; Jiang, X.L.; Cui, L.L.; Qian, Y.M.; Zhuang, J.H.; Gao, L.P.; Xia, T. Functional characterization of a new tea (Camellia sinensis) flavonoid glycosyltransferase. J. Agric. Food Chem. 2017, 65, 2074–2083. [Google Scholar] [CrossRef] [PubMed]
  23. Fang, Z.T.; Song, C.J.; Xu, H.R.; Ye, J.H. Dynamic changes in flavonol glycosides during production of green, yellow, white, oolong and black teas from Camellia sinensis L. (cv. Fudingdabaicha). Int. J. Food Sci. Technol. 2018, 54, 490–498. [Google Scholar] [CrossRef]
  24. Lu, Z.G.; Wang, X.W.; Lin, X.; Mostafa, S.; Zou, H.L.; Wang, L.; Jin, B. Plant anthocyanins: Classification, biosynthesis, regulation, bioactivity, and health benefits. Plant Physiol. Biochem. 2024, 217, 109268. [Google Scholar] [CrossRef]
  25. Terahara, N.; Takeda, Y.; Nesumi, A.; Honda, T. Anthocyanins from red flower tea (Benibana-cha), Camellia sinensis. Phytochemistry 2001, 56, 359–361. [Google Scholar] [CrossRef]
  26. Saito, T.; Honma, D.; Tagashira, M.; Kanda, T.; Nesumi, A.; Maeda-Yamamoto, M. Anthocyanins from new red leaf tea ‘Sunrouge’. J. Agric. Food Chem. 2011, 59, 4779–4782. [Google Scholar] [CrossRef]
  27. Jiang, L.H.; Shen, X.J.; Shoji, T.; Kanda, T.; Zhou, J.C.; Zhao, L.M. Characterization and activity of anthocyanins in zijuan tea (Camellia sinensis var. Kitamura). J. Agric. Food Chem. 2013, 61, 3306–3310. [Google Scholar] [CrossRef]
  28. Lim, E.K.; Ashford, D.A.; Hou, B.; Jackson, R.G.; Bowles, D.J. Arabidopsis glycosyltransferases as biocatalysts in fermentation for regioselective synthesis of diverse quercetin glucosides. Biotechnol. Bioeng. 2004, 87, 623–631. [Google Scholar] [CrossRef]
  29. Zhao, J.; Dixon, R.A. The ‘ins’ and ‘outs’ of flavonoid transport. Trends Plant Sci. 2010, 15, 72–80. [Google Scholar] [CrossRef]
  30. Lai, Y.S.; Li, S.; Tang, Q.; Li, H.X.; Chen, S.X.; Li, P.W.; Xu, J.Y.; Xu, Y.; Guo, X. The dark-purple tea cultivar ‘Ziyan’ accumulates a large amount of delphinidin-related anthocyanins. J. Agric. Food Chem. 2016, 64, 2719–2726. [Google Scholar] [CrossRef]
  31. Achnine, L.; Huhman, D.V.; Farag, M.A.; Sumner, L.W.; Blount, J.W.; Dixon, R.A. Genomics-based selection and functional characterization of triterpene glycosyltransferases from the model legume Medicago truncatula. Plant J. 2005, 41, 875–887. [Google Scholar] [CrossRef] [PubMed]
  32. Devaiah, S.P.; Owens, D.K.; Sibhatu, M.B.; Sarkar, T.R.; Strong, C.L.; Mallampalli, V.K.; Asiago, J.; Cooke, J.; Kiser, S.; Lin, Z.; et al. Identification, recombinant expression, and biochemical analysis of putative secondary product glucosyltransferases from citrus paradisi. J. Agric. Food Chem. 2016, 64, 1957–1969. [Google Scholar] [CrossRef] [PubMed]
  33. Su, X.; Wang, W.; Xia, T.; Gao, L.; Shen, G.; Pang, Y. Characterization of a heat responsive UDP: Flavonoid glucosyltransferase gene in tea plant (Camellia sinensis). PLoS ONE 2018, 13, e0207212. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, X.; Dai, X.; Gao, L.; Guo, L.; Zhuang, J.; Liu, Y.; Ma, X.; Wang, R.; Xia, T.; Wang, Y. Functional analysis of an uridine diphosphate glycosyltransferase involved in the biosynthesis of polyphenolic glucoside in tea plants (Camellia sinensis). J. Agric. Food Chem. 2017, 65, 10993–11001. [Google Scholar] [CrossRef] [PubMed]
  35. Jiang, X.L.; Shi, Y.F.; Dai, X.L.; Zhuang, J.H.; Fu, Z.P.; Zhao, X.Q.; Liu, Y.J.; Gao, L.P.; Xia, T. Four flavonoid glycosyltransferases present in tea overexpressed in model plants Arabidopsis thaliana and Nicotiana tabacum for functional identification. J. Chromatography B Anal. Technol. Biomed. Life Sci. 2018, 1100, 148–157. [Google Scholar] [CrossRef]
  36. Li, P.; Xia, E.; Fu, J.; Xu, Y.; Zhao, X.; Tong, W.; Tang, Q.; Tadege, M.; Fernie, A.R.; Zhao, J. Diverse roles of MYB transcription factors in regulating secondary metabolite biosynthesis, shoot development, and stress responses in tea plants (Camellia sinensis). Plant J. 2022, 110, 1144–1165. [Google Scholar] [CrossRef]
  37. Xu, J.; Li, J.; Liu, Y.; Zheng, P.; Liu, S.; Sun, B. A genus-specific R2R3 MYB transcription factor, CsMYB34, regulates galloylated catechin biosynthesis in Camellia sinensis. Plant Physiol. Biochem. 2024, 219, 109401. [Google Scholar] [CrossRef]
  38. Li, P.; Fu, J.; Xu, Y.; Shen, Y.; Zhang, Y.; Ye, Z.; Tong, W.; Zeng, X.; Yang, J.; Tang, D.; et al. CsMYB1 integrates the regulation of trichome development and catechins biosynthesis in tea plant domestication. New Phytol. 2022, 234, 902–917. [Google Scholar] [CrossRef]
  39. Zhao, X.; Zeng, X.; Lin, N.; Yu, S.; Fernie, A.R.; Zhao, J. CsbZIP1-CsMYB12 mediates the production of bitter-tasting flavonols in tea plants (Camellia sinensis) through a coordinated activator-repressor network. Hortic. Res. 2021, 8, 110. [Google Scholar] [CrossRef]
  40. Hang, Y.; Wang, J.; Xiao, Y.; Wu, Y.; Li, N.; Ding, C.; Hao, X.; Yu, Y.; Wang, L.; Wang, X. CsWRKY12 interacts with CsVQ4L to promote the accumulation of galloylated catechins in tender leaves of tea plants. Plant J. 2024, 120, 2861–2873. [Google Scholar]
  41. Luo, Y.; Yu, S.; Li, J.; Li, Q.; Wang, K.; Huang, J.; Liu, Z. Molecular characterization of WRKY transcription factors that act as negative regulators of O-methylated catechin biosynthesis in tea plants (Camellia sinensis L.). J. Agric. Food Chem. 2018, 66, 11234–11243. [Google Scholar] [CrossRef]
  42. Zhu, J.; Yan, X.; Liu, S.; Xia, X.; An, Y.; Xu, Q.; Zhao, S.; Liu, L.; Guo, R.; Zhang, Z.; et al. Alternative splicing of CsJAZ1 negatively regulates flavan-3-ol biosynthesis in tea plants. Plant J. 2022, 110, 243–261. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, G.; Huang, K.; Ke, J.; Chen, C.; Bao, G.H.; Wan, X. Novel Camellia sinensis O-methyltransferase regulated by CsMADSL1 specifically methylates EGCG in cultivar “GZMe4”. J. Agric. Food Chem. 2023, 71, 6706–6716. [Google Scholar] [CrossRef] [PubMed]
  44. Yang, Y.; He, Z.; Bing, Q.; Duan, X.; Chen, S.; Zeng, M.; Liu, X. Two Dof transcription factors promote flavonoid synthesis in kumquat fruit by activating C-glucosyltransferase. Plant Sci. 2022, 318, 111234. [Google Scholar] [CrossRef]
  45. Skirycz, A.; Jozefczuk, S.; Stobiecki, M.; Muth, D.; Zanor, M.I.; Witt, I.; Mueller-Roeber, B. Transcription factor AtDOF4;2 affects phenylpropanoid metabolism in Arabidopsis thaliana. New Phytol. 2007, 175, 425–438. [Google Scholar] [CrossRef]
  46. Liu, M.Y.; Jiao, Z.; Lou, H.; Tang, D.; Tian, X.; Zhou, B.W.; Ruan, J.; Fernie, A.R.; Zhang, Q. A glutamine synthetase-Dof transcription factor module regulates nitrogen remobilization from source to sink tissues in tea plants. Plant Physiol. 2024, 197, kiae644. [Google Scholar] [CrossRef]
  47. Liu, W.M.; Liu, S.Y.; Zhang, K.Y.; Xie, M.W.; Sun, H.W.; Huang, X.Q.; Zhang, L.X.; Li, M. Chlorophyllase is transcriptionally regulated by CsMYB308/CsDOF3 in young leaves of tea plant. Hortic. Plant J. 2023, 9, 1162–1176. [Google Scholar] [CrossRef]
  48. Li, H.; Huang, W.; Liu, Z.W.; Wang, Y.X.; Zhuang, J. Transcriptome-based analysis of Dof family transcription factors and their responses to abiotic stress in tea plant (Camellia sinensis). Int. J. Genom. 2016, 2016, 5614142. [Google Scholar] [CrossRef]
  49. Zhang, K.K.; Lin, C.Y.; Chen, B.Y.; Lin, Y.E.; Su, H.F.; Du, Y.Y.; Zhang, H.; Zhou, H.; Ji, R.Q.; Zhang, L.Y. A light responsive transcription factor CsbHLH89 positively regulates anthocyanidin synthesis in tea (Camellia sinensis). Sci. Hortic. 2024, 327, 112784. [Google Scholar] [CrossRef]
  50. Chen, X.J.; Wang, P.J.; Zheng, Y.C.; Gu, M.Y.; Lin, X.Y.; Wang, S.Y.; Jin, S.; Ye, N.X. Comparison of metabolome and transcriptome of flavonoid biosynthesis pathway in a purple-leaf tea germplasm Jinmingzao and a green-leaf tea germplasm Huangdan reveals Their Relationship with genetic mechanisms of color formation. Int. J. Mol. Sci. 2020, 21, 4167. [Google Scholar] [CrossRef]
  51. Pang, Y.; Abeysinghe, I.S.; He, J.; He, X.; Huhman, D.; Mewan, K.M.; Sumner, L.W.; Yun, J.; Dixon, R.A. Functional characterization of proanthocyanidin pathway enzymes from tea and their application for metabolic engineering. Plant Physiol. 2013, 161, 1103–1116. [Google Scholar] [CrossRef]
  52. Wang, P.Q.; Zhang, L.J.; Jiang, X.L.; Dai, X.L.; Xu, L.J.; Li, T.; Xing, D.W.; Li, Y.Z.; Li, M.Z.; Gao, L.P.; et al. Evolutionary and functional characterization of leucoanthocyanidin reductases from Camellia sinensis. Planta 2018, 247, 139–154. [Google Scholar] [CrossRef] [PubMed]
  53. Li, Y.; Zhang, C.; Ma, C.; Chen, L.; Yao, M. Transcriptome and biochemical analyses of a chlorophyll-deficient bud mutant of tea plant (Camellia sinensis). Int. J. Mol. Sci. 2023, 24, 15070. [Google Scholar] [CrossRef] [PubMed]
  54. Cui, L.L.; Yao, S.B.; Dai, X.L.; Yin, Q.G.; Liu, Y.J.; Jiang, X.L.; Wu, Y.H.; Qian, Y.M.; Pang, Y.Z.; Gao, L.P.; et al. Identification of UDP-glycosyltransferases involved in the biosynthesis of astringent taste compounds in tea (camellia sinensis). J. Exp. Bot. 2016, 67, 2285–2297. [Google Scholar] [CrossRef] [PubMed]
  55. Mei, X.; Zhang, K.K.; Lin, Y.E.; Su, H.F.; Lin, C.Y.; Chen, B.Y.; Yang, H.J.; Zhang, L.Y. Metabolic and transcriptomic profiling reveals etiolated mechanism in Huangyu tea (Camellia sinensis) leaves. Int. J. Mol. Sci. 2022, 23, 15044. [Google Scholar] [CrossRef]
  56. Cai, J.; Lv, L.; Zeng, X.F.; Zhang, F.; Chen, Y.L.; Tian, W.L.; Li, J.R.; Li, X.Y.; Li, Y. Integrative analysis of metabolomics and transcriptomics reveals molecular mechanisms of anthocyanin metabolism in the Zikui tea plant (Camellia sinensis cv. Zikui). Int. J. Mol. Sci. 2022, 23, 4780. [Google Scholar] [CrossRef]
  57. Zhang, X.Y.; Li, L.Y.; He, Y.Q.; Lang, Z.L.; Zhao, Y.; Tao, H.; Li, Q.S.; Hong, G.J. The CsHSFA-CsJAZ6 module-mediated high temperature regulates flavonoid metabolism in Camellia sinensis. Plant Cell Environ. 2023, 46, 2401–2418. [Google Scholar] [CrossRef]
  58. Xia, E.H.; Tong, W.; Hou, Y.; An, Y.L.; Chen, L.B.; Wu, Q.; Liu, Y.L.; Yu, J.; Li, F.D.; Li, P.H.; et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into its genome evolution and adaptation. Mol. Plant 2020, 13, 1013–1026. [Google Scholar] [CrossRef]
  59. Tian, F.; Yang, D.C.; Meng, Y.Q.; Jin, J.P.; Gao, G. PlantRegMap: Charting functional regulatory maps in plants. Nucleic Acids Res. 2020, 48, D1104–D1113. [Google Scholar] [CrossRef]
Figure 1. Phenotype and metabolite profiling of tea cultivars with different color phenotypes. (A) Leaf phenotypes of CO, LF, and ZJ. (B) Classification of total identified metabolites. (C) Heatmap displaying accumulation of all metabolites in CO, LF, and ZJ leaves. (D) PCA for samples based on all metabolites. (E) Heatmap displaying changes in 13 metabolite categories.
Figure 1. Phenotype and metabolite profiling of tea cultivars with different color phenotypes. (A) Leaf phenotypes of CO, LF, and ZJ. (B) Classification of total identified metabolites. (C) Heatmap displaying accumulation of all metabolites in CO, LF, and ZJ leaves. (D) PCA for samples based on all metabolites. (E) Heatmap displaying changes in 13 metabolite categories.
Plants 14 00625 g001
Figure 2. Phenotype and metabolite profiling of tea cultivars with different color phenotypes. (A) Leaf phenotypes of CO, LF, and ZJ. (B) Classification of total identified metabolites.
Figure 2. Phenotype and metabolite profiling of tea cultivars with different color phenotypes. (A) Leaf phenotypes of CO, LF, and ZJ. (B) Classification of total identified metabolites.
Plants 14 00625 g002
Figure 3. The accumulation profiles of 47 differentially accumulated flavonoids.
Figure 3. The accumulation profiles of 47 differentially accumulated flavonoids.
Plants 14 00625 g003
Figure 4. Differentially accumulated flavonoids in CO, LF, and ZJ leaves. (A) Biosynthetic pathways of anthocyanidins, flavan 3-ols, and flavonol glycosides. (B) Differentially accumulated anthocyanidins, flavan 3-ols, and flavonol glycosides.
Figure 4. Differentially accumulated flavonoids in CO, LF, and ZJ leaves. (A) Biosynthetic pathways of anthocyanidins, flavan 3-ols, and flavonol glycosides. (B) Differentially accumulated anthocyanidins, flavan 3-ols, and flavonol glycosides.
Plants 14 00625 g004
Figure 5. RNA-seq analysis of CO, LF, and ZJ leaves. (A) PCA for samples based on all genes. (B) Numbers of differentially expressed genes (DEGs) in different comparisons. (C) Expression patterns and KEGG enrichment terms of DEGs in different comparisons.
Figure 5. RNA-seq analysis of CO, LF, and ZJ leaves. (A) PCA for samples based on all genes. (B) Numbers of differentially expressed genes (DEGs) in different comparisons. (C) Expression patterns and KEGG enrichment terms of DEGs in different comparisons.
Plants 14 00625 g005
Figure 6. DEGs related to anthocyanidins, flavan 3-ols, and flavonol glycosides. (A) Clustering heatmap of DEGs related to anthocyanidins, flavan 3-ols, and flavonol glycosides. (B) Pearson correlation coefficients between selected DEGs and anthocyanidins, flavan 3-ols, and flavonol glycosides.
Figure 6. DEGs related to anthocyanidins, flavan 3-ols, and flavonol glycosides. (A) Clustering heatmap of DEGs related to anthocyanidins, flavan 3-ols, and flavonol glycosides. (B) Pearson correlation coefficients between selected DEGs and anthocyanidins, flavan 3-ols, and flavonol glycosides.
Plants 14 00625 g006
Figure 7. Identification of potential transcription factors (TFs) involved in biosynthesis of anthocyanidins, flavan 3-ols, and flavonol glycosides. (A) Cis-elements present in promoters of selected genes. (B) Expression of differentially expressed MYB and Dof genes. (C) Potential regulatory networks for anthocyanidin, flavan 3-ols, and flavonol glycoside biosynthesis in CO, LF, and ZJ.
Figure 7. Identification of potential transcription factors (TFs) involved in biosynthesis of anthocyanidins, flavan 3-ols, and flavonol glycosides. (A) Cis-elements present in promoters of selected genes. (B) Expression of differentially expressed MYB and Dof genes. (C) Potential regulatory networks for anthocyanidin, flavan 3-ols, and flavonol glycoside biosynthesis in CO, LF, and ZJ.
Plants 14 00625 g007
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

Shan, R.; Zhang, Y.; You, X.; Kong, X.; Zhang, Y.; Li, X.; Wang, L.; Wang, X.; Chen, C. Revealing the Molecular Regulatory Mechanism of Flavonoid Accumulation in Tender Leaves of Tea Plants by Transcriptomic and Metabolomic Analyses. Plants 2025, 14, 625. https://doi.org/10.3390/plants14040625

AMA Style

Shan R, Zhang Y, You X, Kong X, Zhang Y, Li X, Wang L, Wang X, Chen C. Revealing the Molecular Regulatory Mechanism of Flavonoid Accumulation in Tender Leaves of Tea Plants by Transcriptomic and Metabolomic Analyses. Plants. 2025; 14(4):625. https://doi.org/10.3390/plants14040625

Chicago/Turabian Style

Shan, Ruiyang, Yongheng Zhang, Xiaomei You, Xiangrui Kong, Yazhen Zhang, Xinlei Li, Lu Wang, Xinchao Wang, and Changsong Chen. 2025. "Revealing the Molecular Regulatory Mechanism of Flavonoid Accumulation in Tender Leaves of Tea Plants by Transcriptomic and Metabolomic Analyses" Plants 14, no. 4: 625. https://doi.org/10.3390/plants14040625

APA Style

Shan, R., Zhang, Y., You, X., Kong, X., Zhang, Y., Li, X., Wang, L., Wang, X., & Chen, C. (2025). Revealing the Molecular Regulatory Mechanism of Flavonoid Accumulation in Tender Leaves of Tea Plants by Transcriptomic and Metabolomic Analyses. Plants, 14(4), 625. https://doi.org/10.3390/plants14040625

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