In field conditions, especially in arid and semi-arid areas, tea plants are often simultaneously exposed to various abiotic stresses such as cold and drought, which have profound effects on leaf senescence process and tea quality. However, most studies of gene expression in stress responses focus on a single inciting agent, and the confounding effect of multiple stresses on crop quality and leaf senescence remain unearthed. Here, global transcriptome profiles of tea leaves under separately cold and drought stress were compared with their combination using RNA-Seq technology. This revealed that tea plants shared a large overlap in unigenes displayed “similar” (26%) expression pattern and avoid antagonistic responses (lowest level of “prioritized” mode: 0%) to exhibit very congruent responses to co-occurring cold and drought stress; 31.5% differential expressed genes and 38% of the transcriptome changes in response to combined stresses were unpredictable from cold or drought single-case studies. We also identified 319 candidate genes for enhancing plant resistance to combined stress. We then investigated the combined effect of cold and drought on tea quality and leaf senescence. Our results showed that drought-induced leaf senescence were severely delayed by (i) modulation of a number of senescence-associated genes and cold responsive genes, (ii) enhancement of antioxidant capacity, (iii) attenuation of lipid degradation, (iv) maintenance of cell wall and photosynthetic system, (v) alteration of senescence-induced sugar effect/sensitivity, as well as (vi) regulation of secondary metabolism pathways that significantly influence the quality of tea during combined stress. Therefore, care should be taken when utilizing a set of stresses to try and maximize leaf longevity and tea quality.

Tea [Camellia sinensis (L) O. Kuntze, Theaceae] is one of the most popular non-alcoholic beverages worldwide. Cold stress is one of the most severe abiotic stresses that limit tea plants’ growth, survival and geographical distribution. However, the genetic regulatory network and signaling pathways involved in cold stress responses in tea plants remain unearthed. Using RNA-Seq, DGE and sRNA-Seq technologies, we performed an integrative analysis of miRNA and mRNA expression profiling and their regulatory network of tea plants under chilling (4℃) and freezing (-5℃) stress. Differentially expressed (DE) miRNA and mRNA profiles were obtained based on fold change analysis, miRNAs and target mRNAs were found to show both coherent and incoherent relationships in the regulatory network. Furthermore, we compared several key pathways (e.g., ‘Photosynthesis’), GO terms (e.g., ‘response to karrikin’) and transcriptional factors (TFs, e.g., DREB1b/CBF1) which were identified as involved in the early chilling and/or freezing response of tea plants. Intriguingly, we found that karrikins, a new group of plant growth regulators, and β-primeverosidase (BPR), a key enzyme functionally relevant with the formation of tea aroma might play an important role in both early chilling and freezing response of tea plants. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-Seq and sRNA-Seq analysis. This is the first study to simultaneously profile the expression patterns of both miRNAs and mRNAs on a genome-wide scale to elucidate the molecular mechanisms of early responses of tea plants to cold stress. In addition to gaining a deeper insight into the cold resistant characteristics of tea plants, we provide a good case study to analyse mRNA/miRNA expression and profiling of non-model plant species using next-generation sequencing technology.

The geographic origin of agri-food products contributes greatly to their quality and market value. Here, we developed a robust method combining metabolomics and machine learning (ML) to authenticate the geographic origin of Wuyi rock tea, a premium oolong tea. The volatiles of 333 tea samples (174 from the core region and 159 from the non-core region) were profiled using gas chromatography time-of-flight mass spectrometry and a series of ML algorithms were tested. Wuyi rock tea from the two regions featured distinct aroma profiles. Multilayer Perceptron achieved the best performance with an average accuracy of 92.7% on the training data using 176 volatile features. The model was benchmarked with two independent test sets, showing over 90% accuracy. Gradient Boosting algorithm yielded the best accuracy (89.6%) when using only 30 volatile features. The proposed methodology holds great promise for its broader applications in identifying the geographic origins of other valuable agri-food products.

Tea plant (Camellia sinensis (L) O. Kuntze) respond to herbivore attack through large changes in defense related metabolism and gene expression. Ectropis oblique (Prout) is one of the most devastating insects that feed on tea leaves and tender buds, which can cause severe production loss and deteriorate the quality of tea. To elucidate the biochemicals and molecular mechanism of defense against tea geometrid (TG), transcriptome and metabolome of TG interaction with susceptible (SG) and resistance (RG) tea genotypes were analyzed by using UPLC-Q-TOF-MS, GC-MS, and RNA-seq technologies. This revealed that jasmonic acid was highly induced in RG, following a plethora of secondary metabolites involved in defense against TG could be induced by jasmonic acid signaling pathway. However, the constitutively present of salicylic acid in SG might be a suppressor of jasmonate signaling and thus misdirect tea plants against TG. Furthermore, flavonoids and terpenoids biosynthesis pathways were highly activated in RG to constitute the chemical barrier on TG feeding behavior. In contrast, fructose and theanine, which can act as feeding stimulants were observed to highly accumulate in SG. Being present in the major hub, 39 transcription factors or protein kinases among putative candidates were identified as master regulators from protein-protein interaction network analysis. Together, the current study provides a comprehensive gene expression and metabolite profiles, which can shed new insights into the molecular mechanism of tea defense against TG. The candidate genes and specific metabolites identified in the present study can serve as a valuable resource for unraveling the possible defense mechanism of plants against various biotic stresses.