Single-cell transcriptome atlas reveals developmental trajectories and a novel metabolic pathway of catechin esters in tea leaves.

The tea plant is an economically important woody beverage crop. The unique taste of tea is evoked by certain metabolites, especially catechin esters, whereas their precise formation mechanism in different cell types remains unclear. Here, a fast protoplast isolation method was established and the transcriptional profiles of 16 977 single cells from 1st and 3rd leaves were investigated. We first identified 79 marker genes based on six isolated tissues and constructed a transcriptome atlas, mapped developmental trajectories and further delineated the distribution of different cell types during leaf differentiation and genes associated with cell fate transformation. Interestingly, eight differently expressed genes were found to co-exist at four branch points. Genes involved in the biosynthesis of certain metabolites showed cell- and development-specific characteristics. An unexpected catechin ester glycosyltransferase was characterized for the first time in plants by a gene co-expression network in mesophyll cells. Thus, the first single-cell transcriptional landscape in woody crop leave was reported and a novel metabolism pathway of catechin esters in plants was discovered.

Menghai Huangye, a novel albino tea germplasm with high theanine content and a high catechin index.

Yunnan Province has a very wide diversity of tea germplasm resources. A variety of special tea germplasms with outstanding traits have been discovered, including tea germplasms with high anthocyanin content and low caffeine content. Albino tea cultivars generally have higher contents of theanine that contribute to the umami taste, and the quality of tea brewed from it is higher. The catechin index (CI), the ratio of dihydroxylated catechins (DIC) to trihydroxylated catechins (TRIC), is a crucial index of suitability for processing tea. In this study, the albino tea plant Menghai Huangye (MHHY) with yellow leaves was identified. Analysis of the biochemical components revealed that MHHY was enriched in theanine and the total catechins (TC) were lower than Yunkang 10 (YK10). In addition, the CI value of MHHY was extremely significantly higher than that of YK10. Metabolic profile of catechins and the related gene expression profile analysis found that the coordinated expression of the key branch genes F3'H and F3'5'Ha for the synthesis of DIC and TRIC in tea plant was closely related to the high CI and low TC of MHHY. Further analysis of the F3'H promoter showed that a 284-bp deletion mutation was present in the F3'H promoter of MHHY, containing the binding sites of the transcriptional repressor MYB4 involved in flavonoid metabolism, which might be an important reason for the up-regulated expression of F3'H in MHHY. Overall, this study provides a theoretical basis for understanding the characteristics of albino tea germplasm resources and efficiently utilizing high-CI tea germplasm resources.

Genome-Wide Identification and Expression Patterns of the C2H2-Zinc Finger Gene Family Related to Stress Responses and Catechins Accumulation in Camellia sinensis [L.] O. Kuntze.

The C2H2-zinc finger protein (C2H2-ZFP) is essential for the regulation of plant development and widely responsive to diverse stresses including drought, cold and salt stress, further affecting the late flavonoid accumulation in higher plants. Tea is known as a popular beverage worldwide and its quality is greatly dependent on the physiological status and growing environment of the tea plant. To date, the understanding of C2H2-ZFP gene family in Camellia sinensis [L.] O. Kuntze is not yet available. In the present study, 134 CsC2H2-ZFP genes were identified and randomly distributed on 15 chromosomes. The CsC2H2-ZFP gene family was classified into four clades and gene structures and motif compositions of CsC2H2-ZFPs were similar within the same clade. Segmental duplication and negative selection were the main forces driving the expansion of the CsC2H2-ZFP gene family. Expression patterns suggested that CsC2H2-ZFPs were responsive to different stresses including drought, salt, cold and methyl jasmonate (MeJA) treatment. Specially, several C2H2-ZFPs showed a significant correlation with the catechins content and responded to the MeJA treatment, which might contribute to the tea quality and specialized astringent taste. This study will lay the foundations for further research of C2H2-type zinc finger proteins on the stress responses and quality-related metabolites accumulation in C. sinensis.

Multiple-genotypes transcriptional analysis revealed candidates genes and nucleotide variants for improvement of quality characteristics in tea (Camellia sinensis (L.) O. Kuntze).

Tea quality is a polygenic trait that exhibits tremendous genetic variability due to accumulation of array of secondary metabolites. To elucidate global molecular insights controlling quality attributes, metabolite profiling and transcriptome sequencing of twelve diverse tea cultivars was performed in tea shoots harvested during quality season. RP-HPLC-DAD analysis of quality parameters revealed significant difference in catechins, theanine and caffeine contents. Transcriptome sequencing resulted into 50,107 non-redundant transcripts with functional annotations of 81.6% (40,847) of the transcripts. Interestingly, 2872 differentially expressed transcripts exhibited significant enrichment in 38 pathways (FDR ≤ 0.05) including secondary metabolism, amino acid and carbon metabolism. Thirty-eight key candidates reportedly involved in biosynthesis of fatty acid derived volatiles, volatile terpenes, glycoside hydrolysis and key quality related pathways (flavonoid, caffeine and theanine-biosynthesis) were highly expressed in catechins-rich tea cultivars. Furthermore, enrichment of candidates involved in flavonoid biosynthesis, transcriptional regulation, volatile terpene and biosynthesis of fatty acid derived volatile in Protein-Protein Interactome network revealed well-coordinated regulation of quality characteristics in tea. Additionally, ascertainment of 23,649 non-synonymous SNPs and validation of candidate SNPs present in quality related genes suggests their potential utility in genome-wide mapping and marker development for expediting breeding of elite compound-rich tea cultivars.

Transcriptomic analysis of tea plant (Camellia sinensis) revealed the co-expression network of 4111 paralogous genes and biosynthesis of quality-related key metabolites under multiple stresses.

The tea plant is an essential economic plant in many countries. However, its growing season renders them vulnerable to stresses. To understand the transcriptomic influences of these stresses on tea plants, we sequenced and analyzed the transcriptomes under drought, high-temperature, and pest. Paralogs were identified by comparing 14 evolutionarily close genomes. The differentially expressed paralog (DEPs) genes were analyzed regarding single or multiple stresses, and 1075 of the 4111 DEPs were commonly found in all the stresses. The co-expression network of the DEPs and TFs indicated that genes of catechin biosynthesis were associated with most transcription factors specific to each stress. The genes playing a significant role in the late response to drought and pest stress mainly functioned in the early response to high-temperature. This study revealed the relationship between stress and regulation of QRM synthesis and the role of QRMs in response to these (a)biotic stresses.

Integrative Transcriptomic and Metabolic Analyses Provide Insights into the Role of Trichomes in Tea Plant (Camellia Sinensis).

Trichomes, which develop from epidermal cells, are regarded as one of the key features that are involved in the evaluation of tea quality and tea germplasm resources. The metabolites from trichomes have been well characterized in tea products. However, little is known regarding the metabolites in fresh tea trichomes and the molecular differences in trichomes and tea leaves per se. In this study, we developed a method to collect trichomes from tea plant tender shoots, and their main secondary metabolites, including catechins, caffeine, amino acids, and aroma compounds, were determined. We found that the majority of these compounds were significantly less abundant in trichomes than in tea leaves. RNA-Seq was used to investigate the differences in the molecular regulatory mechanism between trichomes and leaves to gain further insight into the differences in trichomes and tea leaves. In total, 52.96 Gb of clean data were generated, and 6560 differentially expressed genes (DEGs), including 4471 upregulated and 2089 downregulated genes, were identified in the trichomes vs. leaves comparison. Notably, the structural genes of the major metabolite biosynthesis pathways, transcription factors, and other key DEGs were identified and comparatively analyzed between trichomes and leaves, while trichome-specific genes were also identified. Our results provide new insights into the differences between tea trichomes and leaves at the metabolic and transcriptomic levels, and open up new doors to further recognize and re-evaluate the role of trichomes in tea quality formation and tea plant growth and development.

Preferential assimilation of NH4+ over NO3- in tea plant associated with genes involved in nitrogen transportation, utilization and catechins biosynthesis.

Physiological effects of ammonium (NH4+) and nitrate (NO3-) on tea have confirmed that tea plants prefer NH4+ as the dominant nitrogen (N) source. To investigate the possible explanations for this preference, studies of 15NH4+ and 15NO3- assimilation using hydroponically grown tea plants were conducted. During the time course of 15NH4+ and 15NO3- assimilation, the absorption of 15N from 15NH4+ was more rapid than that from 15NO3-, as there was a more efficient expression pattern of NH4+ transporters compared with that of NO3- transporters. 15NH4+-fed tea plants accumulated more 15N than 15NO3- fed plants, which was demonstrated by that genes related to primary N assimilation, like CsNR, CsNiR, CsGDH and CsGOGAT, were more affected by 15NH4+ than 15NO3-. Markedly higher NH4+ concentrations were observed in 15NH4+-fed tea roots in comparison with NO3- treatment, whereas tea plants maintained a balanced concentration of NH4+ in tea leaves under both these two N forms. This maintenance was achieved through the increased expression of genes involved in theanine biosynthesis and the inhibition of genes related to catechins derived from phenylpropanoid pathway. The current results suggest that efficient NH4+ transportation, assimilation, and reutilization enables tea plant as an ammonium preferring plant species.

Discerning between Two Tuscany (Italy) Ancient Apple cultivars, 'Rotella' and 'Casciana', through Polyphenolic Fingerprint and Molecular Markers.

Ancient apple cultivars usually have higher nutraceutical value than commercial ones, but in most cases their variability in pomological traits does not allow us to discriminate among them. Fruit of two Tuscany ancient apple cultivars, 'Casciana' and 'Rotella', picked from eight different orchards (four for each cultivar) were analyzed for their pomological traits, organoleptic qualities, polyphenolic profile and antiradical activity. The effectiveness of a polyphenol-based cluster analysis was compared to molecular markers (internal transcribed spacers, ITS1 and ITS2) to unequivocally discern the two apples. 'Casciana' and 'Rotella' fruit had a higher nutraceutical value than some commercial cultivars, in terms of phenolic abundance, profile and total antiradical activity. Although pedo-climatic conditions of different orchards influenced the phenolic profile of both apples, the polyphenolic discriminant analysis clearly separated the two cultivars, principally due to higher amounts of procyanidin B2, procyanidin B3 and p-coumaroylquinic acid in 'Casciana' than in 'Rotella' fruit. These three polyphenols can be used proficiently as biochemical markers for distinguishing the two apples when pomological traits cannot. Conversely, ITS1 and ITS2 polymorphism did not allow us to distinguish 'Casciana' from 'Rotella' fruit. Overall, the use of polyphenolic fingerprint might represent a valid tool to ensure the traceability of products with a high economic value.

Characterization of the transcriptional regulator CsbHLH62 that negatively regulates EGCG3"Me biosynthesis in Camellia sinensis.

Epigallocatechin-3-O-(3-O-methyl) gallate (EGCG3"Me) in tea (Camellia sinensis (L.) O. Kuntze) is a major source of O-methylated catechin and renowned for a wide range of health effects. However, the transcriptional regulation mechanisms of EGCG3"Me biosynthesis remain unclear. In the present work, the basic Helix-Loop-Helix (bHLH) transcription factor, designated as CsbHLH62, belonging to GBOF group of bHLH families, was isolated and characterized from Camellia sinensis. CsbHLH62 contains an Open Reading Frame of 1662 bp and encodes a polypeptide of 553 amino acids. Subcellular location and transcriptional activity analysis showed it as a nucleus protein and possessed transcriptional inhibition activity. Furthermore, the expression of CsbHLH62 was decreased during EGCG3"Me accumulation. More importantly, E-box motifs (5'-CANNTG-3') were found in the promoters of CCoAOMT, CsLAR, and CsDFR, and further transient expression assays showed that CsbHLH62 repressed the transcription of CCoAOMT, CsLAR, and CsDFR. Collectively, these results suggest that CsbHLH62 acts as a transcriptional repressor that might be negatively affecting the accumulation of EGCG3"Me. These findings provide novel insights into the regulatory mechanism of EGCG3"Me biosynthesis, which might help to breed high EGCG3"Me-content tea plants.

The acaricidal activity and mechanism of eugenol on Psoroptes cuniculi.

In this study, the acaricidal effect of eugenol was measured and its mechanism of action investigated. The results showed that eugenol possessed the effect of killing Psoroptes cuniculi, and could regulate the mRNA expression of glutathione S-transferase (GST), catechinic acid (Ca) and thioredoxin (Trx). PPAR, NF-kappa B, TNF, Rap 1 and Ras signaling pathways might be the main pathways that involved into the process of killing mites. These findings suggested that eugenol could be developed into a new kind of acaricide, and further expand current knowledge on the mechanisms of eugenol for killing Psoroptes cuniculi of eugenol.

Identification of key genes involved in catechin metabolism in tea seedlings based on transcriptomic and HPLC analysis.

Tea is a non-alcoholic beverage with many benefits to human health and thereby widely consumed in the world. It contains plenty of secondary metabolites and tea catechins are the characteristic compounds. To further elucidate the biosynthetic and regulatory mechanisms of catechins in tea, high performance liquid chromatography (HPLC) and transcriptome analysis were performed in tea seedlings of different growth stages. A combined method of differential expression and correlation analysis was then conducted. The results showed that the order of total catechin (TC) contents was leaves > stems > roots, irrespective of growth stages. For transcriptome analysis, a total of 355.81 million clean reads were generated and mapped to the referencing tea genome. Further real time PCR analysis of 18 selected genes confirmed RNA-Seq results. A total of 7 structural genes and 35 transcription factors (TFs) were identified to be significantly correlated with TC changes. Among them, three TFs homologous to ANL2, WRKY44 and AtMYB113 might play key roles in catechin regulation. The de novo transcriptome data of different organs in tea seedlings provided new insights into the biosynthetic and metabolic pathways of catechins.

Molecular Evidence for Catechin Synthesis and Accumulation in Tea Buds (Camellia sinensis).

Early spring buds of the Camellia sinensis variety Shuchazao were separated into two parts, including the shoot tip (ST) and non-expanded young leaves (YL), in which the synthesis and accumulation of catechins in the two parts were assessed by high-performance liquid chromatography (HPLC), p-dimethylaminocinnamaldehyde (DMACA) staining, quantitative real-time polymerase chain reaction (qRT-PCR), and in situ hybridization. HPLC showed that (-)-epigallocatechin-3-gallate (EGCG) and (-)-epicatechin-3-gallate (ECG) amounts in YL were increased significantly by 74.0 and 71.8%, respectively. The results of DMACA staining indicated that catechins in buds accumulated mainly in mesophyll cells and the bud shaft of YL. Meanwhile, qRT-PCR demonstrated that the relative expression levels of genes related to flavonoid metabolism, including CsPAL1, CsC4H1, CsC4H2, CsCHS2, CsF3'5'H1, CsDFR1, CsDFR2, and CsANR1, were significantly higher in YL than in the ST. In situ hybridization revealed that CsDFR1, CsDFR2, CsLAR, and CsANR1 were expressed in leaf primordia and YL but not in the apical meristem. These findings highlight the synthesis and accumulation patterns of catechins in different parts of the ST in C. sinensis, providing a theoretical basis for the assessment of synthesis, accumulation, and transfer patterns of catechins in tea plants.

Transcriptional profiling of catechins biosynthesis genes during tea plant leaf development.

MAIN CONCLUSION: A total of 299,113 unigenes were generated and 15,817 DEGs were identified. We identified candidate genes associated with the regulation of catechins biosynthesis during leaf development in tea plant. The tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most economically significant crops worldwide because of its positive effects on human health. The health benefits of tea are mainly attributed to catechins, which are the predominant polyphenols that accumulate in tea. Catechins are products of the phenylpropanoid and flavonoid biosynthetic pathways. Although catechins were identified in tea leaves long ago, the molecular mechanisms regulating catechins biosynthesis remain unclear. To identify candidate genes involved in catechins biosynthesis, we analyzed the transcriptomes of tea leaves during five different leaf stages of development using RNA-seq. Approximately 809 million high-quality reads were obtained, trimmed, and assembled into 299,113 unigenes with an average length of 565 bp. A total of 15,817 unigenes were differentially expressed during the different stages of leaf development. These differentially expressed genes were enriched in a variety of processes such as the regulation of the cell cycle, starch and sucrose metabolism, photosynthesis, phenylpropanoid biosynthesis, phenylalanine metabolism, and flavonoid biosynthesis. Based on their annotations, 51 of these differentially expressed unigenes are involved in phenylpropanoid and flavonoid biosynthesis. Furthermore, transcription factors such as MYB, bHLH and MADS, which may involve in the regulation of catechins biosynthesis, were identified through co-expression analysis of transcription factors and structural genes. Real-time PCR analysis of candidate genes indicated a good correlation with the transcriptome data. These findings increase our understanding of the molecular mechanisms regulating catechins biosynthesis in the tea plant.

Development of CAPS markers to identify Indian tea (Camellia sinensis) clones with high catechin content.

Tea leaves are rich in plant secondary phenolics, especially flavonoids. Catechins are considered to be the most valuable flavonoids, and the catechin content in tea is an important trait for determining its quality. We have developed cleaved amplified polymorphic sequence (CAPS)-based markers for evaluating total catechin content that target two important secondary metabolite pathway genes, PAL (phenylalanine ammonia-lyase) and CHS (chalcone synthase). Catechin content levels in the tea samples tested ranged from 9 to 33 mg/mg. The CAPS technique identified clones with the homozygous profile PRc1, which has relatively lower catechin content than clones with the heterozygous profile PRc2. A significant difference (t = 16.85) in the level of catechin content was also detected between heterozygotes and homozygotes in the tea seed stock TS379. We found a polynomial relationship between the marker developed for CHS2 and catechin content in these tea samples with R2 = 0.9788. Moreover, PAL has less of a relationship with catechin content. Therefore, we recommend tea clones with heterozygous CAPS profiles for the gene CHS2+RcaI for the further improvement in these clones.

Transcriptome Analysis Reveals Key Flavonoid 3'-Hydroxylase and Flavonoid 3',5'-Hydroxylase Genes in Affecting the Ratio of Dihydroxylated to Trihydroxylated Catechins in Camellia sinensis.

The ratio of dihydroxylated to trihydroxylated catechins (RDTC) is an important indicator of tea quality and biochemical marker for the study of genetic diversity. It is reported to be under genetic control but the underlying mechanism is not well understood. Flavonoid 3'-hydroxylase (F3'H) and flavonoid 3',5'-hydroxylase (F3'5'H) are key enzymes involved in the formation of dihydroxylated and trihydroxylated catechins. The transcriptome and HPLC analysis of tea samples from Longjing43 and Zhonghuang2 under control and shading treatment were performed to assess the F3'H and F3'5'H genes that might affect RDTC. A total of 74.7 million reads of mRNA seq (2×101bp) data were generated. After de novo assembly, 109,909 unigenes were obtained, and 39,982 of them were annotated using 7 public databases. Four key F3'H and F3'5'H genes (including CsF3'5'H1, CsF3'H1, CsF3'H2 and CsF3'H3) were identified to be closely correlated with RDTC. Shading treatment had little effect on RDTC, which was attributed to the stable expression of these key F3'H and F3'5'H genes. The correlation of the coexpression of four key genes and RDTC was further confirmed among 13 tea varieties by real time PCR and HPLC analysis. The coexpression of three F3'H genes and a F3'5'H gene may play a key role in affecting RDTC in Camellia sinensis. The current results may establish valuable foundation for further research about the mechanism controlling catechin composition in tea.

Transcriptomic and phytochemical analysis of the biosynthesis of characteristic constituents in tea (Camellia sinensis) compared with oil tea (Camellia oleifera).

BACKGROUND: Tea plants (Camellia sinensis) are used to produce one of the most important beverages worldwide. The nutritional value and healthful properties of tea are closely related to the large amounts of three major characteristic constituents including polyphenols (mainly catechins), theanine and caffeine. Although oil tea (Camellia oleifera) belongs to the genus Camellia, this plant lacks these three characteristic constituents. Comparative analysis of tea and oil tea via RNA-Seq would help uncover the genetic components underlying the biosynthesis of characteristic metabolites in tea. RESULTS: We found that 3,787 and 3,359 bud genes, as well as 4,042 and 3,302 leaf genes, were up-regulated in tea and oil tea, respectively. High-performance liquid chromatography (HPLC) analysis revealed high levels of all types of catechins, theanine and caffeine in tea compared to those in oil tea. Activation of the genes involved in the biosynthesis of these characteristic compounds was detected by RNA-Seq analysis. In particular, genes encoding enzymes involved in flavonoid, theanine and caffeine pathways exhibited considerably different expression levels in tea compared to oil tea, which were also confirmed by quantitative RT-PCR (qRT-PCR). CONCLUSION: We assembled 81,826 and 78,863 unigenes for tea and oil tea, respectively, based on their differences at the transcriptomic level. A potential connection was observed between gene expression and content variation for catechins, theanine and caffeine in tea and oil tea. The results demonstrated that the metabolism was activated during the accumulation of characteristic metabolites in tea, which were present at low levels in oil tea. From the molecular biological perspective, our comparison of the transcriptomes and related metabolites revealed differential regulatory mechanisms underlying secondary metabolic pathways in tea versus oil tea.

De novo assembly and transcriptome characterization: novel insights into catechins biosynthesis in Camellia sinensis.

BACKGROUND: Tea is a popular natural non-alcoholic beverage consumed worldwide due to its bioactive ingredients, particularly catechins (flavan-3-ols). Catechins not only contribute to tea quality but also serve important functions in the anti-stress regulation of secondary metabolic pathways. However, the percentages of various catechins are different among tea plant [Camellia sinensis (L.) O. Kuntze] cultivars. This study aimed to elucidate the biosynthetic mechanism of catechins. Transcriptomes from leaf tissues of four tea plant cultivars, 'Yunnanshilixiang', 'Chawansanhao', 'Ruchengmaoyecha', and 'Anjibaicha', were sequenced using the high-throughput sequencing platform Illumina HiSeq™ 2000. De novo assemble were also performed. Catechins contents were measured through reversed-phase high-performance liquid chromatography (RP-HPLC), and the biosynthetic pathway was also surveyed. RESULTS: We constructed a unified unigene database. A total of 146,342 pairs of putative orthologs from the four tea plant cultivars, 'Yunnanshilixiang', 'Chawansanhao', 'Ruchengmaoyecha', and 'Anjibaicha' were generated. Approximately 68,890 unigenes (47.1%) were aligned to the sequences of seven public databases with a cut-off E-value of 1E-5. A total of 217 differentially expressed genes were found through RPKM values, and 150 unigenes were assigned to the flavonoid biosynthetic pathway using the integrated function annotation. The (-)-EGC and (-)-EC contents were significantly lower and the (+)-GC and (+)-C contents were abnormally higher in 'Ruchengmaoyecha' than in 'Yunnanshilixiang', 'Chawansanhao', and 'Anjibaicha'. The proportion of catechins was confirmed by selecting critical genes (ANS, ANR, and LAR) for qRT-PCR analysis. CONCLUSIONS: This study provided a global survey of transcriptomes from four tea plant cultivars and serves as an available resource of genetic diversity. The analyses of transcriptome profiles and physiological indicators not only identified the putative genes involved in the flavonoid biosynthetic pathway but also provided some novel insights for the mechanisms of catechins biosynthesis.

Engineering the production of major catechins by Escherichia coli carrying metabolite genes of Camellia sinensis.

A mimicked biosynthetic pathway of catechin metabolite genes from C. sinensis, consisting of flavanone 3 hydroxylase (F3H), dihydroflavonol reductase (DFR), and leucoanthocyanidin reductase (LCR), was designed and arranged in two sets of constructs: (a) single promoter in front of F3H and ribosome-binding sequences both in front of DFR and LCR; (b) three different promoters with each in the front of the three genes and ribosome-binding sequences at appropriate positions. Recombinant E. coli BL (DE3) harbouring the constructs were cultivated for 65 h at 26 °C in M9 medium consisting of 40 g/L glucose, 1 mM IPTG, and 3 mM eriodictyol. Compounds produced were extracted in ethyl acetate in alkaline conditions after 1 h at room temperature and identified by HPLC. Two of the four major catechins, namely, (-)-epicatechin (0.01) and (-)-epicatechin gallate (0.36 mg/L), and two other types ((+)-catechin hydrate (0.13 mg/L) and (-)-catechin gallate (0.04 mg/L)) were successfully produced.

Differential gene expression in tea (Camellia sinensis L.) calli with different morphologies and catechin contents.

Tea (Camellia sinensis) is a commercially important crop that contains valuable secondary metabolites. To understand the molecular regulation of secondary metabolism in tea, we selected and analyzed two cell lines of tea callus (Yunjing63Y and Yunjing63X) that showed different morphological characteristics and catechin contents. Yunjing63Y callus was yellow and tight, while yunjing63X callus was white and loose. HPLC analyses showed that Yunjing63Y contained 3.71 times higher levels of catechins than Yunjing63X. Using cDNA amplified fragment-length polymorphism (cDNA-AFLP) we identified 68 genes that were differentially expressed between the two lines. Of the 68 differentially expressed ESTs, 40 showed higher expressions in Yunjing63Y and 28 showed higher expressions in Yunjing63X. BLASTX comparisons classified these ESTs into seven functional groups; phenylpropanoid metabolism (2.9%), UDPG-dependent glucosyl transferase (8.8%), transcription factors (11.8%), transporters (13.2%), signal transduction (19.1%), other metabolism (26.5%), and unknown (17.7%). We used qRT-PCR to validate the expression of genes and ESTs, and found that genes associated with flavan-3-ols biosynthesis and metabolism were expressed at higher levels in Yunjing63Y than in Yunjing63X. In addition, the expression of ESTs associated with flavonoid biosynthesis, regulation and transport were higher in Yunjing63Y than in Yunjing63X. The full-length cDNA of a EST coding for a putative MYB transcription factor was amplified using rapid amplification of cDNA ends (RACE). The resulting 1270 bp long cDNA, named CsMYB1, contained a 933-bp ORF encoding a 310-amino acid protein with a predicted molecular weight of 105.27 kDa and a predicted isoelectric point of 4.85 and showed highest homology to plant MYBs likely involved in stress signaling.