Search

The unique zigzag-patterned tea plant is a rare germplasm resource. However, the molecular mechanism behind the formation of zigzag stems remains unclear. To address this, a BC1 genetic population of tea plants with zigzag stems was studied using histological observation and bulked segregant RNA-seq. The analysis revealed 1494 differentially expressed genes (DEGs) between the upright and zigzag stem groups. These DEGs may regulate the transduction and biosynthesis of plant hormones, and the effects on the phenylpropane biosynthesis pathways may cause the accumulation of lignin. Tissue sections further supported this finding, showing differences in cell wall thickness between upright and curved stems, potentially due to lignin accumulation. Additionally, 262 single-nucleotide polymorphisms (SNPs) across 38 genes were identified as key SNPs, and 5 genes related to zigzag stems were identified through homologous gene function annotation. Mutations in these genes may impact auxin distribution and content, resulting in the asymmetric development of vascular bundles in curved stems. In summary, we identified the key genes associated with the tortuous phenotype by using BSR-seq on a BC1 population to minimize genetic background noise.

Among all types of tea, black tea is produced in the largest amount worldwide, and its consumption is still increasing. Enzymatic fermentation is considered majorly contribute to quality formation of black tea, and the information about the roles of bacterial community in black tea processing is scarce. This study aimed to analyze the dynamic changes in composition, structure, and function of microbial communities during black tea processing and reveal the roles of bacterial community in black tea processing. Results showed that the genera Sphingomonas and Variovorax were dominant throughout the processing of black tea. Prediction function analysis of bacterial community showed that the mean proportions of glucuronoarabinoxylan endo - 1,4 - beta - xylanase, aminopeptidase B, phosphoserine phosphatase, homoserine O-acetyltransferase, glycolysis related enzymes, pyruvate dehydrogenase, tricarboxylic acid cycle related enzymes, and glyoxylate bypass were significantly elevated in the rolling and fermentation stages. The contents of amino acids, soluble sugar, theaflavins, thearubigins, and theabrownins increased greatly during the rolling and fermentation processes. Redundancy and Pearson's correlation analyses showed that the relative abundance of bacteria was closely related to the contents of water extract, tea polyphenols, epigallocatechin, epicatechin gallate, catechin gallate, thearubigins, theaflavins, and theabrownins. Overall, the findings provided new insights into the variation of bacterial community during black tea processing and improved our understanding of the core functional bacteria involved in black tea processing., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

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.

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 (Y f vs G f ) and parent plants (Y p vs G p ), respectively. The significantly smaller number of DEGs in Y f versus G f than in Y p versus G p 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 F 1 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.

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.

The effects of blue (BL) and green light (GL) treatment during the dark period were examined in Camellia sinensis as a first step to understanding the spectral effects of artificial BL and GL on plant secondary metabolism and light signaling interactions. BL could induce the expression of CRY2/3, SPAs, HY5, and R2R3-MYBs to promote the accumulation of anthocyanins and catechins in tea plants. GL, on the other hand, could stimulate the accumulation of several functional substances (e.g., procyanidin B2/B3 and l-ascorbate) and temper these BL responses via down-regulation of  CRY2/3 and PHOT2. Furthermore, the molecular events that triggered by BL and GL signals were partly overlapped with abiotic/biotic stress responses. We indicate the possibility of a targeted use of BL and GL to regulate the amount of functional metabolites to enhance tea quality and taste, and to potentially trigger defense mechanisms of tea plants.

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.

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.

Tea plant (Camellia sinensis (L) O. Kuntze) respond to herbivore attack through large changes in defense related metabolism and gene expression. Ectropis oblique (Prout) is one of the most devastating insects that feed on tea leaves and tender buds, which can cause severe production loss and deteriorate the quality of tea. To elucidate the biochemicals and molecular mechanism of defense against tea geometrid (TG), transcriptome and metabolome of TG interaction with susceptible (SG) and resistance (RG) tea genotypes were analyzed by using UPLC-Q-TOF-MS, GC-MS, and RNA-seq technologies. This revealed that jasmonic acid was highly induced in RG, following a plethora of secondary metabolites involved in defense against TG could be induced by jasmonic acid signaling pathway. However, the constitutively present of salicylic acid in SG might be a suppressor of jasmonate signaling and thus misdirect tea plants against TG. Furthermore, flavonoids and terpenoids biosynthesis pathways were highly activated in RG to constitute the chemical barrier on TG feeding behavior. In contrast, fructose and theanine, which can act as feeding stimulants were observed to highly accumulate in SG. Being present in the major hub, 39 transcription factors or protein kinases among putative candidates were identified as master regulators from protein-protein interaction network analysis. Together, the current study provides a comprehensive gene expression and metabolite profiles, which can shed new insights into the molecular mechanism of tea defense against TG. The candidate genes and specific metabolites identified in the present study can serve as a valuable resource for unraveling the possible defense mechanism of plants against various biotic stresses., Competing Interests: The authors have declared that no competing interests exist.

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.

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.

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.

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., (Copyright © 2016 Elsevier Masson SAS. All rights reserved.)

Tea caffeine synthase 1 (TCS1) is an enzyme that catalyzes the methylation of N-3 and N-1 and considered to be the most critical enzyme in the caffeine biosynthetic pathway of tea plant. This study shows that TCS1 has six types of allelic variations, namely, TCS1a, TCS1b, TCS1c, TCS1d, TCS1e, and TCS1f, with a 252 bp insertion/deletion mutation in the 5'-untranslated region. Among tea plant and its related species, TCS1a is the predominant allele, and TCS1b-f are the rare alleles that mainly appear in few wild germplasms. The full-length cDNA sequences of three new alleles, namely, TCS1d, TCS1e, and TCS1f, were isolated from specific germplasms, and all of recombinant proteins have higher caffeine synthase (CS, EC 2.1.1.160) activity than theobromine synthase (TS, EC 2.1.1.159). Amino acid residue 269 is responsible for the difference in TCS activity and substrate recognition, which was demonstrated by using site-directed mutagenesis experiments. Furthermore, natural variations in TCS1 change the transcription levels. There are two molecular mechanisms controlling the caffeine biosynthesis in low-caffeine-accumulating tea germplasms, i.e., TCS1 allele with low transcription level or its encoded protein with only TS activity. Allelic variations of TCS1 play a crucial role in caffeine biosynthesis. Taken together, our work provides valuable foundation for a comprehensive understanding of the mechanism of caffeine biosynthesis in section Thea plants and useful guidance for effective breeding., (Copyright © 2016 Elsevier Masson SAS. All rights reserved.)

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.

'Anji Baicha' is an albino tea cultivar with white shoots at low air temperature and green shoots at high air temperature in early spring. The metabolite contents in the shoots dynamically vary with the color changes and with shoot development. To investigate the metabolomic variation during the albescent and re-greening stages, gas chromatography-mass spectrometry combined with multivariate analysis were applied to analyze the metabolite profiles in the different color stages during the development of 'Anji Baicha' leaves. The metabolite profiles of three albescent stages, including the yellow-green stage, the early albescent stage, and the late albescent stage, as well as the re-greening stage were distinguished using principal component analysis, revealing that the distinct developmental stages were likely responsible for the observed metabolic differences. Furthermore, a group classification and pairwise discrimination was revealed among the three albescent stages and re-greening stage by partial least squares discriminant analysis. A total of 65 differential metabolites were identified with a variable influence on projection greater than 1. The main differential metabolic pathways of the albescent stages compared with the re-greening stage included carbon fixation in photosynthetic organisms and the phenylpropanoid and flavonoid biosynthesis pathways. Compared with the re-greening stage, the carbohydrate and amino acid metabolic pathways were disturbed during the albescent stages. During the albescent stages, the sugar (fructofuranose), sugar derivative (glucose-1-phosphate) and epicatechin concentrations decreased, whereas the amino acid (mainly glycine, serine, tryptophan, citrulline, glutamine, proline, and valine) concentrations increased. These results reveal the changes in metabolic profiling that occur during the color changes associated with the development of the albino tea plant leaves.

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.

To understand tea germplasms better and to use them effectively for production and breeding, the catechin content of 403 accessions of representative tea germplasms collected from various locations in China were studied using HPLC. The catechin content of these tea germplasms varied from 56.6 to 231.9 mg/g and averaged 154.5 ± 18.1 mg/g. One germplasm with low total catechin (TC) content (<60 mg/g) and three with high TC (>200 mg/g) contents were found. Averages of the TC content of the three varieties of Camellia sinensis (L.) O. Kuntze, namely, sinensis, assamica, and pubilimba, were 152.9 ± 16.2 mg/g, 162.8 ± 22.3 mg/g, and 165.1 ± 21.3 mg/g, respectively. The TC content of the sinensis variety was significantly lower (P < 0.05) than that of the other two varieties. The assamica variety had the highest levels of (-)-epicatechin gallate (ECG), and (-)-epicatechin (EC), whereas the pubilimba variety had the highest levels of (-)-epigallocatechin gallate (EGCG), (+)-gallocatechin (GC), (+)-catechin (C), and (-)-gallocatechin gallate (GCG). Factor analysis indicated that GC, C, GCG, catechin index, and ECG greatly influenced the classification. The TC content of germplasms collected from the various provinces showed significant differences (P < 0.05). Tea germplasms of the southern provinces had higher degrees of variation in TC.

Catechins are the most important bioactive compounds in tea, and have been demonstrated to possess a wide variety of pharmacological activities. To characterize quantitative trait loci (QTLs) for catechins content in the tender shoots of tea plant, we constructed a moderately saturated genetic map using 406 simple sequence repeat (SSR) markers, based on a pseudo-testcross population of 183 individuals derived from an intraspecific cross of two Camellia sinensis varieties with diverse catechins composition. The map consisted of fifteen linkage groups (LGs), corresponding to the haploid chromosome number of tea plant (2n = 2x = 30). The total map length was 1,143.5 cM, with an average locus spacing of 2.9 cM. A total of 25 QTLs associated with catechins content were identified over two measurement years. Of these, nine stable QTLs were validated across years, and clustered into four main chromosome regions on LG03, LG11, LG12 and LG15. The population variability explained by each QTL was predominantly at moderate-to-high levels and ranged from 2.4% to 71.0%, with an average of 17.7%. The total number of QTL for each trait varied from four to eight, while the total population variability explained by all QTLs for a trait ranged between 38.4% and 79.7%. This is the first report on the identification of QTL for catechins content in tea plant. The results of this study provide a foundation for further cloning and functional characterization of catechin QTLs for utilization in improvement of tea plant.

Background: Tea is the most popular non-alcoholic health beverage in the world. The tea plant (Camellia sinensis (L.) O. Kuntze) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in tea plants. To elucidate the molecular mechanisms of cold acclimation, we employed RNA-Seq and digital gene expression (DGE) technologies to the study of genome-wide expression profiles during cold acclimation in tea plants., Results: Using the Illumina sequencing platform, we obtained approximately 57.35 million RNA-Seq reads. These reads were assembled into 216,831 transcripts, with an average length of 356 bp and an N50 of 529 bp. In total, 1,770 differentially expressed transcripts were identified, of which 1,168 were up-regulated and 602 down-regulated. These include a group of cold sensor or signal transduction genes, cold-responsive transcription factor genes, plasma membrane stabilization related genes, osmosensing-responsive genes, and detoxification enzyme genes. DGE and quantitative RT-PCR analysis further confirmed the results from RNA-Seq analysis. Pathway analysis indicated that the "carbohydrate metabolism pathway" and the "calcium signaling pathway" might play a vital role in tea plants' responses to cold stress., Conclusions: Our study presents a global survey of transcriptome profiles of tea plants in response to low, non-freezing temperatures and yields insights into the molecular mechanisms of tea plants during the cold acclimation process. It could also serve as a valuable resource for relevant research on cold-tolerance and help to explore the cold-related genes in improving the understanding of low-temperature tolerance and plant-environment interactions.