Tea plant genomics: achievements, challenges and perspectives

A genus-specific R2R3 MYB transcription factor, CsMYB34, regulates galloylated catechin biosynthesis in Camellia sinensis

Galloylated catechins are the dominant polyphenols in Camellia sinensis (L.) O. Kuntze. The mechanisms responsible for accumulation of these specialized metabolites in tea plants remains unclear. This paper presents an extended member of subgroup 5 of transcription factors R2R3-MYB, CsMYB34, as a critical gene specifically regulating galloylated catechin biosynthesis. CsMYB34 has a TT2-type motif [VIRTKATRCSKVFIP]. Its transcription levels were positively correlated with galloylated catechin content in 19 tea varieties, with correlation coefficients ≥0.79. Suppression of CsMYB34 expression caused a significant decrease in galloylated catechin content, as well as reduced expression levels of the key galloylated catechin biosynthesis gene CsSCPL4. Yeast one-hybrid (Y1H), electrophoretic mobile shift assay (EMSA) and dual-luciferase reporter system (DLR) showed that CsMYB34 interacts directly with the promoter region of CsSCPL4, thereby upregulating its transcription. This research indicates that the CsMYB34 transcription factor selectively modulates the biosynthetic pathway of galloylated catechins, thereby offering a plausible rationale for the observed elevated levels of these compounds in tea leaves.

Distribution characteristics, source analysis and ecological risk assessment of PAHs in tea garden soil in China

In this study, we collected 177 soil samples from major tea-producing areas in China, systematically investigated the spatial distribution characteristics of the polycyclic aromatic hydrocarbons (PAHs) in the soil of these tea plantations and discussed the environmental factors influencing of the PAHs in tea garden soil. The feature ratio method and source analysis methods were used to determine the PAHs source in tea garden soil, and the potential risk of PAHs in tea garden soil was also evaluated. The results showed that the concentrations of the 16 PAHs in 177 samples ranged from 6.21 to 4068.91 ng g-1, with an average of 257.00 ng g-1. The majority of PHAs in tea garden soils predominantly contained a 5-6 ring pattern, and the highest content was indeno (1,2,3-cd) pyrene (InP, 23%) and benzo (b) fluoranthrene (BbF, 16%). In addition, 10.16% of the PAHs in tea plantation soils contained a 2-3-ring pattern, with naphthalene (NAP) having the highest content. PAH source in Chinese tea garden soil was predominantly mixed combustion, such as incomplete biomass combustion, petroleum combustion, coal combustion and wood combustion. The PAHs distribution was mainly affected by the industrial structure, geographical location of tea plantation, climatic conditions, soil properties and other factors in different regions. According to the Dutch Maliszewska-Kordybach grading standard, 79% of the soil samples from Chinese tea plantations were classified as unpolluted, 13% as mildly polluted, and 2% and 6% as moderately and severely polluted, respectively. Although the PAH pollution in tea plantations was generally low, BaP and InP pose significant ecological risk in some areas. Therefore, strategies such as effective guidelines and environmentally friendly technologies, must be developed to reduce the risk of PAH pollution in tea plantation soils.

Genetic diversity analysis and core collection construction of tea plant from the Yunnan Province of China using ddRAD sequencing

Tea plants are economically important woody plants that originated in southwestern China. The Yunnan Province in China is regarded as the central origin of tea plants owing to the abundance and diversity of the tea germplasm present in this region. However, there is a lack of knowledge regarding the genetic diversity and evolutionary relationships among tea plants in this region. Here, a total of 266,397 high-quality single nucleotide polymorphisms were obtained using double digest restriction-site associated DNA sequencing from 468 tea accessions collected from Yunnan. The phylogenetic relationship and population structure of the 468 tea accessions revealed remarkable inter- and intraspecific introgression across different sect. Thea species. We found that Camellia taliensis (W. W. Smith) Melchior is the main genetic donor that greatly contributed to the domestication of C. sinensis (L.) O. Kuntze, as evidenced by the frequent generation of genetic intermediates derived from the hybridization between C. sinensis and C. taliensis in Lincang, Xishuangbanna, and Pu'er. Notable genetic differences were observed across the populations from distinct regions of Yunnan. Additionally, the Lancang River may have attenuated the genetic interflow between Lincang and Pu'er. Finally, a core collection of 50 tea accessions was constructed based on the genetic diversity of the 468 tea accessions. Our results provide novel insights into the evolutionary and domestication history of tea plants in Yunnan Province and can facilitate the development of strategies for the conservation, breeding, and utilization of the core collection.

New perspectives of post-GWAS analyses: From markers to causal genes for more precise crop breeding

Crop breeding advancement is hindered by the imperfection of methods to reveal genes underlying key traits. Genome-wide Association Study (GWAS) is one such method, identifying genomic regions linked to phenotypes. Post-GWAS analyses predict candidate genes and assist in causative mutation (CM) recognition. Here, we assess post-GWAS approaches, address limitations in omics data integration and stress the importance of evaluating associated variants within a broader context of publicly available datasets. Recent advances in bioinformatics tools and genomic strategies for CM identification and allelic variation exploration are reviewed. We discuss the role of markers and marker panel development for more precise breeding. Finally, we highlight the perspectives and challenges of GWAS-based CM prediction for complex quantitative traits.

CsWRKY29, a key transcription factor in tea plant for freezing tolerance, ABA sensitivity, and sugar metabolism

Tea plants (Camellia sinensis L.) are prone to spring frosts, leading to substantial economic damage. WRKY transcription factors are key in plant abiotic stress responses, yet the role of CsWRKY29 in freezing tolerance is unclear. In this study, quantitative real-time PCR (qRT-PCR) and transient green fluorescent protein assay revealed that CsWRKY29 localizes to the nucleus and its expression is induced by cold and abscisic acid (ABA). CsWRKY29 overexpression in Arabidopsis enhanced freezing tolerance, reduced electrolyte leakage, increased soluble sugars, and boosted superoxide dismutase activity, with upregulated COR genes. These lines also showed heightened ABA and glucose sensitivity. Cold treatment of CsWRKY29-overexpressing lines upregulated AtABI5, AtHXK1, and AtSUS4 compared to wild type, and yeast one-hybrid assays confirmed CsWRKY29 binding to the W-box in the CsABI5 promoter. Furthermore, the application of virus-induced gene silencing (VIGS) technology to reduce CsWRKY29 expression in tea plants revealed a significant decrease in the transcript levels of CsCBFs, CsABI5, CsHXK1, and CsSUS4 in the silenced plants. In summary, our findings indicate that CsWRKY29 may serve as a critical transcription factor that contributes to freezing tolerance, ABA responsiveness, and sugar metabolism within tea plants.

Deciphering the anthocyanin metabolism gene network in tea plant (Camellia sinensis) through structural equation modeling

BACKGROUND

Tea is an important cash crop that significantly contributes to rural development, poverty reduction and food security in many developing countries. It provides livelihoods for millions of smallholder producers and aids their economic stability. Anthocyanins in tea leaves provides excellent commercial quality and germplasm exploration potential. These compounds give tea leaves vibrant colors and increase health benefits. The current understanding of the synergistic regulation mechanisms responsible for color changes in purple tea, attributed to anthocyanin degradation, remains unclear.

RESULTS

In this study, we have identified 30 gene families within the genome that are associated to with anthocyanin metabolism from tea. These gene families play distinct roles in the biosynthesis of anthocyanin including the formation of the core, structure, modification of the molecular framework, facilitation of transport process, regulation of gene expression, breakdown pathways, sugar transportation and iron ion respectively. Subsequently, we investigated the synergistic mechanisms of anthocyanin metabolism related gene families within tea leaves using structural equation modeling. The results showed that sugar transport positively affects anthocyanin transportation, and promotes anthocyanin degradation during leaf pigmentation, whereas, it inhibits anthocyanin degradation during the fading of leaf color. Further, Iron ions facilitate the degradation of anthocyanins during their deposition and conversely, impede this degradation process during digestion. These finding suggests that tea plants may regulate the synthesis and degradation of anthocyanins through sugar transport and iron ions ensure healthy levels and vibrant colors.

CONCLUSIONS

Our study contributes valuable information into the dynamic equilibrium anthocyanin mechanism and sheds light on complex regulatory mechanisms that govern the synthesis, transport and degradation of these pigments. These insights could be further used to develop strategies for enhancing anthocyanins content in unique tea germplasm to aid tea industry in producing new tea products with increased health benefits and aesthetic appeals.

In Vitro vs. In Vivo Transcriptomic Approach Revealed Core Pathways of Nitrogen Deficiency Response in Tea Plant (Camellia sinensis (L.) Kuntze)

For the first time, we used an in vitro vs. in vivo experimental design to reveal core pathways under nitrogen deficiency (ND) in an evergreen tree crop. These pathways were related to lignin biosynthesis, cell redox homeostasis, the defense response to fungus, the response to Karrikin, amino acid transmembrane transport, the extracellular region, the cellular protein catabolic process, and aspartic-type endopeptidase activity. In addition, the mitogen-activated protein kinase pathway and ATP synthase (ATP)-binding cassette transporters were significantly upregulated under nitrogen deficiency in vitro and in vivo. Most of the MAPK downstream genes were related to calcium signaling (818 genes) rather than hormone signaling (157 genes). Moreover, the hormone signaling pathway predominantly contained auxin- and abscisic acid-related genes, indicating the crucial role of these hormones in ND response. Overall, 45 transcription factors were upregulated in both experiments, 5 WRKYs, 3 NACs, 2 MYBs, 2 ERFs, HD-Zip, RLP12, bHLH25, RADIALIS-like, and others, suggesting their ND regulation is independent from the presence of a root system. Gene network reconstruction displayed that these transcription factors participate in response to fungus/chitin, suggesting that nitrogen response and pathogen response have common regulation. The upregulation of lignin biosynthesis genes, cytochrome genes, and strigalactone response genes was much more pronounced under in vitro ND as compared to in vivo ND. Several cell wall-related genes were closely associated with cytochromes, indicating their important role in flavanols biosynthesis in tea plant. These results clarify the signaling mechanisms and regulation of the response to nitrogen deficiency in evergreen tree crops.

Metabolite profiling and transcriptome analyses reveal defense regulatory network against pink tea mite invasion in tea plant

BACKGROUND

The tea plant Camellia sinensis (L.) O. Kuntze is a perennial crop, invaded by diversity of insect pest species, and pink tea mite is one of the most devastating pests for sustainable tea production. However, molecular mechanism of defense responses against pink tea mites in tea is still unknown. In this study, metabolomics and transcriptome profiles of susceptible and resistant tea varieties were compared before and after pink tea mite infestation.

RESULTS

Metabolomics analysis revealed that abundance levels of polyphenol-related compounds changed significantly before and after infestation. At the transcript level, nearly 8 GB of clean reads were obtained from each sequenced library, and a comparison of infested plants of resistant and susceptible tea varieties revealed 9402 genes with significant differential expression. An array of genes enriched in plant pathogen interaction and biosynthetic pathways of phenylpropanoids showed significant differential regulation in response to pink tea mite invasion. In particular, the functional network linkage of disease resistant proteins, phenylalanine ammonia lyase, flavanone -3-hydroxylase, hydroxycinnamoyl-CoA shikimate transferase, brassinosteroid-6-oxidase 1, and gibberellin 2 beta-dioxygenase induced dynamic defense signals to suppress prolonged pink tea mite attacks. Further integrated analyses identified a complex network of transcripts and metabolites interlinked with precursors of various flavonoids that are likely modulate resistance against to pink tea mite.

CONCLUSIONS

Our results characterized the profiles of insect induced metabolic and transcript reprogramming and identified a defense regulatory network that can potentially be used to fend off pink tea mites damage.

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

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 and Characterization of Top-Quality Japanese Green Tea Cultivar 'Seimei'

Japanese green tea, an essential beverage in Japanese culture, is characterized by the initial steaming of freshly harvested leaves during production. This process efficiently inactivates endogenous enzymes such as polyphenol oxidases, resulting in the production of sencha, gyokuro and matcha that preserves the vibrant green color of young leaves. Although genome sequences of several tea cultivars and germplasms have been published, no reference genome sequences are available for Japanese green tea cultivars. Here, we constructed a reference genome sequence of the cultivar 'Seimei', which is used to produce high-quality Japanese green tea. Using the PacBio HiFi and Hi-C technologies for chromosome-scale genome assembly, we obtained 15 chromosome sequences with a total genome size of 3.1 Gb and an N50 of 214.9 Mb. By analyzing the genomic diversity of 23 Japanese tea cultivars and lines, including the leading green tea cultivars 'Yabukita' and 'Saemidori', it was revealed that several candidate genes could be related to the characteristics of Japanese green tea. The reference genome of 'Seimei' and information on genomic diversity of Japanese green tea cultivars should provide crucial information for effective breeding of such cultivars in the future.

Plant-derived citronellol can significantly disrupt cell wall integrity maintenance of Colletotrichum camelliae

Anthracnose, a fungal disease, commonly infects tea plants and severely impacts the yield and quality of tea. One method for controlling anthracnose is the application of citronellol, a plant extract that exhibits broad-spectrum antimicrobial activity. Herein, the physiological and biochemical mechanism by which citronellol controls anthracnose caused by Colletotrichum camelliae was investigated. Citronellol exhibited excellent antifungal activity based on direct and indirect mycelial growth inhibition assays, with EC50 values of 76.88 mg/L and 29.79 μL/L air, respectively. Citronellol also exhibited good control effects on C. camelliae in semi-isolated leaf experiments. Optical and scanning electron microscopy revealed that citronellol caused C. camelliae mycelia to thin, fracture, fold and deform. Transmission electron microscopy revealed that the mycelial cell walls collapsed inward and separated, and the organelles became blurred after treatment with citronellol. The sensitivity of C. camelliae to calcofluor white staining was significantly enhanced by citronellol, while PI staining showed minimal fluorescence, and the relative conductivity of mycelia were not significantly different. Under citronellol treatment, the expression levels of β-1,3-glucanase, chitin synthase, and chitin deacetylase-related genes were significantly decreased, while the expression levels of chitinase genes were increased, leading to lower chitinase activity and increased β-1,3-glucanase activity. Therefore, citronellol disrupted the cell wall integrity of C. camelliae and inhibited normal mycelial growth.

LC-MS and GC-MS based metabolomics analysis revealed the impact of tea trichomes on the chemical and flavor characteristics of white tea

To explore the influence of tea trichomes on the quality of white tea, liquid chromatography quadrupole time-of-flight mass spectrometry (LC-Q-TOF-MS), and headspace solid phase microextraction gas chromatography-mass spectrometry (HS-SPME-GC-MS) were used to identify non-volatile and volatile compounds white tea without trichomes (WTwt) and pure trichomes (PT). It was found that the bitter and astringent compounds, caffeine (CAF), epigallocatechin gallate (EGCG), epicatechin gallate (ECG) and flavonol glycosides, were mainly enriched in the WTwt, with 16.3-fold, 47.1-fold and 28.7-fold decrease in CAF and EGCG and ECG, respectively, and the content of these compounds in PT were lower than the taste thresholds. In PT, kaempferol-3-O-(p-coumaroyl)-glucoside and kaempferol-3-O-(di-p-coumaroyl)-glucoside were non-volatile marker compounds, and decanal was significant aroma contributor with rOAV = 250.86. Moreover, the compounds in trichomes mainly contributed to the fruity and floral aroma of white tea, among which benzyl alcohol, (E)-geranylacetone, decanal, dodecanal and 6-methyl-5-hepten-2-one were the crucial aroma components, which were 2.1, 1.7, 1.8, 1.4 and 2.2 times as much as the WTwt in the PT, respectively. In conclusion, trichomes can improve the quality of white tea by reducing the bitterness and astringency, increasing the umami, as well as enhancing the fruity and floral aromas.

Efficient vegetation indices for phenotyping of abiotic stress tolerance in tea plant (Camellia sinensis (L.) Kuntze)

Early non-destructive detection of stress effect is crucial for efficient breeding strategies and germplasm characterization. Recently developed hyperspectral technologies allow to perform fast real-time phenotyping through reflectance-based vegetation indices. However, efficiency of these vegetation indices has to be validated for each crop in different environment. The aim of this study was to reveal efficient vegetation indices for phenotyping of abiotic stress (cold, freezing and nitrogen deficiency) response in tea plant. Among 31 studied VIs, few indices were efficient to distinguish tolerant and susceptible tea plants under abiotic stress: ZMI (Zarco-Tejada & Miller Index), VREI1,2,3 (Vogelmann Red Edge Indices), RENDVI (Red Edge Normalized Difference Vegetation Index), CTR1 and CTR2 (Carter Indices). Most of these indices are calculated based on reflectance in near-infrared area at 705-760 nm, indicating this range as promising for tea germplasm characterization under abiotic stresses. Tolerant tea plants showed the following values under freezing: ZMI ≥1.90, VREI1 ≥ 1.40, RENDVI ≥0.38, Ctr1 ≤ 1.74. The leaf N-content was positively correlated (Pearson's) with the following indices ZMI, VREI1, RENDVI, while negatively correlated with CTR, and VREI2,3. These results will be useful for tea germplasm management, genomics and breeding research aimed at abiotic stress tolerance of tea plant.

The Ability of Different Tea Tree Germplasm Resources in South China to Aggregate Rhizosphere Soil Characteristic Fungi Affects Tea Quality

It is generally recognized that the quality differences in plant germplasm resources are genetically determined, and that only a good "pedigree" can have good quality. Ecological memory of plants and rhizosphere soil fungi provides a new perspective to understand this phenomenon. Here, we selected 45 tea tree germplasm resources and analyzed the rhizosphere soil fungi, nutrient content and tea quality. We found that the ecological memory of tea trees for soil fungi led to the recruitment and aggregation of dominant fungal populations that were similar across tea tree varieties, differing only in the number of fungi. We performed continuous simulation and validation to identify four characteristic fungal genera that determined the quality differences. Further analysis showed that the greater the recruitment and aggregation of Saitozyma and Archaeorhizomyces by tea trees, the greater the rejection of Chaetomium and Trechispora, the higher the available nutrient content in the soil and the better the tea quality. In summary, our study presents a new perspective, showing that ecological memory between tea trees and rhizosphere soil fungi leads to differences in plants' ability to recruit and aggregate characteristic fungi, which is one of the most important determinants of tea quality. The artificial inoculation of rhizosphere fungi may reconstruct the ecological memory of tea trees and substantially improve their quality.

Recruitment and Aggregation Capacity of Tea Trees to Rhizosphere Soil Characteristic Bacteria Affects the Quality of Tea Leaves

There are obvious differences in quality between different varieties of the same plant, and it is not clear whether they can be effectively distinguished from each other from a bacterial point of view. In this study, 44 tea tree varieties (Camellia sinensis) were used to analyze the rhizosphere soil bacterial community using high-throughput sequencing technology, and five types of machine deep learning were used for modeling to obtain characteristic microorganisms that can effectively differentiate different varieties, and validation was performed. The relationship between characteristic microorganisms, soil nutrient transformation, and tea quality formation was further analyzed. It was found that 44 tea tree varieties were classified into two groups (group A and group B) and the characteristic bacteria that distinguished them came from 23 genera. Secondly, the content of rhizosphere soil available nutrients (available nitrogen, available phosphorus, and available potassium) and tea quality indexes (tea polyphenols, theanine, and caffeine) was significantly higher in group A than in group B. The classification result based on both was consistent with the above bacteria. This study provides a new insight and research methodology into the main reasons for the formation of quality differences among different varieties of the same plant.

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

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

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

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

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

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

TPIA2: an updated tea plant information archive for Camellia genomics

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

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

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

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

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

Introduction

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

Methods

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

Results

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

Discussion

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Background

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

Results

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

Discussion

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

Supplementary Information

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

Temporal metabolite responsiveness of microbiota in the tea plant phyllosphere promotes continuous suppression of fungal pathogens

INTRODUCTION

A broad spectrum of rhizosphere bacteria and fungi were shown to play a central role for health, fitness and productivity of their host plants. However, implications of host metabolism on microbiota assembly in the phyllosphere and potential consequences for holobiont functioning were sparsely addressed. Previous observations indicated that tea plants might reduce disease occurrence in various forests located in their proximity; the underlying mechanisms and potential implications of the phyllosphere microbiota remained elusive.

OBJECTIVES

This study aimed atdeciphering microbiome assembly in the tea plant phyllosphere throughout shoot development as well as elucidating potential implications of host metabolites in this process. The main focus was to explore hidden interconnections between the homeostasis of the phyllosphere microbiome and resistance to fungal pathogens.

METHODS

Profiling of host metabolites and microbiome analyses based on high-throughput sequencing were integrated to identify drivers of microbiome assembly throughout shoot development in the phyllosphere of tea plants. This was complemented by tracking of beneficial microorganisms in all compartments of the plant. Synthetic assemblages (SynAss), bioassays and field surveys were implemented to verify functioning of the phyllosphere microbiota.

RESULTS

Theophylline and epigallocatechin gallate, two prevalent metabolites at the early and late shoot development stage respectively, were identified as the main drivers of microbial community assembly. Flavobacterium and Myriangium were distinct microbial responders at the early stage, while Parabacteroides and Mortierella were more enriched at the late stage. Reconstructed, stage-specific SynAss suppressed various tree phytopathogens by 13.0%-69.3% in vitro and reduced disease incidence by 8.24%-41.3% in vivo.

CONCLUSION

The findings indicate that a functional phyllosphere microbiota was assembled along with development-specific metabolites in tea plants, which continuously suppressed prevalent fungal pathogens. The insights gained into the temporally resolved metabolite response of the tea plant microbiota could provide novel solutions for disease management.