Dynamic DNA Methylation Regulates Season-Dependent Secondary Metabolism in the New Shoots of Tea Plants.

Plant secondary metabolites are critical quality-conferring compositions of plant-derived beverages, medicines, and industrial materials. The accumulations of secondary metabolites are highly variable among seasons; however, the underlying regulatory mechanism remains unclear, especially in epigenetic regulation. Here, we used tea plants to explore an important epigenetic mark DNA methylation (5mC)-mediated regulation of plant secondary metabolism in different seasons. Multiple omics analyses were performed on spring and summer new shoots. The results showed that flavonoids and theanine metabolism dominated in the metabolic response to seasons in the new shoots. In summer new shoots, the genes encoding DNA methyltransferases and demethylases were up-regulated, and the global CG and CHG methylation reduced and CHH methylation increased. 5mC methylation in promoter and gene body regions influenced the seasonal response of gene expression; the amplitude of 5mC methylation was highly correlated with that of gene transcriptions. These differentially methylated genes included those encoding enzymes and transcription factors which play important roles in flavonoid and theanine metabolic pathways. The regulatory role of 5mC methylation was further verified by applying a DNA methylation inhibitor. These findings highlight that dynamic DNA methylation plays an important role in seasonal-dependent secondary metabolism and provide new insights for improving tea quality.

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.

Analysis of Functional Single-Nucleotide Polymorphisms (SNPs) and Leaf Quality in Tea Collection under Nitrogen-Deficient Conditions.

This study discusses the genetic mutations that have a significant association with economically important traits that would benefit tea breeders. The purpose of this study was to analyze the leaf quality and SNPs in quality-related genes in the tea plant collection of 20 mutant genotypes growing without nitrogen fertilizers. Leaf N-content, catechins, L-theanine, and caffeine contents were analyzed in dry leaves via HPLC. Additionally, the photochemical yield, electron transport efficiency, and non-photochemical quenching were analyzed using PAM-fluorimetry. The next generation pooled amplicon-sequencing approach was used for SNPs-calling in 30 key genes related to N metabolism and leaf quality. The leaf N content varied significantly among genotypes (p ≤ 0.05) from 2.3 to 3.7% of dry mass. The caffeine content varied from 0.7 to 11.7 mg g-1, and the L-theanine content varied from 0.2 to 5.8 mg g-1 dry leaf mass. Significant positive correlations were detected between the nitrogen content and biochemical parameters such as theanine, caffeine, and most of the catechins. However, significant negative correlations were observed between the photosynthetic parameters (Y, ETR, Fv/Fm) and several biochemical compounds, including rutin, Quercetin-3-O-glucoside, Kaempferol-3-O-rutinoside, Kaempferol-3-O-glucoside, Theaflavin-3'-gallate, gallic acid. From our SNP-analysis, three SNPs in WRKY57 were detected in all genotypes with a low N content. Moreover, 29 SNPs with a high or moderate effect were specific for #316 (high N-content, high quality) or #507 (low N-content, low quality). The use of a linear regression model revealed 16 significant associations; theaflavin, L-theanine, and ECG were associated with several SNPs of the following genes: ANSa, DFRa, GDH2, 4CL, AlaAT1, MYB4, LHT1, F3'5'Hb, UFGTa. Among them, seven SNPs of moderate effect led to changes in the amino acid contents in the final proteins of the following genes: ANSa, GDH2, 4Cl, F3'5'Hb, UFGTa. These results will be useful for further evaluations of the important SNPs and will help to provide a better understanding of the mechanisms of nitrogen uptake efficiency in tree crops.

The complete chloroplast genome sequence of Camellia atrothea (Ericales: Theaceae).

Camellia atrothea H.T. Chang, H.S. Wang & B.H. Chen 1983 is a native Camellia species in China, which has significant economic and breeding values. Here, we assembled and analyzed the complete chloroplast genome sequence of C. atrothea based on the whole genome sequencing data from Illumina NovaSeq6000 platform. Results show that the complete chloroplast genome of C. atrothea is 157,099 bp in length, and comprises a large single-copy (LSC, 86,643 bp) region, a small single-copy (SSC, 18,276 bp) region, and a pair of inverted repeats (IRs, 26,090 bp) region. Annotation of the C. atrothea chloroplast genome predicts a total of 133 genes, including 88 protein-coding genes, 8 ribosomal RNA genes, and 37 transfer RNA genes. The overall GC content accounts for 37.30% of the genome, which is clearly consistent with that of other previously published Camellia species from section Thea. Phylogenetic analyses of the whole chloroplast genome sequences of 14 Camellia species using Actinidia chinensis as outgroup reveals that C. atrothea has a close relationship with C. leptophylla.

Recent progress and perspectives on physiological and molecular mechanisms underlying cold tolerance of tea plants.

Tea is one of the most consumed and widely planted beverage plant worldwide, which contains many important economic, healthy, and cultural values. Low temperature inflicts serious damage to tea yields and quality. To cope with cold stress, tea plants have evolved a cascade of physiological and molecular mechanisms to rescue the metabolic disorders in plant cells caused by the cold stress; this includes physiological, biochemical changes and molecular regulation of genes and associated pathways. Understanding the physiological and molecular mechanisms underlying how tea plants perceive and respond to cold stress is of great significance to breed new varieties with improved quality and stress resistance. In this review, we summarized the putative cold signal sensors and molecular regulation of the CBF cascade pathway in cold acclimation. We also broadly reviewed the functions and potential regulation networks of 128 cold-responsive gene families of tea plants reported in the literature, including those particularly regulated by light, phytohormone, and glycometabolism. We discussed exogenous treatments, including ABA, MeJA, melatonin, GABA, spermidine and airborne nerolidol that have been reported as effective ways to improve cold resistance in tea plants. We also present perspectives and possible challenges for functional genomic studies on cold tolerance of tea plants in the future.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant.

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants.

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Diverse roles of MYB transcription factors in regulating secondary metabolite biosynthesis, shoot development, and stress responses in tea plants (Camellia sinensis).

Tea (Camellia sinensis) is concocted from tea plant shoot tips that produce catechins, caffeine, theanine, and terpenoids, which collectively determine the rich flavors and health benefits of the infusion. However, little is known about the integrated regulation of shoot tip development and characteristic secondary metabolite biosynthesis in tea plants. Here, we demonstrate that MYB transcription factors (TFs) play key and yet diverse roles in regulating leaf and stem development, secondary metabolite biosynthesis, and environmental stress responses in tea plants. By integrating transcriptomic and metabolic profiling data in different tissues at a series of developmental stages or under various stress conditions, alongside biochemical and genetic analyses, we predicted the MYB TFs involved in regulating shoot development (CsMYB2, 98, 107, and 221), epidermal cell initiation (CsMYB184, 41, 139, and 219), stomatal initiation (CsMYB113 and 153), and the biosynthesis of flavonoids (including catechins, anthocyanins, and flavonols; CsMYB8 and 99), caffeine (CsMYB85 and 86), theanine (CsMYB9 and 49), carotenoids (CsMYB110), mono-/sesquiterpenoid volatiles (CsMYB68, 147, 148, and 193), lignin (CsMYB164 and 192), and indolic compounds (CsMYB139, 162, and 198), as well as the MYB TFs that are likely involved in hormone signaling-mediated environmental stress and defense responses. We characterized the functions of some key MYBs in regulating flavonoid and carotenoid biosynthesis for tea quality and flavor. This study provides a cross-family analysis of MYBs in tea alongside new insights into the coordinated regulation of tea plant shoot development and secondary metabolism, paving the way towards understanding of tea quality trait formation and genetic improvement of quality tea plants.

CsMYB1 integrates the regulation of trichome development and catechins biosynthesis in tea plant domestication.

Tea trichomes synthesize numerous specialized metabolites to protect plants from environmental stresses and contribute to tea flavours, but little is known about the regulation of trichome development. Here, we showed that CsMYB1 is involved in the regulation of trichome formation and galloylated cis-catechins biosynthesis in tea plants. The variations in CsMYB1 expression levels are closely correlated with trichome indexes and galloylated cis-catechins contents in tea plant populations. Genome resequencing showed that CsMYB1 may be selected in modern tea cultivars, since a 192-bp insertion in CsMYB1 promoter was found exclusively in modern tea cultivars but not in the glabrous wild tea Camellia taliensis. Several enhancers in the 192-bp insertion increased CsMYB1 transcription in modern tea cultivars that coincided with their higher galloylated cis-catechins contents and trichome indexes. Biochemical analyses and transgenic data showed that CsMYB1 interacted with CsGL3 and CsWD40 and formed a MYB-bHLH-WD40 (MBW) transcriptional complex to activate the trichome regulator genes CsGL2 and CsCPC, and the galloylated cis-catechins biosynthesis genes anthocyanidin reductase and serine carboxypeptidase-like 1A. CsMYB1 integratively regulated trichome formation and galloylated cis-catechins biosynthesis. Results suggest that CsMYB1, trichome and galloylated cis-catechins are coincidently selected during tea domestication by harsh environments for improved adaption and by breeders for better tea flavours.

Identifying multi-functional bioactive peptide functions using multi-label deep learning.

The bioactive peptide has wide functions, such as lowering blood glucose levels and reducing inflammation. Meanwhile, computational methods such as machine learning are becoming more and more important for peptide functions prediction. Most of the previous studies concentrate on the single-functional bioactive peptides prediction. However, the number of multi-functional peptides is on the increase; therefore, novel computational methods are needed. In this study, we develop a method MLBP (Multi-Label deep learning approach for determining the multi-functionalities of Bioactive Peptides), which can predict multiple functions including anti-cancer, anti-diabetic, anti-hypertensive, anti-inflammatory and anti-microbial simultaneously. MLBP model takes the peptide sequence vector as input to replace the biological and physiochemical features used in other peptides predictors. Using the embedding layer, the dense continuous feature vector is learnt from the sequence vector. Then, we extract convolution features from the feature vector through the convolutional neural network layer and combine with the bidirectional gated recurrent unit layer to improve the prediction performance. The 5-fold cross-validation experiments are conducted on the training dataset, and the results show that Accuracy and Absolute true are 0.695 and 0.685, respectively. On the test dataset, Accuracy and Absolute true of MLBP are 0.709 and 0.697, with 5.0 and 4.7% higher than those of the suboptimum method, respectively. The results indicate MLBP has superior prediction performance on the multi-functional peptides identification. MLBP is available at https://github.com/xialab-ahu/MLBP and http://bioinfo.ahu.edu.cn/MLBP/.

PredAPP: Predicting Anti-Parasitic Peptides with Undersampling and Ensemble Approaches.

Anti-parasitic peptides (APPs) have been regarded as promising therapeutic candidate drugs against parasitic diseases. Due to the fact that the experimental techniques for identifying APPs are expensive and time-consuming, there is an urgent need to develop a computational approach to predict APPs on a large scale. In this study, we provided a computational method, termed PredAPP (Prediction of Anti-Parasitic Peptides) that could effectively identify APPs using an ensemble of well-performed machine learning (ML) classifiers. Firstly, to solve the class imbalance problem, a balanced training dataset was generated by the undersampling method. We found that the balanced dataset based on cluster centroid achieved the best performance. Then, nine groups of features and six ML algorithms were combined to generate 54 classifiers and the output of these classifiers formed 54 feature representations, and in each feature group, we selected the feature representation with best performance for classification. Finally, the selected feature representations were integrated using logistic regression algorithm to construct the prediction model PredAPP. On the independent dataset, PredAPP achieved accuracy and AUC of 0.880 and 0.922, respectively, compared to 0.739 and 0.873 of AMPfun, a state-of-the-art method to predict APPs. The web server of PredAPP is freely accessible at http://predapp.xialab.info and https://github.com/xialab-ahu/PredAPP .

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

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

Identification of Orphan Genes in Unbalanced Datasets Based on Ensemble Learning.

Orphan genes are associated with regulatory patterns, but experimental methods for identifying orphan genes are both time-consuming and expensive. Designing an accurate and robust classification model to detect orphan and non-orphan genes in unbalanced distribution datasets poses a particularly huge challenge. Synthetic minority over-sampling algorithms (SMOTE) are selected in a preliminary step to deal with unbalanced gene datasets. To identify orphan genes in balanced and unbalanced Arabidopsis thaliana gene datasets, SMOTE algorithms were then combined with traditional and advanced ensemble classified algorithms respectively, using Support Vector Machine, Random Forest (RF), AdaBoost (adaptive boosting), GBDT (gradient boosting decision tree), and XGBoost (extreme gradient boosting). After comparing the performance of these ensemble models, SMOTE algorithms with XGBoost achieved an F1 score of 0.94 with the balanced A. thaliana gene datasets, but a lower score with the unbalanced datasets. The proposed ensemble method combines different balanced data algorithms including Borderline SMOTE (BSMOTE), Adaptive Synthetic Sampling (ADSYN), SMOTE-Tomek, and SMOTE-ENN with the XGBoost model separately. The performances of the SMOTE-ENN-XGBoost model, which combined over-sampling and under-sampling algorithms with XGBoost, achieved higher predictive accuracy than the other balanced algorithms with XGBoost models. Thus, SMOTE-ENN-XGBoost provides a theoretical basis for developing evaluation criteria for identifying orphan genes in unbalanced and biological datasets.

Genome-Wide Identification of Seven Polyamine Oxidase Genes in Camellia sinensis (L.) and Their Expression Patterns Under Various Abiotic Stresses.

Polyamines (PAs) in plant play a critical role in growth and development and in response to environmental stress. Polyamine oxidase (PAO) is a flavin adenine dinucleotide dependent enzyme that plays a major role in PA catabolism. For the first time, PAO genes in tea plant were screened for the whole genome-wide and seven CsPAO genes were identified, which were named CsPAO1-7. Phylogenetic tree analysis revealed seven CsPAO protein sequences classed into three groups, including clade I, III, and IV. Compared with other plants, the tea plant lacked clade II members. Genetic structure and tissue specific expression analysis showed that there were significant differences among members of the CsPAO gene family. Among members of the CsPAOs family, CsPAO4 and CsPAO5 contain more introns and are highly expressed in various organizations. CsPAO1, CsPAO4, and CsPAO5 genes were cloned and expressed heterologously to verify theirs function. Heat map showed high response of CsPAO5 to drought stress, while CsPAO1 and CsPAO2 were sensitive to changes in nitrogen nutrition. Furthermore, exogenous abscisic acid (ABA) treatment indicated that the expression of most CsPAO genes in roots and leaves was significantly induced. In the root, Spm content increased significantly, while Put and Spd content decreased, suggesting that ABA has great influence on the biosynthesis of PAs. Anaerobic treatment of picked tea leaves showed that the decomposition of PAs was promoted to a certain extent. The above data help to clarify the role of CsPAO in response abiotic and nitrogen nutritional stresses in tea plants, and provide a reference perspective for the potential influence of PAs on the tea processing quality.

Transcriptome data reveal conserved patterns of fruiting body development and response to heat stress in the mushroom-forming fungus Flammulina filiformis.

Mushroom-forming fungi are complex multicellular organisms that form the basis of a large industry, yet, our understanding of the mechanisms of mushroom development and its responses to various stresses remains limited. The winter mushroom (Flammulina filiformis) is cultivated at a large commercial scale in East Asia and is a species with a preference for low temperatures. This study investigated fruiting body development in F. filiformis by comparing transcriptomes of 4 developmental stages, and compared the developmental genes to a 200-genome dataset to identify conserved genes involved in fruiting body development, and examined the response of heat sensitive and -resistant strains to heat stress. Our data revealed widely conserved genes involved in primordium development of F. filiformis, many of which originated before the emergence of the Agaricomycetes, indicating co-option for complex multicellularity during evolution. We also revealed several notable fruiting-specific genes, including the genes with conserved stipe-specific expression patterns and the others which related to sexual development, water absorption, basidium formation and sporulation, among others. Comparative analysis revealed that heat stress induced more genes in the heat resistant strain (M1) than in the heat sensitive one (XR). Of particular importance are the hsp70, hsp90 and fes1 genes, which may facilitate the adjustment to heat stress in the early stages of fruiting body development. These data highlighted novel genes involved in complex multicellular development in fungi and aid further studies on gene function and efforts to improve the productivity and heat tolerance in mushroom-forming fungi.

Improved hybrid de novo genome assembly and annotation of African wild rice, Oryza longistaminata, from Illumina and PacBio sequencing reads.

African wild rice Oryza longistaminata, one of the eight AA- genome species in the genus Oryza, possesses highly valued traits, such as the rhizomatousness for perennial rice breeding, strong tolerance to biotic and abiotic stresses, and high biomass production on poor soils. To obtain the high-quality reference genome for O. longistaminata we employed a hybrid assembly approach through incorporating Illumina and PacBio sequencing datasets. The final genome assembly comprised only 107 scaffolds and was approximately ∼363.5 Mb, representing ∼92.7% of the estimated African wild rice genome (∼392 Mb). The N50 lengths of the assembled contigs and scaffolds were ∼46.49 Kb and ∼6.83 Mb, indicating ∼3.72-fold and ∼18.8-fold improvement in length compared to the earlier released assembly (∼12.5 Kb and 364 Kb, respectively). Aided with Hi-C data and syntenic relationship with O. sativa, these assembled scaffolds were anchored into 12 pseudo-chromosomes. Genome annotation and comparative genomic analysis reveal that lineage-specific expansion of gene families that respond to biotic- and abiotic stresses are of great potential for mining novel alleles to overcome major diseases and abiotic adaptation in rice breeding programs. This reference genome of African wild rice will greatly enlarge the existing database of rice genome resources and unquestionably form a solid base to understand genomic basis underlying highly valued phenotypic traits and search for novel gene sources in O. longistaminata for the future rice breeding programs.

Genome-Wide Characterization of the C-repeat Binding Factor (CBF) Gene Family Involved in the Response to Abiotic Stresses in Tea Plant (Camellia sinensis).

C-repeat (CRT)/dehydration responsive element (DRE)-binding factor CBFs, a small family of genes encoding transcriptional activators, play important roles in plant cold tolerance. In this study, a comprehensive genome-wide analysis was carried out to identify and characterize the functional dynamics of CsCBFs in tea plant (Camellia sinensis). A total of 6 CBF genes were obtained from the tea plant genome and named CBF1-6. All of the CsCBFs had an AP2/ERF DNA-binding domain and nuclear localization signal (NLS) sequence. CsCBF-eGFP fusion and DAPI staining analysis confirmed the nuclear localization of the CsCBFs. Transactivation assays showed that the CsCBFs, except CsCBF1, had transcriptional activity. CsCBF expression was differentially induced by cold, heat, PEG, salinity, ABA, GA, MeJA, and SA stresses. In particular, the CsCBF genes were significantly induced by cold treatments. To further characterize the functions of CsCBF genes, we overexpressed the CsCBF3 gene in Arabidopsis thaliana plants. The resulting transgenic plants showed increased cold tolerance compared with the wild-type Arabidopsis plant. The enhanced cold tolerance of the transgenic plants was potentially achieved through an ABA-independent pathway. This study will help to increase our understanding of CsCBF genes and their contributions to stress tolerance in tea plants.

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

BACKGROUND: Tea plant is one of the most important non-alcoholic beverage crops worldwide. While potassium (K+) is an essential macronutrient and greatly affects the growth and development of plants, the molecular mechanism underlying K+ uptake and transport in tea plant root, especially under limited-K+ conditions, is still poorly understood. In plants, HAK/KUP/KT family members play a crucial role in K+ acquisition and translocation, growth and development, and response to stresses. Nevertheless, the biological functions of these genes in tea plant are still in mystery, especially their roles in K+ uptake and stress responses. RESULTS: In this study, a total of 21 non-redundant HAK/KUP/KT genes (designated as CsHAKs) were identified in tea plant. Phylogenetic and structural analysis classified the CsHAKs into four clusters (I, II, III, IV), containing 4, 8, 4 and 5 genes, respectively. Three major categories of cis-acting elements were found in the promoter regions of CsHAKs. Tissue-specific expression analysis indicated extremely low expression levels in various tissues of cluster I CsHAKs with the exception of a high root expression of CsHAK4 and CsHAK5, a constitutive expression of clusters II and III CsHAKs, and a moderate cluster IV CsHAKs expression. Remarkably, the transcript levels of CsHAKs in roots were significantly induced or suppressed after exposure to K+ deficiency, salt and drought stresses, and phytohormones treatments. Also notably, CsHAK7 was highly expressed in all tissues and was further induced under various stress conditions. Therefore, functional characterization of CsHAK7 was performed, and the results demostrated that CsHAK7 locates on plasma membrane and plays a key role in K+ transport in yeast. Taken together, the results provide promising candidate CsHAKs for further functional studies and contribute to the molecular breeding for new tea plants varieties with highly efficient utilization of K+. CONCLUSION: This study demonstrated the first genome-wide analysis of CsHAK family genes of tea plant and provides a foundation for understanding the classification and functions of the CsHAKs in tea plants.