Ye-Yun Li

Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds.

BackgroundTea is one of the most popular non-alcoholic beverages worldwide. However, the tea plant, Camellia sinensis, is difficult to culture in vitro, to transform, and has a large genome, rendering little genomic information available. Recent advances in large-scale RNA sequencing (RNA-seq) provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.ResultsUsing high-throughput Illumina RNA-seq, the transcriptome from poly (A)+ RNA of C. sinensis was analyzed at an unprecedented depth (2.59 gigabase pairs). Approximate 34.5 million reads were obtained, trimmed, and assembled into 127,094 unigenes, with an average length of 355 bp and an N50 of 506 bp, which consisted of 788 contig clusters and 126,306 singletons. This number of unigenes was 10-fold higher than existing C. sinensis sequences deposited in GenBank (as of August 2010). Sequence similarity analyses against six public databases (Uniprot, NR and COGs at NCBI, Pfam, InterPro and KEGG) found 55,088 unigenes that could be annotated with gene descriptions, conserved protein domains, or gene ontology terms. Some of the unigenes were assigned to putative metabolic pathways. Targeted searches using these annotations identified the majority of genes associated with several primary metabolic pathways and natural product pathways that are important to tea quality, such as flavonoid, theanine and caffeine biosynthesis pathways. Novel candidate genes of these secondary pathways were discovered. Comparisons with four previously prepared cDNA libraries revealed that this transcriptome dataset has both a high degree of consistency with previous EST data and an approximate 20 times increase in coverage. Thirteen unigenes related to theanine and flavonoid synthesis were validated. Their expression patterns in different organs of the tea plant were analyzed by RT-PCR and quantitative real time PCR (qRT-PCR).ConclusionsAn extensive transcriptome dataset has been obtained from the deep sequencing of tea plant. The coverage of the transcriptome is comprehensive enough to discover all known genes of several major metabolic pathways. This transcriptome dataset can serve as an important public information platform for gene expression, genomics, and functional genomic studies in C. sinensis.

CsICE1 and CsCBF1: two transcription factors involved in cold responses in Camellia sinensis

C-repeat/dehydration-responsive element binding factors (CBFs) can induce the expression of a suite of cold-responsive genes to increase plant cold tolerance, and inducer of CBF expression 1 (ICE1) is a major activator for CBF. In the present study, we isolated the full-length cDNAs of ICE1 and CBF from Camellia sinensis, designated as CsICE1 and CsCBF1, respectively. The deduced protein CsICE1 contains a highly conserved basic helix-loop-helix (bHLH) domain and C-terminal region of ICE1-like proteins. CsCBF1 contains all conserved domains of CBFs in other plant species and can specifically bind to the C-repeat/dehydration-responsive element (CRT/DRE) as confirmed by electrophoretic mobility shift assay. The transcription of CsICE1 had no apparent alteration after chilling treatment (4°C). CsCBF1 expression was not detected in normal temperature (20°C) but was induced immediately and significantly by low temperature (4°C). Our results suggest that ICE1–CBF cold-response pathway is conserved in tea plants. CsICE1 and CsCBF1, two components of this pathway, play roles in cold responses in tea plants.

Molecular cloning, functional analysis of three cinnamyl alcohol dehydrogenase (CAD) genes in the leaves of tea plant, Camellia sinensis.

Cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.195) is considered to be a key enzyme in lignin biosynthesis, but little was known about CADs in tea plants (Camellia sinensis). A full-length cDNA sequence (CsCAD2) was isolated by suppressive subtractive hybridization (SSH) in Ectropis oblique feeding-induced tea plants, and another two full-length cDNA sequences (CsCAD1 and CsCAD3) were obtained from a transcriptome obtained by deep sequencing. However, they showed only 20-54% identities. Phylogenetic analysis revealed that they belonged to three different families. DNA gel blotting analysis revealed that two copies of CsCAD1 and CsCAD2 genes existed in tea genome, but CsCAD3 likely had only one copy. Recombinant proteins of these CsCADs were produced in Escherichia coli. The activity of purified recombinant CsCAD2 protein was up to 0.43 μmol min(-1) mg(-1). However, the other two recombinant proteins had lower activities, probably due to incomplete refolding. qRT-PCR analysis indicated that while CsCAD3 was strongly up-regulated in tea plants after E. oblique attack and mechanical damage, CsCAD1 and CsCAD2 showed only moderate or no changes in transcript levels. Treatment of defence-related hormones methyl jasmonate (MeJA) and salicylic acid (SA) elevated the expression of CsCAD1 and CsCAD2, but decreased the transcript abundance of CsCAD3. The transcript levels of CsCAD2 did not change after applying abscisic acid (ABA), whereas CsCAD1 and CsCAD3 were induced. These results suggested that these three CsCAD genes in tea plants may play a role in defense against insects and pathogens and adaptation to abiotic stresses and these genes likely have divergant functions.

Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality

Significance A high-quality genome assembly of Camellia sinensis var. sinensis facilitates genomic, transcriptomic, and metabolomic analyses of the quality traits that make tea one of the world’s most-consumed beverages. The specific gene family members critical for biosynthesis of key tea metabolites, monomeric galloylated catechins and theanine, are indicated and found to have evolved specifically for these functions in the tea plant lineage. Two whole-genome duplications, critical to gene family evolution for these two metabolites, are identified and dated, but are shown to account for less amplification than subsequent paralogous duplications. These studies lay the foundation for future research to understand and utilize the genes that determine tea quality and its diversity within tea germplasm. Tea, one of the world’s most important beverage crops, provides numerous secondary metabolites that account for its rich taste and health benefits. Here we present a high-quality sequence of the genome of tea, Camellia sinensis var. sinensis (CSS), using both Illumina and PacBio sequencing technologies. At least 64% of the 3.1-Gb genome assembly consists of repetitive sequences, and the rest yields 33,932 high-confidence predictions of encoded proteins. Divergence between two major lineages, CSS and Camellia sinensis var. assamica (CSA), is calculated to ∼0.38 to 1.54 million years ago (Mya). Analysis of genic collinearity reveals that the tea genome is the product of two rounds of whole-genome duplications (WGDs) that occurred ∼30 to 40 and ∼90 to 100 Mya. We provide evidence that these WGD events, and subsequent paralogous duplications, had major impacts on the copy numbers of secondary metabolite genes, particularly genes critical to producing three key quality compounds: catechins, theanine, and caffeine. Analyses of transcriptome and phytochemistry data show that amplification and transcriptional divergence of genes encoding a large acyltransferase family and leucoanthocyanidin reductases are associated with the characteristic young leaf accumulation of monomeric galloylated catechins in tea, while functional divergence of a single member of the glutamine synthetase gene family yielded theanine synthetase. This genome sequence will facilitate understanding of tea genome evolution and tea metabolite pathways, and will promote germplasm utilization for breeding improved tea varieties.

Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing

Tea is one of the most popular beverages across the world and is made exclusively from cultivars of Camellia sinensis. Many wild relatives of the genus Camellia that are closely related to C. sinensis are native to Southwest China. In this study, we first identified the distinct genetic divergence between C. sinensis and its wild relatives and provided a glimpse into the artificial selection of tea plants at a genome-wide level by analyzing 15,444 genomic SNPs that were identified from 18 cultivated and wild tea accessions using a high-throughput genome-wide restriction site-associated DNA sequencing (RAD-Seq) approach. Six distinct clusters were detected by phylogeny inferrence and principal component and genetic structural analyses, and these clusters corresponded to six Camellia species/varieties. Genetic divergence apparently indicated that C. taliensis var. bangwei is a semi-wild or transient landrace occupying a phylogenetic position between those wild and cultivated tea plants. Cultivated accessions exhibited greater heterozygosity than wild accessions, with the exception of C. taliensis var. bangwei. Thirteen genes with non-synonymous SNPs exhibited strong selective signals that were suggestive of putative artificial selective footprints for tea plants during domestication. The genome-wide SNPs provide a fundamental data resource for assessing genetic relationships, characterizing complex traits, comparing heterozygosity and analyzing putatitve artificial selection in tea plants.

Molecular Cloning and Characterization of Hydroperoxide Lyase Gene in the Leaves of Tea Plant (Camellia sinensis).

Hydroperoxide lyase (HPL, E.C. 4.1.2.) is the major enzyme in the biosynthesis of natural volatile aldehydes and alcohols in plants, however, little was known about HPL in tea plants (Camellia sinensis). A unique cDNA fragment was isolated by suppressive subtractive hybridization (SSH) from a tea plant subjected to herbivory by tea geometrid Ectropis obliqua. This full length cDNA acquired by RACE was 1476 bp and encoded 491 amino acids. DNA and protein BLAST searches showed high homology to HPL sequences from other plants. The His-tag expression vector pET-32a(+)/CsHPL was constructed and transferred into Escherichia coli Rosetta (DE3). The expression product of recombinant CsHPL in E. coli was about 60 kDa. The enzyme activity of CsHPL was 0.20 μmol·min(-1)·mg(-1). Quantitative RT-PCR analysis indicated CsHPL was strongly up-regulated in tea plants after Ectropis obliqua attack, suggesting that it may be an important candidate for defense against insects in tea plants.

Comparative analysis of the response and gene regulation in cold resistant and susceptible tea plants

Cold environment is the main constraint for tea plants (Camellia sinensis) distribution and tea farming. We identified two tea cultivars, called var. sinensis cv. Shuchazao (SCZ) with a high cold-tolerance and var. assamica cv. Yinghong9 (YH9) with low cold-tolerance. To better understand the response mechanism of tea plants under cold stress for improving breeding, we compared physiological and biochemical responses, and associated genes expression in response to 7-day and 14-day cold acclimation, followed by 7-day de-acclimation in these two tea cultivars. We found that the low EL50, low Fv/Fm, and high sucrose and raffinose accumulation are responsible for higher cold tolerance in SCZ comparing with YH9. We then measured the expression of 14 key homologous genes, known as involved in these responses in other plants, for each stages of treatment in both cultivars using RT-qPCR. Our results suggested that the increased expression of CsCBF1 and CsDHNs coupling with the accumulation of sucrose play key roles in conferring higher cold resistance in SCZ. Our findings have revealed key genes regulation responsible for cold resistance, which help to understand the cold-resistant mechanisms and guide breeding in tea plants.

Critical factors determining fluoride concentration in tea leaves produced from Anhui province, China.

This study investigated the fluoride present in tea plants (Camellia sinensis (L.) O. Kuntze) and its relationship to soils, varieties, seasons and tea leaf maturity. The study also explored how different manufacturing processes affect the leaching of fluoride into tea beverages. The fluoride concentration in the tea leaves was significantly correlate to the concentration of water-soluble fluoride in the soil. Different tea varieties accumulated different levels of fluoride, with varieties, Anji baicha having the highest and Nongkang zao having the lowest fluoride concentration. In eight different varieties of tea plant harvested over three tea seasons, fluoride concentration were highest in the summer and lowest in the spring in china. The fluoride concentration in tea leaves was directly related to the maturity of the tea leaves at harvest. Importantly, the tea manufacturing process did not introduced fluoride contamination. The leaching of fluoride was 6.8% and 14.1% higher in black and white tea, respectively, than in fresh tea leaves. The manufacturing step most affecting the leaching of fluoride into tea beverage was withering used in white, black and oolong tea rather than rolling or fermentation. The exposure and associated health risks for fluoride concentration in infusions of 115 commercially available teas from Chinese tea markets was determined. The fluoride concentration ranged from 5.0 to 306.0mgkg(-1), with an average of 81.7mgkg(-1). The hazard quotient (HQ) of these teas indicated that there was no risk of fluorosis from drinking tea, based on statistical analysis by Monte Carlo simulation.

Genome-Wide Analysis of the Biosynthesis and Deactivation of Gibberellin-Dioxygenases Gene Family in Camellia sinensis (L.) O. Kuntze

Gibberellins (GAs), a class of diterpenoid phytohormones, play a key role in regulating diverse processes throughout the life cycle of plants. Bioactive GA levels are rapidly regulated by Gibberellin-dioxygenases (GAox), which are involved in the biosynthesis and deactivation of gibberellin. In this manuscript, a comprehensive genome-wide analysis was carried out to find all GAox in Camellia sinensis. For the first time in a tea plant, 14 CsGAox genes, containing two domains, DIOX_N (PF14226) and 2OG-FeII_Oxy, were identified (PF03171). These genes all belong to 2-oxoglutarate-dependent dioxygenases (2-ODD), including four CsGA20ox (EC: 1.14.11.12), three CsGA3ox (EC: 1.14.11.15), and seven CsGA2ox (EC: 1.14.11.13). According to the phylogenetic classification as in Arabidopsis, the CsGAox genes spanned five subgroups. Each CsGAox shows tissue-specific expression patterns, although these vary greatly. Some candidate genes, which may play an important role in response to external abiotic stresses, have been identified with regards to patterns, such as CsGA20ox2, CsGA3ox2, CsGA3ox3, CsGA2ox1, CsGA2ox2, and CsGA2ox4. The bioactive GA levels may be closely related to the GA20ox, GA3ox and GA2ox genes. In addition, the candidate genes could be used as marker genes for abiotic stress resistance breeding in tea plants.

Differential expression of the ascorbate oxidase multigene family of Camellia sinensis in response to stress

ABSTRACT Apoplastic ascorbate oxidase (AO) plays a major role in cell growth. Although AO genes have been studied in depth, some articles have mistakenly identified AO homologues as AO genes. Overall, the divergence between AO genes and AO homologues has not been explored. Meanwhile, there is little information concerning AO and the AO homologue with respect to Camellia sinensis. In the present study, one CsAO homologue and three CsAOs were confirmed by RT-PCR amplification, cloning and sequencing. Multicopper oxidase type 1 (PF00394), type 2 (PF07731) or type 3 (PF07732) domains and one transmembrane helix were the key domains for each member of the cupredoxin family. The CsAOs with their counterparts from seven dicotyledonous plants and three monocotyledonous plants were used to build phylogenetic tree and compare the deduced polypeptides. CsAO may be strongly expressed in the stretch expanded tissues, including bud and root. The abiotic stress-induced expression pattern of the CsAO homologue (CsAO2) is similar to those of CsAO1, CsAO3 and CsAO4. A new and very large group of AO homologues, which may function as AO genes, was present in both dicotyledonous and monocotyledonous plants. Our study may help in identifying stress-responsive AO genes of plants.