Tea plant genomics: achievements, challenges and perspectives

Simultaneous estimation of total phenolic and alkaloid contents in the tea samples by utilizing the catechin and caffeine oxidation signals through the square-wave voltammetry technique

This work outlines the simultaneous estimation of the total phenolic and alkaloid contents in the tea samples by using catechin (C) and caffeine (CAF) oxidation signals at a non-modified boron-doped diamond (BDD) electrode. Two irreversible oxidation peaks, about + 1.03 (for C) and + 1.45 V (for CAF) vs Ag/AgCl in acetate buffer solution at pH 4.7, were seen in the cyclic voltammetric profile of the binary mixtures of C and CAF. In optimal conditions and utilizing the square-wave mode, the BDD electrode allows for simultaneous quantification of C and CAF within the concentration ranges of 5.0-100.0 µg mL-1 (1.72 × 10-5 - 3.45 × 10-3 mol/L) and 1.0-50.0 µg mL-1 (5.15 × 10-6 - 2.57 × 10-4 mol/L) respectively. The corresponding detection limits are 1.22 µg mL-1 (4.21 × 10-6 mol/L) for C and 0.11 µg mL-1 (5.66 × 10-7 mol/L) for CAF. Other phenolic compounds (like tannic acid, gallic acid, epicatechin, and epigallocatechin gallate) and other alkaloids (theophylline and theobromine) present in tea samples were examined for selectivity assessment. Ultimately, the applicability of the proposed approach was demonstrated by estimating the total phenolic and alkaloid contents in the black and green tea samples, expressed as C and CAF equivalents. The results obtained were contrasted against those acquired using UV-Vis spectrometry.

Physiological and Transcriptome Analyses Reveal the Protective Effect of Exogenous Trehalose in Response to Heat Stress in Tea Plant (Camellia sinensis)

It is well known that application of exogenous trehalose can enhance the heat resistance of plants. To investigate the underlying molecular mechanisms by which exogenous trehalose induces heat resistance in C. sinensis, a combination of physiological and transcriptome analyses was conducted. The findings revealed a significant increase in the activity of superoxide dismutase (SOD) and peroxidase (POD) upon treatment with 5.0 mM trehalose at different time points. Moreover, the contents of proline (PRO), endogenous trehalose, and soluble sugar exhibited a significant increase, while malondialdehyde (MDA) content decreased following treatment with 5.0 mM trehalose under 24 h high-temperature stress (38 °C/29 °C, 12 h/12 h). RNA-seq analysis demonstrated that the differentially expressed genes (DEGs) were significantly enriched in the MAPK pathway, plant hormone signal transduction, phenylpropanoid biosynthesis, flavone and flavonol biosynthesis, flavonoid biosynthesis, and the galactose metabolism pathway. The capability to scavenge free radicals was enhanced, and the expression of a heat shock factor gene (HSFB2B) and two heat shock protein genes (HSP18.1 and HSP26.5) were upregulated in the tea plant. Consequently, it was concluded that exogenous trehalose contributes to alleviating heat stress in C. sinensis. Furthermore, it regulates the expression of genes involved in diverse pathways crucial for C. sinensis under heat-stress conditions. These findings provide novel insights into the molecular mechanisms underlying the alleviation of heat stress in C. sinensis with trehalose.

A method for mining condition-specific co-expressed genes in Camellia sinensis based on k-means clustering

BACKGROUND

As one of the world's most important beverage crops, tea plants (Camellia sinensis) are renowned for their unique flavors and numerous beneficial secondary metabolites, attracting researchers to investigate the formation of tea quality. With the increasing availability of transcriptome data on tea plants in public databases, conducting large-scale co-expression analyses has become feasible to meet the demand for functional characterization of tea plant genes. However, as the multidimensional noise increases, larger-scale co-expression analyses are not always effective. Analyzing a subset of samples generated by effectively downsampling and reorganizing the global sample set often leads to more accurate results in co-expression analysis. Meanwhile, global-based co-expression analyses are more likely to overlook condition-specific gene interactions, which may be more important and worthy of exploration and research.

RESULTS

Here, we employed the k-means clustering method to organize and classify the global samples of tea plants, resulting in clustered samples. Metadata annotations were then performed on these clustered samples to determine the "conditions" represented by each cluster. Subsequently, we conducted gene co-expression network analysis (WGCNA) separately on the global samples and the clustered samples, resulting in global modules and cluster-specific modules. Comparative analyses of global modules and cluster-specific modules have demonstrated that cluster-specific modules exhibit higher accuracy in co-expression analysis. To measure the degree of condition specificity of genes within condition-specific clusters, we introduced the correlation difference value (CDV). By incorporating the CDV into co-expression analyses, we can assess the condition specificity of genes. This approach proved instrumental in identifying a series of high CDV transcription factor encoding genes upregulated during sustained cold treatment in Camellia sinensis leaves and buds, and pinpointing a pair of genes that participate in the antioxidant defense system of tea plants under sustained cold stress.

CONCLUSIONS

To summarize, downsampling and reorganizing the sample set improved the accuracy of co-expression analysis. Cluster-specific modules were more accurate in capturing condition-specific gene interactions. The introduction of CDV allowed for the assessment of condition specificity in gene co-expression analyses. Using this approach, we identified a series of high CDV transcription factor encoding genes related to sustained cold stress in Camellia sinensis. This study highlights the importance of considering condition specificity in co-expression analysis and provides insights into the regulation of the cold stress in Camellia sinensis.

Fingerprinting black tea: When spectroscopy meets machine learning a novel workflow for geographical origin identification

Food fraud, along with many challenges to the integrity and sustainability, threatens the prosperity of businesses and society as a whole. Tea is the second most commonly consumed non-alcoholic beverage globally. Challenges to tea authenticity require the development of highly efficient and rapid solutions to improve supply chain transparency. This study has produced an innovative workflow for black tea geographical indications (GI) discrimination based on non-targeted spectroscopic fingerprinting techniques. A total of 360 samples originating from nine GI regions worldwide were analysed by Fourier Transform Infrared (FTIR) and Near Infrared spectroscopy. Machine learning algorithms (k-nearest neighbours and support vector machine models) applied to the test data greatly improved the GI identification achieving 100% accuracy using FTIR. This workflow will provide a low-cost and user-friendly solution for on-site and real-time determination of black tea geographical origin along supply chains.

Characterization of the complete chloroplast genome sequence of Camellia tetracocca (Theaceae)

Camellia tetracocca H.T. Chang 1981 is an important wild relative of cultivated tea plants. Its leaves are widely used by local people to make tea, showing great economic and breeding values. We here report the complete chloroplast genome of C. tetracocca using Illumina sequencing technology. The complete chloroplast genome of C. tetracocca is 157,026 bp in length, and structurally contains a pair of inverted repeat regions (IRa and IRb, 26,052 bp) separated by a large single-copy region (LSC, 86,669 bp) and a small single-copy region (SSC, 18,253 bp). It is composed of 131 predicted genes, including 86 protein-coding genes, 37 transfer RNA genes, and eight ribosomal RNA genes. The overall GC content is 37.31%. Phylogenetic analysis among four Camellia species and 11 other close species reveals a close relationship between C. tetracocca and C. gymnogyna.

Agrobacterium Transformation of Tea Plants (Camellia sinensis (L.) KUNTZE): A Small Experiment with Great Prospects

Tea has historically been one of the most popular beverages, and it is currently an economically significant crop cultivated in over 50 countries. The Northwestern Caucasus is one of the northernmost regions for industrial tea cultivation worldwide. The domestication of the tea plant in this region took approximately 150 years, during which plantations spreading from the Ozurgeti region in northern Georgia to the southern city of Maykop in Russia. Consequently, tea plantations in the Northern Caucasus can serve as a source of unique genotypes with exceptional cold tolerance. Tea plants are known to be recalcitrant to Agrobacterium-mediated transfection. Research into optimal transfection and regeneration methodologies, as well as the identification of tea varieties with enhanced transformation efficiency, is an advanced strategy for improving tea plant culture. The aim of this study was to search for the optimal Agrobacterium tumefaciens-mediated transfection protocol for the Kolkhida tea variety. As a result of optimizing the transfection medium with potassium phosphate buffer at the stages of pre-inoculation, inoculation and co-cultivation, the restoration of normal morphology and improvement in the attachment of Agrobacterium cells to the surface of tea explants were observed by scanning electron microscopy. And an effective method of high-efficiency Agrobacteria tumefaciens-mediated transfection of the best local tea cultivar, Kolkhida, was demonstrated for the first time.

Multi-Omics Research Accelerates the Clarification of the Formation Mechanism and the Influence of Leaf Color Variation in Tea (Camellia sinensis) Plants

Tea is a popular beverage with characteristic functional and flavor qualities, known to be rich in bioactive metabolites such as tea polyphenols and theanine. Recently, tea varieties with variations in leaf color have been widely used in agriculture production due to their potential advantages in terms of tea quality. Numerous studies have used genome, transcriptome, metabolome, proteome, and lipidome methods to uncover the causes of leaf color variations and investigate their impacts on the accumulation of crucial bioactive metabolites in tea plants. Through a comprehensive review of various omics investigations, we note that decreased expression levels of critical genes in the biosynthesis of chlorophyll and carotenoids, activated chlorophyll degradation, and an impaired photosynthetic chain function are related to the chlorina phenotype in tea plants. For purple-leaf tea, increased expression levels of late biosynthetic genes in the flavonoid synthesis pathway and anthocyanin transport genes are the major and common causes of purple coloration. We have also summarized the influence of leaf color variation on amino acid, polyphenol, and lipid contents and put forward possible causes of these metabolic changes. Finally, this review further proposes the research demands in this field in the future.

Complete chloroplast genomes of 13 species of sect. Tuberculata Chang (Camellia L.): genomic features, comparative analysis, and phylogenetic relationships

BACKGROUND

Sect. Tuberculata belongs to Camellia, and its members are characterized by a wrinkled pericarp and united filaments. All the plants in this group, which are endemic to China, are highly valuable for exploring the evolution of Camellia and have great potential for use as an oil source. However, due to the complex and diverse phenotypes of these species and the difficulty of investigating them in the field, their complex evolutionary history and interspecific definitions have remained largely unelucidated.

RESULTS

Therefore, we newly sequenced and annotated 12 chloroplast (cp) genomes and retrieved the published cp genome of Camellia anlungensis Chang in sect. Tuberculata. In this study, comparative analysis of the cp genomes of the thirteen sect. Tuberculata species revealed a typical quadripartite structure characterized by a total sequence length ranging from 156,587 bp to 157,068 bp. The cp.genome arrangement is highly conserved and moderately differentiated. A total of 130 to 136 genes specific to the three types were identified by annotation, including protein-coding genes (coding sequences (CDSs)) (87-91), tRNA genes (35-37), and rRNA genes (8). The total observed frequency ranged from 23,045 (C. lipingensis) to 26,557 (C. anlungensis). IR region boundaries were analyzed to show that the ycf1 gene of C. anlungensis is located in the IRb region, while the remaining species are present only in the IRa region. Sequence variation in the SSC region is greater than that in the IR region, and most protein-coding genes have high codon preferences. Comparative analyses revealed six hotspot regions (tRNA-Thr(GGT)-psbD, psbE-petL, ycf15-tRNA-Leu(CAA), ndhF-rpl32, ndhD, and trnL(CAA)-ycf15) in the cp genomes that could serve as potential molecular markers. In addition, the results of phylogenetic tree construction based on the cp genomes showed that the thirteen sect. Tuberculata species formed a monophyletic group and were divided into two evolutionarily independent clades, confirming the independence of the section.

CONCLUSIONS

In summary, we obtained the cp genomes of thirteen sect. Tuberculata plants and performed the first comparative analysis of this group. These results will help us better characterize the plants in this section, deepen our understanding of their genetic characteristics and phylogenetic relationships, and lay the theoretical foundation for their accurate classification, elucidation of their evolutionary changes, and rational development and utilization of this section in the future.

Comparative Phylogenetic Analysis of Ancient Korean Tea "Hadong Cheon-Nyeon Cha (Camellia sinensis var. sinensis)" Using Complete Chloroplast Genome Sequences

Wild teas are valuable genetic resources for studying evolution and breeding. Here, we report the complete chloroplast genome of the ancient Korean tea 'Hadong Cheon-nyeon Cha' (C. sinensis var. sinensis), which is known as the oldest tea tree in Korea. This study determined seven Camellia sinensis var. sinenesis, including Hadong Cheon-nyeon Cha (HCNC) chloroplast genome sequences, using Illumina sequencing technology via de novo assembly. The chloroplast genome sizes ranged from 157,019 to 157,114 bp and were organized into quadripartite regions with the typical chloroplast genomes. Further, differences in SNPs and InDels were detected across the seven chloroplast genomes through variance analysis. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These genomic resources provide evolutionary insight into Korean tea plant cultivars and lay the foundation for a better understanding of the ancient Korean tea plant HCNC.

TPIA2: an updated tea plant information archive for Camellia genomics

The genus Camellia consists of about 200 species, which include many economically important species widely used for making tea, ornamental flowers and edible oil. Here, we present an updated tea plant information archive for Camellia genomics (TPIA2; http://tpia.teaplants.cn) by integrating more novel large-scale genomic, transcriptomic, metabolic and genetic variation datasets as well as a variety of useful tools. Specifically, TPIA2 hosts all currently available and well assembled 10 Camellia genomes and their comprehensive annotations from three major sections of Camellia. A collection of 15 million SNPs and 950 950 small indels from large-scale genome resequencing of 350 diverse tea accessions were newly incorporated, followed by the implementation of a novel 'Variation' module to facilitate data retrieval and analysis of the functionally annotated variome. Moreover, 116 Camellia transcriptomes were newly assembled and added, leading to a significant extension of expression profiles of Camellia genes to 13 developmental stages and eight abiotic/biotic treatments. An updated 'Expression' function has also been implemented to provide a comprehensive gene expression atlas for Camellia. Two novel analytic tools (e.g. Gene ID Convert and Population Genetic Analysis) were specifically designed to facilitate the data exchange and population genomics in Camellia. Collectively, TPIA2 provides diverse updated valuable genomic resources and powerful functions, and will continue to be an important gateway for functional genomics and population genetic studies in Camellia.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

Genome-Wide Association Studies of Biluochun Tea Plant Populations in Dongting Mountain and Comprehensive Identification of Candidate Genes Associated with Core Agronomic Traits by Four Analysis Models

The elite germplasm resources are key to the beautiful appearance and pleasant flavor of Biluochun tea. We collected and measured the agronomic traits of 95 tea plants to reveal the trait diversity and breeding value of Biluochun tea plant populations. The results revealed that the agronomic traits of Biluochun tea plant populations were diverse and had high breeding value. Additionally, we resequenced these tea plant populations to reveal genetic diversity, population structure, and selection pressure. The Biluochun tea plant populations contained two groups and were least affected by natural selection based on the results of population structure and selection pressure. More importantly, four non-synonymous single nucleotide polymorphisms (nsSNPs) and candidate genes associated with (-)-gallocatechin gallate (GCG), (-)-gallocatechin (GC), and caffeine (CAF) were detected using at least two GWAS models. The results will promote the development and application of molecular markers and the utilization of elite germplasm from Biluochun populations.

Complete Chloroplast Genomes of Saussurea katochaete, Saussurea superba, and Saussurea stella: Genome Structures and Comparative and Phylogenetic Analyses

Saussurea plants are widely distributed in Asia and Europe; however, their complex phylogenetic relationships have led to many difficulties in phylogenetic studies and interspecific identification. In this study, we assembled, annotated, and analyzed the chloroplast genomes of three Saussurea plants: Saussurea katochaete, Saussurea superba, and Saussurea stella. The results showed that the full-length sequences of the three Saussurea plants were 152,561 bp, 151,452 bp, and 152,293 bp, respectively, which represent the typical quadripartite structure, and the genomes were relatively conserved. The gene annotation results showed that the chloroplast genomes of S. katochaete, S. superba, and S. stella were annotated with 128, 124, and 127 unique genes, respectively, which included 83, 80, and 83 protein-coding genes (PCGs), respectively, 37, 36, and 36 tRNA genes, respectively, and 8 rRNA genes. Moreover, 46, 45, and 43 SSR loci, respectively, and nine highly variable regions (rpl32-trnL-UAG, rpl32, ndhF-rpl32, ycf1, trnC-GCA-petN, trnC-GCA, rpcL, psbE-petL, and rpl16-trnG-UUG) were identified and could be used as potential molecular markers for population identification and phylogenetic study of Saussurea plants. Phylogenetic analyses strongly support the sisterhood of S. katochaete with S. superba and S. stella, and are all clustered with S. depsagensis, S. inversa, S. medusa, and S. gossipihora, of which S. gossipiphora is most closely related. Additionally, the phylogenetic results indicate a high frequency of differentiation among different species of Saussurea plants, and many different species or genera are morphologically very different from each other, which may be related to certain genetic material in the chloroplasts. This study provides an important reference for the identification of Saussurea plants and studies their evolution and phylogenetics.

An Integrated Analysis of microRNAs and the Transcriptome Reveals the Molecular Mechanisms Underlying the Regulation of Leaf Development in Xinyang Maojian Green Tea (Camellia sinensis)

Xinyang Maojian (XYMJ) tea is one of the world's most popular green teas; the development of new sprouts directly affects the yield and quality of tea products, especially for XYMJ, which has hairy tips. Here, we used transcriptome and small RNA sequencing to identify mRNAs and miRNAs, respectively, involved in regulating leaf development in different plant tissues (bud, leaf, and stem). We identified a total of 381 conserved miRNAs. Given that no genomic data for XYMJ green tea are available, we compared the sequencing data for XYMJ green tea with genomic data from a closely related species (Tieguanyin) and the Camellia sinensis var. sinensis database; we identified a total of 506 and 485 novel miRNAs, respectively. We also identified 11 sequence-identical novel miRNAs in the tissues of XYMJ tea plants. Correlation analyses revealed 97 miRNA-mRNA pairs involved in leaf growth and development; the csn-miR319-2/csnTCP2 and miR159-csnMYB modules were found to be involved in leaf development in XYMJ green tea. Quantitative real-time PCR was used to validate the expression levels of the miRNAs and mRNAs. The miRNAs and target genes identified in this study might shed new light on the molecular mechanisms underlying the regulation of leaf development in tea plants.

A geraniol synthase regulates plant defense via alternative splicing in tea plants

Geraniol is an important contributor to the pleasant floral scent of tea products and one of the most abundant aroma compounds in tea plants; however, its biosynthesis and physiological function in response to stress in tea plants remain unclear. The proteins encoded by the full-length terpene synthase (CsTPS1) and its alternative splicing isoform (CsTPS1-AS) could catalyze the formation of geraniol when GPP was used as a substrate in vitro, whereas the expression of CsTPS1-AS was only significantly induced by Colletotrichum gloeosporioides and Neopestalotiopsis sp. infection. Silencing of CsTPS1 and CsTPS1-AS resulted in a significant decrease of geraniol content in tea plants. The geraniol content and disease resistance of tea plants were compared when CsTPS1 and CsTPS1-AS were silenced. Down-regulation of the expression of CsTPS1-AS reduced the accumulation of geraniol, and the silenced tea plants exhibited greater susceptibility to pathogen infection than control plants. However, there was no significant difference observed in the geraniol content and pathogen resistance between CsTPS1-silenced plants and control plants in the tea plants infected with two pathogens. Further analysis showed that silencing of CsTPS1-AS led to a decrease in the expression of the defense-related genes PR1 and PR2 and SA pathway-related genes in tea plants, which increased the susceptibility of tea plants to pathogens infections. Both in vitro and in vivo results indicated that CsTPS1 is involved in the regulation of geraniol formation and plant defense via alternative splicing in tea plants. The results of this study provide new insights into geraniol biosynthesis and highlight the role of monoterpene synthases in modulating plant disease resistance via alternative splicing.

Analysis of the effect of different withering methods on tea quality based on transcriptomics and metabolomics

Withering is very important to the quality of Wuyi rock tea. In this study, transcriptomics and metabolomics were used to analyze the effects of different withering methods on tea quality formation. The results showed that sunlight withering (SW) was most beneficial in increasing the gene expression of ubiquinone and other terpenoid-quinone biosynthesis (ko00130), pyruvate metabolism (ko00620), starch and sucrose metabolism (ko00500), and tryptophan metabolism (ko00380) pathways, and increasing the content of nucleotides and derivatives, terpenoids, organic acids and lipids, thus enhancing the mellowness, fresh and brisk taste and aroma of tea. Withering trough withering (WW) was most beneficial in increasing the gene expression of glutathione metabolism (ko00480), phenylpropanoid biosynthesis (ko00940) pathways, increasing the content of phenolic acids and flavonoids, thus enhancing tea bitterness. A comprehensive evaluation of the metabolite content and taste characteristics of tea leaves showed SW to be the best quality and charcoal fire withering (FW) to be the worst quality. This study provided an important basis for guiding the processing of Wuyi rock tea with different flavors.

Polyphenols in Plants: Structure, Biosynthesis, Abiotic Stress Regulation, and Practical Applications (Review)

Phenolic compounds or polyphenols are among the most common compounds of secondary metabolism in plants. Their biosynthesis is characteristic of all plant cells and is carried out with the participation of the shikimate and acetate-malonate pathways. In this case, polyphenols of various structures are formed, such as phenylpropanoids, flavonoids, and various oligomeric and polymeric compounds of phenolic nature. Their number already exceeds 10,000. The diversity of phenolics affects their biological activity and functional role. Most of their representatives are characterized by interaction with reactive oxygen species, which manifests itself not only in plants but also in the human body, where they enter through food chains. Having a high biological activity, phenolic compounds are successfully used as medicines and nutritional supplements for the health of the population. The accumulation and biosynthesis of polyphenols in plants depend on many factors, including physiological-biochemical, molecular-genetic, and environmental factors. In the review, we present the latest literature data on the structure of various classes of phenolic compounds, their antioxidant activity, and their biosynthesis, including their molecular genetic aspects (genes and transfactors). Since plants grow with significant environmental changes on the planet, their response to the action of abiotic factors (light, UV radiation, temperature, and heavy metals) at the level of accumulation and composition of these secondary metabolites, as well as their metabolic regulation, is considered. Information is given about plant polyphenols as important and necessary components of functional nutrition and pharmaceutically valuable substances for the health of the population. Proposals on promising areas of research and development in the field of plant polyphenols are presented.

A review on matcha: Chemical composition, health benefits, with insights on its quality control by applying chemometrics and multi-omics

This review discussed the origin, manufacturing process, chemical composition, factors affecting quality and health benefits of matcha (Camellia sinensis), and the application of chemometrics and multi-omics in the science of matcha. The discussion primarily distinguishes between matcha and regular green tea with processing and compositional factors, and demonstrates beneficial health effects of consuming matcha. Preferred Reporting Items for Systematic Reviews and Meta-Analyses was adopted to search for relevant information in this review. Boolean operators were incorporated to explore related sources in various databases. Notably, climate, cultivar, maturity of tea leaves, grinding process and brewing temperature impact on the overall quality of matcha. Besides, sufficient shading prior to harvesting significantly increases the contents of theanine and chlorophyll in the tea leaves. Furthermore, the ground whole tea leaf powder delivers matcha with the greatest benefits to the consumers. The health promoting benefits of matcha are mainly contributed by its micro-nutrients and the antioxidative phytochemicals, specifically epigallocatechin-gallate, theanine and caffeine. Collectively, the chemical composition of matcha affected its quality and health benefits significantly. To this end, more studies are required to elucidate the biological mechanisms of these compounds for human health. Chemometrics and multi-omics technologies are useful to fill up the research gaps identified in this review.

Genome-wide identification of the PYL gene family of tea plants (Camellia sinensis) revealed its expression profiles under different stress and tissues

BACKGROUND

PYL (Pyrabactin resistance 1-like) protein is a receptor of abscisic acid (ABA), which plays an important role in ABA signaling and influences plant growth and development and stress response. However, studies on PYL gene family in tea plants have not been reported.

RESULTS

In this study, we identified 20 PYL genes from the reference genome of tea plant ('Shuchazao'). Phylogeny analysis indicated that PYLs from tea and other plant species were clustered into seven groups. The promoter region of PYL genes contains a large number of cis-elements related to hormones and stresses. A large number of PYL genes responding to stress were found by analyzing the expression levels of abiotic stress and biotic stress transcriptome data. For example, CSS0047272.1 were up-regulated by drought stress, and CSS0027597.1 could respond to both anthracnose disease and geometrid feeding treatments. In addition, 10 PYL genes related to growth and development were verified by RT-qPCR and their tissue expression characteristics were revealed.

CONCLUSIONS

Our results provided a comprehensive characteristic of the PYL gene family in tea plants and provided an important clue for further exploring its functions in the growth and development, and resistance to stress of tea plants.

Field Control Effect and Initial Mechanism: A Study of Isobavachalcone against Blister Blight Disease

Blister blight (BB) disease is caused by the obligate biotrophic fungal pathogen Exobasidium vexans Massee and seriously affects the yield and quality of Camellia sinensis. The use of chemical pesticides on tea leaves substantially increases the toxic risks of tea consumption. Botanic fungicide isobavachalcone (IBC) has the potential to control fungal diseases on many crops but has not been used on tea plants. In this study, the field control effects of IBC were evaluated by comparison and in combination with natural elicitor chitosan oligosaccharides (COSs) and the chemical pesticide pyraclostrobin (Py), and the preliminary action mode of IBC was also investigated. The bioassay results for IBC or its combination with COSs showed a remarkable control effect against BB (61.72% and 70.46%). IBC, like COSs, could improve the disease resistance of tea plants by enhancing the activity of tea-plant-related defense enzymes, including polyphenol oxidase (PPO), catalase (CAT), phenylalanine aminolase (PAL), peroxidase (POD), superoxide dismutase (SOD), β-1,3-glucanase (Glu), and chitinase enzymes. The fungal community structure and diversity of the diseased tea leaves were examined using Illumina MiSeq sequencing of the internal transcribed spacer (ITS) region of the ribosomal rDNA genes. It was obvious that IBC could significantly alter the species' richness and the diversity of the fungal community in affected plant sites. This study broadens the application range of IBC and provides an important strategy for the control of BB disease.

Establishment of an efficient transformation system and its application in regulatory mechanism analysis of biological macromolecules in tea plants

Tea (Camellia sinensis), one of the most important beverage crops originated from China and is now cultivated worldwide, provides numerous secondary metabolites that account for its health benefits and rich flavor. However, the lack of an efficient and reliable genetic transformation system has seriously hindered the gene function investigation and precise breeding of C. sinensis. In this study, we established a highly efficient, labor-saving, and cost-effective Agrobacterium rhizogenes-mediated hairy roots genetic transformation system for C. sinensis, which can be used for gene overexpression and genome editing. The established transformation system was simple to operate, bypassing tissue culture and antibiotic screening, and only took two months to complete. We used this system to conduct function analysis of transcription factor CsMYB73 and found that CsMYB73 negatively regulates L-theanine synthesis in tea plant. Additionally, callus formation was successfully induced using transgenic roots, and the transgenic callus exhibited normal chlorophyll production, enabling the study of the corresponding biological functions. Furthermore, this genetic transformation system was effective for multiple C. sinensis varieties and other woody plant species. By overcoming technical obstacles such as low efficiency, long experimental periods, and high costs, this genetic transformation will be a valuable tool for routine gene investigation and precise breeding in tea plants.

Light Intensity Regulates Low-Temperature Adaptability of Tea Plant through ROS Stress and Developmental Programs

Low-temperature stress limits global tea planting areas and production efficiency. Light is another essential ecological factor that acts in conjunction with temperature in the plant life cycle. However, it is unclear whether the differential light environment affects the low temperature adaptability of tea plant (Camellia sect. Thea). In this study, tea plant materials in three groups of light intensity treatments showed differentiated characteristics for low-temperature adaptability. Strong light (ST, 240 μmol·m^-2·s^-1) caused the degradation of chlorophyll and a decrease in peroxidase (POD), superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), and polyphenol oxidase (PPO) activities, as well as an increase in soluble sugar, soluble protein, malondialdehyde (MDA), and relative conductivity in tea leaves. In contrast, antioxidant enzyme activities, chlorophyll content, and relative conductivity were highest in weak light (WT, 15 μmol·m^-2·s^-1). Damage was observed in both ST and WT materials relative to moderate light intensity (MT, 160 μmol·m^-2·s^-1) in a frost resistance test. Chlorophyll degradation in strong light was a behavior that prevented photodamage, and the maximum photosynthetic quantum yield of PS II (Fv/Fm) decreased with increasing light intensity. This suggests that the browning that occurs on the leaf surface of ST materials through frost may have been stressed by the previous increase in reactive oxygen species (ROS). Frost intolerance of WT materials is mainly related to delayed tissue development and tenderness holding. Interestingly, transcriptome sequencing revealed that stronger light favors starch biosynthesis, while cellulose biosynthesis is enhanced in weaker light. It showed that light intensity mediated the form of carbon fixation in tea plant, and this was associated with low-temperature adaptability.

Genome-wide identification of UDP-glycosyltransferases in the tea plant (Camellia sinensis) and their biochemical and physiological functions

Tea (Camellia sinensis) has been an immensely important commercially grown crop for decades. This is due to the presence of essential nutrients and plant secondary metabolites that exhibit beneficial health effects. UDP-glycosyltransferases (UGTs) play an important role in the diversity of such secondary metabolites by catalysing the transfer of an activated sugar donor to acceptor molecules, and thereby creating a huge variety of glycoconjugates. Only in recent years, thanks to the sequencing of the tea plant genome, have there been increased efforts to characterise the UGTs in C. sinensis to gain an understanding of their physiological role and biotechnological potential. Based on the conserved plant secondary product glycosyltransferase (PSPG) motif and the catalytically active histidine in the active site, UGTs of family 1 in C. sinensis are identified here, and shown to cluster into 21 groups in a phylogenetic tree. Building on this, our current understanding of recently characterised C. sinensis UGTs (CsUGTs) is highlighted and a discussion on future perspectives made.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of wild tea plant in Guizhou Plateau

BACKGROUND

Tea, the second largest consumer beverage in the world after water, is widely cultivated in tropical and subtropical areas. However, the effect of environmental factors on the distribution of wild tea plants is unclear.

RESULTS

A total of 159 wild tea plants were collected from different altitudes and geological types of the Guizhou Plateau. Using the genotyping-by-sequencing method, 98,241 high-quality single nucleotide polymorphisms were identified. Genetic diversity, population structure analysis, principal component analysis, phylogenetic analysis, and linkage disequilibrium were performed. The genetic diversity of the wild tea plant population from the Silicate Rock Classes of Camellia gymnogyna was higher than that from the Carbonate Rock Classes of Camellia tachangensis. In addition, the genetic diversity of wild tea plants from the second altitude gradient was significantly higher than that of wild tea plants from the third and first altitude gradients. Two inferred pure groups (GP01 and GP02) and one inferred admixture group (GP03) were identified by population structure analysis and were verified by principal component and phylogenetic analyses. The highest differentiation coefficients were determined for GP01 vs. GP02, while the lowest differentiation coefficients were determined for GP01 vs. GP03.

CONCLUSIONS

This study revealed the genetic diversity and geographical distribution characteristics of wild tea plants in the Guizhou Plateau. There are significant differences in genetic diversity and evolutionary direction between Camellia tachangensis with Carbonate Rock Classes at the first altitude gradient and Camellia gymnogyna with Silicate Rock Classes at the third altitude gradient. Geological environment, soil mineral element content, soil pH, and altitude markedly contributed to the genetic differentiation between Camellia tachangensis and Camellia gymnogyna.

Long-term cold, freezing and drought: overlapping and specific regulatory mechanisms and signal transduction in tea plant (Camellia sinensis (L.) Kuntze)

Introduction

Low temperatures and drought are two main environmental constraints reducing the yield and geographical distribution of horticultural crops worldwide. Understanding the genetic crosstalk between stress responses has potential importance for crop improvement.

Methods

In this study, Illumina RNA-seq and Pac-Bio genome resequencing were used to annotate genes and analyze transcriptome dynamics in tea plants under long-term cold, freezing, and drought.

Results

The highest number of differentially expressed genes (DEGs) was identified under long-term cold (7,896) and freezing (7,915), with 3,532 and 3,780 upregulated genes, respectively. The lowest number of DEGs was observed under 3-day drought (47) and 9-day drought (220), with five and 112 genes upregulated, respectively. The recovery after the cold had 6.5 times greater DEG numbers as compared to the drought recovery. Only 17.9% of cold-induced genes were upregulated by drought. In total, 1,492 transcription factor genes related to 57 families were identified. However, only 20 transcription factor genes were commonly upregulated by cold, freezing, and drought. Among the 232 common upregulated DEGs, most were related to signal transduction, cell wall remodeling, and lipid metabolism. Co-expression analysis and network reconstruction showed 19 genes with the highest co-expression connectivity: seven genes are related to cell wall remodeling (GATL7, UXS4, PRP-F1, 4CL, UEL-1, UDP-Arap, and TBL32), four genes are related to calcium-signaling (PXL1, Strap, CRT, and CIPK6), three genes are related to photo-perception (GIL1, CHUP1, and DnaJ11), two genes are related to hormone signaling (TTL3 and GID1C-like), two genes are involved in ROS signaling (ERO1 and CXE11), and one gene is related to the phenylpropanoid pathway (GALT6).

Discussion

Based on our results, several important overlapping mechanisms of long-term stress responses include cell wall remodeling through lignin biosynthesis, o-acetylation of polysaccharides, pectin biosynthesis and branching, and xyloglucan and arabinogalactan biosynthesis. This study provides new insight into long-term stress responses in woody crops, and a set of new target candidate genes were identified for molecular breeding aimed at tolerance to abiotic stresses.

Understanding the Origin and Evolution of Tea (Camellia sinensis [L.]): Genomic Advances in Tea

Tea, which is processed by the tender shoots or leaves of tea plant (Camellia sinensis), is one of the most popular nonalcoholic beverages in the world and has numerous health benefits for humans. Along with new progress in biotechnologies, the refined chromosome-scale reference tea genomes have been achieved, which facilitates great promise for the understanding of fundamental genomic architecture and evolution of the tea plants. Here, we summarize recent achievements in genome sequencing in tea plants and review the new progress in origin and evolution of tea plants by population sequencing analysis. Understanding the genomic characterization of tea plants is import to improve tea quality and accelerate breeding in tea plants.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94 C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Topics and trends in fresh tea (Camellia sinensis) leaf research: A comprehensive bibliometric study

Tea plant (Camellia sinensis) is a widely cultivated cash crop and tea is a favorite functional food in the world. Fresh tea leaves (FTLs) play a critical role in bridging the two fields closely related to tea cultivation and tea processing, those are, tea plant biology and tea biochemistry. To provide a comprehensive overview of the development stages, authorship collaboration, research topics, and hotspots and their temporal evolution trends in the field of FTLs research, we conducted a bibliometric analysis, based on 971 publications on FTLs-related research published during 2001-2021 from Web of Science Core Collection. CiteSpace, R package Bibliometrix, and VOSviewer were employed in this research. The results revealed that the development history can be roughly divided into three stages, namely initial stage, slow development stage and rapid development stage. Journal of Agricultural & Food Chemistry published most articles in this field, while Frontiers in Plant Science held the highest total citations and h-index. The most influential country, institution, and author in this field was identified as China, the Chinese Academy of Agricultural Sciences, and Xiaochun Wan, respectively. FTLs-related research can be categorized into three main topics: the regulation mechanism of key genes, the metabolism and features of essential compounds, and tea plants' growth and stress responses. The most concerning hotspots are the application of advanced technologies, essential metabolites, leaf color variants, and effective cultivation treatments. There has been a shift from basic biochemical and enzymatic studies to studies of molecular mechanisms that depend on multi-omics technologies. We also discussed the future development in this field. This study provides a comprehensive summary of the research field, making it easier for researchers to be informed about its development history, status, and trends.

Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis

Tea plants are widely grown all over the world because they are an important economic crop. The purity and authenticity of tea varieties are frequent problems in the conservation and promotion of germplasm resources in recent years, which has brought considerable inconvenience and uncertainty to the selection of parental lines for breeding and the research and cultivation of superior varieties. However, the development of core SNP markers can quickly and accurately identify the germplasm, which plays an important role in germplasm identification and the genetic relationship analysis of tea plants. In this study, based on 179,970 SNP loci from the whole genome of the tea plant, all of 142 cultivars were clearly divided into three groups: Assam type (CSA), Chinese type (CSS), and transitional type. Most CSA cultivars are from Yunnan Province, which confirms that Yunnan Province is the primary center of CSA origin and domestication. Most CSS cultivars are distributed in east China; therefore, we deduced that east China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. Moreover, 45 core markers were screened using strict criteria to 179,970 SNP loci, and we analyzed 117 well-Known tea cultivars in China with 45 core SNP markers. The results were as follows: (1) In total, 117 tea cultivars were distinguished by eight markers, which were selected to construct the DNA fingerprint, and the remaining markers were used as standby markers for germplasm identification. (2) Ten pairs of parent and offspring relationships were confirmed or identified, and among them, seven pairs were well-established pedigree relationships; the other three pairs were newly identified. In this study, the east of China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. The 45 core SNP markers were developed, which provide a scientific basis at the molecular level to identify the superior tea germplasm, undertake genetic relationship analysis, and benefit subsequent breeding work.

Host genetic determinants drive compartment-specific assembly of tea plant microbiomes

Diverse host factors drive microbial variation in plant-associated environments, whereas their genetic mechanisms remain largely unexplored. To address this, we coupled the analyses of plant genetics and microbiomes in this study. Using 100 tea plant (Camellia sinensis) cultivars, the microbiomes of rhizosphere, root endosphere and phyllosphere showed clear compartment-specific assembly, whereas the subpopulation differentiation of tea cultivars exhibited small effects on microbial variation in each compartment. Through microbiome genome-wide association studies, we examined the interactions between tea genetic loci and microbial variation. Notably, genes related to the cell wall and carbon catabolism were heavily linked to root endosphere microbial composition, whereas genes related to the metabolism of metal ions and small organic molecules were overrepresented in association with rhizosphere microbial composition. Moreover, a set of tea genetic variants, including the cytoskeleton-related formin homology interacting protein 1 gene, were strongly associated with the β-diversity of phyllosphere microbiomes, implying their interactions with the overall structure of microbial communities. Our results create a catalogue of tea genetic determinants interacting with microbiomes and reveal the compartment-specific microbiome assembly driven by host genetics.

Chromosome-scale genome assembly of Camellia sinensis combined with multi-omics provides insights into its responses to infestation with green leafhoppers

The tea plant (Camellia sinensis) is an important economic crop, which is becoming increasingly popular worldwide, and is now planted in more than 50 countries. Tea green leafhopper is one of the major pests in tea plantations, which can significantly reduce the yield and quality of tea during the growth of plant. In this study, we report a genome assembly for DuyunMaojian tea plants using a combination of Oxford Nanopore Technology PromethION^™ with high-throughput chromosome conformation capture technology and used multi-omics to study how the tea plant responds to infestation with tea green leafhoppers. The final genome was 3.08 Gb. A total of 2.97 Gb of the genome was mapped to 15 pseudo-chromosomes, and 2.79 Gb of them could confirm the order and direction. The contig N50, scaffold N50 and GC content were 723.7 kb, 207.72 Mb and 38.54%, respectively. There were 2.67 Gb (86.77%) repetitive sequences, 34,896 protein-coding genes, 104 miRNAs, 261 rRNA, 669 tRNA, and 6,502 pseudogenes. A comparative genomics analysis showed that DuyunMaojian was the most closely related to Shuchazao and Yunkang 10, followed by DASZ and tea-oil tree. The multi-omics results indicated that phenylpropanoid biosynthesis, α-linolenic acid metabolism, flavonoid biosynthesis and 50 differentially expressed genes, particularly peroxidase, played important roles in response to infestation with tea green leafhoppers (Empoasca vitis Göthe). This study on the tea tree is highly significant for its role in illustrating the evolution of its genome and discovering how the tea plant responds to infestation with tea green leafhoppers will contribute to a theoretical foundation to breed tea plants resistant to insects that will ultimately result in an increase in the yield and quality of tea.

A scheduling route planning algorithm based on the dynamic genetic algorithm with ant colony binary iterative optimization for unmanned aerial vehicle spraying in multiple tea fields

The complex environments and weak infrastructure constructions of hilly mountainous areas complicate the effective path planning for plant protection operations. Therefore, with the aim of improving the current status of complicated tea plant protections in hills and slopes, an unmanned aerial vehicle (UAV) multi-tea field plant protection route planning algorithm is developed in this paper and integrated with a full-coverage spraying route method for a single region. By optimizing the crossover and mutation operators of the genetic algorithm (GA), the crossover and mutation probabilities are automatically adjusted with the individual fitness and a dynamic genetic algorithm (DGA) is proposed. The iteration period and reinforcement concepts are then introduced in the pheromone update rule of the ant colony optimization (ACO) to improve the convergence accuracy and global optimization capability, and an ant colony binary iteration optimization (ACBIO) is proposed. Serial fusion is subsequently employed on the two algorithms to optimize the route planning for multi-regional operations. Simulation tests reveal that the dynamic genetic algorithm with ant colony binary iterative optimization (DGA-ACBIO) proposed in this study shortens the optimal flight range by 715.8 m, 428.3 m, 589 m, and 287.6 m compared to the dynamic genetic algorithm, ant colony binary iterative algorithm, artificial fish swarm algorithm (AFSA) and particle swarm optimization (PSO), respectively, for multiple tea field scheduling route planning. Moreover, the search time is reduced by more than half compared to other bionic algorithms. The proposed algorithm maintains advantages in performance and stability when solving standard traveling salesman problems with more complex objectives, as well as the planning accuracy and search speed. In this paper, the research on the planning algorithm of plant protection route for multi-tea field scheduling helps to shorten the inter-regional scheduling range and thus reduces the cost of plant protection.

MTAGCN: predicting miRNA-target associations in Camellia sinensis var. assamica through graph convolution neural network

Background

MircoRNAs (miRNAs) play a central role in diverse biological processes of Camellia sinensis var.assamica (CSA) through their associations with target mRNAs, including CSA growth, development and stress response. However, although the experiment methods of CSA miRNA-target identifications are costly and time-consuming, few computational methods have been developed to tackle the CSA miRNA-target association prediction problem.

Results

In this paper, we constructed a heterogeneous network for CSA miRNA and targets by integrating rich biological information, including a miRNA similarity network, a target similarity network, and a miRNA-target association network. We then proposed a deep learning framework of graph convolution networks with layer attention mechanism, named MTAGCN. In particular, MTAGCN uses the attention mechanism to combine embeddings of multiple graph convolution layers, employing the integrated embedding to score the unobserved CSA miRNA-target associations.

Discussion

Comprehensive experiment results on two tasks (balanced task and unbalanced task) demonstrated that our proposed model achieved better performance than the classic machine learning and existing graph convolution network-based methods. The analysis of these results could offer valuable information for understanding complex CSA miRNA-target association mechanisms and would make a contribution to precision plant breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04819-3.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Genome Assembly and Analysis of the Flavonoid and Phenylpropanoid Biosynthetic Pathways in Fingerroot Ginger (Boesenbergia rotunda)

Boesenbergia rotunda (Zingiberaceae), is a high-value culinary and ethno-medicinal plant of Southeast Asia. The rhizomes of this herb have a high flavanone and chalcone content. Here we report the genome analysis of B. rotunda together with a complete genome sequence as a hybrid assembly. B. rotunda has an estimated genome size of 2.4 Gb which is assembled as 27,491 contigs with an N50 size of 12.386 Mb. The highly heterozygous genome encodes 71,072 protein-coding genes and has a 72% repeat content, with class I TEs occupying ~67% of the assembled genome. Fluorescence in situ hybridization of the 18 chromosome pairs at the metaphase showed six sites of 45S rDNA and two sites of 5S rDNA. An SSR analysis identified 238,441 gSSRs and 4604 EST-SSRs with 49 SSR markers common among related species. Genome-wide methylation percentages ranged from 73% CpG, 36% CHG and 34% CHH in the leaf to 53% CpG, 18% CHG and 25% CHH in the embryogenic callus. Panduratin A biosynthetic unigenes were most highly expressed in the watery callus. B rotunda has a relatively large genome with a high heterozygosity and TE content. This assembly and data (PRJNA71294) comprise a source for further research on the functional genomics of B. rotunda, the evolution of the ginger plant family and the potential genetic selection or improvement of gingers.

Genome-Wide Identification and Characterization of GARP Transcription Factor Gene Family Members Reveal Their Diverse Functions in Tea Plant (Camellia sinensis)

Golden2, ARR-B, Psr1 (GARP) proteins are plant-specific transcription factors that play vital and diverse roles in plants. However, systematic research on the GARP gene family in plants, including tea plant (Camellia sinensis), is scarce. In this study, a total of 69 GARP genes were identified and characterized from the tea plant genome based on the B-motif sequence signature. The CsGARP genes were clustered into five subfamilies: PHR1/PHL1, KAN, NIGT1/HRS1/HHO, GLK and ARR-B subfamilies. The phylogenetic relationships, gene structures, chromosomal locations, conserved motifs and regulatory cis-acting elements of the CsGARP family members were comprehensively analyzed. The expansion of CsGARP genes occurred via whole-genome duplication/segmental duplication, proximal duplication, and dispersed duplication under purifying selective pressure. The expression patterns of the CsGARP genes were systematically explored from various perspectives: in different tissues during different seasons; in different leaf color stages of tea plant; under aluminum treatment and nitrogen treatment; and in response to abiotic stresses such as cold, drought and salt and to biotic stress caused by Acaphylla theae. The results demonstrate that CsGARP family genes are ubiquitously expressed and play crucial roles in the regulation of growth and development of tea plant and the responses to environmental stimuli. Collectively, these results not only provide valuable information for further functional investigations of CsGARPs in tea plant but also contribute to broadening our knowledge of the functional diversity of GARP family genes in plants.

A survey of transcriptome complexity using full-length isoform sequencing in the tea plant Camellia sinensis

Tea is one of the most popular beverages and its leaves are rich in catechins, contributing to the diverse flavor as well as beneficial for human health. However, the study of the post-transcriptional regulatory mechanism affecting the synthesis of catechins remains insufficient. Here, we sequenced the transcriptome using PacBio sequencing technology and obtained 63,111 full-length high-quality isoforms, including 1302 potential novel genes and 583 highly reliable fusion transcripts. We also identified 1204 lncRNAs with high quality, containing 188 known and 1016 novel lncRNAs. In addition, 311 mis-annotated genes were corrected based on the high-quality Isoseq reads. A large number of alternative splicing (AS) events (3784) and alternative polyadenylation (APA) genes (18,714) were analyzed, accounting for 8.84% and 43.7% of the total annotated genes, respectively. We also found that 2884 genes containing AS and APA features exhibited higher expression levels than other genes. These genes are mainly involved in amino acid biosynthesis, carbon fixation in photosynthetic organisms, phenylalanine, tyrosine, tryptophan biosynthesis, and pyruvate metabolism, suggesting that they play an essential role in the catechins content of tea polyphenols. Our results further improved the level of genome annotation and indicated that post-transcriptional regulation plays a crucial part in synthesizing catechins.

Effects of Different Shading Treatments on the Biomass and Transcriptome Profiles of Tea Leaves (Camellia sinensis L.) and the Regulatory Effect on Phytohormone Biosynthesis

Our previous study showed that colored net shading treatments had comparable effects on the reduction of bitter and astringent compounds such as flavonol glycosides in tea leaves, compared with black net shading treatment, whereas the effects on the biomass and phytohormones are still unclear. In this study, we investigated the phytohormone and transcriptome profiles of tea leaves under different shading treatments, using black, blue, and red nets with the same shade percentages. The bud density, fresh weight of 100 buds, and yield under blue net shading treatments were greatly elevated by 2.00-fold, 1.24-fold, and 2.48-fold, compared with black net shading treatment, while their effects on flavonoid composition were comparable with black net shading treatment. The transcriptome profiles of different shade net-treated samples were well resolved and discriminated from control. The KEGG result indicated that the pathways of phenylpropanoid biosynthesis, MAPK signaling pathways, and plant hormone signal transduction were differentially regulated by different shading treatments. The co-expression analysis showed that the contents of salicylic acid and melatonin were closely correlated with certain light signal perception and signaling genes (p < 0.05), and UVR8, PHYE, CRY1, PHYB, PHOT2, and HY5 had more close interactions with phytohormone biosynthetic genes (p < 0.05). Our results suggest that different shading treatments can mediate the growth of tea plants, which could be attributed to the regulatory effect on phytohormones levels, providing an instruction for the production of summer/autumn tea and matcha.

Regulation of biosynthesis of the main flavor-contributing metabolites in tea plant (Camellia sinensis): A review

In the process of adapting to the environment, tea plants (Camellia sinensis) endow tea with unique flavor and health functions, which should be attributed to secondary metabolites, including catechins, L-theanine, caffeine and terpene volatiles. Since the content of these flavor-contributing metabolites are mainly determined by the growth of tea plant, it is very important to understand their alteration and regulation mechanisms. In the present work, we first summarize the distribution, change characteristics of the main flavor-contributing metabolites in different cultivars, organs and under environmental stresses of tea plant. Subsequently, we discuss the regulating mechanisms involved in the biosynthesis of these metabolites based on the existing evidence. Finally, we propose the remarks and perspectives on the future study relating flavor-contributing metabolites. This review would contribute to the acceleration of research on the characteristic secondary metabolites and the breeding programs in tea plants.

Comparative Analysis of Purine Alkaloids and Main Quality Components of the Three Camellia Species in China

Tea (Camelliasinensis var. sinensis) is a widely consumed caffeine-containing beverage, however the Camellia genus also includes other species, which are consumed as tea in their local growing regions. Presently, HPLC analysis assessed 126 unique Camellia germplasms belonging to three Camellia species, C. sinensis var. pubilimba Chang (Csp), C. gymnogyna Chang (CgC) and C. crassicolumna Chang (CcC). Theobromine was the predominant purine alkaloid in all species, representing over 90% of purine alkaloids in Csp and CgC, and 50% in CcC. Significant variability existed in purine alkaloid patterns both between and within species, and some germplasms possessed highly unique alkaloid profiles. Sensory evaluation and quality composition analysis of green tea products produced from the three Camellia species suggested their unsuitability for use in tea production due to their unpalatable flavor. The results of this study revealed the differences in purine alkaloids and main quality components between Camellia species and tea, which contributed to understand why tea, rather than other Camellia species, has become a popular beverage in the world after long-term artificial selection. In addition, unique alkaloid profiles suggest usefulness of these germplasm resources in future breeding of decaffeinated tea plant varieties and alkaloid metabolism research.

Tea and tea drinking: China's outstanding contributions to the mankind

Background

Tea trees originated in southwest China 60 million or 70 million years ago. Written records show that Chinese ancestors had begun drinking tea over 3000 years ago. Nowadays, with the aging of populations worldwide and more people suffering from non-communicable diseases or poor health, tea beverages have become an inexpensive and fine complementary and alternative medicine (CAM) therapy. At present, there are 3 billion people who like to drink tea in the world, but few of them actually understand tea, especially on its development process and the spiritual and cultural connotations.

Methods

We searched PubMed, Google Scholar, Web of Science, CNKI, and other relevant platforms with the key word “tea”, and reviewed and analyzed tea-related literatures and pictures in the past 40 years about tea’s history, culture, customs, experimental studies, and markets.

Results

China is the hometown of tea, tea trees, tea drinking, and tea culture. China has the oldest wild and planted tea trees in the world, fossil of a tea leaf from 35,400,000 years ago, and abundant tea-related literatures and art works. Moreover, tea may be the first Chinese herbal medicine (CHM) used by Chinese people in ancient times. Tea drinking has many benefits to our physical health via its antioxidant, anti-inflammatory, immuno-regulatory, anticancer, cardiovascular-protective, anti-diabetic, and anti-obesity activities. At the moment, COVID-19 is wreaking havoc across the globe and causing severe damages to people’s health and lives. Tea has anti-COVID-19 functions via the enhancement of the innate immune response and inhibition of viral growth. Besides, drinking tea can allow people to acquire a peaceful, relaxed, refreshed and cheerful enjoyment, and even longevity. According to the meridian theory of traditional Chinese medicine, different kinds of tea can activate different meridian systems in the human body. At present, black tea (fermented tea) and green tea (non-fermented tea) are the most popular in the world. Black tea accounts for over 90% of all teas sold in western countries. The world’s top-grade black teas include Qi Men black in China, Darjeeling and Assam black tea in India, and Uva black tea in Sri Lanka. However, all top ten famous green teas in the world are produced in China, and Xi Hu Long Jing tea is the most famous among all green teas. More than 700 different kinds of components and 27 mineral elements can be found in tea. Tea polyphenols and theaflavin/thearubigins are considered to be the major bioactive components of black tea and green tea, respectively. Overly strong or overheated tea liquid should be avoided when drinking tea.

Conclusions

Today, CAM provides an array of treatment modalities for the health promotion in both developed and developing countries all over the world. Tea drinking, a simple herb-based CAM therapy, has become a popular man-made non-alcoholic beverage widely consumed worldwide, and it can improve the growth of economy as well. Tea can improve our physical and mental health and promote the harmonious development of society through its chemical and cultural elements.

Genetic Diversity and Genome Size Variability in the Russian Genebank Collection of Tea Plant [Camellia sinensis (L). O. Kuntze]

The tea collection of the FRC SSC RAS (Sochi, Maykop in Russia) represents one of the northernmost germplasm comprising a number of locally derived cultivars and ɣ-irradiation mutants. The latter are often characterized by larger genome size, which may lead to better adaptation to biotic and abiotic stress. Such genotypes may be a valuable genetic resource for better adaptability to extreme environmental conditions, which could enable tea cultivation outside global growing regions. Microsatellite markers are often the best choice for genetic diversity analysis in genebank collections. However, their use in polyploid species is questionable because simple sequence repeat (SSR) allele dosage cannot be readily determined. Therefore, the efficiency of SSR and start codon targeted (SCoT) markers was investigated using 43 selected cultivars from the Russian genebank collection derived from mutant breeding and clonal selection. Previously, the increase in genome size was confirmed in 18 mutants within this collection. Despite the presence of polyploid tea genotypes, our study revealed higher efficiency of SSR markers than SCoT markers. Subsequent SSR analysis of the 106 genotypes in the Russian genebank collection revealed three distinct genetic clusters after STRUCTURE analysis. Greater genetic variation was observed within genetic clusters than between clusters, indicating low genetic variation between collections. Nevertheless, the northernmost tea collection exhibited a greater genetic distance from the other two clusters than they did from each other. Close genetic relationships were found between many cultivars with particularly large leaves and mutant forms. Pearson's correlation analysis revealed a significant, moderate correlation between genome size and leaf area size. Our study shows that microsatellite fingerprinting is useful to estimate the genetic diversity and genetic background of tea germplasm in Russia despite polyploid tea accessions. Thus, the results of our study contribute to the development of future tea germplasm conservation strategies and modern tea breeding programs.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.

The reference genome of camellia chekiangoleosa provides insights into camellia evolution and tea oil biosynthesis

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.

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

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

A first-generation haplotype map (HapMap-1) of tea (Camellia sinensis L. O. Kuntz)

MOTIVATION

Tea is a cross-pollinated woody perennial plant, which is why, application of conventional breeding is limited for its genetic improvement. However, lack of the genome-wide high-density SNP markers and genome-wide haplotype information has greatly hampered the utilization of tea genetic resources toward fast-track tea breeding programs. To address this challenge, we have generated a first-generation haplotype map of tea (Tea HapMap-1). Out-crossing and highly heterozygous nature of tea plants, make them more complicated for DNA-level variant discovery.

RESULTS

In this study, whole genome re-sequencing data of 369 tea genotypes were used to generate 2,334,564 biallelic SNPs and 1,447,985 InDels. Around 2928.04 million paired-end reads were used with an average mapping depth of ∼0.31× per accession. Identified polymorphic sites in this study will be useful in mapping the genomic regions responsible for important traits of tea. These resources lay the foundation for future research to understand the genetic diversity within tea germplasm and utilize genes that determine tea quality. This will further facilitate the understanding of tea genome evolution and tea metabolite pathways thus, offers an effective germplasm utilization for breeding the tea varieties.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pan-transcriptome assembly combined with multiple association analysis provides new insights into the regulatory network of specialized metabolites in the tea plant Camellia sinensis

Specialized metabolites not only play important roles in biotic and abiotic stress adaptation of tea plants (Camellia sinensis (L.) O. Kuntze) but also contribute to the unique flavor of tea, the most important nonalcoholic beverage. However, the molecular networks and major genes that regulate specialized metabolites in tea plants are not well understood. Here, we constructed a population-level pan-transcriptome of the tea plant leaf using second-leaf transcriptome data from 134 accessions to investigate global expression differences in the population, expression presence or absence variations (ePAVs), and differentially expressed genes (DEGs) between pure Camellia sinensis var. assamica (CSA) and pure Camellia sinensis var. sinensis (CSS) accessions. Next, we used a genome-wide association study, a quantitative trait transcript study, and a transcriptome-wide association study to integrate genotypes, accumulation levels of specialized metabolites, and expression levels of pan-transcriptome genes to identify candidate regulatory genes for flavor-related metabolites and to construct a regulatory network for specialized metabolites in tea plants. The pan-transcriptome contains 30 482 expressed genes, 4940 and 5506 of which were newly annotated from a de novo transcriptome assembly without a reference and a genome reference-based assembly, respectively. DEGs and ePAVs indicated that CSA and CSS were clearly differentiated at the population transcriptome level, and they were closely related to abiotic tolerance and secondary metabolite synthesis phenotypes of CSA and CSS based on gene annotations. The regulatory network contained 212 specialized metabolites, 3843 candidate genes, and 3407 eQTLs, highlighting many pleiotropic candidate genes, candidate gene-rich eQTLs, and potential regulators of specialized metabolites. These included important transcription factors in the AP2/ERF-ERF, MYB, WD40, and bHLH families. CsTGY14G0001296, an ortholog of AtANS, appeared to be directly related to variation in proanthocyanins in the tea plant population, and the CsTGY11G0002074 gene encoding F3'5'H was found to contribute to the biased distribution of catechins between pure CSAs and pure CSSs. Together, these results provide a new understanding of the metabolite diversity in tea plants and offer new insights for more effective breeding of better-flavored tea varieties.