Non-Volatile Metabolic Profiling and Regulatory Network Analysis in Fresh Shoots of Tea Plant and Its Wild Relatives.

There are numerous non-volatile metabolites in the fresh shoots of tea plants. However, we know little about the complex relationship between the content of these metabolites and their gene expression levels. In investigating this, this study involved non-volatile metabolites from 68 accessions of tea plants that were detected and identified using untargeted metabolomics. The tea accessions were divided into three groups from the results of a principal component analysis based on the relative content of the metabolites. There were differences in variability between the primary and secondary metabolites. Furthermore, correlations among genes, gene metabolites, and metabolites were conducted based on Pearson's correlation coefficient (PCC) values. This study offered several significant insights into the co-current network of genes and metabolites in the global genetic background. Thus, the study is useful for providing insights into the regulatory relationship of the genetic basis for predominant metabolites in fresh tea shoots.

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.

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.

Identification and distribution of a single nucleotide polymorphism responsible for the catechin content in tea plants.

Catechins are the predominant products in tea plants and have essential functions for both plants and humans. Several genes encoding the enzymes regulating catechin biosynthesis have been identified, and the identification of single nucleotide polymorphisms (SNPs) resulting in nonsynonymous mutations within these genes can be used to establish a functional link to catechin content. Therefore, the transcriptomes of two parents and four filial offspring were sequenced using next-generation sequencing technology and aligned to the reference genome to enable SNP mining. Subsequently, 176 tea plant accessions were genotyped based on candidate SNPs using kompetitive allele-specific polymerase chain reaction (KASP). The catechin contents of these samples were characterized by high-performance liquid chromatography (HPLC), and analysis of variance (ANOVA) was subsequently performed to determine the relationship between genotypes and catechin content. As a result of these efforts, a SNP within the chalcone synthase (CHS) gene was shown to be functionally associated with catechin content. Furthermore, the geographical and interspecific distribution of this SNP was investigated. Collectively, these results will contribute to the early evaluation of tea plants and serve as a rapid tool for accelerating targeted efforts in tea breeding.

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.

Carthamin Yellow Protects the Heart Against Ischemia/Reperfusion Injury With Reduced Reactive Oxygen Species Release and Inflammatory Response.

Carthamin yellow (CY) is a flavonoid compound isolated from safflower, which is widely used clinically in China. It has various pharmacological effects including promoting blood circulation to remove blood stasis and alleviating pain. Ischemic heart disease is one of the main culprits of illness and death. Here, in this study, ex vivo and in vivo models were used to investigate whether CY reduces ischemia/reperfusion injury. In vitro experiments further verify and explain the potential mechanisms of CY cardioprotective function. Isolated hearts from male rats with or without CY pretreatment before ischemia which underwent 30-minute ischemia followed by 60-minute reperfusion showed that CY pretreatment significantly reduced the infarct size and lactate dehydrogenase release. The in vivo experiments also indicated CY preadministration (i.v.) reduced infarct size and improved the heart function, which was impaired by myocardial ischemia/reperfusion injury. The in vitro model on myocardial cell also showed that CY reduced ischemia/reperfusion injury by reducing the lactate dehydrogenase and reactive oxygen species (ROS) releasing. Eliminating ROS with N-acetylcysteine or preinject CY into rat jugular vein reduces the expression of IL-6, TNF-a, and, especially, IL-1b in an in vivo I/R model. Also, CY pretreatment strongly reduces ischemia/reperfusion-induced NLRP3 up-expression and caspase-1 activation. Our results indicated CY reduced ischemia-reperfusion injury when administered before reperfusion. The reduction in injury is accompanied by a reduced ROS release and decreased inflammatory response.

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.

The complete chloroplast genome of an albino tea, Camellia sinensis cultivar 'Baiye 1'.

For obtaining the sequence and phylogenetic position of Camellia sinensis cultivar 'Baiye1', the complete chloroplast genome was determined. This chloroplast genome is 156,691 bp in length with overall GC content of 37.3%. It was comprised of a large single-copy (LSC) region of 86,585bp, a small single-copy (SSC) region of 18,276bp, and two inverted repeat (IR) regions of 26,083 bp. It contains 87 protein-coding, 8 rRNA, and 35 tRNA genes. Phylogenetic analysis showed 'Baiye1' and C. sinensis cv. 'Longjing43' were clustered into a group. These results may contribute to the further understanding of the albino phenotype and genetic evolution.

The complete chloroplast genome sequence of Camellia tachangensis F. C. Zhang (Theaceae).

To understand genetic background and phylogenetic position of Camellia tachangensis, we determined its complete chloroplast genome sequence which is 157,026 bp in length with overall GC content of 36.7%. It has four sub regions: a large single-copy (LSC) region (86,669 bp) and a small single-copy (SSC) region (18,253 bp) are separated by two inverted repeat (IR) regions (26,052 bp each). A total of 129 genes were annotated, containing 86 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Phylogenetic trees showed C. tachangensis clustered with Camellia gymnogyna and Camellia taliensis and separated from Camellia sinensis and its two varieties, Camellia sinensis var. assamica and Camellia sinensis var. pubilimba.

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.

Expression and Histopathological Significance of Disabled-2 in Aldosterone-Producing Adenoma.

The current pathological diagnosis of aldosterone-producing adenoma (APA) is challenging because no histological markers of aldosterone production are available in routine practice. A previous study demonstrated that Disabled-2 (DAB2) is a specific marker of the zona glomerulosa (ZG) in rodents. The aim of the present study was to investigate the significance of immunohistochemical staining to detect DAB2 in the adrenal tissue of patients with APA. We investigated the expression of DAB2 in 36 adrenal glands with APA, 23 adrenal glands with cortisol-producing adenoma (CPA), and 33 adrenal glands with non-functioning adenoma (NFA). Immunohistochemical staining was performed using anti-DAB2 antibodies on paraffin-embedded sections. We analysed the expression of DAB2 semi-quantitatively by scoring staining intensity, and assessed the correlation of this information with the clinical findings. DAB2 mRNA expression in adenoma tissues was evaluated by RT-PCR. DAB2 was highly expressed in the ZG in normal human adrenal glands. DAB2 expression was heterogeneous in APA, with spotted, strong staining noted in most samples (25 of 36 APA). CPA and NFA also exhibited extensive low or moderate DAB2 expression. DAB2 mRNA was significantly increased and positively correlated with CYP11B2 in APA (p<0.05). In APA, the DAB2 score adjusted for tumour volume was positively correlated with plasma aldosterone (p<0.05). Patients with low or moderate DAB2 staining more frequently exhibited high blood pressure and were diagnosed at a younger age compared with patients with high DAB2 staining. The present study clearly demonstrates that DAB2 is a specific marker of the ZG in normal human adrenal glands but that DAB2 immunostaining is not sufficiently powerful for histopathological diagnosis of APA. DAB2 might be involved in excessive aldosterone biosynthesis and correlate with specific clinical characteristics of APA patients.

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.