Identification of Candidate Anthocyanin-Related Genes by Transcriptomic Analysis of 'Furongli' Plum (Prunus salicina Lindl.) during Fruit Ripening Using RNA-Seq

Insights into tissue-specific anthocyanin accumulation in Japanese plum (Prunus salicina L.) fruits: A comparative study of three cultivars

In the present study, three matured Japanese plum cultivars with different colored peel and flesh were selected to mine the key transcription factors regulating anthocyanin formation in tissues. Results showed that PsMYB10 was correlated with structural genes C4H, F3H, and ANS. PsMYB6 could positively regulate C4H (r = 0.732) and accumulated anthocyanins in Sanhua plum's flesh. Sanhua plum has the highest phenolic and anthocyanin contents (10.24 ± 0.37 gallic acid equivalent mg g-1 dry weight (DW) and 68.95 ± 1.03 μg g-1 DW), resulting itself superior biological activity as 367.1 ± 42.9 Trolox equivalent mg g-1 DW in oxygen radical absorbance capacity value and 72.79 ± 4.34 quercetin equivalent mg g-1 DW in cellular antioxidant activity value. The present work provides new insights into the regulatory mechanism of tissue-specific anthocyanin biosynthesis, confirming the pivotal role of anthocyanins in the biological activity of plums, providing essential support for the development of horticultural products enriched with anthocyanins.

Polyphenol Composition, Antioxidant Capacity and Xanthine Oxidase Inhibition Mechanism of Furong Plum Fruits at Different Maturity Stages

An experiment was conducted on the polyphenol content, flavonoid content, anthocyanin content, and antioxidant capacity of Furong plum (Prunus salicina Lindl. cv. "furong") at different maturity stages to determine the most suitable maturity stage. The inhibition of plum polyphenols on xanthine oxidase (XOD) was measured, and its kinetics were studied to reveal the inhibitory mechanism. The experimental results showed that the polyphenol, flavonoid and anthocyanin contents of plums at the ripe stage were the highest, reaching 320.46 mg GAE/100 g FW, 204.21 mg/100 g FW, and 66.24 mg/100 g FW, respectively, in comparison those of the plums at the immature and mid-ripe stages. The antioxidant capacity of the ripe plums was stronger than it was during the other stages of the plums growth. Among them, the total polyphenols of the ripe plums exhibited the strongest antioxidant capacity (IC50 values against DPPH and hydroxyl radicals were 28.19 ± 0.67 μg/mL and 198.16 ± 7.55 μg/mL, respectively), which was between the antioxidant capacity of the free polyphenols and bound polyphenols. The major phenolic monomer compounds of plum polyphenols were flavan-3-ols (epicatechin, catechin, proanthocyanidin, and procyanidin B2), flavonols (myricetin), and phenolic acids (chlorogenic acid, ferulic acid, and protocatechuic acid). Additionally, plum polyphenols exhibited a strong inhibitory effect on XOD, with an IC50 value of 77.64 μg/mL. The inhibition kinetics showed that plum polyphenols are mixed-type inhibitors that inhibit XOD activity and that the inhibition process is reversible. The calculated values of Ki and α were 16.53 mmol/L and 0.26, respectively.

Function verification of a chlorophyll a/b binding protein gene through a newly established tobacco rattle virus-induced gene silencing system in Kandelia obovata

As an important mangrove species, Kandelia obovata plays an irreplaceable role in the coastal ecosystem. However, due to a lack of genetic technology, there is limited research on its functional genes. As such, establishing an efficient and rapid functional verification system is particularly important. In this study,tobacco rattle virus (TRV) and the phytoene desaturase gene KoPDS were used as the vector and target gene, respectively, to establish a virus-induced gene silencing system (VIGS) in K. obovata. Besides, the system was also used to verify the role of a Chlorophyll a/b binding protein (Cab) gene KoCAB in leaf carbon sequestration of K. obovata. RNA-Seq and qRT-PCR showed that the highest gene-silencing efficiency could reach 90% after 10 days of inoculation and maintain above 80% after 15 days, which was achieved with resuspension buffer at pH 5.8 and Agrobacterium culture at OD600 of 0.4-0.6. Taken together, the TRV-mediated VIGS system established herein is the first genetic analysis tool for mangroves, which may greatly impel functional genomics studies in mangrove plants.

Genome-Wide Identification and Expression Profiling Analysis of SWEET Family Genes Involved in Fruit Development in Plum (Prunus salicina Lindl)

SWEETs (sugars will eventually be exported transporters) play a vital role in longer-distance sugar transportation, and thus control carbon flow and energy metabolism in plants. SWEET genes have been identified in various plant species, but their functions in fruit development remain uncharacterized. Here, we isolated 15 putative PsSWEETs from the Prunus salicina genome. For further analysis, comprehensive bioinformatics methods were applied to determine the gene structure, chromosome distribution, phylogeny, cis-acting regulatory elements, and expression profiles of PsSWEETs. qRT-PCR analysis suggested that these SWEETs might have diverse functions in the development of plum fruit. The relative expression levels of PsSWEET1 and PsSWEET9 were obviously higher in ripened fruit than the ones in other developmental stages, suggesting their possible roles in the transport and accumulation of sugars in plum fruit. Positive correlations were found between the expression level of PsSWEET3/10/13 and the content of sucrose, and the expression level of PsSWEET2 and the content of fructose, respectively, during the development of 'Furongli' fruit, suggesting their possible roles in the accumulation of sucrose and fructose. The current study investigated the initial genomic characterization and expression patterns of the SWEET gene family in plum, which could provide a foundation for the further understanding of the functional analysis of the SWEET gene family.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Genome–Transcriptome Transition Approaches to Characterize Anthocyanin Biosynthesis Pathway Genes in Blue, Black and Purple Wheat

Colored wheat has gained enormous attention from the scientific community, but the information available on the anthocyanin biosynthetic genes is very minimal. The study involved their genome-wide identification, in silico characterization and differential expression analysis among purple, blue, black and white wheat lines. The recently released wheat genome mining putatively identified eight structural genes in the anthocyanin biosynthesis pathway with a total of 1194 isoforms. Genes showed distinct exon architecture, domain profile, regulatory elements, chromosome emplacement, tissue localization, phylogeny and synteny, indicative of their unique function. RNA sequencing of developing seeds from colored (black, blue and purple) and white wheats identified differential expressions in 97 isoforms. The F3H on group two chromosomes and F3′5′H on 1D chromosomes could be significant influencers in purple and blue color development, respectively. Apart from a role in anthocyanin biosynthesis, these putative structural genes also played an important role in light, drought, low temperature and other defense responses. The information can assist in targeted anthocyanin production in the wheat seed endosperm.

Identification of Candidate Genes Associated with Pulp Color by Transcriptomic Analysis of 'Huaxiu' Plum (Prunus salicina Lindl.) during Fruit-Ripening

The plum (Prunus salicina Lindl.) is one of the traditional and economically important stone fruit trees in China. Anthocyanins are important pigments in plums. However, little is known about the molecular mechanisms underlying anthocyanin accumulation in plum fruits, which has hindered research on the molecular mechanism of its utilization. Our research shows that the chlorophyll content was gradually decreased and the contents of anthocyanin and flavonoid increased during the coloring process of the pulp in 'Huaxiu' plums (P. salicina). Then, the RNA-Seq technique was used to analyze the transcriptome of pulp color changes with three different stages (yellow, orange, and red) in the 'Huaxiu' plum (P. salicina). A total of 57,119 unigenes with a mean length of 953 bp were generated, and 61.6% of them were annotated to public databases. The Gene Ontology (GO) database assigned 21,438 unigenes with biological process, cellular components, and molecular function. In addition, 32,146 unigenes were clustered into 25 categories for functional classification by the COG database, and 7595 unigenes were mapped to 128 KEGG pathways by the KEGG pathway database. Of these, 1095 (YS-versus-OS), 4947 (YS-versus-RS), and 3414 (OS-versus-RS) genes were significantly expressed differentially between two coloration stages. The GO and KEGG pathway enrichment analysis revealed that 20 and 1 differentially expressed genes (DEG) are involved in flavonoid biosynthesis and anthocyanin biosynthesis, respectively. Finally, we mainly identified three structural genes as candidate genes. The transcriptome information in this study provide a basis for further studies of pulp colors in plum and contribute to our understanding of the molecular mechanisms underlying anthocyanin biosynthesis in pulp.

An LTR retrotransposon in the promoter of a PsMYB10.2 gene associated with the regulation of fruit flesh color in Japanese plum

Japanese plums exhibit wide diversity of fruit coloration. The red to black hues are caused by the accumulation of anthocyanins, while their absence results in yellow, orange or green fruits. In Prunus, MYB10 genes are determinants for anthocyanin accumulation. In peach, QTLs for red plant organ traits map in an LG3 region with three MYB10 copies (PpMYB10.1, PpMYB10.2 and PpMYB10.3). In Japanese plum the gene copy number in this region differs with respect to peach: there are at least three copies of PsMYB10.1, with the expression of one of them (PsMYB10.1a) correlating with fruit skin color. The objective of this study was to determine a possible role of LG3-PsMYB10 genes in the natural variability of the flesh color trait and to develop a molecular marker for marker-assisted selection (MAS). We explored the variability within the LG3-PsMYB10 region using long-range sequences obtained in previous studies through CRISPR-Cas9 enrichment sequencing. We found that the PsMYB10.2 gene was only expressed in red flesh fruits. Its role in promoting anthocyanin biosynthesis was validated by transient overexpression in Japanese plum fruits. The analysis of long-range sequences identified an LTR retrotransposon in the promoter of the expressed PsMYB10.2 gene that explained the trait in 93.1% of the 145 individuals analyzed. We hypothesize that the LTR retrotransposon may promote the PsMYB10.2 expression and activate the anthocyanin biosynthesis pathway. We propose for the first time the PsMYB10.2 gene as candidate for the flesh color natural variation in Japanese plum and provide a molecular marker for MAS.

Color Stability Enhancement and Antioxidation Improvement of Sanhua Plum Wine under Circulating Ultrasound

Anthocyanins contribute to the attractive color of fruit wine, and their excessive degradation is deleterious to quality, especially for wine with an inherently low anthocyanin content, such as Sanhua plum wine. Ultrasonic treatment is well recognized for wine color maintenance. In the present study, fresh Sanhua plum wine was ultrasonic-treated and aged in barrels for three months. Our results demonstrate that ultrasonic treatment at 28 and 40 kHz improves color performance, as expressed by an increase in a*, b*, and C* values and color intensity, which is highly related to copigmentation. This successful conservation was attributed to the inactivation of polyphenol oxidase and the corresponding reduction in anthocyanin degradation. Finally, the increased antioxidative ability was verified due to the hydrogen donating ability of the surviving anthocyanins. This study indicates the reliability of ultrasonic treatment for providing superior colorfastness during Sanhua plum wine aging, which is also of great potential in processing different fruit wines.

Monitoring Apricot (Prunus armeniaca L.) Ripening Progression through Candidate Gene Expression Analysis

This study aimed at the monitoring of the apricot (Prunus armeniaca L.) ripening progression through the expression analysis of 25 genes related to fruit quality traits in nine cultivars with great differences in fruit color and ripening date. The level of pigment compounds, such as anthocyanins and carotenoids, is a key factor in food taste, and is responsible for the reddish blush color or orange skin and flesh color in apricot fruit, which are desirable quality traits in apricot breeding programs. The construction of multiple linear regression models to predict anthocyanins and carotenoids content from gene expression allows us to evaluate which genes have the strongest influence over fruit color, as these candidate genes are key during biosynthetic pathways or gene expression regulation, and are responsible for the final fruit phenotype. We propose the gene CHS as the main predictor for anthocyanins content, CCD4 and ZDS for carotenoids content, and LOX2 and MADS-box for the beginning and end of the ripening process in apricot fruit. All these genes could be applied as RNA markers to monitoring the ripening stage and estimate the anthocyanins and carotenoids content in apricot fruit during the ripening process.

Identification of a Strong Anthocyanin Activator, VbMYBA, From Berries of Vaccinium bracteatum Thunb

Wufanshu (Vaccinium bracteatum Thunb.), which is a wild member of the genus Vaccinium, accumulates high concentration of anthocyanin in its berries. In this study, the accumulated anthocyanins and their derivatives in Wufanshu berries were identified through UHPLC-MS/MS analysis. Candidate anthocyanin biosynthetic genes were identified from the transcriptome of Wufanshu berries. qRT-PCR analyses showed that the expression of anthocyanin structural genes correlated with anthocyanin accumulation in berries. The R2R3-MYB, VbMYBA, which is a homolog of anthocyanin promoting R2R3-MYBs from other Vaccinium species, was also identified. Transient expression of VbMYBA in Nicotiana tabacum leaves confirmed its role as an anthocyanin regulator, and produced a higher anthocyanin concentration when compared with blueberry VcMYBA expression. Dual-luciferase assays further showed that VbMYBA can activate the DFR and UFGT promoters from other Vaccinium species. VbMYBA has an additional 23 aa at the N terminus compared with blueberry VcMYBA, but this was shown not to affect the ability to regulate anthocyanins. Taken together, our results provide important information on the molecular mechanisms responsible for the high anthocyanin content in Wufanshu berries.

Metabolome and transcriptome profiling of Theobroma cacao provides insights into the molecular basis of pod color variation

The Theobroma cacao presents a wide diversity in pod color among different cultivars. Although flavonoid biosynthesis has been studied in many plants, molecular mechanisms governing the diversity of coloration in cacao pods are largely unknown. The flavonoid metabolite profiles and flavonoid biosynthetic gene expression in the pod exocarps of light green pod 'TAS 410' (GW), green pod 'TAS 166' (GF), and mauve pod 'TAS 168' (PF) were determined. Changes in flavonoid metabolites, particularly the anthocyanins (cyanidin 3-O-galactoside, cyanidin 3-O-glucoside, and cyanidin O-syringic acid) were significantly up-accumulated in the mauve phenotype (PF) compared to the light green or green phenotypes, endowing the pod color change from light green or green to mauve. Consistently, the PF phenotype showed different expression patterns of flavonoid biosynthetic structural genes in comparison with GW/GF phenotypes. The expression level of LAR and ANR in GW/GF was significantly higher than PF, while the expression level of UFGT in GW/GF was lower than PF. These genes likely generated more anthocyanins in the exocarps samples of PF than that of GW/GF. Simultaneously, colorless flavan-3-ols (catechin, epicatechin and proanthocyanidin) content in the exocarp samples of PF was lower than GW/GF. Additionally, MYB (gene18079) and bHLH (gene5045 and gene21575) may participate in the regulation of the pod color. This study sheds light on the molecular basis of cacao pod color variation, which will contribute to breeding cacao varieties with enhanced flavonoid profiles for nutritional applications.

Transcriptome analysis and identification of genes associated with floral transition and fruit development in rabbiteye blueberry (Vaccinium ashei)

Flowering and fruit set are important traits affecting fruit quality and yield in rabbiteye blueberry (Vaccinium ashei). Intense efforts have been made to elucidate the influence of vernalization and phytohormones on flowering, but the molecular mechanisms of flowering and fruit set remain unclear. To unravel these mechanisms, we performed transcriptome analysis to explore blueberry transcripts from flowering to early fruit stage. We divided flowering and fruit set into flower bud (S2), initial flower (S3), bloom flower (S4), pad fruit (S5), and cup fruit (S6) based on phenotype and identified 1,344, 69, 658, and 189 unique differentially expressed genes (DEGs) in comparisons of S3/S2, S4/S3, S5/S4, and S6/S5, respectively. There were obviously more DEGs in S3/S2 and S5/S4 than in S4/S3, and S6/S5, suggesting that S3/S2 and S5/S4 represent major transitions from buds to fruit in blueberry. GO and KEGG enrichment analysis indicated these DEGs were mostly enriched in phytohormone biosynthesis and signaling, transporter proteins, photosynthesis, anthocyanins biosynthesis, disease resistance protein and transcription factor categories, in addition, transcript levels of phytohormones and transporters changed greatly throughout the flowering and fruit set process. Gibberellic acid and jasmonic acid mainly acted on the early stage of flowering development like expression of the florigen gene FT, while the expression of auxin response factor genes increased almost throughout the process from bud to fruit development. Transporter proteins were mainly associated with minerals during the early flowering development stage and sugars during the early fruit stage. At the early fruit stage, anthocyanins started to accumulate, and the fruit was susceptible to diseases such as fungal infection. Expression of the transcription factor MYB86 was up-regulated during initial fruit development, which may promote anthocyanin accumulation. These results will aid future studies exploring the molecular mechanism underlying flowering and fruit set of rabbiteye blueberry.

Transcriptome-based analysis reveals that the biosynthesis of anthocyanins is more active than that of flavonols and proanthocyanins in the colorful flowers of Lagerstroemia indica

Mechanisms associated with the control of flower color in crape myrtle varieties have yet to be sufficiently elucidated, which has tended to hamper the use of modern molecular and genetic strategies in the breeding programs for this plant. The whole transcriptome of four L. indica varieties characterized by different flower colors (white, light purple, deep purplish pink, and strong red) was sequenced, and we performed bioinformatic, quantitative PCR, and co-expression analyses of R2R3 MYB transcription factor and anthocyanin/flavonol pathway genes. We obtained a total of 49,980 transcripts with full-length coding sequences. Both transcriptome and qPCR analyses revealed that anthocyanin/flavonol pathway genes were differentially expressed among the four different flowers types, with the expression of LiPAL, LiCHS, LiCHI, LiDFR, LiANS/LDOX, and LiUFGT being induced in colorful flowers, whereas that of LiF3´5´H, LiFLS, and LiLAR was found to be inhibited. Base on phylogenetic analysis, seven R2R3 MYB transcriptional factors were identified as putative regulators of flower color. The molecular characteristics and co-expression patterns indicated that these MYBs differentially modulate their target genes, with two probably acting as activators, three as repressors, and one contributing to the regulation of vacuolar pH. The findings of this study indicate that the anthocyanin biosynthesis is more active than the flavonol and proanthocyanin in the colorful flowers. These observations provide new genomic information on L. indica and contribute gene resources for the flower color-targeted breeding of crape myrtle.

Activation of PsMYB10.2 Transcription Causes Anthocyanin Accumulation in Flesh of the Red-Fleshed Mutant of 'Sanyueli' (Prunus salicina Lindl.)

Plum is one of the most important stone fruits in the world and anthocyanin-rich plums are increasingly popular due to their health-promoting potential. In this study, we investigated the mechanisms of anthocyanin accumulation in the flesh of the red-fleshed mutant of the yellow-fleshed plum 'Sanyueli'. RNA-Seq and qRT-PCR showed that anthocyanin biosynthetic genes and the transcription factor PsMYB10.2 were upregulated in the flesh of the mutant. Functional testing in tobacco leaves indicated that PsMYB10.2 was an anthocyanin pathway activator and can activate the promoter of the anthocyanin biosynthetic genes PsUFGT and PsGST. The role of PsMYB10.2 in anthocyanin accumulation in the flesh of plum was further confirmed by virus-induced gene silencing. These results provide information for further elucidating the underlying mechanisms of anthocyanin accumulation in the flesh of plum and for the breeding of new red-fleshed plum cultivars.

De novo transcriptome sequencing and anthocyanin metabolite analysis reveals leaf color of Acer pseudosieboldianum in autumn

Background

Leaf color is an important ornamental trait of colored-leaf plants. The change of leaf color is closely related to the synthesis and accumulation of anthocyanins in leaves. Acer pseudosieboldianum is a colored-leaf tree native to Northeastern China, however, there was less knowledge in Acer about anthocyanins biosynthesis and many steps of the pathway remain unknown to date.

Results

Anthocyanins metabolite and transcript profiling were conducted using HPLC and ESI-MS/MS system and high-throughput RNA sequencing respectively. The results demonstrated that five anthocyanins were detected in this experiment. It is worth mentioning that Peonidin O-hexoside and Cyanidin 3, 5-O-diglucoside were abundant, especially Cyanidin 3, 5-O-diglucoside displayed significant differences in content change at two periods, meaning it may be play an important role for the final color. Transcriptome identification showed that a total of 67.47 Gb of clean data were obtained from our sequencing results. Functional annotation of unigenes, including comparison with COG and GO databases, yielded 35,316 unigene annotations. 16,521 differentially expressed genes were identified from a statistical analysis of differentially gene expression. The genes related to leaf color formation including PAL, ANS, DFR, F3H were selected. Also, we screened out the regulatory genes such as MYB, bHLH and WD40. Combined with the detection of metabolites, the gene pathways related to anthocyanin synthesis were analyzed.

Conclusions

Cyanidin 3, 5-O-diglucoside played an important role for the final color. The genes related to leaf color formation including PAL, ANS, DFR, F3H and regulatory genes such as MYB, bHLH and WD40 were selected. This study enriched the available transcriptome information for A. pseudosieboldianum and identified a series of differentially expressed genes related to leaf color, which provides valuable information for further study on the genetic mechanism of leaf color expression in A. pseudosieboldianum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07715-x.

Functional Annotation of a Full-Length Transcriptome and Identification of Genes Associated with Flower Development in Rhododendronsimsii (Ericaceae)

Rhododendronsimsii is one of the top ten famous flowers in China. Due to its historical value and high aesthetic, it is widely popular among Chinese people. Various colors are important breeding objectives in Rhododendron L. The understanding of the molecular mechanism of flower color formation can provide a theoretical basis for the improvement of flower color in Rhododendron L. To generate the R.simsii transcriptome, PacBio sequencing technology has been used. A total of 833,137 full-length non-chimeric reads were obtained and 726,846 high-quality full-length transcripts were found. Moreover, 40,556 total open reading frames were obtained; of which 36,018 were complete. In gene annotation analyses, 39,411, 18,565, 16,102 and 17,450 transcriptions were allocated to GO, Nr, KEGG and COG databases, correspondingly. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with Protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT) and Coding Non Coding Index (CNCI) databases and observed 6170, 2265, 4084 and 1240 lncRNAs, respectively. Based on the results, most genes were enriched in the flavonoid biosynthetic pathway. The eight key genes on the anthocyanin biosynthetic pathway were further selected and analyzed by qRT-PCR. The F3'H and ANS showed an upward trend in the developmental stages of R. simsii. The highest expression of F3'5'H and FLS in the petal color formation of R. simsii was observed. This research provided a huge number of full-length transcripts, which will help to proceed genetic analyses of R.simsii. native, which is a semi-deciduous shrub.

Developmental transcriptome profiling uncovered carbon signaling genes associated with almond fruit drop

Almond is one of the most featured nut crops owing to its high nutritional value. However, due to three different waves of flower and fruitlet drop, fruit drop is a major concern for growers. In this study, we carried out a time-course transcriptome analysis to investigate gene expression differences between normal and abnormal fruitlet development. By de novo assembly analysis, we identified 33,577 unigenes and provided their functional annotations. In total, we identified 7,469 differentially expressed genes and observed the most apparent difference between normal and abnormal fruits at 12 and 17 days after flowering. Their biological functions were enriched in carbon metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction. RT-qPCR validated the expression pattern of 14 representative genes, including glycosyltransferase like family 2, MYB39, IAA13, gibberellin-regulated protein 11-like and POD44, which confirmed the reliability of our transcriptome data. This study provides an insight into the association between abnormal fruit development and carbohydrate signaling from the early developmental stages and could be served as useful information for understanding the regulatory mechanisms related to almond fruit drop.

Molecular Bases of Fruit Quality in Prunus Species: An Integrated Genomic, Transcriptomic, and Metabolic Review with a Breeding Perspective

In plants, fruit ripening is a coordinated developmental process that requires the change in expression of hundreds to thousands of genes to modify many biochemical and physiological signal cascades such as carbohydrate and organic acid metabolism, cell wall restructuring, ethylene production, stress response, and organoleptic compound formation. In Prunus species (including peaches, apricots, plums, and cherries), fruit ripening leads to the breakdown of complex carbohydrates into sugars, fruit firmness reductions (softening by cell wall degradation and cuticle properties alteration), color changes (loss of green color by chlorophylls degradation and increase in non-photosynthetic pigments like anthocyanins and carotenoids), acidity decreases, and aroma increases (the production and release of organic volatile compounds). Actually, the level of information of molecular events at the transcriptional, biochemical, hormonal, and metabolite levels underlying ripening in Prunus fruits has increased considerably. However, we still poorly understand the molecular switch that occurs during the transition from unripe to ripe fruits. The objective of this review was to analyze of the molecular bases of fruit quality in Prunus species through an integrated metabolic, genomic, transcriptomic, and epigenetic approach to better understand the molecular switch involved in the ripening process with important consequences from a breeding point of view.

Stone Fruit as Biofactories of Phytochemicals With Potential Roles in Human Nutrition and Health

Phytochemicals or secondary metabolites present in fruit are key components contributing to sensory attributes like aroma, taste, and color. In addition, these compounds improve human nutrition and health. Stone fruits are an important source of an array of secondary metabolites that may reduce the risk of different diseases. The first part of this review is dedicated to the description of the main secondary organic compounds found in plants which include (a) phenolic compounds, (b) terpenoids/isoprenoids, and (c) nitrogen or sulfur containing compounds, and their principal biosynthetic pathways and their regulation in stone fruit. Then, the type and levels of bioactive compounds in different stone fruits of the Rosaceae family such as peach (Prunus persica), plum (P. domestica, P. salicina and P. cerasifera), sweet cherries (P. avium), almond kernels (P. dulcis, syn. P. amygdalus), and apricot (P. armeniaca) are presented. The last part of this review encompasses pre- and postharvest treatments affecting the phytochemical composition in stone fruit. Appropriate management of these factors during pre- and postharvest handling, along with further characterization of phytochemicals and the regulation of their synthesis in different cultivars, could help to increase the levels of these compounds, leading to the future improvement of stone fruit not only to enhance organoleptic characteristics but also to benefit human health.

Genomic Dissection of a Wild Region in a Superior Solanum pennellii Introgression Sub-Line with High Ascorbic Acid Accumulation in Tomato Fruit

The Solanum pennellii introgression lines (ILs) have been exploited to map quantitative trait loci (QTLs) and identify favorable alleles that could improve fruit quality traits in tomato varieties. Over the past few years, ILs exhibiting increased content of ascorbic acid in the fruit have been selected, among which the sub-line R182. The aims of this work were to identify the genes of the wild donor S. pennellii harbored by the sub-line and to detect genes controlling ascorbic acid accumulation by using genomics tools. A Genotyping-By-Sequencing (GBS) approach confirmed that no wild introgressions were present in the sub-line besides one region on chromosome 7. By using a dense single nucleotide polymorphism (SNP) map obtained by RNA sequencing (RNA-Seq), the wild region of the sub-line was finely identified; thus, defining 39 wild genes that replaced 33 genes of the ILs genetic background (cv. M82). The differentially expressed genes mapping in the region and the variants detected among the cultivated and the wild alleles evidenced the potential role of the novel genes present in the wild region. Interestingly, one upregulated gene, annotated as a major facilitator superfamily protein, showed a novel structure in R182, with respect to the parental lines. These genes will be further investigated using gene editing strategies.

Transcriptional reprogramming strategies and miRNA-mediated regulation networks of Taxus media induced into callus cells from tissues

Background

Taxus cells are a potential sustainable and environment-friendly source of taxol, but they have low survival ratios and slow grow rates. Despite these limitations, Taxus callus cells induced through 6 months of culture contain more taxol than their parent tissues. In this work, we utilized 6-month-old Taxus media calli to investigate their regulatory mechanisms of taxol biosynthesis by applying multiomics technologies. Our results provide insights into the adaptation strategies of T. media by transcriptional reprogramming when induced into calli from parent tissues.

Results

Seven out of 12 known taxol, most of flavonoid and phenylpropanoid biosynthesis genes were significantly upregulated in callus cells relative to that in the parent tissue, thus indicating that secondary metabolism is significantly strengthened. The expression of genes involved in pathways metabolizing biological materials, such as amino acids and sugars, also dramatically increased because all nutrients are supplied from the medium. The expression level of 94.1% genes involved in photosynthesis significantly decreased. These results reveal that callus cells undergo transcriptional reprogramming and transition into heterotrophs. Interestingly, common defense and immune activities, such as “plant–pathogen interaction” and salicylic acid- and jasmonic acid-signaling transduction, were repressed in calli. Thus, it’s an intelligent adaption strategy to use secondary metabolites as a cost-effective defense system. MiRNA- and degradome-sequencing results showed the involvement of a precise regulatory network in the miRNA-mediated transcriptional reprogramming of calli. MiRNAs act as direct regulators to enhance the metabolism of biological substances and repress defense activities. Given that only 17 genes of secondary metabolite biosynthesis were effectively regulated, miRNAs are likely to play intermediate roles in the biosynthesis of secondary metabolites by regulating transcriptional factors (TFs), such as ERF, WRKY, and SPL.

Conclusion

Our results suggest that increasing the biosynthesis of taxol and other secondary metabolites is an active regulatory measure of calli to adapt to heterotrophic culture, and this alteration mainly involved direct and indirect miRNA-induced transcriptional reprogramming. These results expand our understanding of the relationships among the metabolism of biological substances, the biosynthesis of secondary metabolites, and defense systems. They also provide a series of candidate miRNAs and transcription factors for taxol biosynthesis.

Transcriptomic Analysis of Ficus carica Peels with a Focus on the Key Genes for Anthocyanin Biosynthesis

Fig (Ficus carica L.), a deciduous fruit tree of the Moraceae, provides ingredients for human health such as anthocyanins. However, little information is available on its molecular structure. In this study, the fig peels in the yellow (Y) and red (R) stages were used for transcriptomic analyses. Comparing the R with the Y stage, we obtained 6224 differentially expressed genes, specifically, anthocyanin-related genes including five CHS, three CHI, three DFR, three ANS, two UFGT and seven R2R3-MYB genes. Furthermore, three anthocyanin biosynthetic genes, i.e., FcCHS1, FcCHI1 and FcDFR1, and two R2R3-MYB genes, i.e., FcMYB21 and FcMYB123, were cloned; sequences analysis and their molecular characteristics indicated their important roles in fig anthocyanin biosynthesis. Heterologous expression of FcMYB21 and FcMYB123 significantly promoted anthocyanin accumulation in both apple fruits and calli, further suggesting their regulatory roles in fig coloration. These findings provide novel insights into the molecular mechanisms behind fig anthocyanin biosynthesis and coloration, facilitating the genetic improvement of high-anthocyanin cultivars and other horticultural traits in fig fruits.

Analysis of Metabolites and Gene Expression Changes Relative to Apricot (Prunus armeniaca L.) Fruit Quality During Development and Ripening

Apricot (Prunus armeniaca L.) is a valuable worldwide agronomical crop, with a delicious fruit highlighted as a functional food with both nutritional and bioactive properties, remarkably beneficial to human health. Apricot fruit ripening is a coordinated developmental process which requires change in the expression of hundreds to thousands of genes to modify many biochemical and physiological processes arising from quality characteristics in ripe fruit. In addition, enhancing fruit and nutraceutical quality is one of the central objectives to be improved in the new varieties developed by breeding programs. In this study we analyzed the contents of main metabolites linked to the nutraceutical value of apricot fruits, together with the most important pomological characteristics and biochemical contents of fruit during the ripening process in two contrasted apricot genotypes. Additionally, the gene expression changes were analyzed using RNA-Seq and real time qPCR. Results showed that genes with differential expression in the biosynthetic pathways, such as phenylpropanoids, flavonoids, starch and sucrose and carotenoid metabolism, could be possible candidates as molecular markers of fruit quality characteristics for fruit color and soluble solid content. The gene involves in carotenoid metabolism carotenoid cleavage dioxygenase 4, and the gene sucrose synthase in starch and sucrose metabolism were identified as candidate genes in the ripening process for white skin ground color and flesh color and high soluble sugar content. The application of these candidate genes on marker-assisted selection in apricot breeding programs may contribute to the early selection of high-quality fruit genotypes with suitable nutraceutical values.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts

BACKGROUND

Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species.

RESULTS

The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration.

CONCLUSIONS

In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Alteration of flesh color and enhancement of bioactive substances via the stimulation of anthocyanin biosynthesis in 'Friar' plum fruit by low temperature and the removal

'Friar' plum (Prunus salicina Lindl.) fruit were transferred to shelf life (25 °C) following different storage periods at low (0 °C) or intermediate (5 °C) temperature. The earliest flesh reddening appeared in plums during shelf life removed after 28 d at 0 °C and 14 d at 5 °C, respectively, in comparison with turning yellow in plums remained at 25 °C immediately after harvest. The flesh reddening developed rapidly thereafter. Anthocyanins, in particular, cyanidin 3-O-glucoside, significantly accumulated in the reddening tissue, and activities of enzymes associated with the phenylpropanoid pathway were considerably activated after the removal. The removal elicited extremely high ethylene production in plums, which might mediate the activation of the anthocyanin biosynthesis in response to cold stress signal. The results provided a potential approach for postharvest regulation of flesh color and accumulation of bioactive substances in plums, which could lead to attractive attributes and health-promoting effects on consumers.

Transcriptome profiling reveals key genes related to astringency during cucumber fruit development

The goal of this study was to provide quantitative data on the catechin contents and underlying molecular regulatory mechanisms in cucumber during fruit development. The dynamic changes in the total catechin contents and RNA-seq-based transcriptome profiling of the flesh and peel of the cucumber cultivar 'YanBai', which is strongly astringent, were examined at three key developmental stages 3, 6 and 9 days post-pollination. The total catechin content decreased as cucumber fruit developed and was significantly lower in the flesh than in the peel. In total, 5092 and 4004 genes were found to be differently expressed in the peel and flesh, respectively. Based on a functional annotation, eight structural genes encode enzymes involved in the catechin biosynthesis pathway. Three genes encoding 4-coumarate-CoA ligases, two genes encoding chalcone isomerases, two genes encoding dihydroflavonol-4-reductase and one gene each encoding a phenylalanine ammonia-lyase, flavanone 3-hydroxylase and cinnamate 4-hydroxylase were identified as affecting the catechin content of cucumber. The transcriptome data also revealed the significance of transcription factors, including WD40-repeat proteins, MYB and bHLH, in regulating catechin biosynthesis. These findings help increase our understanding of the molecular mechanisms controlling catechin biosynthesis and astringency development in cucumber fruit.

A bHLH transcription factor TsMYC2 is associated with the blue grain character in triticale (Triticum × Secale)

RNA-Seq was employed to compare the transcriptome differences between the triticale lines and to identify the key gene responsible for the blue aleurone trait. The accumulation of anthocyanins in the aleurone of triticale results in the formation of the blue-grained trait, but the identity of the genes associated with anthocyanin biosynthesis in the aleurone has not yet been reported. In this manuscript, RNA-Seq was employed to compare the transcriptome differences between the triticale lines HM13 (blue aleurone) and HM5 (white aleurone), and to identify the key genes responsible for the blue aleurone trait. There were 32,406 differentially expressed genes between HM13 and HM5. Seventy-three unigenes were homologous to the structural genes related to anthocyanin biosynthesis, and the average transcript level of the structural genes was higher in HM13 than in HM5, so that quantitative differences between the two lines in transcription rates could be the cause of the blue aleurone. The MYB and bHLH transcription factors had two homologous unigenes, but contained only one differentially expressed unigene each. The relative transcript level of bHLH Unigene5672_All (TsMYC2) in HM13 was 42.71 times that in HM5, while the relative transcript level of the MYB transcription factor Unigene12228_All in HM13 was 2.20 times that in HM5. qPCR experiments determined the relative transcript level of TsMYC2 in developing grain, with the expression of TsMYC2 in grain being the highest compared with that in root, stem or leaf tissue. TsMYC2 was homologous to the bHLH transcription factor regulating anthocyanin biosynthesis and contained three entire functional domains: bHLH-MYC_N, HLH and ACT-like, which were important for exercising regulation of anthocyanin biosynthesis as a bHLH transcription factor. Transient expression of ZmC1 and TsMYC2 could induce anthocyanin biosynthesis in white wheat coleoptile cells, demonstrating that TsMYC2 was a functional bHLH transcription factor. These results indicated that TsMYC2 was associated with the blue aleurone trait and could prove to be a valuable gene with which to breed new triticale cultivars with the blue aleurone trait.

Changes in secondary metabolites, organic acids and soluble sugars during the development of plum fruit cv. 'Furongli' (Prunus salicina Lindl)

BACKGROUND

Organic acids, sugars and pigments are key components that determine the taste and flavor of plum fruit. However, metabolism of organic acid and sugar is not fully understood during the development of plum fruit cv. 'Furongli'.

RESULTS

Mature fruit of 'Furongli' has the highest content of anthocyanins and the lowest content of total phenol compounds and flavonoids. Malate is the predominant organic acid anion in 'Furongli' fruit, followed by citrate and isocitrate. Glucose was the predominant sugar form, followed by fructose and sucrose. Correlation analysis indicated that malate content increased with increasing phosphoenolpyruvate carboxylase (PEPC) activity and decreasing nicotinamide adenine dinucleotide-malate dehydrogenase (NAD-MDH) activity. Citrate and isocitrate content increased with increasing PEPC and aconitase (ACO) activities, respectively. Both acid invertase and neutral invertase had higher activities at the early stage than later stage of fruit development. Fructose content decreased with increasing phosphoglucoisomerase (PGI) activity, whereas glucose content increased with decreasing hexokinase (HK) activity.

CONCLUSION

Dynamics in organic acid anions were not solely controlled by a single enzyme but regulated by the integrated activity of enzymes such as nicotinamide adenine dinucleotide phosphate-malic enzyme (NADP-ME), NAD-ME, PEPC, ACO and NADP-isocitrate dehydrogenase. Sugar metabolism enzymes such as PGI, invertase and HK may play vital roles in the regulation of individual sugar metabolic processes. © 2018 Society of Chemical Industry.

New insights in the control of antioxidants accumulation in tomato by transcriptomic analyses of genotypes exhibiting contrasting levels of fruit metabolites

Background

Tomato is an economically important crop with fruits that are a significant source of bioactive compounds such as ascorbic acid and phenolics. Nowadays, the majority of the enzymes of the biosynthetic pathways and of the structural genes controlling the production and the accumulation of antioxidants in plants are known; however, the mechanisms that regulate the expression of these genes are yet to be investigated. Here, we analyzed the transcriptomic changes occurring during ripening in the fruits of two tomato cultivars (E1 and E115), characterized by a different accumulation of antioxidants, in order to identify candidate genes potentially involved in the biosynthesis of ascorbic acid and phenylpropanoids.

Results

RNA sequencing analyses allowed identifying several structural and regulator genes putatively involved in ascorbate and phenylpropanoids biosynthesis in tomato fruits. Furthermore, transcription factors that may control antioxidants biosynthesis were identified through a weighted gene co-expression network analysis (WGCNA). Results obtained by RNA-seq and WGCNA analyses were further confirmed by RT-qPCR carried out at different ripening stages on ten cultivated tomato genotypes that accumulate different amount of bioactive compounds in the fruit. These analyses allowed us to identify one pectin methylesterase, which may affect the release of pectin-derived D-Galacturonic acid as metabolic precursor of ascorbate biosynthesis. Results reported in the present work allowed also identifying one L-ascorbate oxidase, which may favor the accumulation of reduced ascorbate in tomato fruits. Finally, the pivotal role of the enzymes chalcone synthases (CHS) in controlling the accumulation of phenolic compounds in cultivated tomato genotypes and the transcriptional control of the CHS genes exerted by Myb12 were confirmed.

Conclusions

By using transcriptomic analyses, candidate genes encoding transcription factors and structural genes were identified that may be involved in the accumulation of ascorbic acid and phenylpropanoids in tomato fruits of cultivated genotypes. These analyses provided novel insights into the molecular mechanisms controlling antioxidants accumulation in ripening tomato fruits. The structural genes and regulators here identified could also be used as efficient genetic markers for selecting high antioxidants tomato cultivars.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5428-4) contains supplementary material, which is available to authorized users.

Prunus genetics and applications after de novo genome sequencing: achievements and prospects

Prior to the availability of whole-genome sequences, our understanding of the structural and functional aspects of Prunus tree genomes was limited mostly to molecular genetic mapping of important traits and development of EST resources. With public release of the peach genome and others that followed, significant advances in our knowledge of Prunus genomes and the genetic underpinnings of important traits ensued. In this review, we highlight key achievements in Prunus genetics and breeding driven by the availability of these whole-genome sequences. Within the structural and evolutionary contexts, we summarize: (1) the current status of Prunus whole-genome sequences; (2) preliminary and ongoing work on the sequence structure and diversity of the genomes; (3) the analyses of Prunus genome evolution driven by natural and man-made selection; and (4) provide insight into haploblocking genomes as a means to define genome-scale patterns of evolution that can be leveraged for trait selection in pedigree-based Prunus tree breeding programs worldwide. Functionally, we summarize recent and ongoing work that leverages whole-genome sequences to identify and characterize genes controlling 22 agronomically important Prunus traits. These include phenology, fruit quality, allergens, disease resistance, tree architecture, and self-incompatibility. Translationally, we explore the application of sequence-based marker-assisted breeding technologies and other sequence-guided biotechnological approaches for Prunus crop improvement. Finally, we present the current status of publically available Prunus genomics and genetics data housed mainly in the Genome Database for Rosaceae (GDR) and its updated functionalities for future bioinformatics-based Prunus genetics and genomics inquiry.

Transcriptomic analysis of Lycium ruthenicum Murr. during fruit ripening provides insight into structural and regulatory genes in the anthocyanin biosynthetic pathway

Fruit development in Lycium ruthenicum Murr. involves a succession of physiological and biochemical changes reflecting the transcriptional modulation of thousands of genes. Although recent studies have investigated the dynamic transcriptomic responses during fruit ripening in L. ruthenicum, most have been limited in scope, and thus systematic data representing the structural genes and transcription factors involved in anthocyanin biosynthesis are lacking. In this study, the transcriptomes of three ripening stages associated with anthocyanin accumulation, including S1 (green ripeness stage), S2 (skin color change) and S3 (complete ripeness stage) in L. ruthenicum were investigated using Illumina sequencing. Of a total of 43,573 assembled unigenes, 12,734 were differentially expressed during fruit ripening in L. ruthenicum. Twenty-five significantly differentially expressed structural genes (including PAL, C4H, 4CL, CHS, CHI, F3H, F3'H, F3'5'H, DFR, ANS and UFGT) were identified that might be associated with anthocyanin biosynthesis. Additionally, several transcription factors, including MYB, bHLH, WD40, NAC, WRKY, bZIP and MADS, were correlated with the structural genes, implying their important interaction with anthocyanin biosynthesis-related genes. Our findings provide insight into anthocyanin biosynthesis and regulation patterns in L. ruthenicum and offer a systematic basis for elucidating the molecular mechanisms governing anthocyanin biosynthesis in L. ruthenicum.

Identification of candidate genes involved in anthocyanin accumulation in the peel of jaboticaba (Myrciaria cauliflora) fruits by transcriptomic analysis

Jaboticaba is a grape-like fruit that accumulates high levels of anthocyanins in the peel and is proposed as a good source of functional pigments. However, the molecular mechanisms underlying anthocyanin accumulation in jaboticaba peel remains to be elucidated. In this study, we employed RNA-seq technique to compare the transcriptomic differences between green-colored and black-colored jaboticaba peels. Over 5 million high-quality reads were assembled into 62,190 unigenes with an average length of 737 bp, 29,320 (47.15%) of them were annotated by public databases. 2152 unigenes were found to be differentially expressed (830 upregulated and 1322 downregulated). Gene ontology analysis and pathway enrichment annotation revealed that 18 differentially expressed genes encode phenylalanine ammonialyase, 4-coumaroyl:CoA-ligase, chalcone synthase, flavanone 3-hydroxylase, flavonoid 3'-hydroxylase, anthocyanidin synthase, UDP-glucose: flavonoid 3-O-glucosyltransferase, glutathione S-transferase, Cytochrome b5 were associated with anthocyanin biosynthesis. Additionally, 54 differentially expressed transcription factors were identified. Furthermore, the expression of genes involved in biosynthesis and signal transduction of ethylene and abscisic acid were negatively and positively correlated with that of anthocyanin pathway genes and anthocyanin accumulation, respectively. Quantitative reverse transcription PCR analysis of candidate genes showed trends similar to those in the RNA-seq analysis. McMYB, a homolog of AtMYB113, induced anthocyanin accumulation in tobacco leaves when co-infiltrated PsbHLH3. These results will contribute to further understanding of the molecular mechanisms regulating anthocyanin accumulation in jaboticaba peel.

Transcriptional profiling of catechins biosynthesis genes during tea plant leaf development

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.