Characterization of two O-methyltransferases involved in the biosynthesis of O-methylated catechins in tea plant.

Tea is known for having a high catechin content, with the main component being (-)-epigallocatechin gallate (EGCG), which has significant bioactivities, including potential anti-cancer and anti-inflammatory activity. The poor intestinal stability and permeability of EGCG, however, undermine these health-improving benefits. O-methylated EGCG derivatives, found in a few tea cultivars in low levels, have attracted considerable interest due to their increased bioavailability. Here, we identify two O-methyltransferases from tea plant: CsFAOMT1 that has a specific O-methyltransferase activity on the 3''-position of EGCG to generate EGCG3''Me, and CsFAOMT2 that predominantly catalyzes the formation of EGCG4″Me. In different tea tissues and germplasms, the transcript levels of CsFAOMT1 and CsFAOMT2 are strongly correlated with the amounts of EGCG3''Me and EGCG4''Me, respectively. Furthermore, the crystal structures of CsFAOMT1 and CsFAOMT2 reveal the key residues necessary for 3''- and 4''-O-methylation. These findings may provide guidance for the future development of tea cultivars with high O-methylated catechin content.

Deeply functional identification of TCS1 alleles provides efficient technical paths for low-caffeine breeding of tea plants.

Caffeine is an important functional component in tea, which has the effect of excitement and nerve stimulation, but excessive intake can cause insomnia and dysphoria. Therefore, the production of tea with low-caffeine content can meet the consumption needs of certain people. Here, in addition to the previous alleles of the tea caffeine synthase (TCS1) gene, a new allele (TCS1h) from tea germplasms was identified. Results of in vitro activity analysis showed that TCS1h had both theobromine synthase (TS) and caffeine synthase (CS) activities. Site-directed mutagenesis experiments of TCS1a, TCS1c, and TCS1h demonstrated that apart from the 225th amino acid residue, the 269th amino acid also determined the CS activity. GUS histochemical analysis and dual-luciferase assay indicated the low promoter activity of TCS1e and TCS1f. In parallel, insertion and deletion mutations in large fragments of alleles and experiments of site-directed mutagenesis identified a key cis-acting element (G-box). Furthermore, it was found that the contents of purine alkaloids were related to the expression of corresponding functional genes and alleles, and the absence or presence and level of gene expression determined the content of purine alkaloids in tea plants to a certain extent. In summary, we concluded TCS1 alleles into three types with different functions and proposed a strategy to effectively enhance low-caffeine tea germplasms in breeding practices. This research provided an applicable technical avenue for accelerating the cultivation of specific low-caffeine tea plants.

A novel TcS allele conferring the high-theacrine and low-caffeine traits and having potential use in tea plant breeding.

Theacrine (1,3,7,9-tetramethyluric acid) is a natural product with remarkable pharmacological activities such as antidepressant, sedative and hypnotic activities, while caffeine (1,3,7-trimethylxanthine) has certain side effects to special populations. Hence, breeding tea plants with high theacrine and low caffeine will increase tea health benefits and promote consumption. In this study, we construct an F1 population by crossing 'Zhongcha 302' (theacrine-free) and a tea germplasm 'Ruyuan Kucha' (RY, theacrine-rich) to identify the causal gene for accumulating theacrine. The results showed that the content of theacrine was highly negatively correlated with caffeine (R2 > 0.9). Bulked segregant RNA sequencing analysis, molecular markers and gene expression analysis indicated that the theacrine synthase (TcS) gene was the candidate gene. The TcS was located in the nucleus and cytoplasm, and the theacrine can be detected in stably genetic transformed tobacco by feeding the substrate 1,3,7-trimethyluric acid. Moreover, an in vitro enzyme activity experiment revealed that the 241st amino acid residue was the key residue. Besides, we amplified the promoter region in several tea accessions with varied theacrine levels, and found a 234-bp deletion and a 271-bp insertion in RY. Both GUS histochemical analysis and dual-luciferase assay showed that TcS promoter activity in RY was relatively high. Lastly, we developed a molecular marker that is co-segregate with high-theacrine individuals in RY's offspring. These results demonstrate that the novel TcS allele in RY results in the high-theacrine and low-caffeine traits and the developed functional marker will facilitate the breeding of characteristic tea plants.

Quantitative trait loci mapping for free amino acid content using an albino population and SNP markers provides insight into the genetic improvement of tea plants.

Free amino acids are one of the main chemical components in tea, and they contribute to the pleasant flavor, function, and quality of tea, notably the level of theanine. Here, a high-density genetic map was constructed to characterize quantitative trait loci (QTLs) for free amino acid content. A total of 2688 polymorphic SNP markers were obtained using genotyping-by-sequencing (GBS) based on 198 individuals derived from a pseudotestcross population of "Longjing 43" × "Baijiguan", which are elite and albino tea cultivars, respectively. The 1846.32 cM high-density map with an average interval of 0.69 cM was successfully divided into 15 linkage groups (LGs) ranging from 93.41 cM to 171.28 cM. Furthermore, a total of 4 QTLs related to free amino acid content (theanine, glutamate, glutamine, aspartic acid and arginine) identified over two years were mapped to LG03, LG06, LG11 and LG14. The phenotypic variation explained by these QTLs ranged from 11.8% to 23.7%, with an LOD score from 3.56 to 7.7. Furthermore, several important amino acid metabolic pathways were enriched based on the upregulated differentially expressed genes (DEGs) among the offspring. These results will be essential for fine mapping genes involved in amino acid pathways and diversity, thereby providing a promising avenue for the genetic improvement of tea plants.

Transcriptomic and Metabolomic Analyses Provide Insights Into an Aberrant Tissue of Tea Plant (Camellia sinensis).

Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops with multiple mutants. Recently, we found a special tea germplasm that has an aberrant tissue on its branches. To figure out whether this aberrant tissue is associated with floral bud (FB) or dormant bud (DB), we performed tissue section, transcriptome sequencing, and metabolomic analysis of these tissues. Longitudinal sections indicated the aberrant tissue internal structure was more like a special bud (SB), but was similar to that of DB. Transcriptome data analysis showed that the number of heterozygous and homozygous SNPs was significantly different in the aberrant tissue compared with FB and DB. Further, by aligning the unmapped sequences of the aberrant tissue to the Non-Redundant Protein Sequences (NR) database, we observed that 36.13% of unmapped sequences were insect sequences, which suggested that the aberrant tissue might be a variation of dormant bud tissue influenced by the interaction of tea plants and insects or pathogens. Metabolomic analysis showed that the differentially expressed metabolites (DEMs) between the aberrant tissue and DB were significantly enriched in the metabolic pathways of biosynthesis of plant hormones and biosynthesis of phenylpropanoids. Subsequently, we analyzed the differentially expressed genes (DEGs) in the above mentioned two tissues, and the results indicated that photosynthetic capacity in the aberrant tissue was reduced, whereas the ethylene, salicylic acid and jasmonic acid signaling pathways were activated. We speculated that exogenous infection induced programmed cell death (PCD) and increased the lignin content in dormant buds of tea plants, leading to the formation of this aberrant tissue. This study advanced our understanding of the interaction between plants and insects or pathogens, providing important clues about biotic stress factors and key genes that lead to mutations and formation of the aberrant tissue.

Identifying Conserved Functional Gene Modules Underlying the Dynamic Regulation of Tea Plant Development and Secondary Metabolism.

Tea plants adjust development and metabolism by integrating environmental and endogenous signals in complex but poorly defined gene networks. Here, we present an integrative analysis framework for the identification of conserved modules controlling important agronomic traits using a comprehensive collection of RNA-seq datasets in Camellia plants including 189 samples. In total, 212 secondary metabolism-, 182 stress response-, and 182 tissue development-related coexpressed modules were revealed. Functional modules (e.g., drought response, theobromine biosynthesis, and new shoot development-related modules) and potential regulators that were highly conserved across diverse genetic backgrounds and/or environmental conditions were then identified by cross-experiment comparisons and consensus clustering. Moreover, we investigate the preservation of gene networks between Camellia sinensis and other Camellia species. This revealed that the coexpression patterns of several recently evolved modules related to secondary metabolism and environmental adaptation were rewired and showed higher connectivity in tea plants. These conserved modules are excellent candidates for modeling the core mechanism of tea plant development and secondary metabolism and should serve as a great resource for hypothesis generation and tea quality improvement.

Repressed Gene Expression of Photosynthetic Antenna Proteins Associated with Yellow Leaf Variation as Revealed by Bulked Segregant RNA-seq in Tea Plant Camellia sinensis.

The young leaves and shoots of albino tea cultivars are usually characterized as having a yellow or pale color, high amino acid, and low catechin. Increasing attention has been paid to albino tea cultivars in recent years because their tea generally shows high umami and reduced astringency. However, the genetic mechanism of yellow-leaf variation in albino tea cultivar has not been elucidated clearly. In this study, bulked segregant RNA-seq (BSR-seq) was performed on bulked yellow- and green-leaf hybrid progenies from a leaf color variation population. A total of 359 and 1134 differentially expressed genes (DEGs) were identified in the yellow and green hybrid bulked groups (Yf vs Gf) and parent plants (Yp vs Gp), respectively. The significantly smaller number of DEGs in Yf versus Gf than in Yp versus Gp indicated that individual differences could be reduced within the same hybrid progeny. Analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes revealed that the photosynthetic antenna protein was most significantly enriched in either the bulked groups or their parents. Interaction was found among light-harvesting chlorophyll a/b -binding proteins (LHC), heat shock proteins (HSPs), and enzymes involved in cuticle formation. Combined with the transcriptomic expression profile, results showed that the repressed genes encoding LHC were closely linked to aberrant chloroplast development in yellow-leaf tea plants. Furthermore, the photoprotection and light stress response possessed by genes involved in HSP protein interaction and cuticle formation were discussed. The expression profile of DEGs was verified via quantitative real-time PCR analysis of the bulked samples and other F1 individuals. In summary, using BSR-seq on a hybrid population eliminated certain disturbing effects of genetic background and individual discrepancy, thereby helping this study to intensively focus on the key genes controlling leaf color variation in yellow-leaf tea plants.

Baiyacha, a wild tea plant naturally occurring high contents of theacrine and 3″-methyl-epigallocatechin gallate from Fujian, China.

Baiyacha (BYC) is a kind of wild tea plant growing and utilizing in the remote mountain area of Fujian province, Southeastern China. However, scientific studies on this plant remain limited. Our results showed that BYC exhibits the typical morphological characteristics of Camellia gymnogyna Chang, a closely related species of C. sinensis (L.) O. Kuntze, which was not found in Fujian before. Chemical profiling revealed that parts of BYC plants are rich in purine alkaloids and catechins, especially featuring high levels of theacrine and 3″-methyl-epigallocatechin gallate (EGCG3″Me), chemical compounds with multiple biological activities that are rarely observed in regular tea plants. The contents of EGCG3″Me and theacrine in BYC both increased with the leaf maturity of tea shoots, whereas the caffeine content decreased significantly. The obtained results provide abundant information about the morphology and chemical compounds of BYC and may be used for tea production, breeding, and scientific research in the future.

The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant.

Tea is one of the most popular nonalcoholic beverages due to its characteristic secondary metabolites with numerous health benefits. Although two draft genomes of tea plant (Camellia sinensis) have been published recently, the lack of chromosome-scale assembly hampers the understanding of the fundamental genomic architecture of tea plant and potential improvement. Here, we performed a genome-wide chromosome conformation capture technique (Hi-C) to obtain a chromosome-scale assembly based on the draft genome of C. sinensis var. sinensis and successfully ordered 2984.7 Mb (94.7%) scaffolds into 15 chromosomes. The scaffold N50 of the improved genome was 218.1 Mb, ~157-fold higher than that of the draft genome. Collinearity comparison of genome sequences and two genetic maps validated the high contiguity and accuracy of the chromosome-scale assembly. We clarified that only one Camellia recent tetraploidization event (CRT, 58.9-61.7 million years ago (Mya)) occurred after the core-eudicot common hexaploidization event (146.6-152.7 Mya). Meanwhile, 9243 genes (28.6%) occurred in tandem duplication, and most of these expanded after the CRT event. These gene duplicates increased functionally divergent genes that play important roles in tea-specific biosynthesis or stress response. Sixty-four catechin- and caffeine-related quantitative trait loci (QTLs) were anchored to chromosome assembly. Of these, two catechin-related QTL hotspots were derived from the CRT event, which illustrated that polyploidy has played a dramatic role in the diversification of tea germplasms. The availability of a chromosome-scale genome of tea plant holds great promise for the understanding of genome evolution and the discovery of novel genes contributing to agronomically beneficial traits in future breeding programs.

Gene Coexpression Networks Reveal Key Drivers of Flavonoid Variation in Eleven Tea Cultivars (Camellia sinensis).

Following the recent completion of the draft genome sequence of the tea plant, high-throughput decoding of gene function, especially for those involved in complex secondary metabolic pathways, has become a major challenge. Here, we profiled the metabolome and transcriptome of 11 tea cultivars, and then illustrated a weighted gene coexpression network analysis (WGCNA)-based system biological strategy to interpret metabolomic flux, predict gene functions, and mine key regulators involved in the flavonoid biosynthesis pathway. We constructed a multilayered regulatory network, which integrated the gene coexpression relationship with the microRNA target and promoter cis-regulatory element information. This allowed us to reveal new uncharacterized TFs (e.g., MADSs, WRKYs, and SBPs) and microRNAs (including 17 conserved and 15 novel microRNAs) that are potentially implicated in different steps of the catechin biosynthesis. Furthermore, we applied metabolic-signature-based association method to capture additional key regulators involved in catechin pathway. This provides important clues for the functional characterization of five SCPL1A acyltransferase family members, which might be implicated in the production balance of anthocyanins, galloylated catechins, and proanthocyanins. Application of an "omics"-based system biology strategy should facilitate germplasm utilization and provide valuable resources for tea quality improvement.

Isolation and Characterization of CsWRKY7, a Subgroup IId WRKY Transcription Factor from Camellia sinensis, Linked to Development in Arabidopsis.

WRKY transcription factors (TFs) containing one or two WRKY domains are a class of plant TFs that respond to diverse abiotic stresses and are associated with developmental processes. However, little has been known about the function of WRKY gene in tea plant. In this study, a subgroup IId WRKY gene CsWRKY7 was isolated from Camellia sinensis, which displayed amino acid sequence homology with Arabidopsis AtWRKY7 and AtWRKY15. Subcellular localization prediction indicated that CsWRKY7 localized to nucleus. Cis-acting elements detected in the promotor region of CsWRKY7 are mainly involved in plant response to environmental stress and growth. Consistently, expression analysis showed that CsWRKY7 transcripts responded to NaCl, mannitol, PEG, and diverse hormones treatments. Additionally, CsWRKY7 exhibited a higher accumulation both in old leaves and roots compared to bud. Seed germination and root growth assay indicated that overexpressed CsWRKY7 in transgenic Arabidopsis was not sensitive to NaCl, mannitol, PEG, and low concentration of ABA treatments. CsWRKY7 overexpressing Arabidopsis showed a late-flowering phenotype under normal conditions compared to wild type. Furthermore, gene expression analysis showed that the transcription levels of the flowering time integrator gene FLOWERING LOCUS T (FT) and the floral meristem identity genes APETALA1 (AP1) and LEAFY (LFY) were lower in WRKY7-OE than in the WT. Taken together, these findings indicate that CsWRKY7 TF may participate in plant growth. This study provides a potential strategy to breed late-blooming tea cultivar.

Quantitative Trait Loci Mapping for Theobromine and Caffeine Contents in Tea Plant ( Camellia sinensis).

Understanding the genetic basis of theobromine and caffeine accumulation in the tea plant is important due to their contribution to tea flavor. Quantitative trait loci (QTL) analyses were carried out to identify genetic variants associated with theobromine and caffeine contents and ratio using a pseudo-testcross population derived from an intervarietal cross between two varieties of Camellia sinensis. A total of 10 QTL controlling caffeine content (CAF), theobromine content (TBR), sum of caffeine and theobromine (SCT), and caffeine-to-theobromine ratio (CTR) were identified over four measurement years. The major QTL controlling CAF, qCAF1, was mapped onto LG01 and validated across years, explaining an average of 20.1% of the phenotypic variance. The other QTL were detected in 1 or 2 years, and of them there were four, two, and three for TBR, SCT, and CTR, respectively. The present results provide valuable information for further fine mapping and cloning functional genes and for genetic improvement in tea plant.

Hongyacha, a Naturally Caffeine-Free Tea Plant from Fujian, China.

Hongyacha (HYC) is a type of new wild tea plant discovered in Fujian Province, China. This tea is helpful to the healing or prevention of disease in its original growing area. However, research on this tea is limited. Our results showed that HYC displayed obvious differences in its morphological characteristics compared with Cocoa tea ( Camellia ptilophylla Chang), a famous caffeine-free tea plant in China. Theobromine and trans-catechins, but not caffeine and cis-catechins, were the dominant purine alkaloids and catechins detected in HYC. HYC might contain abundant gallocatechin-(4 → 8)-gallocatechin gallate, 1,3,4,6-tetra- O-galloyl-ß-d-glucopyranose, and (-)-gallocatechin-3,5-di- O-gallate, which were not detected in regular tea. We also found that the TCS1 of HYC was distinct, and the responding recombinant protein exhibited only theobromine synthase activity. The obtained results showed that HYC is a new kind of caffeine-free tea plant and may be used for scientific protection and efficient utilization in the future.

A Novel F3' 5' H Allele with 14 bp Deletion Is Associated with High Catechin Index Trait of Wild Tea Plants and Has Potential Use in Enhancing Tea Quality.

Catechins are important chemical components determining the quality of tea. The catechin index (CI, ratio of dihydroxylated catechin (DIC)/trihydroxylated catechin (TRIC)) in the green leaf has a major influence on the amounts of theaflavins in black tea. In this work, the major catechin profiles of wild tea plants originating from Guizhou Province with high CI trait were investigated. We identified a novel flavonoid 3',5' hydroxylase gene ( F3' 5' H) allele with a 14 bp deletion in the upstream regulation region and developed an insertion/deletion (InDel) marker accordingly. The 14 bp deletion in the novel  F3' 5' H allele was associated with low F3' 5' H mRNA expression, thereby resulting in low TRIC content and high CI value. The allelic variant in the novel F3' 5' H allele associated with high CI values and DIC contents was confirmed by the introgression lines derived from a distant cross population. The novel F3' 5' H allele in wild tea plants is a valuable gene resource, which could be applied to breeding improvement on tea quality.

Comprehensive Dissection of Metabolic Changes in Albino and Green Tea Cultivars.

Albino tea cultivars are special mutants of tea plants with white or yellow leaf color. In this study, three albino tea cultivars, including 'Anji Baicha', 'Huangjinya', and 'Baijiguan', and two green tea cultivars, 'Longjing 43' and 'Fuding Dabaicha', were applied to metabolite profiling by gas chromatography-mass spectrometry and ultraperformance liquid chromatography-mass spectrometry. Multivariate analyses revealed significantly different metabolite phenotypes in leaves among albino cultivars and green cultivars. The differential metabolite-related pathways included galactose metabolism, tryptophan metabolism, phenylpropanoid biosynthesis, and flavonoid biosynthesis. For the young leaves of albino cultivars, the sugar (sorbitol and erythrose) and amino acid (mainly proline, isoleucine, ornithine, aspartic acid, threonine, and valine) concentrations increased, whereas gallocatechin and epigallocatechin gallate concentrations decreased. These results reveal the divergence in metabolic profiling between tea plant cultivars with different leaf colors. With the development of leaves, the concentrations of flavonoids increased largely in the older leaves of albino cultivars.

Functional natural allelic variants of flavonoid 3',5'-hydroxylase gene governing catechin traits in tea plant and its relatives.

MAIN CONCLUSION: Functional allelic variants of the flavonoid 3',5'-hydroxylase (F3'5'H) gene provides new information of F3'5'H function of tea plant and its relatives. This insight may serve as the foundation upon which to advance molecular breeding in the tea plant. Catechins are the active components of tea that determine its quality and health attributes. This study established the first integrated genomic strategy for deciphering the genetic basis of catechin traits of tea plant. With the RNA-sequencing analysis of bulked segregants representing the tails of a F1 population segregated for total catechin content, we identified a flavonoid 3',5'-hydroxylase (F3'5'H) gene. F3'5'H had one copy in the genomic DNA of tea plant. Among 202 tea accessions, we identified 120 single nucleotide polymorphisms (SNPs) at F3'5'H locus. Seventeen significant marker-trait associations were identified by association mapping in multiple environments, which were involved in 10 SNP markers, and the traits including the ratio of di/tri-hydroxylated catechins and catechin contents. The associated individual and combination of SNPs explained 4.5-25.2 and 53.0-63.0% phenotypic variations, respectively. In the F1 population (validation population), the catechin trait variation percentages explained by F3'5'H diplotype were 6.9-74.3%. The genotype effects of ten functional SNPs in the F1 population were all consistent with the association population. Furthermore, the function of SNP-711/-655 within F3'5'H was validated by gene expression analysis. Altogether, our work indicated functional SNP allelic variants within F3'5'H governing the ratio of di/tri-hydroxylated catechins and catechin contents. The strong catechin-associated SNPs identified in this study can be used for future marker-assisted selection to improve tea quality.

Biochemical and transcriptomic analyses reveal different metabolite biosynthesis profiles among three color and developmental stages in 'Anji Baicha' (Camellia sinensis).

BACKGROUND: The new shoots of the albino tea cultivar 'Anji Baicha' are yellow or white at low temperatures and turn green as the environmental temperatures increase during the early spring. 'Anji Baicha' metabolite profiles exhibit considerable variability over three color and developmental stages, especially regarding the carotenoid, chlorophyll, and theanine concentrations. Previous studies focused on physiological characteristics, gene expression differences, and variations in metabolite abundances in albino tea plant leaves at specific growth stages. However, the molecular mechanisms regulating metabolite biosynthesis in various color and developmental stages in albino tea leaves have not been fully characterized. RESULTS: We used RNA-sequencing to analyze 'Anji Baicha' leaves at the yellow-green, albescent, and re-greening stages. The leaf transcriptomes differed considerably among the three stages. Functional classifications based on Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that differentially expressed unigenes were mainly related to metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, and carbon fixation in photosynthetic organisms. Chemical analyses revealed higher ß-carotene and theanine levels, but lower chlorophyll a levels, in the albescent stage than in the green stage. Furthermore, unigenes involved in carotenoid, chlorophyll, and theanine biosyntheses were identified, and the expression patterns of the differentially expressed unigenes in these biosynthesis pathways were characterized. Through co-expression analyses, we identified the key genes in these pathways. These genes may be responsible for the metabolite biosynthesis differences among the different leaf color and developmental stages of 'Anji Baicha' tea plants. CONCLUSIONS: Our study presents the results of transcriptomic and biochemical analyses of 'Anji Baicha' tea plants at various stages. The distinct transcriptome profiles for each color and developmental stage enabled us to identify changes to biosynthesis pathways and revealed the contributions of such variations to the albino phenotype of tea plants. Furthermore, comparisons of the transcriptomes and related metabolites helped clarify the molecular regulatory mechanisms underlying the secondary metabolic pathways in different stages.

Association mapping of caffeine content with TCS1 in tea plant and its related species.

Caffeine is the most abundant purine alkaloid in majority of tea plant and its related species. This purine alkaloid contributes to the important flavor and health attributes of tea. Tea caffeine synthase 1 (TCS1, EC 2.1.1.159/2.1.1.160) gene plays a crucial role in caffeine biosynthesis. The objective of this study was to investigate the genetic relationship between the TCS1 and caffeine content of tea plant and its related species using association mapping. We identified 87 single-nucleotide polymorphisms (SNPs, π = 0.00447) by resequencing the TCS1 locus of 44 tea accessions. Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r(2) < 0.1, within 1000 bp). Two cleaved amplified polymorphism sequence (CAPS) markers were developed from sequence variations (SNP4318 and SNP6252). By association mapping, we identified SNP4318 associated with caffeine content in four environments, explaining 4.0%-7.7% of the phenotypic variance. We also validated the significant marker-trait associations in site-directed mutagenesis experiments. Examination of allelic variation and linkage disequilibrium by a candidate-gene-based approach can help to decipher the genetic basis of caffeine biosynthesis. Moreover, the SNP marker identified in this study can potentially be applied for future marker-assisted selection to improve tea quality.

Transcriptomic Analysis of Tea Plant Responding to Drought Stress and Recovery.

Tea plant (Camellia sinensis) is an economically important beverage crop. Drought stress (DS) seriously limits the growth and development of tea plant, thus affecting crop yield and quality. To elucidate the molecular mechanisms of tea plant responding to DS, we performed transcriptomic analysis of tea plant during the three stages [control (CK) and during DS, and recovery (RC) after DS] using RNA sequencing (RNA-Seq). Totally 378.08 million high-quality trimmed reads were obtained and assembled into 59,674 unigenes, which were extensively annotated. There were 5,955 differentially expressed genes (DEGs) among the three stages. Among them, 3,948 and 1,673 DEGs were up-regulated under DS and RC, respectively. RNA-Seq data were further confirmed by qRT-PCR analysis. Genes involved in abscisic acid (ABA), ethylene, and jasmonic acid biosynthesis and signaling were generally up-regulated under DS and down-regulated during RC. Tea plant potentially used an exchange pathway for biosynthesis of indole-3-acetic acid (IAA) and salicylic acid under DS. IAA signaling was possibly decreased under DS but increased after RC. Genes encoding enzymes involved in cytokinin synthesis were up-regulated under DS, but down-regulated during RC. It seemed probable that cytokinin signaling was slightly enhanced under DS. In total, 762 and 950 protein kinases belonging to 26 families were differentially expressed during DS and RC, respectively. Overall, 547 and 604 transcription factor (TF) genes belonging to 58 families were induced in the DS vs. CK and RC vs. DS libraries, respectively. Most members of the 12 TF families were up-regulated under DS. Under DS, genes related to starch synthesis were down-regulated, while those related to starch decomposition were up-regulated. Mannitol, trehalose and sucrose synthesis-related genes were up-regulated under DS. Proline was probably mainly biosynthesized from glutamate under DS and RC. The mechanism by which ABA regulated stomatal movement under DS and RC was partly clarified. These results document the global and novel responses of tea plant during DS and RC. These data will serve as a valuable resource for drought-tolerance research and will be useful for breeding drought-resistant tea cultivars.