The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant

The Serine Acetyltransferase (SAT) Gene Family in Tea Plant (Camellia sinensis): Identification, Classification and Expression Analysis under Salt Stress

Cysteine plays a pivotal role in the sulfur metabolism network of plants, intimately influencing the conversion rate of organic sulfur and the plant's capacity to withstand abiotic stresses. In tea plants, the serine acetyltransferase (SAT) genes emerge as a crucial regulator of cysteine metabolism, albeit with a notable lack of comprehensive research. Utilizing Hidden Markov Models, we identified seven CssSATs genes within the tea plant genome. The results of the bioinformatics analysis indicate that these genes exhibit an average molecular weight of 33.22 kD and cluster into three distinct groups. Regarding gene structure, CssSAT1 stands out with ten exons, significantly more than its family members. In the promoter regions, cis-acting elements associated with environmental responsiveness and hormone induction predominate, accounting for 34.4% and 53.1%, respectively. Transcriptome data revealed intricate expression dynamics of CssSATs under various stress conditions (e.g., PEG, NaCl, Cold, MeJA) and their tissue-specific expression patterns in tea plants. Notably, qRT-PCR analysis indicated that under salt stress, CssSAT1 and CssSAT3 expression levels markedly increased, whereas CssSAT2 displayed a downregulatory trend. Furthermore, we cloned CssSAT1-CssSAT3 genes and constructed corresponding prokaryotic expression vectors. The resultant recombinant proteins, upon induction, significantly enhanced the NaCl tolerance of Escherichia coli BL21, suggesting the potential application of CssSATs in bolstering plant stress resistance. These findings have enriched our comprehension of the multifaceted roles played by CssSATs genes in stress tolerance mechanisms, laying a theoretical groundwork for future scientific endeavors and research pursuits.

Chromosome-scale genome assembly of oil-tea tree Camellia crapnelliana

Camellia crapnelliana Tutch., belonging to the Theaceae family, is an excellent landscape tree species with high ornamental values. It is particularly an important woody oil-bearing plant species with high ecological, economic, and medicinal values. Here, we first report the chromosome-scale reference genome of C. crapnelliana with integrated technologies of SMRT, Hi-C and Illumina sequencing platforms. The genome assembly had a total length of ~2.94 Gb with contig N50 of ~67.5 Mb, and ~96.34% of contigs were assigned to 15 chromosomes. In total, we predicted 37,390 protein-coding genes, ~99.00% of which could be functionally annotated. The chromosome-scale genome of C. crapnelliana will become valuable resources for understanding the genetic basis of the fatty acid biosynthesis, and greatly facilitate the exploration and conservation of C. crapnelliana.

New insights into evolution and functional diversification of Camellia sinensis LRR-RLKs

Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subgroup of receptor-like kinases (RLKs) in plants. While some LRR-RLK members play a role in regulating various plant growth processes related to morphogenesis, disease resistance, and stress response, the functions of most LRR-RLK genes remain unclear. In this study, we identified 397 LRR-RLK genes from the genome of Camellia sinensis and categorized them into 16 subfamilies. Approximately 62% of CsLRR-RLK genes are situated in regions resulting from segmental duplications, suggesting that the expansion of CsLRR-RLK genes is due to segmental duplications. Analysis of gene expression patterns revealed differential expression of CsLRR-RLK genes across different tissues and in response to stress. Furthermore, we demonstrated that CssEMS1 localizes to the cell membrane and can complement Arabidopsis ems1 mutant. This study is the initial in-depth evolutionary examination of LRR-RLKs in tea and provides a basis for future investigations into their functionality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01458-1.

Circadian rhythm response and its effect on photosynthetic characteristics of the Lhcb family genes in tea plant

BACKGROUND

The circadian clock, also known as the circadian rhythm, is responsible for predicting daily and seasonal changes in the environment, and adjusting various physiological and developmental processes to the appropriate times during plant growth and development. The circadian clock controls the expression of the Lhcb gene, which encodes the chlorophyll a/b binding protein. However, the roles of the Lhcb gene in tea plant remain unclear.

RESULTS

In this study, a total of 16 CsLhcb genes were identified based on the tea plant genome, which were distributed on 8 chromosomes of the tea plant. The promoter regions of CsLhcb genes have a variety of cis-acting elements including hormonal, abiotic stress responses and light response elements. The CsLhcb family genes are involved in the light response process in tea plant. The photosynthetic parameter of tea leaves showed rhythmic changes during the two photoperiod periods (48 h). Stomata are basically open during the day and closed at night. Real-time quantitative PCR results showed that most of the CsLhcb family genes were highly expressed during the day, but were less expressed at night.

CONCLUSIONS

Results indicated that CsLhcb genes were involved in the circadian clock process of tea plant, it also provided potential references for further understanding of the function of CsLhcb gene family in tea plant.

Bulked Segregant RNA-Seq Reveals Different Gene Expression Patterns and Mutant Genes Associated with the Zigzag Pattern of Tea Plants (Camellia sinensis)

The unique zigzag-patterned tea plant is a rare germplasm resource. However, the molecular mechanism behind the formation of zigzag stems remains unclear. To address this, a BC1 genetic population of tea plants with zigzag stems was studied using histological observation and bulked segregant RNA-seq. The analysis revealed 1494 differentially expressed genes (DEGs) between the upright and zigzag stem groups. These DEGs may regulate the transduction and biosynthesis of plant hormones, and the effects on the phenylpropane biosynthesis pathways may cause the accumulation of lignin. Tissue sections further supported this finding, showing differences in cell wall thickness between upright and curved stems, potentially due to lignin accumulation. Additionally, 262 single-nucleotide polymorphisms (SNPs) across 38 genes were identified as key SNPs, and 5 genes related to zigzag stems were identified through homologous gene function annotation. Mutations in these genes may impact auxin distribution and content, resulting in the asymmetric development of vascular bundles in curved stems. In summary, we identified the key genes associated with the tortuous phenotype by using BSR-seq on a BC1 population to minimize genetic background noise.

Genome-wide identification and comparative analysis of YABBY transcription factors in oil tea and tea tree

The plant-specific transcription factor gene family, YABBY, plays an important role in plant development and stress response. Although YABBY genes have been identified in numerous species, a comprehensive characterization of YABBYs in tea tree and oil tea has been lacking. In this study, ten and three YABBY genes were identified in Camellia sinensis and C. oleifera, respectively. YABBY proteins could be divided into five subfamilies. Most YABBY genes in the same clade had similar structures and conserved motifs. Protein evolutionary analysis revealed that FIL/YAB3 displayed high conservation in all positions, followed by INO, YAB2, YAB5, and CRC. Specific site analysis suggested that the YABBY family was polyphyletic during the evolution. Compared to C. oleifera, two segmentally duplicated gene pairs were formed in C. sinensis during recent WGD events generated 30.69 and 45.08 Mya, respectively. Cis-acting element indicated that most YABBY genes contain box4, ARE, and MYB elements. A total of 120 SSR loci were found within CsYABBYs, consisting of six types, while 48 SSR loci were identified within CoYABBY, consisting of three types. Transcriptome results revealed that CRC and INO clades were specifically expressed in floral organs. The expression of CsYABBY10 and CsYABBY5 was significantly up-regulated under drought and salt treatments, respectively, as confirmed by qRT-PCR. CoYABBY genes were more susceptible to salt stress, as CoYABBY3 increased by about 15-fold. Furthermore, functional differentiation may have occurred in duplicated genes. These discoveries provide important information for further research on YABBYs in tea tree and oil tea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03940-9.

Depicting the genetic and metabolic panorama of chemical diversity in the tea plant

As a frequently consumed beverage worldwide, tea is rich in naturally important bioactive metabolites. Combining genetic, metabolomic and biochemical methodologies, here, we present a comprehensive study to dissect the chemical diversity in tea plant. A total of 2837 metabolites were identified at high-resolution with 1098 of them being structurally annotated and 63 of them were structurally identified. Metabolite-based genome-wide association mapping identified 6199 and 7823 metabolic quantitative trait loci (mQTL) for 971 and 1254 compounds in young leaves (YL) and the third leaves (TL), respectively. The major mQTL (i.e., P < 1.05 × 10-5, and phenotypic variation explained (PVE) > 25%) were further interrogated. Through extensive annotation of the tea metabolome as well as network-based analysis, this study broadens the understanding of tea metabolism and lays a solid foundation for revealing the natural variations in the chemical composition of the tea plant. Interestingly, we found that galloylations, rather than hydroxylations or glycosylations, were the largest class of conversions within the tea metabolome. The prevalence of galloylations in tea is unusual, as hydroxylations and glycosylations are typically the most prominent conversions of plant specialized metabolism. The biosynthetic pathway of flavonoids, which are one of the most featured metabolites in tea plant, was further refined with the identified metabolites. And we demonstrated the further mining and interpretation of our GWAS results by verifying two identified mQTL (including functional candidate genes CsUGTa, CsUGTb, and CsCCoAOMT) and completing the flavonoid biosynthetic pathway of the tea plant.

Decoding the complete organelle genomic architecture of Stewartia gemmata: an early-diverging species in Theaceae

BACKGROUND

Theaceae, comprising 300 + species, holds significance in biodiversity, economics, and culture, notably including the globally consumed tea plant. Stewartia gemmata, a species of the earliest diverging tribe Stewartieae, is critical to offer insights into Theaceae's origin and evolutionary history.

RESULT

We sequenced the complete organelle genomes of Stewartia gemmata using short/long reads sequencing technologies. The chloroplast genome (158,406 bp) exhibited a quadripartite structure including the large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs); 114 genes encoded 80 proteins, 30 tRNAs, and four rRNAs. The mitochondrial genome (681,203 bp) exhibited alternative conformations alongside a monocyclic structure: 61 genes encoding 38 proteins, 20 tRNAs, three rRNAs, and RNA editing-impacting genes, including ATP6, RPL16, COX2, NAD4L, NAD5, NAD7, and RPS1. Comparative analyses revealed frequent recombination events and apparent rRNA gene gains and losses in the mitochondrial genome of Theaceae. In organelle genomes, the protein-coding genes exhibited a strong A/U bias at codon endings; ENC-GC3 analysis implies selection-driven codon bias. Transposable elements might facilitate interorganelle sequence transfer. Phylogenetic analysis confirmed Stewartieae's early divergence within Theaceae, shedding light on organelle genome characteristics and evolution in Theaceae.

CONCLUSIONS

We studied the detailed characterization of organelle genomes, including genome structure, composition, and repeated sequences, along with the identification of lateral gene transfer (LGT) events and complexities. The discovery of a large number of repetitive sequences and simple sequence repeats (SSRs) has led to new insights into molecular phylogenetic markers. Decoding the Stewartia gemmata organellar genome provides valuable genomic resources for further studies in tea plant phylogenomics and evolutionary biology.

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.

Global transcriptome analysis reveals fungal disease responsive core gene regulatory landscape in tea

Fungal infections are the inevitable limiting factor for productivity of tea. Transcriptome reprogramming recruits multiple regulatory pathways during pathogen infection. A comprehensive meta-analysis was performed utilizing previously reported, well-replicated transcriptomic datasets from seven fungal diseases of tea. The study identified a cumulative set of 18,517 differentially expressed genes (DEGs) in tea, implicated in several functional clusters, including the MAPK signaling pathway, transcriptional regulation, and the biosynthesis of phenylpropanoids. Gene set enrichment analyses under each pathogen stress elucidated that DEGs were involved in ethylene metabolism, secondary metabolism, receptor kinase activity, and various reactive oxygen species detoxification enzyme activities. Expressional fold change of combined datasets highlighting 2258 meta-DEGs shared a common transcriptomic response upon fungal stress in tea. Pervasive duplication events caused biotic stress-responsive core DEGs to appear in multiple copies throughout the tea genome. The co-expression network of meta-DEGs in multiple modules demonstrated the coordination of appropriate pathways, most of which involved cell wall organization. The functional coordination was controlled by a number of hub genes and miRNAs, leading to pathogenic resistance or susceptibility. This first-of-its-kind meta-analysis of host-pathogen interaction generated consensus candidate loci as molecular signatures, which can be associated with future resistance breeding programs in tea.

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.

Genome and whole-genome resequencing of Cinnamomum camphora elucidate its dominance in subtropical urban landscapes

BACKGROUND

Lauraceae is well known for its significant phylogenetic position as well as important economic and ornamental value; however, most evergreen species in Lauraceae are restricted to tropical regions. In contrast, camphor tree (Cinnamomum camphora) is the most dominant evergreen broadleaved tree in subtropical urban landscapes.

RESULTS

Here, we present a high-quality reference genome of C. camphora and conduct comparative genomics between C. camphora and C. kanehirae. Our findings demonstrated the significance of key genes in circadian rhythms and phenylpropanoid metabolism in enhancing cold response, and terpene synthases (TPSs) improved defence response with tandem duplication and gene cluster formation in C. camphora. Additionally, the first comprehensive catalogue of C. camphora based on whole-genome resequencing of 75 accessions was constructed, which confirmed the crucial roles of the above pathways and revealed candidate genes under selection in more popular C. camphora, and indicated that enhancing environmental adaptation is the primary force driving C. camphora breeding and dominance.

CONCLUSIONS

These results decipher the dominance of C. camphora in subtropical urban landscapes and provide abundant genomic resources for enlarging the application scopes of evergreen broadleaved trees.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

The lineage-specific evolution of the oleosin family in Theaceae

Oleosins play essential roles in stabilization of lipid droplets (LDs) and seed oil production. However, evolution of this gene family has not been reported in Theaceae, a large plant family that contains many important tea and oil tea species. In this study, a total of 65 oleosin genes were identified in nine genome-sequenced Theaceae species. Among these genomes, the gene number of oleosin showed significant difference, with Camellia sinensis var. sinensis cv. Shuchazao and Camellia lanceoleosa displayed more oleosin numbers than other species. Phylogenetic analyses revealed that Theaceae oleosin genes were classified into three clades (U, SL, SH) respectively. Proteins within the same clade had similar gene structure and motif composition. Segmental duplication was the primary driving force for the evolution of oleosin genes in Shuchazao (SCZ), Huangdan (HD), C.lanceoleosa (Cla), and wild tea (DASZ). Synteny analysis showed that most oleosin genes displayed inter-species synteny among tea and oil tea species. Expression analysis demonstrated that oleosin genes were specifically expressed in seed and kernel of Huangdan (HD) and C.lanceoleosa. Moreover, expression divergence was observed in paralogous pairs and ∼1-2 oleosin genes in each clade have become activate. This study leads to a comprehensive understanding of evolution of oleosin family in Theaceae, and provides a rich resource to further address the functions of oleosin in tea and oil tea species.

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.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

A Tea Plant (Camellia sinensis) FLOWERING LOCUS C-like Gene, CsFLC1, Is Correlated to Bud Dormancy and Triggers Early Flowering in Arabidopsis

Flowering and bud dormancy are crucial stages in the life cycle of perennial angiosperms in temperate climates. MADS-box family genes are involved in many plant growth and development processes. Here, we identified three MADS-box genes in tea plant belonging to the FLOWERING LOCUS C (CsFLC) family. We monitored CsFLC1 transcription throughout the year and found that CsFLC1 was expressed at a higher level during the winter bud dormancy and flowering phases. To clarify the function of CsFLC1, we developed transgenic Arabidopsis thaliana plants heterologously expressing 35S::CsFLC1. These lines bolted and bloomed earlier than the WT (Col-0), and the seed germination rate was inversely proportional to the increased CsFLC1 expression level. The RNA-seq of 35S::CsFLC1 transgenic Arabidopsis showed that many genes responding to ageing, flower development and leaf senescence were affected, and phytohormone-related pathways were especially enriched. According to the results of hormone content detection and RNA transcript level analysis, CsFLC1 controls flowering time possibly by regulating SOC1, AGL42, SEP3 and AP3 and hormone signaling, accumulation and metabolism. This is the first time a study has identified FLC-like genes and characterized CsFLC1 in tea plant. Our results suggest that CsFLC1 might play dual roles in flowering and winter bud dormancy and provide new insight into the molecular mechanisms of FLC in tea plants as well as other plant species.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

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.

United States tea: A synopsis of ongoing tea research and solutions to United States tea production issues

Tea is a steeped beverage made from the leaves of Camellia sinensis. Globally, this healthy, caffeine-containing drink is one of the most widely consumed beverages. At least 50 countries produce tea and most of the production information and tea research is derived from international sources. Here, we discuss information related to tea production, genetics, and chemistry as well as production issues that affect or are likely to affect emerging tea production and research in the United States. With this review, we relay current knowledge on tea production, threats to tea production, and solutions to production problems to inform this emerging market in the United States.

Phylogenomics Resolves the Phylogeny of Theaceae by Using Low-Copy and Multi-Copy Nuclear Gene Makers and Uncovers a Fast Radiation Event Contributing to Tea Plants Diversity

Simple Summary

The Theaceae includes more than 300 species of great morphological diversity and has immense economic, cultural, and ornamental values. However, the evolutionary history of this family remains elusive. We integrated 91 genomes and transcriptome datasets of Theaceae and successfully resolved the phylogeny of Theaceae including relatives of cultivated tea plants from both extensive low-copy and multi-copy nuclear gene markers. Bayes-based molecular dating revealed that the ancestor of the tea family originated slightly earlier than the K-Pg boundary (Mass extinction events including the extinction of dinosaurs) with early diversification of three tribes associated with the Early Eocene Climatic Optimum. Further speciation analysis suggested a sole significant diversification shift rate in the common ancestor of Camellia associated with the Mid-Miocene Climatic Optimum. Collectively, polyploidy events, and key morphological innovation characters, such as pericarp with seed coat hardening, could possibly contribute to the Theaceae species diversity.

Abstract

Tea is one of the three most popular nonalcoholic beverages globally and has extremely high economic and cultural value. Currently, the classification, taxonomy, and evolutionary history of the tea family are largely elusive, including phylogeny, divergence, speciation, and diversity. For understanding the evolutionary history and dynamics of species diversity in Theaceae, a robust phylogenetic framework based on 1785 low-copy and 79,103 multi-copy nuclear genes from 91 tea plant genomes and transcriptome datasets had been reconstructed. Our results maximumly supported that the tribes Stewartieae and Gordonieae are successive sister groups to the tribe Theeae from both coalescent and super matrix ML tree analyses. Moreover, in the most evolved tribe, Theeae, the monophyletic genera Pyrenaria, Apterosperma, and Polyspora are the successive sister groups of Camellia. We also yield a well-resolved relationship of Camellia, which contains the vast majority of Theaceae species richness. Molecular dating suggests that Theaceae originated in the late L-Cretaceous, with subsequent early radiation under the Early Eocene Climatic Optimal (EECO) for the three tribes. A diversification rate shift was detected in the common ancestors of Camellia with subsequent acceleration in speciation rate under the climate optimum in the early Miocene. These results provide a phylogenetic framework and new insights into factors that likely have contributed to the survival of Theaceae, especially a successful radiation event of genus Camellia members to subtropic/tropic regions. These novel findings will facilitate the efficient conservation and utilization of germplasm resources for breeding cultivated tea and oil-tea. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Theaceae.

Genome-Wide Analysis of the TCP Gene Family and Their Expression Pattern Analysis in Tea Plant (Camellia sinensis)

TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors TEOSINTE BRANCHED1/CYCLOIDEA/PCF have been suggested to control the cell growth and proliferation in meristems and lateral organs. A total of 37 CsTCP genes were identified and divided into two classes, class I (PCF, group 1) and class II (CIN CYC/TB1, groups 2, and 3). The residues of TEOSINTE BRANCHED1/CYCLOIDEA/PCF of Camellia sinensis (Tea plant) (CsTCP) proteins between class I and class II were definitely different in the loop, helix I, and helix II regions; however, eighteen conserved tandem was found in bHLH. There are a large number of CsTCP homologous gene pairs in three groups. Additionally, most CsTCP proteins have obvious differences in motif composition. The results illuminated that CsTCP proteins in different groups are supposed to have complementary functions, whereas those in the same class seem to display function redundancies. There is no relationship between the number of CsTCP gene members and genome size, and the CsTCP gene family has only expanded since the divergence of monocots and eudicots. WGD/segmental duplication played a vital role in the expansion of the CsTCP gene family in tea plant, and the CsTCP gene family has expanded a lot. Most CsTCP genes of group 1 are more widely and non-specifically expressed, and the CsTCP genes of group 2 are mainly expressed in buds, flowers, and leaves. Most genes of group 1 and some genes of group 2 were up-/downregulated in varying degrees under different stress, CsTCP genes of group 3 basically do not respond to stress. TCP genes involved in abiotic stress response mostly belong to PCF group. Some CsTCP genes may have the same function as the homologous genes in Arabidopsis, but there is functional differentiation.

Genome-wide identification of PME gene family and expression of candidate genes associated with aluminum tolerance in tea plant (Camellia sinensis)

BACKGROUND

The major aluminum (Al) detoxication mechanism of tea plant (Camellia sinensis), as an Al hyperaccumulator plant, is the fixation of almost 70% of Al in the cell walls. Pectin is the primary constituent of cell walls, a degree of methylation of pectin polysaccharides regulated by the pectin methylesterase (PME) genes can greatly affect the Al binding capacity. The knowledge on PME gene family in tea plant is still poor.

RESULTS

We identified 66 (CsPME1-CsPME66) PME genes from C. sinensis genome. We studied their protein characterization, conserved motifs, gene structure, systematic evolution and gene expression under Al treatments, to establish a basis for in-depth research on the function of PMEs in tea plant. Gene structures analysis revealed that the majority of PME genes had 2-4 exons. Phylogenetic results pointed out that the PME genes from the same species displayed comparatively high sequence consistency and genetic similarity. Selective pressure investigation suggested that the Ka/Ks value for homologous genes of PME family was less than one. The expression of CsPMEs under three Al concentration treatments was tissue specific, eight PME genes in leaves and 15 in roots displayed a trend similar to of the Al contents and PME activities under Al concentration treatments, indicating that the degree of pectin de-esterification regulated by PME was crucial for Al tolerance of tea plant.

CONCLUSIONS

Sixty-six CsPME genes were identified for the first time in tea plant. The genome-wide identification, classification, evolutionary and transcription analyses of the PME gene family provided a new direction for further research on the function of PME gene in Al tolerance of tea plant.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Genome-Wide Investigation of the MiR166 Family Provides New Insights into Its Involvement in the Drought Stress Responses of Tea Plants (Camellia sinensis (L.) O. Kuntze)

MicroRNA166 (miR166) is a highly conserved plant miRNA that plays a crucial role in plant growth and the resistance to various abiotic stresses. However, the miR166s in tea (Camellia sinensis (L.) O. Kuntze) have not been comprehensively identified and analyzed. This study identified 30 mature miR166s and twelve pre-miR166s in tea plants. An evolutionary analysis revealed that csn-miR166s originating from the 3′ arm of their precursors were more conserved than the csn-miR166s derived from the 5′ arm of their precursors. The twelve pre-miR166s in tea were divided into two groups, with csn-MIR166 Scaffold364-2 separated from the other precursors. The Mfold-based predictions indicated that the twelve csn-MIR166s formed typical and stable structures comprising a stem-loop hairpin, with minimum free energy ranging from −110.90 to −71.80 kcal/mol. An analysis of the CsMIR166 promoters detected diverse cis-acting elements, including those related to light responses, biosynthesis and metabolism, abiotic stress defenses, and hormone responses. There was no one-to-one relationship between the csn-miR166s and their targets, but most csn-miR166s targeted HD-Zip III genes. Physiological characterization of tea plants under drought stress showed that leaf water content proportionally decreased with the aggravation of drought stress. In contrast, tea leaves’ malondialdehyde (MDA) content proportionally increased. Moreover, the cleavage site of the ATHB-15-like transcript was identified according to a modified 5′ RNA ligase-mediated rapid amplification of cDNA ends. The RT-qPCR data indicated that the transcription of nine csn-miR166s was negatively correlated with their target gene.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

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.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Genome-wide identification and characterization of functionally relevant microsatellite markers from transcription factor genes of Tea (Camellia sinensis (L.) O. Kuntze)

Tea, being one of the most popular beverages requires large set of molecular markers for genetic improvement of quality, yield and stress tolerance. Identification of functionally relevant microsatellite or simple sequence repeat (SSR) marker resources from regulatory “Transcription factor (TF) genes” can be potential targets to expedite molecular breeding efforts. In current study, 2776 transcripts encoding TFs harbouring 3687 SSR loci yielding 1843 flanking markers were identified from traits specific transcriptome resource of 20 popular tea cultivars. Of these, 689 functionally relevant SSR markers were successfully validated and assigned to 15 chromosomes (Chr) of CSS genome. Interestingly, 589 polymorphic markers including 403 core-set of TF-SSR markers amplified 2864 alleles in key TF families (bHLH, WRKY, MYB-related, C2H2, ERF, C3H, NAC, FAR1, MYB and G2-like). Their significant network interactions with key genes corresponding to aroma, quality and stress tolerance suggests their potential implications in traits dissection. Furthermore, single amino acid repeat reiteration in CDS revealed presence of favoured and hydrophobic amino acids. Successful deployment of markers for genetic diversity characterization of 135 popular tea cultivars and segregation in bi-parental population suggests their wider utility in high-throughput genotyping studies in tea.

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.

Combining novel technologies with interdisciplinary basic research to enhance horticultural crops

Horticultural crops mainly include fruits, vegetables, ornamental trees and flowers, and tea trees (Melaleuca alternifolia). They produce a variety of nutrients for the daily human diet in addition to the nutrition provided by staple crops, and some of them additionally possess ornamental and medicinal features. As such, horticultural crops make unique and important contributions to both food security and a colorful lifestyle. Under the current climate change scenario, the growing population and limited arable land means that agriculture, and especially horticulture, has been facing unprecedented challenges to meet the diverse demands of human daily life. Breeding horticultural crops with high quality and adaptability and establishing an effective system that combines cultivation, post-harvest handling, and sales becomes increasingly imperative for horticultural production. This review discusses characteristic and recent research highlights in horticultural crops, focusing on the breeding of quality traits and the mechanisms that underpin them. It additionally addresses challenges and potential solutions in horticultural production and post-harvest practices. Finally, we provide a prospective as to how emerging technologies can be implemented alongside interdisciplinary basic research to enhance our understanding and exploitation of horticultural crops.

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.

Feasible strategies for studying the involvement of DNA methylation and histone acetylation in the stress-induced formation of quality-related metabolites in tea (Camellia sinensis)

Tea plants are subjected to multiple stresses during growth, development, and postharvest processing, which affects levels of secondary metabolites in leaves and influences tea functional properties and quality. Most studies on secondary metabolism in tea have focused on gene, protein, and metabolite levels, whereas upstream regulatory mechanisms remain unclear. In this review, we exemplify DNA methylation and histone acetylation, summarize the important regulatory effects that epigenetic modifications have on plant secondary metabolism, and discuss feasible research strategies to elucidate the underlying specific epigenetic mechanisms of secondary metabolism regulation in tea. This information will help researchers investigate the epigenetic regulation of secondary metabolism in tea, providing key epigenetic data that can be used for future tea genetic breeding.

The Laccase Gene Family Mediate Multi-Perspective Trade-Offs during Tea Plant (Camellia sinensis) Development and Defense Processes

Laccase (LAC) plays important roles in different plant development and defense processes. In this study, we identified laccase genes (CsLACs) in Camellia sinensis cv ‘Longjing43′ cultivars, which were classified into six subclades. The expression patterns of CsLACs displayed significant spatiotemporal variations across different tissues and developmental stages. Most members in subclades II, IV and subclade I exhibited contrasting expression patterns during leaf development, consistent with a trade-off model for preferential expression in the early and late developmental stages. The extensive transcriptional changes of CsLACs under different phytohormone and herbivore treatment were observed and compared, with the expression of most genes in subclades I, II and III being downregulated but genes in subclades IV, V and VI being upregulated, suggesting a growth and defense trade-off model between these subclades. Taken together, our research reveal that CsLACs mediate multi-perspective trade-offs during tea plant development and defense processes and are involved in herbivore resistance in tea plants. More in-depth research of CsLACs upstream regulation and downstream targets mediating herbivore defense should be conducted in the future.

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.

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

Chemical composition and antifungal activity of Zanthoxylum riedelianum stem bark essential oil

The phytopathogenic fungi Sclerotinia sclerotiorum and Rhizopus stolonifer infect different crop plants. Essential oils have been used as an alternative to chemical control methods. Therefore, the objective of this study was to analyze the essential oil chemical composition of the stem bark of Zanthoxylum riedelianum (ZREO) and evaluate its antifungal potential. The ZREO obtained by hydrodistillation was analyzed using gas chromatography coupled to mass spectrometry and evaluated for in vitro antifungal activity. The major components identified in the chemical analysis were E-nerolidol (67.21%), α-selinene (14.94%), and β-selinene (7.41%). The antifungal potential of ZREO against S. sclerotiorum and R. stolonifer was detected at all concentrations evaluated in a dose-dependent manner. The best results were against R. stolonifer, with the concentration of 150 µLmL^-1 inhibiting more than 80% of mycelial growth. On the basis of the current knowledge, this study describes for the first time the chemical composition and their antifungal activity.

Biochemical characterization of specific Alanine Decarboxylase (AlaDC) and its ancestral enzyme Serine Decarboxylase (SDC) in tea plants (Camellia sinensis)

Background

Alanine decarboxylase (AlaDC), specifically present in tea plants, is crucial for theanine biosynthesis. Serine decarboxylase (SDC), found in many plants, is a protein most closely related to AlaDC. To investigate whether the new gene AlaDC originate from gene SDC and to determine the biochemical properties of the two proteins from Camellia sinensis, the sequences of CsAlaDC and CsSDC were analyzed and the two proteins were over-expressed, purified, and characterized.

Results

The results showed that exon-intron structures of AlaDC and SDC were quite similar and the protein sequences, encoded by the two genes, shared a high similarity of 85.1%, revealing that new gene AlaDC originated from SDC by gene duplication. CsAlaDC and CsSDC catalyzed the decarboxylation of alanine and serine, respectively. CsAlaDC and CsSDC exhibited the optimal activities at 45 °C (pH 8.0) and 40 °C (pH 7.0), respectively. CsAlaDC was stable under 30 °C (pH 7.0) and CsSDC was stable under 40 °C (pH 6.0–8.0). The activities of the two enzymes were greatly enhanced by the presence of pyridoxal-5′-phosphate. The specific activity of CsSDC (30,488 IU/mg) was 8.8-fold higher than that of CsAlaDC (3467 IU/mg).

Conclusions

Comparing to CsAlaDC, its ancestral enzyme CsSDC exhibited a higher specific activity and a better thermal and pH stability, indicating that CsSDC acquired the optimized function after a longer evolutionary period. The biochemical properties of CsAlaDC might offer reference for theanine industrial production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00674-x.

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.

How does tea (Camellia sinensis) produce specialized metabolites which determine its unique quality and function: a review

Tea (Camellia sinensis) is both a plant and a foodstuff. Many bioactive compounds, which are present in the final tea product and related to its quality or functional properties, are produced during the tea manufacturing process. However, the characteristic secondary metabolites, which give tea its unique qualities and are beneficial to human health, are produced mainly in the leaves during the process of plant growth. Therefore, it is important to understand how tea leaves produce these specialized metabolites. In this review, we first compare the common metabolites and specialized metabolites in tea, coffee, cocoa, and grape and discuss the occurrence of characteristic secondary metabolites in tea. Progress in research into the formation of these characteristic secondary metabolites in tea is summarized, including establishing a biological database and genetic transformation system, and the biosynthesis of characteristic secondary metabolites. Finally, speculation on future research into the characteristic secondary metabolites of tea is provided from the viewpoints of the origin, resources, cultivation, and processing of tea. This review provides an important reference for future research on the specialized metabolites of tea in terms of its characteristics.

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.

Metabolite signatures of diverse Camellia sinensis tea populations

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

Identifying Conserved Functional Gene Modules Underlying the Dynamic Regulation of Tea Plant Development and Secondary Metabolism

Tea plants adjust development and metabolism by integrating environmental and endogenous signals in complex but poorly defined gene networks. Here, we present an integrative analysis framework for the identification of conserved modules controlling important agronomic traits using a comprehensive collection of RNA-seq datasets in Camellia plants including 189 samples. In total, 212 secondary metabolism-, 182 stress response-, and 182 tissue development-related coexpressed modules were revealed. Functional modules (e.g., drought response, theobromine biosynthesis, and new shoot development-related modules) and potential regulators that were highly conserved across diverse genetic backgrounds and/or environmental conditions were then identified by cross-experiment comparisons and consensus clustering. Moreover, we investigate the preservation of gene networks between Camellia sinensis and other Camellia species. This revealed that the coexpression patterns of several recently evolved modules related to secondary metabolism and environmental adaptation were rewired and showed higher connectivity in tea plants. These conserved modules are excellent candidates for modeling the core mechanism of tea plant development and secondary metabolism and should serve as a great resource for hypothesis generation and tea quality improvement.