Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis).

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

R2R3-MYB transcription factor family in tea plant (Camellia sinensis): Genome-wide characterization, phylogeny, chromosome location, structure and expression patterns.

MYB transcription factors play essential roles in many biological processes and environmental stimuli. However, the functions of the MYB transcription factor family in tea plants have not been elucidated. Here, a total of 122 CsR2R3-MYB genes were identified from the chromosome level genome of tea plant (Camellia sinensis). The CsR2R3-MYB genes were phylogenetically classified into 25 groups. Results from the structure analysis of the gene, conserved motifs, and chromosomal distribution supported the relative conservation of the R2R3-MYB genes family in the tea plant. Synteny analysis indicated that 122, 34, and 112 CsR2R3-MYB genes were orthologous to Arabidopsis thaliana, Oryza sativa and C. sinensis var. 'huangdan' (HD), respectively. Tissue-specific expression showed that all CsR2R3-MYB genes had different expression patterns in the tea plant tissues, indicating that these genes may perform diverse functions. The expression patterns of representative R2R3-MYB genes and the regulatory network of the main anthocyanin components were analyzed, which suggested that CsMYB17 may played a key role in the regulation of cya-3-O-gal, del-3-O-gal, cya-3-O-glu and pel-3-O-glu. Results from the qRT-PCR validation of selected genes suggested that CsR2R3-MYB genes were induced in response to drought, cold, GA, and ABA treatments. Overall, this study provides comprehensive and systematic information for research on the function of R2R3-MYB genes in tea plants.

Lipidomics analysis unravels changes from flavor precursors in different processing treatments of purple-leaf tea.

BACKGROUND: Lipids are one of the most important bioactive compounds, affecting the character and quality of tea. However, the contribution of lipids to tea productions is still elusive. Here, we systematically identified the lipid profiles of green, oolong, and black teas in purple-leaf tea (Jinmingzao, JMZ) and green-leaf tea (Huangdan, HD), respectively. RESULTS: The lipids analysis showed regular accumulation in tea products with different manufacturing processes, among which the fatty acids, glycerolipids, glycerophospholipids, and sphingolipids contribute to the quality characteristics of tea products, including typical fatty acyl (FA), monogalactosyldiacylglycerol (MGDG), digalactosyldiacylglycerols (DGDG), and phosphatidylcholine (PC). Compared tea materials with products, levels of fatty acids were up-regulated, while glycerolipids and glycerophospholipids were down-regulated in tea products. FA 18:3, FA 16:0, MGDG 36:6, DGDG 36:6, PC 34:3, and PC 36:6 were the negative contributors to green tea flavor formation of purple-leaf tea. The pathway analysis of significant lipids in materials and products of purple-leaf tea were enriched linolenic acid metabolism pathway and glycerolipid metabolism. CONCLUSION: This study provides insights into the lipid metabolism profiles of different tea leaf colors, and found that fatty acids are essential precursors of black tea flavor formation. © 2021 Society of Chemical Industry.

Identification of Co-Expressed Genes Related to Theacrine Synthesis in Tea Flowers at Different Developmental Stages.

Jiaocheng kucha is the first reported tea germplasm resource which contains theacrine founded in Fujian Province. Currently, the anabolic mechanism of theacrine within tea leaves is clear, but there are few studies focused on its flowers. In order to further explore the mechanism of theacrine synthesis and related genes in flowers, current study applied Jiaocheng kucha flowers (JC) as test materials and Fuding Dabaicha flowers (FD) as control materials to make transcriptome sequencing, and determination of purine alkaloid content in three different developmental periods (flower bud stage, whitening stage and full opening stage). The results showed that the flower in all stages of JC contained theacrine. The theacrine in the flower bud stage was significantly higher than in the other stages. The differentially expressed genes (DEGs) at three different developmental stages were screened from the transcriptome data, and were in a total of 5642, 8640 and 8465. These DEGs related to the synthesis of theacrine were primarily annotated to the pathways of purine alkaloids. Among them, the number of DEGs in xanthine synthesis pathway was the largest and upregulated in JC, while it was the smallest in caffeine synthesis pathway and downregulated in JC. Further weighted gene co-expression network (WGCNA) indicated that ADSL (CsTGY03G0002327), ADSL (CsTGY09G0001824) and UAZ (CsTGY06G0002694) may be a hub gene for the regulation of theacrine metabolism in JC. Our results will contribute to the identification of candidate genes related to the synthesis of theacrine in tea flowers, and explore the molecular mechanism of theacrine synthesis in JC at different developmental stages.

Transcriptome and metabolite analyses provide insights into zigzag-shaped stem formation in tea plants (Camellia sinensis).

BACKGROUND: Shoot orientation is important for plant architecture formation, and zigzag-shaped shoots are a special trait found in many plants. Zigzag-shaped shoots have been selected and thoroughly studied in Arabidopsis; however, the regulatory mechanism underlying zigzag-shaped shoot development in other plants, especially woody plants, is largely unknown. RESULTS: In this study, tea plants with zigzag-shaped shoots, namely, Qiqu (QQ) and Lianyuanqiqu (LYQQ), were investigated and compared with the erect-shoot tea plant Meizhan (MZ) in an attempt to reveal the regulation of zigzag-shaped shoot formation. Tissue section observation showed that the cell arrangement and shape of zigzag-shaped stems were aberrant compared with those of normal shoots. Moreover, a total of 2175 differentially expressed genes (DEGs) were identified from the zigzag-shaped shoots of the tea plants QQ and LYQQ compared to the shoots of MZ using transcriptome sequencing, and the DEGs involved in the "Plant-pathogen interaction", "Phenylpropanoid biosynthesis", "Flavonoid biosynthesis" and "Linoleic acid metabolism" pathways were significantly enriched. Additionally, the DEGs associated with cell expansion, vesicular trafficking, phytohormones, and transcription factors were identified and analysed. Metabolomic analysis showed that 13 metabolites overlapped and were significantly changed in the shoots of QQ and LYQQ compared to MZ. CONCLUSIONS: Our results suggest that zigzag-shaped shoot formation might be associated with the gravitropism response and polar auxin transport in tea plants. This study provides a valuable foundation for further understanding the regulation of plant architecture formation and for the cultivation and application of horticultural plants in the future.

Comparative transcriptomic analysis of resistant and susceptible tea cultivars in response to Empoasca onukii (Matsuda) damage.

MAIN CONCLUSION: Transcriptomic studies in resistant and susceptible tea cultivars have been performed to reveal the different defense molecular mechanisms of tea after E. onukii feeding. The molecular mechanism by which tea plants respond to small green leafhopper Empoasca onukii (Matsuda) damage is unclear. Using the resistant tea plant cultivar Juyan (JY) and the susceptible tea plant cultivar Enbiao (EB) as materials, this study performed RNA-seq on tea leaf samples collected at three time points (6 h, 12 h, 24 h) during exposure of the plants to leafhopper to reveal the molecular mechanisms that are activated in susceptible and resistant tea plant cultivars in response to leafhopper damage. The numbers of DEGs in the susceptible tea cultivar during early (6 h) and late (24 h) stages of leafhopper induction were higher than those in the resistant cultivar at the same time points. The stress responses to leafhopper were most intense at 12 h in both tea cultivars. Pathway enrichment analysis showed that most up-regulated DEGs and their related metabolic pathways were similar in the two tea cultivars. However, during the early stage of leafhopper induction (6 h), jasmonic acid (JA)-related genes were significantly up-regulated in the resistant cultivar. The terpenoid biosynthetic pathway and the α-linolenic acid metabolic pathway were activated earlier in the resistant cultivar and remained activated until the late stage of leafhopper damage. Our results confirmed that after leafhopper damage, the resistant tea cultivar activated its defense responses earlier than the susceptible cultivar, and these defense responses were mainly related to terpenoid metabolism and JA biosynthetic pathway. The results provide important clues for further studies on resistance strategy of tea plants to pest.

Rapid and direct identification of the origin of white tea with proton transfer reaction time-of-flight mass spectrometry.

RATIONALE: White tea has become very popular in recent years, but there has been no scientific identification of white tea from different origins. For product authentication and valorization, every kind of white tea must be marked with an indication of its origin. METHODS: Volatile profiles of white tea leaf samples from their main origins in China (Fuding City, Zhenghe City and Jianyang City) were analyzed using proton transfer reaction time-of-flight mass spectrometry (PTR-TOFMS). Tentative identifications of the volatile organic compounds (VOCs) were obtained by PTR-TOFMS of the headspace. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were performed to evaluate the differences among the various origins. RESULTS: Teas from different origins were shown to have characteristic VOCs and profiles. Thus, white teas from different origins could be separated by characterizing the volatile emissions from the dry tea leaves. The ability of the two classification models to use the volatile fingerprints in origin discrimination was investigated. CONCLUSIONS: Two classification models (PCA and OPLS-DA) were applied to the PTR-TOFMS data obtained from the VOCs of various white teas. The classification models were shown to be useful in identifying the origin of white tea samples, providing a reference for white tea identification.

Comparison of Metabolome and Transcriptome of Flavonoid Biosynthesis Pathway in a Purple-Leaf Tea Germplasm Jinmingzao and a Green-Leaf Tea Germplasm Huangdan reveals Their Relationship with Genetic Mechanisms of Color Formation.

Purple-leaf tea is a phenotype with unique color because of its high anthocyanin content. The special flavor of purple-leaf tea is highly different from that of green-leaf tea, and its main ingredient is also of economic value. To probe the genetic mechanism of the phenotypic characteristics of tea leaf color, we conducted widely targeted metabolic and transcriptomic profiling. The metabolites in the flavonoid biosynthetic pathway of purple- and green-leaf tea were compared, and results showed that phenolic compounds, including phenolic acids, flavonoids, and tannins, accumulated in purple-leaf tea. The high expression of genes related to flavonoid biosynthesis (e.g., PAL and LAR) exhibits the specific expression of biosynthesis and the accumulation of these metabolites. Our result also shows that two CsUFGTs were positively related to the accumulation of anthocyanin. Moreover, genes encoding transcription factors that regulate flavonoids were identified by coexpression analysis. These results may help to identify the metabolic factors that influence leaf color differentiation and provide reference for future research on leaf color biology and the genetic improvement of tea.

Exploration of the Effects of Different Blue LED Light Intensities on Flavonoid and Lipid Metabolism in Tea Plants via Transcriptomics and Metabolomics.

Blue light extensively regulates multiple physiological processes and secondary metabolism of plants. Although blue light quantity (fluence rate) is important for plant life, few studies have focused on the effects of different blue light intensity on plant secondary metabolism regulation, including tea plants. Here, we performed transcriptomic and metabolomic analyses of young tea shoots (one bud and two leaves) under three levels of supplemental blue light, including low-intensity blue light (LBL, 50 µmol m-2 s-1), medium-intensity blue light (MBL, 100 µmol m-2 s-1), and high-intensity blue light (HBL, 200 µmol m-2 s-1). The total number of differentially expressed genes (DEGs) in LBL, MBL and HBL was 1, 7 and 1097, respectively, indicating that high-intensity blue light comprehensively affects the transcription of tea plants. These DEGs were primarily annotated to the pathways of photosynthesis, lipid metabolism and flavonoid synthesis. In addition, the most abundant transcription factor (TF) families in DEGs were bHLH and MYB, which have been shown to be widely involved in the regulation of plant flavonoids. The significantly changed metabolites that we detected contained 15 lipids and 6 flavonoid components. Further weighted gene co-expression network analysis (WGCNA) indicated that CsMYB (TEA001045) may be a hub gene for the regulation of lipid and flavonoid metabolism by blue light. Our results may help to establish a foundation for future research investigating the regulation of woody plants by blue light.

Comparison of Volatiles in Different Jasmine Tea Grade Samples Using Electronic Nose and Automatic Thermal Desorption-Gas Chromatography-Mass Spectrometry Followed by Multivariate Statistical Analysis.

Chinese jasmine tea is a type of flower-scented tea, which is produced by mixing green tea with the Jasminum sambac flower repeatedly. Both the total amount and composition of volatiles absorbed from the Jasminum sambac flower are mostly responsible for its sensory quality grade. This study aims to compare volatile organic compound (VOC) differences in authoritative jasmine tea grade samples. Automatic thermal desorption-gas-chromatography-mass spectrometry (ATD-GC-MS) and electronic nose (E-nose), followed by multivariate data analysis is conducted. Consequently, specific VOCs with a positive or negative correlation to the grades are screened out. Partial least squares-discriminant analysis (PLS-DA) and hierarchical cluster analysis (HCA) show a satisfactory discriminant effect on rank. It is intriguing to find that the E-nose is good at distinguishing the grade difference caused by VOC concentrations but is deficient in identifying essential aromas that attribute to the unique characteristics of excellent grade jasmine tea.

Identification, expression, and putative target gene analysis of nuclear factor-Y (NF-Y) transcription factors in tea plant (Camellia sinensis).

MAIN CONCLUSION: Genome-wide identification and characterization of nuclear factor-Y family in tea plants, and their expression profiles and putative targets provide the basis for further elucidation of their biological functions. The nuclear factor-Y (NF-Y) transcription factors (TFs) are crucial regulators of plant growth and physiology. However, the NF-Y TFs in tea plant (Camellia sinensis) have not yet been elucidated, and its biological functions, especially the putative target genes within the genome range, are still unclear. In this study, we identified 35 CsNF-Y encoding genes in the tea plant genome, including 10 CsNF-YAs, 15 CsNF-YBs and 10 CsNF-YCs. Their conserved domains and motifs, phylogeny, duplication event, gene structure, and promoter were subsequently analyzed. Tissue expression analysis revealed that CsNF-Ys exhibited three distinct expression patterns in eight tea tree tissues, among which CsNF-YAs were moderately expressed. Drought and abscisic acid (ABA) treatment indicated that CsNF-YAs may have a greater impact than other subunit members. Furthermore, through the genome-wide investigation of the presence of the CCAAT box, we found that CsNF-Ys may participate in the development of tea plants by regulating target genes of multiple physiological pathways, including photosynthesis, chlorophyll metabolism, fatty acid biosynthesis, and amino acid metabolism pathways. Our findings will contribute to the functional analysis of NF-Y genes in woody plants and the cultivation of high-quality tea plant cultivars.

Identification of CBF Transcription Factors in Tea Plants and a Survey of Potential CBF Target Genes under Low Temperature.

C-repeat binding factors (CBFs) are key signaling genes that can be rapidly induced by cold and bind to the C-repeat/dehydration-responsive motif (CRT/DRE) in the promoter region of the downstream cold-responsive (COR) genes, which play a vital role in the plant response to low temperature. However, the CBF family in tea plants has not yet been elucidated, and the possible target genes regulated by this family under low temperature are still unclear. In this study, we identified five CsCBF family genes in the tea plant genome and analyzed their phylogenetic tree, conserved domains and motifs, and cis-elements. These results indicate that CsCBF3 may be unique in the CsCBF family. This is further supported by our findings from the low-temperature treatment: all the CsCBF genes except CsCBF3 were significantly induced after treatment at 4 °C. The expression profiles of eight tea plant tissues showed that CsCBFs were mainly expressed in winter mature leaves, roots and fruits. Furthermore, 685 potential target genes were identified by transcriptome data and CRT/DRE element information. These target genes play a functional role under the low temperatures of winter through multiple pathways, including carbohydrate metabolism, lipid metabolism, cell wall modification, circadian rhythm, calcium signaling, transcriptional cascade, and hormone signaling pathways. Our findings will further the understanding of the stress regulatory network of CsCBFs in tea plants.

Transcriptome and Metabolite Profiling Reveal Novel Insights into Volatile Heterosis in the Tea Plant (Camellia Sinensis).

Tea aroma is a key indicator for evaluating tea quality. Although notable success in tea aroma improvement has been achieved with heterosis breeding technology, the molecular basis underlying heterosis remains largely unexplored. Thus, the present report studies the tea plant volatile heterosis using a high-throughput next-generation RNA-seq strategy and gas chromatography-mass spectrometry. Phenotypically, we found higher terpenoid volatile and green leaf volatile contents by gas chromatography-mass spectrometry in the F1 hybrids than in their parental lines. Volatile heterosis was obvious in both F1 hybrids. At the molecular level, the comparative transcriptomics analysis revealed that approximately 41% (9027 of 21,995) of the genes showed non-additive expression, whereas only 7.83% (1723 of 21,995) showed additive expression. Among the non-additive genes, 42.1% showed high parental dominance and 17.6% showed over-dominance. Among different expression genes with high parental dominance and over-dominance expression patterns, KEGG and GO analyses found that plant hormone signal transduction, tea plant physiological process related pathways and most pathways associated with tea tree volatiles were enriched. In addition, we identified multiple genes (CsDXS, CsAATC2, CsSPLA2, etc.) and transcription factors (CsMYB1, CsbHLH79, CsWRKY40, etc.) that played important roles in tea volatile heterosis. Based on transcriptome and metabolite profiling, we conclude that non-additive action plays a major role in tea volatile heterosis. Genes and transcription factors involved in tea volatiles showing over-dominance expression patterns can be considered candidate genes and provide novel clues for breeding high-volatile tea varieties.

Differential expression of gibberellin- and abscisic acid-related genes implies their roles in the bud activity-dormancy transition of tea plants.

KEY MESSAGE: Thirty genes involved in GA and ABA metabolism and signalling were identified, and the expression profiles indicated that they play crucial roles in the bud activity-dormancy transition in tea plants. Gibberellin (GA) and abscisic acid (ABA) are fundamental phytohormones that extensively regulate plant growth and development, especially bud dormancy and sprouting transition in perennial plants. However, there is little information on GA- and ABA-related genes and their expression profiles during the activity-dormancy transition in tea plants. In the present study, 30 genes involved in the metabolism and signalling pathways of GA and ABA were first identified, and their expression patterns in different tissues were assessed. Further evaluation of the expression patterns of selected genes in response to GA3 and ABA application showed that CsGA3ox, CsGA20ox, CsGA2ox, CsZEP and CsNCED transcripts were differentially expressed after exogenous treatment. The expression profiles of the studied genes during winter dormancy and spring sprouting were investigated, and somewhat diverse expression patterns were found for GA- and ABA-related genes. This diversity was associated with the bud activity-dormancy cycle of tea plants. These results indicate that the genes involved in the metabolism and signalling of GA and ABA are important for regulating the bud activity-dormancy transition in tea plants.

Identification and Expression Analyses of SBP-Box Genes Reveal Their Involvement in Abiotic Stress and Hormone Response in Tea Plant (Camellia sinensis).

The SQUAMOSA promoter binding protein (SBP)-box gene family is a plant-specific transcription factor family. This family plays a crucial role in plant growth and development. In this study, 20 SBP-box genes were identified in the tea plant genome and classified into six groups. The genes in each group shared similar exon-intron structures and motif positions. Expression pattern analyses in five different tissues demonstrated that expression in the buds and leaves was higher than that in other tissues. The cis-elements and expression patterns of the CsSBP genes suggested that the CsSBP genes play active roles in abiotic stress responses; these responses may depend on the abscisic acid (ABA), gibberellic acid (GA), and methyl jasmonate (MeJA) signaling pathways. Our work provides a comprehensive understanding of the CsSBP family and will aid in genetically improving tea plants.

Volatiles Emitted at Different Flowering Stages of Jasminum sambac and Expression of Genes Related to α-Farnesene Biosynthesis.

Fresh jasmine flowers have been used to make jasmine teas in China, but there has been no complete information about volatile organic compound emissions in relation to flower developmental stages and no science-based knowledge about which floral stage should be used for the infusion. This study monitored volatile organic compounds emitted from living flowers of Jasminum sambac (L.) Ait. 'Bifoliatum' at five developmental stages and also from excised flowers. Among the compounds identified, α-farnesene, linalool, and benzyl acetate were most abundant. Since α-farnesene is synthesized through the Mevalonate pathway, four genes encoding 3-hydroxy-3-methylglutaryl coenzyme A synthase, 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR), farnesyl pyrophosphate synthase, and terpene synthase were isolated. Their expression patterns in living flowers at the five stages and in excised flowers coincided with the emission patterns of α-farnesene. Application of lovastatin, a HMGR inhibitor, significantly reduced the expression of the genes and greatly decreased the emission of α-farnesene. The sweet scent was diminished from lovastatin-treated flowers as well. These results indicate that α-farnesene is an important compound emitted from jasmine flowers, and its emission patterns suggest that flowers at the opening stage or flower buds 8 h after excision should be used for the infusion of tea leaves.

Production regions discrimination of Huangguanyin oolong tea by using the content of chemical components and rare earth elements.

Oolong tea is one of the most popular tea beverages in China. Tea cultivars, processing technology and origin of production affect the quality and price of oolong teas. To investigate the differences in Huangguanyin oolong tea from different production regions, the chemical components, mineral elements and rare earth elements of Huangguanyin oolong tea produced in Yunxiao (YX) and Wuyishan (WY) were analyzed by using spectrophotometry methods, targeted metabolomics and inductive plasma coupled mass spectrometry (ICP-MS). The results of spectrophotometry methods revealed that there were significant differences in thearubigin, tea polyphenols and water extract between Huangguanyin oolong teas from different production regions. Targeted metabolomics identified a total of 31 chemical components in Huangguanyin oolong teas from the two production regions, of which 14 chemical components were significantly different and contributed to the regional differentiation of Huangguanyin oolong tea. Yunxiao Huangguanyin had relatively higher contents of (-)-Epigallocatechin-3-O-(3-O-methylgallate) (EGCG3″Me), ornithine (Orn) and histidine (His), while Wuyishan Huangguanyin had relatively higher contents of glutamic acid (Glu), γ-aminobutyric acid (GABA), ß-aminobutyric acid (ß-ABA) and other components. Moreover, ICP-MS identified a total of 15 mineral elements and 15 rare earth elements in Huangguanyin oolong tea from the two production regions, of which 15 elements were significantly different between YX and WY, and contributed to the regional differentiation of Huangguanyin oolong tea. K had a relatively higher content in Yunxiao Huangguanyin, while rare earth elements had relatively higher contents in Wuyishan Huangguanyin. The classification results by the production region showed that the discrimination rate of the support vector machine (SVM) model based on the 14 different chemical components reached 88.89%, while the SVM model based on the 15 elements reached 100%. Therefore, we used targeted metabolomics and ICP-MS techniques to screen and explore the chemical components, mineral elements and rare earth elements differences among two production regions, which indicated the feasibility of Huangguanyin oolong tea classification by production regions in the study. The results will provide some reference for the distinction between the two production regions of Huangguanyin oolong tea.

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.

Changes in Non-Volatile and Volatile Metabolites Associated with Heterosis in Tea Plants (Camellia sinensis).

Heterosis or hybrid vigor is extensively used in plant breeding. However, the contribution of metabolites to heterosis is still elusive. Here, we systematically identified the non-volatile and volatile metabolites of two hybrids and their parents in Camellia sinensis. The metabolomics analysis showed prevalent non-additive accumulation in hybrids, among which the non-additive nucleotides, alkaloids, organic acids, and tannins contribute to the positive heterosis of hybrids, including typical inosine, guanosine, adenosine, caffeine, succinic acid, adipic acid, xylonic acid, and gallic acid. The catechins and free amino acids in hybrids showed negative heterosis compared to its maternal cultivar TGY. Furthermore, the significant accumulation of non-additive terpenes combined with the mild heterosis of other types of volatiles contributes to the aroma of tea plant hybrids. The genetics of volatiles from different parents affect the aroma of hybrids processed into oolong tea. The comprehensive heterosis of these non-additive metabolites may play an important role in the formation of desirable breeding traits for hybrids. Our results provide insights into the utilization of heterosis breeding and the regulation of heterosis metabolites in tea plants.