Comprehensive analysis and expression profiles of the AP2/ERF gene family during spring bud break in tea plant (Camellia sinensis).

BACKGROUND: AP2/ERF transcription factors (AP2/ERFs) are important regulators of plant physiological and biochemical metabolism. Evidence suggests that AP2/ERFs may be involved in the regulation of bud break in woody perennials. Green tea is economically vital in China, and its production value is significantly affected by the time of spring bud break of tea plant. However, the relationship between AP2/ERFs in tea plant and spring bud break remains largely unknown. RESULTS: A total of 178 AP2/ERF genes (CsAP2/ERFs) were identified in the genome of tea plant. Based on the phylogenetic analysis, these genes could be classified into five subfamilies. The analysis of gene duplication events demonstrated that whole genome duplication (WGD) or segmental duplication was the primary way of CsAP2/ERFs amplification. According to the result of the Ka/Ks value calculation, purification selection dominated the evolution of CsAP2/ERFs. Furthermore, gene composition and structure analyses of CsAP2/ERFs indicated that different subfamilies contained a variety of gene structures and conserved motifs, potentially resulting in functional differences among five subfamilies. The promoters of CsAP2/ERFs also contained various signal-sensing elements, such as abscisic acid responsive elements, light responsive elements and low temperature responsive elements. The evidence presented here offers a theoretical foundation for the diverse functions of CsAP2/ERFs. Additionally, the expressions of CsAP2/ERFs during spring bud break of tea plant were analyzed by RNA-seq and grouped into clusters A-F according to their expression patterns. The gene expression changes in clusters A and B were more synchronized with the spring bud break of tea plant. Moreover, several potential correlation genes, such as D-type cyclin genes, were screened out through weighted correlation network analysis (WGCNA). Temperature and light treatment experiments individually identified nine candidate CsAP2/ERFs that may be related to the spring bud break of tea plant. CONCLUSIONS: This study provides new evidence for role of the CsAP2/ERFs in the spring bud break of tea plant, establishes a theoretical foundation for analyzing the molecular mechanism of the spring bud break of tea plant, and contributes to the improvement of tea cultivars.

Metabolite signatures of diverse Camellia sinensis tea populations.

The tea plant (Camellia sinensis) presents an excellent system to study evolution and diversification of the numerous classes, types and variable contents of specialized metabolites. Here, we investigate the relationship among C. sinensis phylogenetic groups and specialized metabolites using transcriptomic and metabolomic data on the fresh leaves collected from 136 representative tea accessions in China. We obtain 925,854 high-quality single-nucleotide polymorphisms (SNPs) enabling the refined grouping of the sampled tea accessions into five major clades. Untargeted metabolomic analyses detect 129 and 199 annotated metabolites that are differentially accumulated in different tea groups in positive and negative ionization modes, respectively. Each phylogenetic group contains signature metabolites. In particular, CSA tea accessions are featured with high accumulation of diverse classes of flavonoid compounds, such as flavanols, flavonol mono-/di-glycosides, proanthocyanidin dimers, and phenolic acids. Our results provide insights into the genetic and metabolite diversity and are useful for accelerated tea plant breeding.

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.

Novel insight into theacrine metabolism revealed by transcriptome analysis in bitter tea (Kucha, Camellia sinensis).

Kucha (Camellia sinensis) is a kind of unique wild tea resources in southwest China, containing sizeable amounts of theacrine (1,3,7,9-tetramethyluric acid) and having a special bitter taste both in fresh leaves and made tea. Theacrine has good healthy function locally. But the molecular mechanism of theacrine metabolism in Kucha was still unclear. In order to illuminate the biosynthesis and catabolism of theacrine in Kucha plants, three tea cultivars, C. sinensis 'Shangyou Zhongye' (SY) with low-theacrine, 'Niedu Kucha 2' (ND2) with middle-theacrine and, 'Niedu Kucha 3' (ND3) with high-theacrine, were used for our research. Purine alkaloid analysis and transcriptome of those samples were performed by High Performance Liquid Chromatography (HPLC) and RNA-Seq, respectively. The related gene expression levels of purine alkaloid were correlated with the content of purine alkaloid, and the results of quantitative real-time (qRT) PCR were also confirmed the reliability of transcriptome. Based on the data, we found that theacrine biosynthesis is a relatively complex process, N-methyltransferase (NMT) encoded by TEA024443 may catalyze the methylation at 9-N position in Kucha plant. Our finding will assist to reveal the molecular mechanism of theacrine biosynthesis, and be applied to selection and breeding of Kucha tea cultivars in the future.

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.

Metabolomics for a Millenniums-Old Crop: Tea Plant ( Camellia sinensis).

Tea cultivation and utilization dates back to antiquity. Today it is the most widely consumed beverage on earth due to its pleasant taste and several beneficial health properties attributed to specific metabolites. Metabolomics has a tremendous potential to correlate tea metabolites with taste and health properties in humans. Our review on the current application of metabolomics in the science of tea suggests that metabolomics is a promising frontier in the evaluation of tea quality, identification of functional genes responsible for key metabolites, investigation of their metabolic regulation, and pathway analysis in the tea plant. Furthermore, the challenges, possible solutions, and the prospects of metabolomics in tea science are reviewed.

Regulation of Growth and Flavonoid Formation of Tea Plants ( Camellia sinensis) by Blue and Green Light.

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.

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.

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.

Transcriptome and metabolome analysis reveal candidate genes and biochemicals involved in tea geometrid defense in Camellia sinensis.

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.

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.

Quantitative Succinyl-Proteome Profiling of Camellia sinensis cv. 'Anji Baicha' During Periodic Albinism.

Lysine succinylation is a novel dynamic and evolutionarily conserved post-translational modification (PTM) that regulates various biological processes. 'Anji Baicha' is an albino tea variety that exhibits temperature-based variability of leaf colour and amino acid concentrations. However, the mechanism underlying albinism in 'Anji Baicha' has not been investigated at the level of succinylation. Here, we identify 3530 lysine succinylation sites mapped to 2132 proteins in 'Anji Baicha', representing the first extensive data on the lysine succinylome in the tea plant. Eleven conserved succinylation motifs were enriched among the identified succinylated peptides. The protein-protein interaction maps were visualized using Cytoscape software. Comparison across three typical developmental stages of 'Anji Baicha' revealed that proteins exhibiting differential succinylation levels were primarily involved in photosynthesis, carbon fixation, biosynthesis of amino acids and porphyrin and chlorophyll metabolism, suggesting that these succinylated proteins are involved in 'Anji Baicha' leaf colour variability. These results not only deepen our understanding of the mechanism underlying 'Anji Baicha' albinism and the regulatory role of succinylation in the tea plant but also provide new insight into molecular breeding for leaf colour variety.

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.

Small RNA and degradome profiling reveals important roles for microRNAs and their targets in tea plant response to drought stress.

Tea (Camellia sinensis) is a popular beverage worldwide. Drought stress (DS) is a major constraint on the growth, yield and quality of tea plants. MicroRNAs (miRNAs) play important roles in plant responses to DS. We constructed eight small RNA libraries from the drought-tolerant 'Ningzhou 2' (NZ2) and drought-susceptible 'Zhuyeqi' (ZYQ) cultivars during four stages [control (CK), the fourth day of DS, the eighth day of DS and after recovery (RC)]. A total of 268 conserved and 62 novel miRNAs were identified using small RNA sequencing. In total, 139 (52.9%) and 96 (36.0%) conserved miRNAs were differentially expressed during the four stages (P ≤ 0.05) in NZ2 and ZYQ, respectively. A total of 814 predicted target genes were identified as differentially regulated by 199 miRNAs through degradome sequencing. Among them, 201 and 218 genes were specific to the NZ2 and ZYQ cultivars, respectively, and 395 were common to both cultivars. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed the biological roles of these targets and showed that some of the targets responded to DS in a stress- and cultivar-dependent manner. Correlated expression patterns between miRNA and their targets showed that specific miRNAs target the miRNA effector Argonaute 1 (AGO1), drought signaling-related receptors and enzymes, transcription factors, and other structural and functional proteins. The predicted regulatory networks provide insights into a potential miRNA-mediated regulatory mechanism. These results will contribute to the breeding of drought-tolerant tea plants and to elucidating miRNA regulation in response to drought.

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.

Natural allelic variations of TCS1 play a crucial role in caffeine biosynthesis of tea plant and its related species.

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.

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.

Differential Metabolic Profiles during the Albescent Stages of 'Anji Baicha' (Camellia sinensis).

'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.

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq).

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.