Tea plant genomics: achievements, challenges and perspectives

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.

Identification and expression analysis of PUB genes in tea plant exposed to anthracnose pathogen and drought stresses

The plant U-box (PUB) gene family, one of the major ubiquitin ligase families in plants, plays important roles in multiple cellular processes including environmental stress responses and resistance. The function of U-box genes has been well characterized in Arabidopsis and other plants. However, little is known about the tea plant (Camellia sinensis) PUB genes. Here, 89 U-box proteins were identified from the chromosome-scale referenced genome of tea plant. According to the domain organization and phylogenetic analysis, the tea plant PUB family were classified into ten classes, named Class I to X, respectively. Using previously released stress-related RNA-seq data in tea plant, we identified 34 stress-inducible CsPUB genes. Specifically, eight CsPUB genes were expressed differentially under both anthracnose pathogen and drought stresses. Moreover, six of the eight CsPUBs were upregulated in response to these two stresses. Expression profiling performed by qRT-PCR was consistent with the RNA-seq analysis, and stress-related cis-acting elements were identified in the promoter regions of the six upregulated CsPUB genes. These results strongly implied the putative functions of U-box ligase genes in response to biotic and abiotic stresses in tea plant.

Parentage analysis of tea cultivars in Japan based on simple sequence repeat markers

Tea cultivars have been bred by individual selection of landraces and by crossbreeding, but the validation of the parentage is limited. In this study, we performed parentage analysis of 79 tea cultivars in Japan based on SSR markers to confirm or identify the parent-offspring relationships among them. The effectiveness of nine SSR markers for parentage analysis was validated by comparing them to the existing cleaved amplified polymorphic sequence markers. The former markers were detectable more alleles than the latter. Simulation of parentage analysis of the tea cultivars predicted biparental origins for 12 cultivars ('Houshun', 'Mie ryokuhou no. 1', 'Surugawase', 'Tenmyo', 'Yamanoibuki', 'Harumidori', 'Koushun', 'Minekaori', 'Okumusashi', 'Saemidori', 'Sofu', and 'Toyoka'), in the first five of which candidate parents of yet-to-be-defined pedigree were newly identified. Comparisons of a total of 41 SSR genotypes confirmed the newly-identified parentages of 'Asahi' for 'Tenmyo', 'Rokurou' for 'Houshun', 'Surugawase', and 'Yamanoibuki', and 'Yamatomidori' for 'Mie ryokuhou no. 1'. The maternity of seven cultivars out of the 12 was also confirmed with chloroplast DNA sequences. Uniparental origins were confirmed for 25 cultivars. This parentage analysis has improved our knowledge of tea pedigrees and will aid in the development of new cultivars.

Genetic relationship analysis and molecular fingerprint identification of the tea germplasms from Guangxi Province, China

The tea plant (Camellia sinensis) is an evergreen woody plant with a high economic value. Guangxi Province is adjacent to the origin center of the tea plant in southern China. It has abundant germplasm resources and is a historically important tea-producing province. However, there is little information about the genetic diversity, genetic introgression, and fingerprints of the tea germplasms from Guangxi Province. Here, we constructed a phylogenetic tree of 126 tea accessions from Guangxi Province using 20 SSR markers. This tree classified these tea accessions into three subgroups containing 19, 47, and 60 members, respectively. High genetic similarity was observed among the three subgroups, and the genetic diversity of the populations was ranked as follows: subgroup 3 > subgroup 2 > subgroup 1. Furthermore, we analyzed the genetic relationships among 168 tea accessions from Guangxi Province and neighboring provinces. The results of the population structure analysis were highly consistent with the clustering results, and genetic introgression was observed. We identified six SSRs as the core marker set, because they could sufficiently distinguish between all 126 tea accessions. The results provide a crucial theoretical basis for utilization and protection of tea germplasms from Guangxi Province, and will help improve the breeding and popularization of elite tea cultivars.

Improving the Tea Withering Process Using Ethylene or UV-C

Using a combination of biochemical, transcriptomic, and physiological analyses, we elucidated the mechanisms of physical and chemical withering of tea shoots subjected to UV-C and ethylene treatments. UV-C irradiation (15 kJ m^-2) initiated oxidation of catechins into theaflavins, increasing theaflavin-3-monogallate and theaflavin digallate by 5- and 13.2-4.4-fold, respectively, at the end of withering. Concomitantly, a rapid change to brown/red, an increase in electrolyte leakage, and the upregulation of peroxidases (viz. Px2, Px4, and Px6) and polyphenol oxidases (PPO-1) occurred. Exogenous ethylene significantly increased the metabolic rate (40%) and moisture loss (30%) compared to control during simulated withering (12 h at 25 °C) and upregulated transcripts associated with responses to dehydration and abiotic stress, such as those in the ethylene signaling pathway (viz. EIN4-like, EIN3-FBox1, and ERFs). Incorporating ethylene during withering could shorten the tea manufacturing process, while UV-C could enhance the accumulation of flavor-related compounds.

Recent insights on tea metabolites, their biosynthesis and chemo-preventing effects: A review

Tea manufactured from the cultivated shoots of Camellia sinensis (L.) O. Kuntze is the most commonly consumed nonalcoholic drink around the world. Tea is an agro-based, environmentally sustainable, labor-intensive, job-generating, and export-oriented industry in many countries. Tea includes phenolic compounds, flavonoids, alkaloids, vitamins, enzymes, crude fibers, protein, lipids, and carbohydrates, among other biochemical constituents. This review described the nature of tea metabolites, their biosynthesis and accumulation with response to various factors. The therapeutic application of various metabolites of tea against microbial diseases, cancer, neurological, and other metabolic disorders was also discussed in detail. The seasonal variation, cultivation practices and genetic variability influence tea metabolite synthesis. Tea biochemical constituents, especially polyphenols and its integral part catechin metabolites, are broadly focused on potential applicability for their action against various diseases. In addition to this, tea also contains bioactive flavonoids that possess health-beneficial effects. The catechin fractions, epigallocatechin 3-gallate and epicatechin 3-gallate, are the main components of tea that has strong antioxidant and medicinal properties. The synergistic function of natural tea metabolites with synthetic drugs provides effective protection against various diseases. Furthermore, the application of nanotechnologies enhanced bioavailability, enhancing the therapeutic potential of natural metabolites against numerous diseases and pathogens.

An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis)

Tea, coffee, and cocoa are the three most popular nonalcoholic beverages in the world and have extremely high economic and cultural value. The genomes of four tea plant varieties have recently been sequenced, but there is some debate regarding the characterization of a whole-genome duplication (WGD) event in tea plants. Whether the WGD in the tea plant is shared with other plants in order Ericales and how it contributed to tea plant evolution remained unanswered. Here we re-analyzed the tea plant genome and provided evidence that tea experienced only WGD event after the core-eudicot whole-genome triplication (WGT) event. This WGD was shared by the Polemonioids-Primuloids-Core Ericales (PPC) sections, encompassing at least 17 families in the order Ericales. In addition, our study identified eight pairs of duplicated genes in the catechins biosynthesis pathway, four pairs of duplicated genes in the theanine biosynthesis pathway, and one pair of genes in the caffeine biosynthesis pathway, which were expanded and retained following this WGD. Nearly all these gene pairs were expressed in tea plants, implying the contribution of the WGD. This study shows that in addition to the role of the recent tandem gene duplication in the accumulation of tea flavor-related genes, the WGD may have been another main factor driving the evolution of tea flavor.

Integrating Omics and Gene Editing Tools for Rapid Improvement of Traditional Food Plants for Diversified and Sustainable Food Security

Indigenous communities across the globe, especially in rural areas, consume locally available plants known as Traditional Food Plants (TFPs) for their nutritional and health-related needs. Recent research shows that many TFPs are highly nutritious as they contain health beneficial metabolites, vitamins, mineral elements and other nutrients. Excessive reliance on the mainstream staple crops has its own disadvantages. Traditional food plants are nowadays considered important crops of the future and can act as supplementary foods for the burgeoning global population. They can also act as emergency foods in situations such as COVID-19 and in times of other pandemics. The current situation necessitates locally available alternative nutritious TFPs for sustainable food production. To increase the cultivation or improve the traits in TFPs, it is essential to understand the molecular basis of the genes that regulate some important traits such as nutritional components and resilience to biotic and abiotic stresses. The integrated use of modern omics and gene editing technologies provide great opportunities to better understand the genetic and molecular basis of superior nutrient content, climate-resilient traits and adaptation to local agroclimatic zones. Recently, realizing the importance and benefits of TFPs, scientists have shown interest in the prospection and sequencing of TFPs for their improvements, cultivation and mainstreaming. Integrated omics such as genomics, transcriptomics, proteomics, metabolomics and ionomics are successfully used in plants and have provided a comprehensive understanding of gene-protein-metabolite networks. Combined use of omics and editing tools has led to successful editing of beneficial traits in several TFPs. This suggests that there is ample scope for improvement of TFPs for sustainable food production. In this article, we highlight the importance, scope and progress towards improvement of TFPs for valuable traits by integrated use of omics and gene editing techniques.

TeaAS: a comprehensive database for alternative splicing in tea plants (Camellia sinensis)

Alternative splicing (AS) increases the diversity of transcripts and proteins through the selection of different splice sites and plays an important role in the growth, development and stress tolerance of plants. With the release of the reference genome of the tea plant (Camellia sinensis) and the development of transcriptome sequencing, researchers have reported the existence of AS in tea plants. However, there is a lack of a platform, centered on different RNA-seq datasets, that provides comprehensive information on AS.To facilitate access to information on AS and reveal the molecular function of AS in tea plants, we established the first comprehensive AS database for tea plants (TeaAS, http://www.teaas.cn/index.php ). In this study, 3.96 Tb reads from 66 different RNA-seq datasets were collected to identify AS events. TeaAS supports four methods of retrieval of AS information based on gene ID, gene name, annotation (non-redundant/Kyoto encyclopedia of genes and genomes/gene ontology annotation or chromosomal location) and RNA-seq data. It integrates data pertaining to genome annotation, type of AS event, transcript sequence, and isoforms expression levels from 66 RNA-seq datasets. The AS events resulting from different environmental conditions and that occurring in varied tissue types, and the expression levels of specific transcripts can be clearly identified through this online database. Moreover, it also provides two useful tools, Basic Local Alignment Search Tool and Generic Genome Browser, for sequence alignment and visualization of gene structure.The features of the TeaAS database make it a comprehensive AS bioinformatics platform for researchers, as well as a reference for studying AS events in woody crops. It could also be helpful for revealing the novel biological functions of AS in gene regulation in tea plants.

Genetic Diversity of Ancient Camellia sinensis (L.) O.Kuntze in Sandu County of Guizhou Province in China

The ancient tea plant germplasm is an important resource for breeding new tea plant varieties and has great economic value. However, due to man-made and natural disturbances, it has become endangered. In order to have a better management of the conserved tea plant germplasm, it is a requirement to understand the genetic and phenotypic diversity. The aim of this study was to evaluate the genetic and phenotypic diversity of 145 ancient tea plant germplasm resources from five populations in Sandu County of Guizhou province in China. To explore the population genetics of tea plant, we successfully identified 15 simple sequence repeat (SSR) markers, which were highly polymorphic. Additionally, we applied traditional phenotypic methods to evaluate the tea plant diversity. The results suggested that the genetic and phenotypic diversity were relatively high. A total of 96 alleles were identified, and the mean polymorphic information content (PIC) value was found to be 0.66. The analysis of molecular variance (AMOVA) showed that genetic variation within the populations was greater than among the populations. Overall, our results are the valuable baseline data in developing more efficient management and breeding plans for one of the most popular non-alcoholic beverage crops, the tea plant species.

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.

Complete chloroplast genomes of Achnatherum inebrians and comparative analyses with related species from Poaceae

This article reports the complete chloroplast genome of Achnatherum inebrians, a poisonous herb that is widely distributed in the rangelands of Northern China. The genome is 137 714 bp in total and consists of a large single‐copy (81 758 bp) region and small single‐copy (12 682 bp) region separated by a pair of inverted repeats (21 637 bp). The genome contains 130 genes, including 84 protein‐coding genes, 38 tRNA genes and 8 ribosomal RNA genes, and the guanine + cytosine content is 36.17%. We subsequently performed comparative analysis of complete genomes from A. inebrians and other Poaceae‐related species from GenBank. Thirty‐eight simple sequence repeats were identified, further demonstrating rapid evolution in Poaceae. Finally, the phylogenetic trees of 37 species of Poaceae and 2 species of Amaranthaceae were constructed by using maximum likelihood and Bayesian inference methods, based on the genes of the complete chloroplast genome. We identified hotspots that can be used as molecular markers and barcodes for phylogenetic analysis, as well as for species identification. Phylogenetic analysis indicated that A. inebrians is a member of the genus Stipa rather than Achnatherum.

CsbZIP1-CsMYB12 mediates the production of bitter-tasting flavonols in tea plants (Camellia sinensis) through a coordinated activator-repressor network

Under high light conditions or UV radiation, tea plant leaves produce more flavonols, which contribute to the bitter taste of tea; however, neither the flavonol biosynthesis pathways nor the regulation of their production are well understood. Intriguingly, tea leaf flavonols are enhanced by UV-B but reduced by shading treatment. CsFLS, CsUGT78A14, CsMYB12, and CsbZIP1 were upregulated by UV-B radiation and downregulated by shading. CsMYB12 and CsbZIP1 bound to the promoters of CsFLS and CsUGT78A14, respectively, and activated their expression individually. CsbZIP1 positively regulated CsMYB12 and interacted with CsMYB12, which specifically activated flavonol biosynthesis. Meanwhile, CsPIF3 and two MYB repressor genes, CsMYB4 and CsMYB7, displayed expression patterns opposite to that of CsMYB12. CsMYB4 and CsMYB7 bound to CsFLS and CsUGT78A14 and repressed their CsMYB12-activated expression. While CsbZIP1 and CsMYB12 regulated neither CsMYB4 nor CsMYB7, CsMYB12 interacted with CsbZIP1, CsMYB4, and CsMYB7, but CsbZIP1 did not physically interact with CsMYB4 or CsMYB7. Finally, CsPIF3 bound to and activated CsMYB7 under shading to repress flavonol biosynthesis. These combined results suggest that UV activation and shading repression of flavonol biosynthesis in tea leaves are coordinated through a complex network involving CsbZIP1 and CsPIF3 as positive MYB activators and negative MYB repressors, respectively. The study thus provides insight into the regulatory mechanism underlying the production of bitter-tasting flavonols in tea plants.

Genetic basis of high aroma and stress tolerance in the oolong tea cultivar genome

Tea plants (Camellia sinensis) are commercially cultivated in >60 countries, and their fresh leaves are processed into tea, which is the most widely consumed beverage in the world. Although several chromosome-level tea plant genomes have been published, they collapsed the two haplotypes and ignored a large number of allelic variations that may underlie important biological functions in this species. Here, we present a phased chromosome-scale assembly for an elite oolong tea cultivar, "Huangdan", that is well known for its high levels of aroma. Based on the two sets of haplotype genome data, we identified numerous genetic variations and a substantial proportion of allelic imbalance related to important traits, including aroma- and stress-related alleles. Comparative genomics revealed extensive structural variations as well as expansion of some gene families, such as terpene synthases (TPSs), that likely contribute to the high-aroma characteristics of the backbone parent, underlying the molecular basis for the biosynthesis of aroma-related chemicals in oolong tea. Our results uncovered the genetic basis of special features of this oolong tea cultivar, providing fundamental genomic resources to study evolution and domestication for the economically important tea crop.

Genome-wide identification, evolutionary relationship and expression analysis of AGO, DCL and RDR family genes in tea

Three gene families in plants viz. Argonaute (AGOs), Dicer-like (DCLs) and RNA dependent RNA polymerase (RDRs) constitute the core components of small RNA mediated gene silencing machinery. The present study endeavours to identify members of these gene families in tea and to investigate their expression patterns in different tissues and various stress regimes. Using genome-wide analysis, we have identified 18 AGOs, 5 DCLs and 9 RDRs in tea, and analyzed their phylogenetic relationship with orthologs of Arabidopsis thaliana. Gene expression analysis revealed constitutive expression of CsAGO1 in all the studied tissues and stress conditions, whereas CsAGO10c showed most variable expression among all the genes. CsAGO10c gene was found to be upregulated in tissues undergoing high meristematic activity such as buds and roots, as well as in Exobasidium vexans infected samples. CsRDR2 and two paralogs of CsAGO4, which are known to participate in biogenesis of hc-siRNAs, showed similarities in their expression levels in most of the tea plant tissues. This report provides first ever insight into the important gene families involved in biogenesis of small RNAs in tea. The comprehensive knowledge of these small RNA biogenesis purveyors can be utilized for tea crop improvement aimed at stress tolerance and quality enhancement.

Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites

The characteristic secondary metabolites in tea (theanine, caffeine, and catechins) are important factors contributing to unique tea flavors. However, there has been relatively little research on molecular markers related to these metabolites. Thus, we conducted a genome-wide association analysis of the levels of these tea flavor-related metabolites in three seasons. The theanine, caffeine, and catechin levels in Population 1 comprising 191 tea plant germplasms were examined, which revealed that their heritability exceeded 0.5 in the analyzed seasons, with the following rank order (highest to lowest heritabilities): (+)-catechin > (-)-gallocatechin gallate > caffeine = (-)-epicatechin > (-)-epigallocatechin-3-gallate > theanine > (-)-epigallocatechin > (-)-epicatechin-3-gallate > catechin gallate > (+)-gallocatechin. The SNPs detected by amplified-fragment SNP and methylation sequencing divided Population 1 into three groups and seven subgroups. An association analysis yielded 307 SNP markers related to theanine, caffeine, and catechins that were common to all three seasons. Some of the markers were pleiotropic. The functional annotation of 180 key genes at the SNP loci revealed that FLS, UGT, MYB, and WD40 domain-containing proteins, as well as ATP-binding cassette transporters, may be important for catechin synthesis. KEGG and GO analyses indicated that these genes are associated with metabolic pathways and secondary metabolite biosynthesis. Moreover, in Population 2 (98 tea plant germplasm resources), 30 candidate SNPs were verified, including 17 SNPs that were significantly or extremely significantly associated with specific metabolite levels. These results will provide a foundation for future research on important flavor-related metabolites and may help accelerate the breeding of new tea varieties.

Transcriptional analysis reveals key insights into seasonal induced anthocyanin degradation and leaf color transition in purple tea (Camellia sinensis (L.) O. Kuntze)

Purple-tea, an anthocyanin rich cultivar has recently gained popularity due to its health benefits and captivating leaf appearance. However, the sustainability of purple pigmentation and anthocyanin content during production period is hampered by seasonal variation. To understand seasonal dependent anthocyanin pigmentation in purple tea, global transcriptional and anthocyanin profiling was carried out in tea shoots with two leaves and a bud harvested during in early (reddish purple: S1_RP), main (dark gray purple: S2_GP) and backend flush (moderately olive green: S3_G) seasons. Of the three seasons, maximum accumulation of total anthocyanin content was recorded in S2_GP, while least amount was recorded during S3_G. Reference based transcriptome assembly of 412 million quality reads resulted into 71,349 non-redundant transcripts with 6081 significant differentially expressed genes. Interestingly, key DEGs involved in anthocyanin biosynthesis [PAL, 4CL, F3H, DFR and UGT/UFGT], vacuolar trafficking [ABC, MATE and GST] transcriptional regulation [MYB, NAC, bHLH, WRKY and HMG] and Abscisic acid signaling pathway [PYL and PP2C] were significantly upregulated in S2_GP. Conversely, DEGs associated with anthocyanin degradation [Prx and lac], repressor TFs and key components of auxin and ethylene signaling pathways [ARF, AUX/IAA/SAUR, ETR, ERF, EBF1/2] exhibited significant upregulation in S3_G, correlating positively with reduced anthocyanin content and purple coloration. The present study for the first-time elucidated genome-wide transcriptional insights and hypothesized the involvement of anthocyanin biosynthesis activators/repressor and anthocyanin degrading genes via peroxidases and laccases during seasonal induced leaf color transition in purple tea. Futuristically, key candidate gene(s) identified here can be used for genetic engineering and molecular breeding of seasonal independent anthocyanin-rich tea cultivars.

Understanding the role of miRNAs for improvement of tea quality and stress tolerance

MicroRNAs (miRNAs) are an emerging class of small non-coding RNAs that exhibit important role in regulation of gene expression, mostly through the mechanism of cleavage and/or inhibition of translation of target mRNAs during or after transcription. Although much has been unravelled about the role of miRNAs in diverse biological processes like maintenance of functional integrity of genes and genome, growth and development, metabolism, and adaptive responses towards biotic and abiotic stresses in plants, not much is known on their specific roles in majority of cash crops - an area of investigation with potentially significant and gainful economic implications. Tea (Camellia sinensis) is globally the second most consumed beverage after water and its cultivation has major agro-economic and social ramifications. In recent years, global tea production has been greatly challenged by many biotic and abiotic stress factors and a deeper understanding of molecular processes regulating stress adaptation in this largely under investigated crop stands to significantly facilitate potential crop improvement strategies towards durable stress tolerance. This review endeavours to highlight recent advances in our understanding of the role of miRNAs in regulating stress tolerance traits in tea plant with additional focus on their role in determining tea quality attributes.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Comparative Expression Analysis of Stress-Inducible Candidate Genes in Response to Cold and Drought in Tea Plant [Camellia sinensis (L.) Kuntze]

Cold and drought are two of the most severe threats affecting the growth and productivity of the tea plant, limiting its global spread. Both stresses cause osmotic changes in the cells of the tea plant by decreasing their water potential. To develop cultivars that are tolerant to both stresses, it is essential to understand the genetic responses of tea plant to these two stresses, particularly in terms of the genes involved. In this study, we combined literature data with interspecific transcriptomic analyses (using Arabidopsis thaliana and Solanum lycopersicum) to choose genes related to cold tolerance. We identified 45 stress-inducible candidate genes associated with cold and drought responses in tea plants based on a comprehensive homologous detection method. Of these, nine were newly characterized by us, and 36 had previously been reported. The gene network analysis revealed upregulated expression in ICE1-related cluster of bHLH factors, HSP70/BAM5 connected genes (hexokinases, galactinol synthases, SnRK complex, etc.) indicating their possible co-expression. Using qRT-PCR we revealed that 10 genes were significantly upregulated in response to both cold and drought in tea plant: HSP70, GST, SUS1, DHN1, BMY5, bHLH102, GR-RBP3, ICE1, GOLS1, and GOLS3. SnRK1.2, HXK1/2, bHLH7/43/79/93 were specifically upregulated in cold, while RHL41, CAU1, Hydrolase22 were specifically upregulated in drought. Interestingly, the expression of CIP was higher in the recovery stage of both stresses, indicating its potentially important role in plant recovery after stress. In addition, some genes, such as DHN3, bHLH79, PEI54, SnRK1.2, SnRK1.3, and Hydrolase22, were significantly positively correlated between the cold and drought responses. CBF1, GOLS1, HXK2, and HXK3, by contrast, showed significantly negative correlations between the cold and drought responses. Our results provide valuable information and robust candidate genes for future functional analyses intended to improve the stress tolerance of the tea plant and other species.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.

HAK/KUP/KT family potassium transporter genes are involved in potassium deficiency and stress responses in tea plants (Camellia sinensis L.): expression and functional analysis

BACKGROUND

Tea plant is one of the most important non-alcoholic beverage crops worldwide. While potassium (K+) is an essential macronutrient and greatly affects the growth and development of plants, the molecular mechanism underlying K+ uptake and transport in tea plant root, especially under limited-K+ conditions, is still poorly understood. In plants, HAK/KUP/KT family members play a crucial role in K+ acquisition and translocation, growth and development, and response to stresses. Nevertheless, the biological functions of these genes in tea plant are still in mystery, especially their roles in K+ uptake and stress responses.

RESULTS

In this study, a total of 21 non-redundant HAK/KUP/KT genes (designated as CsHAKs) were identified in tea plant. Phylogenetic and structural analysis classified the CsHAKs into four clusters (I, II, III, IV), containing 4, 8, 4 and 5 genes, respectively. Three major categories of cis-acting elements were found in the promoter regions of CsHAKs. Tissue-specific expression analysis indicated extremely low expression levels in various tissues of cluster I CsHAKs with the exception of a high root expression of CsHAK4 and CsHAK5, a constitutive expression of clusters II and III CsHAKs, and a moderate cluster IV CsHAKs expression. Remarkably, the transcript levels of CsHAKs in roots were significantly induced or suppressed after exposure to K+ deficiency, salt and drought stresses, and phytohormones treatments. Also notably, CsHAK7 was highly expressed in all tissues and was further induced under various stress conditions. Therefore, functional characterization of CsHAK7 was performed, and the results demostrated that CsHAK7 locates on plasma membrane and plays a key role in K+ transport in yeast. Taken together, the results provide promising candidate CsHAKs for further functional studies and contribute to the molecular breeding for new tea plants varieties with highly efficient utilization of K+.

CONCLUSION

This study demonstrated the first genome-wide analysis of CsHAK family genes of tea plant and provides a foundation for understanding the classification and functions of the CsHAKs in tea plants.

The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation

Tea plant is an important economic crop, which is used to produce the world's oldest and most widely consumed tea beverages. Here, we present a high-quality reference genome assembly of the tea plant (Camellia sinensis var. sinensis) consisting of 15 pseudo-chromosomes. LTR retrotransposons (LTR-RTs) account for 70.38% of the genome, and we present evidence that LTR-RTs play critical roles in genome size expansion and the transcriptional diversification of tea plant genes through preferential insertion in promoter regions and introns. Genes, particularly those coding for terpene biosynthesis proteins, associated with tea aroma and stress resistance were significantly amplified through recent tandem duplications and exist as gene clusters in tea plant genome. Phylogenetic analysis of the sequences of 81 tea plant accessions with diverse origins revealed three well-differentiated tea plant populations, supporting the proposition for the southwest origin of the Chinese cultivated tea plant and its later spread to western Asia through introduction. Domestication and modern breeding left significant signatures on hundreds of genes in the tea plant genome, particularly those associated with tea quality and stress resistance. The genomic sequences of the reported reference and resequenced tea plant accessions provide valuable resources for future functional genomics study and molecular breeding of improved cultivars of tea plants.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance

BACKGROUND

Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.).

RESULTS

In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca²⁺. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response.

CONCLUSIONS

A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

Integrative transcriptome, proteome, and microRNA analysis reveals the effects of nitrogen sufficiency and deficiency conditions on theanine metabolism in the tea plant (Camellia sinensis)

Nitrogen (N) is associated with amino acid metabolism in higher plants. Theanine is an important amino acid in tea plants. To explore the relationship between theanine metabolism and N conditions, we examined the differentially expressed genes (DEGs), proteins (DEPs), and microRNAs (DEMs) involved in theanine metabolism in tea plant shoots and roots under N sufficiency and deficiency conditions. Transcriptome, proteome, and microRNA analyses were performed on tea plant shoots and roots under N sufficiency and deficiency conditions. The contents of theanine, expression levels of genes involved in theanine metabolism, contents of proteinogenic amino acids, and activity of enzymes were analyzed. The DEP-DEG correlation pairs and negative DEM-DEG interactions related to theanine metabolism were identified based on correlation analyses. The expression profiles of DEGs and negative DEM-DEG pairs related to theanine biosynthesis were consistent with the sequencing results. Our results suggest that the molecular and physiological mechanism of theanine accumulation is significantly affected by N sufficiency and deficiency conditions. The DEGs, DEPs, and DEMs and the activity of the enzymes involved in theanine biosynthesis might play vital roles in theanine accumulation under N sufficiency and deficiency conditions in the shoots and roots of tea plants.

Exploring plant metabolic genomics: chemical diversity, metabolic complexity in the biosynthesis and transport of specialized metabolites with the tea plant as a model

The diversity and complexity of secondary metabolites in tea plants contribute substantially to the popularity of tea, by determining tea flavors and their numerous health benefits. The most significant characteristics of tea plants are that they concentrate the complex plant secondary metabolites into one leaf: flavonoids, alkaloids, theanine, volatiles, and saponins. Many fundamental questions regarding tea plant secondary metabolism remain unanswered. This includes how tea plants accumulate high levels of monomeric galloylated catechins, unlike the polymerized flavan-3-ols in most other plants, as well as how they are evolved to selectively synthesize theanine and caffeine, and how tea plants properly transport and store these cytotoxic products and then reuse them in defense. Tea plants coordinate many metabolic pathways that simultaneously take place in young tea leaves in response to both developmental and environmental cues. With the available genome sequences of tea plants and high-throughput metabolomic tools as great platforms, it is of particular interest to launch metabolic genomics studies using tea plants as a model system. Plant metabolic genomics are to investigate all aspects of plant secondary metabolism at the genetic, genome, and molecular levels. This includes plant domestication and adaptation, divergence and convergence of secondary metaboloic pathways. The biosynthesis, transport, storage, and transcriptional regulation mechanisms of all metabolites are of core interest in the plant as a whole. This review highlights relevant contexts of metabolic genomics, outstanding questions, and strategies for answering them, with aim to guide future research for genetic improvement of nutrition quality for healthier plant foods.

Revealing Distinctions in Genetic Diversity and Adaptive Evolution Between Two Varieties of Camellia sinensis by Whole-Genome Resequencing

Camellia sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA) are the two most economically important tea varieties. They have different characteristics and geographical distribution. Their genetic diversity and differentiation are unclear. Here, we identified 18,903,625 single nucleotide polymorphisms (SNPs) and 7,314,133 insertion-deletion mutations (indels) by whole-genome resequencing of 30 cultivated and three wild related species. Population structure and phylogenetic tree analyses divided the cultivated accessions into CSS and CSA containing 6,440,419 and 6,176,510 unique variations, respectively. The CSS subgroup possessed higher genetic diversity and was enriched for rare alleles. The CSA subgroup had more non-synonymous mutations and might have experienced a greater degree of balancing selection. The evolution rate (dN/dS) and KEGG enrichment indicated that genes involved in the synthesis and metabolism of flavor substances were positively selected in both CSS and CSA subpopulations. However, there are extensive genome differentiation regions (2959 bins and approximately 148 M in size) between the two subgroups. Compared with CSA (141 selected regions containing 124 genes), the CSS subgroup (830 selected regions containing 687 genes) displayed more selection regions potentially related to environmental adaptability. Fifty-three pairs of polymorphic indel markers were developed. Some markers were located in hormone-related genes with distinct alleles in the two cultivated subgroups. These identified variations and selected regions provide clues for the differentiation and adaptive evolution of tea varieties. The newly developed indel markers will be valuable in further genetic research on tea plants.