Next Article in Journal
Molecular Mechanisms Employed by Neurons to Receive and Transduce Signals Essential for Learning and Memory
Next Article in Special Issue
Tissue-Specific Transcriptome and Metabolome Analysis Reveals the Response Mechanism of Brassica napus to Waterlogging Stress
Previous Article in Journal
Leveraging Natural Killer Cell Innate Immunity against Hematologic Malignancies: From Stem Cell Transplant to Adoptive Transfer and Beyond
Previous Article in Special Issue
Genome-Wide Analysis and Expression Profiling of Glutathione Reductase Gene Family in Oat (Avena sativa) Indicate Their Responses to Abiotic Stress during Seed Imbibition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 1
10.3390/ijms24010202
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification, Characterization and Expression Profiling of the RS Gene Family during the Withering Process of White Tea in the Tea Plant (Camellia sinensis) Reveal the Transcriptional Regulation of CsRS8

by
Tao Wang
1,
Yiqing Wang
1,
Jiamin Zhao
1,
Jiumei Kong
1,
Lingzhi Zhang
1,
Siyu Qi
1,
Jiajia Chen
1,
Zhidan Chen
2,
Wen Zeng
1,* and
Weijiang Sun
1,*
1
College of Horticulture, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
Anxi College of Tea Science, Fujian Agriculture and Forestry University, Quanzhou 362000, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(1), 202; https://doi.org/10.3390/ijms24010202
Submission received: 25 November 2022 / Revised: 17 December 2022 / Accepted: 18 December 2022 / Published: 22 December 2022
(This article belongs to the Special Issue Regulatory Mechanism and Network of Abiotic Stress-Response in Plants)
Browse Figures
Review Reports Versions Notes

Abstract

:
Raffinose synthetase (RS) is a key enzyme in the process of raffinose (Raf) synthesis and is involved in plant development and stress responses through regulating Raf content. As a sweetener, Raf makes an important contribution to the sweet taste of white tea. However, studies on the identification, analysis and transcriptional regulation of CsRSs (Camellia sinensis RS genes) are still lacking. In this study, nine CsRSs were identified from the tea plant (Camellia sinensis) genome database. The CsRSs were classified into five groups in the phylogenetic tree. Expression level analysis showed that the CsRSs varied in different parts of the tea plant. Transcriptome data showed that CsRSs could respond to persistent drought and cold acclimation. Except for CsRS5 and CsRS9, the expression pattern of all CsRSs increased at 12 h and decreased at 30 h during the withering process of white tea, consistent with the change trend of the Raf content. Furthermore, combining yeast one-hybrid assays with expression analysis, we found that CsDBB could potentially regulate the expression of CsRS8. Our results provide a new perspective for further research into the characterization of CsRS genes and the formation of the white tea flavour.

1. Introduction

Raffinose (Raf) is an oligosaccharide of the raffinose family and is widely present in plants [1]. The first step of Raf biosynthesis is to use UDP-galactose (UDP-Gal) and L-myo-inositol as substrates to produce galactinol (Gol, 1-O-α-D-galactopyranosyl-L-myo-inositol) under the catalysis of galactitol synthase (GolS, EC 2.4.1.123) [2]. Next, the galactosyl from myo-inositol galactinol is connected with sucrose (Suc) through α-1 and 6-glycosidic bonds to generate Raf, which is catalyzed by raffinose synthase (RS, EC 2.4.1.82) [3]. Notably, the reaction process is reversible [4]. Current studies show that Raf mostly accumulates in plant seeds, which can promote the process of cytoplasmic vitrification [5]. In addition, Raf can participate in the protection of the cell membrane by inhibiting the crystallization of sucrose during seed dehydration, improving seed storage tolerance and maintaining seed vitality [6,7,8]. Recent studies have found that Raf can also accumulate in leaves and roots as a response mechanism to abiotic stress and protect plant growth. For example, the content of Raf in maize leaves significantly increases under salt stress and high temperature stress [9]; the content of Raf in Arabidopsis shows a positive correlation with its freezing tolerance [10]; the Raf in the cytoplasm is transported to the chloroplast to protect the chloroplast structure in response to cold stress [11]; moreover, as a scavenger of ROS in plant cells, Raf can reduce the damage caused by the ROS oxidative stress reaction [12]. In the food industry, Raf is widely used as a food-grade sweet additive. Although Raf cannot be absorbed by humans, it can be used by beneficial bacteria (Bifidobacterium and Lactobacillus) in the intestine and inhibit the growth of harmful bacteria [13,14,15].
RS is a key enzyme that catalyzes the formation of Raf, and the RS gene family has been identified and characterized in many plants, including Arabidopsis [3], corn [16], lentil [17], cotton [18] and poplar [19]. The expression of RS genes can be induced by plant seed development and stress; GmRS2A and GmRS2B are highly expressed in the middle of soybean seed maturity, and the content of Raf also increases [20]; the seed vigor of maize mutants with ZmRS knockout decreases [21]; the overexpression of AtRS5 in Arabidopsis improves the seed germination rate and its seedlings’ drought resistance [22]. Abiotic stress induces the expression of the RS gene in plant leaves. For example, the expression of the RS gene in rice [23], maize [24] and melon [25] was upregulated under cold stress; under drought stress, high-temperature stress, high-salt stress and other stress treatments, the expression of RS was induced to varying degrees [26,27]; AtRS5 (At5g40390) in Arabidopsis is the only AtRS member involved in Raf accumulation in leaves under abiotic stress [3]. Some transcription factors are involved in regulating the expression of RS, which can affect the resistance by regulating the biosynthesis of Raf. The overexpression of the Arabidopsis heat shock transcription factor AtHsf2 can induce an increase in AtRS2 expression and in the Raf content and produce transgenic plants with greater stress resistance [12]. In addition, the overexpression of BnHsf4a in Arabidopsis can also induce the expression of AtRS2 [3]; ZmHsf2 is a heat shock transcription factor in maize, which can induce the expression of RS and improve the heat tolerance of maize [28]; the overexpression of OsWRKY11 in rice can increase the expression of OsRS and the content of Raf, resulting in improved drought resistance [29]; in addition, transcription factors such as CBF and ERF can also improve the plant’s resistance to stress by regulating the expression of RS [30,31].
Tea is a global consumer beverage, second only to water in consumption, and is processed from leaves of the tea plant(Camellia sinensis) [32]. According to different processing technologies, tea is divided into six categories (white tea, oolong tea, black tea, green tea, yellow tea and black tea) [33]. In recent years, white tea has become popular with consumers because of its health benefits and fresh, sweet taste, and its output and consumption are increasing rapidly [34]. Compared with other types of tea, the processing technology of white tea is simpler, including long-term withering and drying. Among these processes, withering is a slow process of water deficit stress in leaves. During this period, various flavor compounds in the leaves change, resulting in the sweet taste of white tea [35]. Amino acids and soluble sugars contribute to the sweet and umami taste of tea infusions, particularly theanine and maltose [36,37]. However, as an important flavor compound, Raf is also involved in the resistance physiology during withering and is involved in determining the plant’s taste and resistance to abiotic stress, but its changes in this process have not been reported. The RS gene has been identified and characterized in various plants but has rarely been reported in Camellia sinensis until now.
Due to its important role in taste contribution and the abiotic stress response, this paper focuses on the whole genome analysis of the CsRS family in Camellia sinensis and compares it with the RS of other plants, including the tea genome database comparison, phylogenetic tree analysis, chromosome distribution and gene structure. The expression pattern of CsRSs and the change in Raf content during the withering process of white tea were analyzed using transcriptome data from the TPIA platform, qRT-PCR and ultrahigh-performance liquid chromatography (UHPLC). In addition, the regulator of CsRS8 was screened out by yeast one-hybrid assays. These results help us better understand the CsRS genes and the changes in chemical compounds during the withering process of white tea and provide a new perspective for the characterization of RS in Camellia sinensis in the future.

2. Results

2.1. Identification of RS Family Genes in Camellia sinensis

We identified nine RS genes from the Camellia sinensis genome, named based on the chromosome ID, from CsRS1 to CsRS9. The CsRS genes exhibited significant variations in sequence length, protein size and physicochemical properties, with the sequence lengths ranging from 1944 bp to 2445 bp, the protein sizes ranging from 647 aa (CsRS2) to 814 aa (CsRS9), and the molecular weights (MWs) ranging from 71.63 kDa (CsRS2) to 89.28 kDa (CsRS9). The prediction results of subcellular localization showed that eight of the nine CsRSs were located in the chloroplast or cytoplasm, and only one(CsRS2) was located in the nucleus. Table S2 shows the details, including protein ID, protein length, molecular weight, isoelectric point, instability index, aliphatic index, grand average of hydropathicity, and subcellular localization.

2.2. Phylogenetic Analysis of the CsRS Family

Based on the phylogenetic tree, all RS members can be categorized into seven different clusters (I-RS, II-RS, III-RS, IV-RS, V-RS, VI-RS, and VII-RS) (Figure 1A). IV-RS and VII are the largest clusters, both with 11 members, and I-RS has only three members. AtRS5 is a protein induced by abiotic stress and involved in the biosynthesis of Raf [3]. Three CsRS members (CsRS1, CsRS2 and CsRS8) were clustered on cluster IV-RS with AtRS5. None of the CsRSs were classified as III-RS or V-RS.

2.3. Chromosome Location and Synteny Analysis of CsRSs

To confirm the CsRSs genomic distribution, we mapped them to the published tea plant genome (Figure 1B). CsRSs were unevenly distributed on chromosomes, and nine CsRSs were distributed on five chromosomes and three contigs. One CsRS was distributed on chromosome 1, chromosome 7, chromosome 8, chromosome 11, contig 284, contig 655 and contig 1151, and two CsRSs were distributed on chromosome 2. The Ka/Ks ratio was used to examine whether segment duplication was influenced by selection pressure during evolution. The Ka/Ks ratio of CsRS1/CsRS2 was 0.099, the Ka/Ks ratio of CsRS1/CsRS8 was 0.115, and the Ka/Ks ratio of CsRS2/CsRS8 was 0.85. The Ka/Ks ratio for all paired genes was less than 1.

2.4. Gene Structure and Conserved Motifs of CsRSs

To classify the nine CsRSs, a phylogenetic tree was constructed using MEGA X [38] based on the sequences. The nine CsRSs were divided into three classes: CsRS3/CsRS5/CsRS7, CsRS4/CsRS6/CsRS9, and CsRS1/CsRS2/CsRS8 (Figure 2A). The prediction results of the conservative motif showed that all CsRS contained 10 motifs, ranging from 6 to 50 amino acids (Figure 2B). All the identified CsRSs contain RS protein domain (Figure 2C). The number of exons per gene ranged from 4 to 14 (Figure 2D), and the genes in the same cluster showed highly similar exon numbers: 3 members (CsRS3, CsRS5 and CsRS7) in Cluster a contained 12 to 14 exons, 3 members in Cluster b (CsRS4, CsRS6 and CsRS9) contained 13 to 14 exons, and 3 members (CsRS1, CsRS2 and CsRS8) in Cluster c contained 4 to 5 exons. Significantly, CsRS1 and CsRS3 had no introns.

2.5. Secondary and Three-Dimensional Structures of CsRS

The predicted secondary structure of the CsRSs mainly includes random coils, α-helices and extended chains (Table 1). Swiss was used to predict the tertiary structure of the CsRSs, and the template coverage of the predicted CsRS structure was more than 99%. The predicted conformation of CsRS proteins was not identical, as shown in Figure S1A. Ramachandran’s online evaluation of protein structure showed that the amino acid sites within the allowable range were 98.5~98.8%, indicating that the constructed tertiary structure model of CsRSs was accurate and reliable (Figure S1B).

2.6. Cis-Elements in the Promoters of CsRSs

Fifty-six cis-element components were identified, with the number of cis-element components ranging from 87 (CsRS5) to 142 (CsRS8) (Figure 2E). Six types of cis-elements are included, namely light responsive, abiotic and biotic stress, phytohormone responsive, transcription factor recognition and binding sites, tissue, and core. As shown in Figure 2B, TATA and CAAT boxes were found in all nine CsRS promoter regions, and the A-box was specifically found in CsRS8 and CsRS9. Seventeen different types of cis-elements correlated with the light response were found in the promoter region of CsRSs, and most CsRSs had Box 4, G-box and GT1-motif elements. The promoter regions of CsRSs also contained the cis-elements of phytohormone responsiveness, including abscisic acid responsiveness (AAGAA-motif and ABRE), MeJA responsiveness (CGTCA-motif and TGACG-motif), and ethylene responsiveness (ERE). Some cis-elements involved in abiotic/biotic stress responses and tissue-specific expression were also included: anaerobic induction (ARE), which was found in all CsRS promoter regions, and stress responsiveness (STRE), coercion responsiveness (WRE3), injury-related elements (WUN-motif) and meristem expression (CAT-box and CCGTTC-motif), which were irregularly distributed in CsRS promoter regions. Specifically, we found nine types of cis-elements of transcription factor recognition and binding sites. The myb-binding site element was present in CsRS promoter regions, except CsRS2 and CsRS4, myb and protein binding sites (ABRE3, AT-rich element, Box III, CCAAT-box, HD-zip 3, MBS and MES) with sporadic distribution. Only CsRS5 and CsRS9 had the MYB recognition site element.

2.7. Expression Patterns of CsRSs in Different Organs

All CsRS genes were differentially expressed among different tissues (Figure 3A). The transcript levels of CsRS1, CsRS3, CsRS7 and CsRS8 were lower than those of the other members in all tissues. The transcript levels of CsRS2, CsRS5 and CsRS9 were high in mature and old leaves. Interestingly, the transcript levels of CsRS4 and CsRS6 were high in all tissues, especially in roots and flowers, followed by young leaves and fruits.

2.8. Expression Patterns of CsRSs during Continuous Drought Stress and Cold Acclimation

To investigate the expression patterns of CsRSs under cold acclimation and persistent drought, we downloaded and analyzed previous transcriptome data from the TPIA platform [39,40]. As shown in Figure 3B, with the increase of drought stress treatment time, all members of cluster IV-RS (CsRS1, CsRS2, and CsRS8) and VII-RS (CsRS6 and CsRS9) were continuously upregulated; members of cluster II were upregulated at 24 h, and downregulated at 48 h and 72 h. Under cold acclimation (CA1-6 h, CA1-7 d and CA2-7 d), the expression of CsRS2, CsRS5 and CsRS6 were upregulated, and downregulated at DA-7 d.

2.9. Expression Patterns of CsRSs and the Raf Content in the Withering Process of White Tea

During the withering process, the expression levels of CsRSs were upregulated at 12 h and downregulated at 30 h, except CsRS5 and CsRS6 (Figure 3C). There was no significant difference in the expression level of CsRS5 at 0 h and 12 h, but it was significantly upregulated at 30 h. Among all CsRSs, the expression level of CsRS8 peaked at 12 h, was approximately 366.5 times higher than at 0 h, and decreased at 30 h. Correspondingly, we detected the Raf content during the withering process and found that its content increased significantly at 12 hours and decreased at 30 hours, but it was still higher than at 0 h (Figure 4A). Correlation analysis showed that there was a significant positive correlation between Raf content and the expression of CsRS1 and CsRS8, and a strong correlation between Raf content and CsRS2, CsRS3, CsRS4, CsRS6 and CsRS7 was found (Figure S2).

2.10. Screening out Transcription Factors Regulating CsRS8

The cDNA library was transformed into a bait yeast strain and inoculated on SD/-Leu + 500 ng/mL AbA medium for screening, and one transcription factor (CsDBB) was identified, which may directly bind to the CsRS8 promoter (Figure S3A). We predicted the subcellular localization of CsDBB, and both the Wolf and Cell-PLoc results showed that CsDBB localized in the nucleus and had typical transcription factor characteristics. According to the conserved domain analysis, we found that CsDBB contains a double B-box zinc finger protein structure, which is a typical feature of DBB transcription factors (Figure S3B). The qRT-PCR results showed that CsDBB was expressed in different tissues (Figure 4B). During the withering process, the expression level of CsDBB first increased and then decreased (Figure 4C). Furthermore, the correlation analysis showed that CsDBB was positively correlated with CsRS8 at the expression level and positively correlated with the Raf content (Figure S3C).

2.11. The Interaction and Regulation Network of CsRSs

To investigate the regulation of the expression of CsRSs by transcription factors (TF), the PlantTFDB database was used to predict transcription factor binding sites on promoters. A total of 8 TF families (AP2, B3, BBR-BPC, C2H2, Dof, ERF, LBD and MIKC_MADS) were predicted can bind the CsRSs (CsRS2, CsRS5, CsRS7, CsRS8 and CsRS9) promoters, which covers 1, 1, 3, 1, 9, 44, 2 and 1 members. The ERF family had the largest number of binding sites (114), while the AP2 family had the smallest binding sites (2). The largest number and variety of TFs were identified in the promoters of CsRS8, but only a TF family(C2H2) was identified in CsRS9 (Figure 5A).
We obtained the protein interaction regulatory network of CsRSs through the yeast one-hybrid, STRING and PlantTFDB databases. Nine CsRSs were compared to three AtRSs in Arabidopsis; CsRS4 corresponded to AtRFS1, CsRS1 and CsRS2; CsRS8 corresponded to AtRFS5, CsRS3 and CsRS5, CsRS6, and CsRS7; CsRS9 corresponded to AtSIP2; a total of 14 functional genes were predicted to interact with CsRSs; 64 transcription factors were found to possibly regulate the expression of CsRSs (Figure 5B). There are 11 galactosidases and 1 stachyose synthase in 14 functional proteins, and 64 transcription factors mainly belonging to the AP2, Dof, C2H2, LOB and other transcription factor families.

3. Discussion

White tea is mainly produced in China, is popular with consumers all over the world and is characterized by a sweet taste [34]. The manufacturing process of white tea includes two steps, withering and drying, and withering is the key process to form the characteristics of white tea. During the withering process, leaves undergo various stresses, mainly drought. Under water deficit stress, some stress-related flavor compounds and aroma compounds increase significantly, which contribute to the unique taste of white tea [35]. Raf plays a key role in resisting abiotic stress and taste formation [12,41]. In vivo, Raf biosynthesis is catalyzed by RS, but there is a lack of relevant research in Camellia sinensis [3]. Therefore, we sought to identify and characterize Raf in Camellia sinensis and the changes in Raf and its related genes during the withering process of white tea. In this study, we conducted a comprehensive and systematic analysis of the RS gene family of Camellia sinensis and explored the expression level and Raf content during the withering period.

3.1. Evolution of the CsRS Gene Family

RS is a key enzyme involved in Raf biosynthesis, and RS genes have been identified in many plants, including Arabidopsis [3,42], cassava [43], maize [16] and sesame [44]. In this study, we identified nine CsRS genes from Camellia sinensis distributed on five chromosomes and three contigs. Collinearity analysis showed that CsRS1, CsRS2 and CsRS8 had gene replication events, and their Ka/Ks values were all less than 1 (Figure 1B), indicating that CsRSs mainly experienced purification selection after replication [45]. Based on the evolutionary tree, we divided the selected RSs into six clusters, and nine CsRSs were irregularly distributed in five clusters (Figure 1A). Similarly, the RS evolutionary trees constructed in sesame [44] and maize [16] were also divided into six clusters. Previous studies have reported that the raffinose synthase domain (PF05691) exists in both RS and STS (stachyose synthase) proteins, but a Gol-Raf galactosyltransferase/galactosyl hydrolase domain occurs exclusively in STS sequences [3]. In our study, the sequences containing complete RS domains were selected for phylogenetic construction, and none of the CsRSs grouped together with the STS from other species. CsRS1, CsRS2, CsRS8 and AtRS5 from Arabidopsis are clustered in IV-RS, and CsRS1, CsRS2 and CsRS8 also have similar gene structures, indicating that these proteins may have similar functions (Figure 1A). AtRS5 is involved in the biosynthesis of raffinose in leaves under abiotic stress, and it is the only RS gene involved in raffinose biosynthesis in Arabidopsis [3]. Therefore, we speculate that CsRS1, CsRS2 and CsRS8 may also participate in the accumulation of Raf in tea leaves under abiotic stress.

3.2. Potential Role of CsRS in Abiotic Stress and the Withering Process of White Tea

The promoter is a specific DNA sequence located upstream of the gene start codon ATG that affects the expression level of foreign genes [46]. In this study, we identified the homeopathic elements of the CsRS promoter region and further classified them into six categories: light response elements, stress response elements, hormone response elements, transcription factor-related elements and plant growth and development elements (Figure 2E). Our results show that in addition to numerous light response elements, there are many hormone stress response elements, including MeJA, GA, ABA, ethylene and auxin response elements, indicating that hormones may affect the expression of RS. Under ABA treatment, the mRNA level of CsRS in Cucumis sativus was positively correlated with the concentration of ABA treatment groups, and the Raf content also increased [47]. In addition, we found that there are numerous stress response elements, such as drought response (DRE core), pressure response (Box S, STRE, TC rich repeats) and low temperature (LTR). The overexpression of AtRS5 increases the content of raffinose in seeds, thus improving the desiccation tolerance of seeds [22]; the overexpression of CsRS from Cucumis sativus in tobacco can promote the accumulation of Raf in leaves, and cold stress can induce the expression of CsRS in cucumber, thereby improving the content of Raf and cold resistance [47]; RS genes are also involved in cold stress and cold acclimation in sugar beets [25], rice [23] and tea plants [48]. Combined with previous research data [39,40], the analysis of transcriptome data during continuous drought stress and cold acclimation and the relative expression of the withering process of white tea revealed that the expression of CsRSs under different stresses differed. The members of the same cluster have the same expression pattern under continuous drought stress, especially the cluster of IV-RS and VII-RS that were up-regulated continuously; CsRS2, CsRS5 and CsRS6 were induced during cold acclimation. The withering is a slow process of inducing water deficit stress in the leaves; CsRS1, CsRS2, CsRS3, CsRS4, CsRS7, CsRS8 and CsRS9 were up-regulated at 12h and down-regulated at 30h, and the content of Raf has the same trend.
In conclusion, we speculate that CsRSs may play an important role in the accumulation of Raf in the process of tea plant stress. Most notably, CsRS8 contains the most types and quantities of cis-acting elements in response to stress, which may play a more crucial role in the response to stress. The qRT-PCR results showed that most of the CsRSs were highly expressed at 12 h of withering, especially CsRS8 (Figure 3C), which had the most significant change, and its expression was significantly positively correlated with the content of Raf (Figure S2), so we speculate that CsRS8 played a significant role in the accumulation of Raf during the withering process.

3.3. CsDBB Potentially Regulates CsRS8 to Affect Raffinose Accumulation during the Withering of White Tea

Transcription factors are special structural proteins that can regulate plant growth and development and regulate the specific expression of target genes by binding cis-acting elements in the promoter region of target genes [49]. Recent studies have shown that RS is regulated by various transcription factors. The overexpression of the heat shock transcription factors Arabidopsis AtHSFA2 and Brassica napus BnHSFA4a can induce the expression of AtRS2 [12,50]. Maize ZmDREB2A can regulate the expression of RS and prolong the preservation time of seeds [51]. In addition, WRKY, ERF and other types of transcription factors are involved in regulating the expression of RS genes [31,52]. In our study, we predicted that 59 transcription factors belonging to eight TF families (AP2, B3, BCR-BPC, C2H2, Dof, ERF, LBD and MIKC_MADS) could potentially regulate the expression of CsRSs. Cis-element prediction results showed that CsRS8 contained the most types and quantities of elements, indicating that CsRS8 may be regulated by multiple transcription factors (Figure 2E). Besides, six categories of transcription factors were identified, all of which could potentially regulate CsRS8 (Figure 5A), further confirming our speculation. DBB belongs to the subfamily of b-box transcription factors and is named because it contains two b-box structures involved in plant growth and development and abiotic stress response [53]. In our study, a CsDBB family transcription factor, CsDBB (CSS0046793), was screened by yeast one-hybrid assays (Figure S3A). We found that CsDBB was significantly upregulated and then decreased during the withering process (Figure 4B). The expression level of CsDBB was positively correlated with CsRS8, CsDBB and Raf content, according to the correlation analysis results (Figure S3C). Therefore, we speculated that CsDBB might positively regulate CsRS8 during the withering process, affecting the accumulation of Raf. We predicted the protein interaction regulation network of CsRSs and found that there are nine functional proteins of the carbohydrate metabolism pathway that have potential interactions with CsRS (Figure 5). We found several transcription factor families, such as AP2, Dof and C2H2, indicating that there is a complex protein interaction regulation network in the carbohydrate metabolism pathway that regulates sugar components in the stress environment to resist external stress.

4. Materials and Methods

4.1. Identification of CsRSs from Tea Genome

The genome database of tea cultivar “Shuchazao” (Camellia sinensis) was used in this study [54]. The identification of RS genes was carried out with HMMsearch using PF05691 as the Pfam domains, which were downloaded from the Pfam database (http://pfam.xfam.org, accessed on 24 December 2021) [55], and the HMMsearch was performed by SPDE v1.2 [56]. To verify the presence of the significant RS domain, all the sequences were examined using the Pfam database [55] and NCBI CDD database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 24 December 2021) [57]. The protein sequence, coding sequence, promoters and gene annotation were downloaded from the Tea Plant Information Archive (TPIA) (http://tpdb.shengxin.ren/index.html, accessed on 24 December 2021) [39]. The chromosome localization, collinearity and gene duplications events were analyzed using TBtools software [58]. Tools at the ProtParam in ExPASy web (https://web.expasy.org/protparam/, accessed on 24 December 2021) [59] and WoLF PSORT (https://wolfpsort.hgc.jp/, accessed on 24 December 2021) [60] were used to predict the physicochemical characteristics and the subcellular localization of the identified CsRS proteins.

4.2. Phylogenetic Tree Construction and Structural Analysis of CsRSs

To clarify the phylogenetic relationship of the CsRS genes, a multi-sequence alignment was constructed with the MEGA X software [38] with the neighbor-joining (NJ) algorithm with 2000 bootstrap replicates, and the phylogenetic tree was refined by the iTOL tool (https://itol.embl.de/, accessed on 29 January 2022) [61]. The exon/intron structures of the CsRSs were determined by TBtools software [58]. The secondary and three-dimensional structures of the CsRS were obtained by homology modeling using the SOMPA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html, accessed on 25 December 2021) and SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 21 September 2022) [62]. The CsRSs protein tertiary structure model of CsRSs was evaluated using the PROCHECK plugin in the PDBsum online database (http://www.ebi.ac.uk/thorntonsrv/databases/cgibin/pdbsum/GetPage.pl?pdbcode=index.html, accessed on 21 September 2022) [63,64].

4.3. Cis-Element Prediction for CsRSs Promoters

The promoter sequences (2 kb upstream of ATG) of CsRSs were exported by TBtools software [58], then all the promoter sequences were submitted to the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 1 February 2022) [65], and TBtools software [58] was used for visualization.

4.4. Plant Materials and Treatments

Samples for relative expression pattern analysis: a total of eight tissue samples (young leaf, mature leaf, old leaf, young stem, old stem, flower bud, seed, and root) were collected from Camellia sinensis cultivar “Huangdan” (HD) growing in Fujian Agriculture and Forestry University, Fuzhou, China (26°05′ N, 119°18′ E), in September 2020. Camellia sinensis cultivar “Fudingdahao” (FDDH) is a kind of tea variety suitable for white tea that is widely planted in southern China. Processing sample: the young shoots with one bud and two leaves were harvested from the FDDH for making white tea. The fresh leaves were processed according to the standard white tea production procedure, including wilting and drying. During the withering period, sampling occurred at 0 h, 12 h and 30 h (3 times), with 0 h as control, and the whole process was carried out in an air-conditioned room (relative humidity 70 ± 3%, temperature 37 ± 3 °C). All samples were immersed in liquid nitrogen immediately after collection, and then stored at −80 °C. Three replicates were collected for each treatment. All samples for chemical compositions and gene expression analyses were extracted and analyzed in triplicate.

4.5. Analysis of Raffinose Content in the Withering Process of White Tea

The contents of Raf were determined by ultra-high performance liquid chromatography (UHPLC) (Waters, Milford, MA, USA). All samples from the withering process of white tea (0 h, 12 h and 30 h) were dried in a vacuum freeze-dryer. A total of 0.25 g of dried and fully crushed tea samples was added to distilled water (2.5 mL, 100 °C), soaked in a 100 °C water bath for 15 min, stirred 2–3 times (every 5 min), and centrifuged to obtain the supernatant. The supernatant was transferred to a 25 mL volumetric flask, and repeated extraction once. The supernatant obtained from two extractions was combined, 15mL acetonitrile was added, and added distilled water to the total volume of 25 mL. Samples were then filtered through a 0.22 µm membrane before analysis. Chromatographic conditions were set as previous studies with minor modifications and an ACQUITY UPLC BEH Amide column (Waters, Milford, MA, USA) was used to separate Raf [66].

4.6. Analysis of CsRSs Expression Patterns

The expression patterns of CsRSs during cold acclimation and continuous drought stress were analyzed using the previous transcriptome date downloaded from the TPIA platform [39,40]. The fragments per kilobase of exon model per million mapped fragments (FPKM) values of CsRS genes were extracted from the transcriptome date, and the heatmap was drawn using TBtools software [56]. The qRT-PCR results were used to analyze the expression patterns of CsRSs in different parts of the tea plant and during the withering process of white tea. Total RNA was obtained from samples using a RNAprep Pure Plant Plus Kit (Tiangen, Beijing, China) according to manufacturer instructions. RNA samples were reverse-transcribed into cDNA using One-Step gDNA Removal (TransGen, Beijing, China), and the synthesis reaction was carried out following the manufacturer protocol. The qRT-PCR analysis was carried out on qTOWER3G (Analytik Jena AG, Jena, Germany) with a 20 μL reaction mixture, with the following thermal parameters: 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. A melt curve was performed to verify product specificity of the PCR at the end of each reaction. Relative gene transcript levels were calculated by the 2−ΔΔct [67]. Gene expression values were log2 transformed, and heatmaps were drafted by TBtools software [58]. The gene-specific primers are listed in Table S1. Changes in mRNA levels of related genes were normalized to that of CsGAPDH.

4.7. Yeast One-Hybrid Screeing

Combined with RS phylogenetic tree and correlation analysis, the 500 bp promoter sequence of the CsRS8 was inserted into the pAbAi vector to construct the bait vector pAbAi-CsRS8 promoter. Then, the bait vector was integrated into the Y1HGold yeast strain by the PEG/LiAc transformation method to form the Y1H-pAbAi-CsRS8 promoter bait yeast strain. The bait yeast strain was inoculated on SD/-Ura medium and SD/-Ura + 100/150/300/500 ng/ml Aureobasidin A(AbA) medium to determine the minimum amount of AbA to inhibit the basal expression level of the bait yeast strain inhibitory concentration. cDNA library plasmids were purchased from Nanjing Ruiyuan Biotechnology Co., Ltd. The cDNA library was transformed into the Y1H-pAbAi-CsRS8 promoter bait yeast strain and inoculated on SD/-Leu + 500 ng/mL AbA medium for selection. The plasmid was extracted from the positive yeast monoclonal, transferred to Y1H-pAbAi-CsRS8 promoter bait yeast strain, and spotted on SD/-Leu medium and SD/-Leu + 500 ng/mL AbA medium. After 2 days of incubation at 30 °C, the interaction between CsDBB and the bait sequence was observed.

4.8. Prediction of CsRSS Protein Interactions and Regulatory Networks

STRING (https://cn.string-db.org/, accessed on 1 October 2022) and PlantTFDB (http://planttfdb.gao-lab.org/, accessed on 5 October 2022) were used to predict the protein interaction regulatory network of CsRS. In the STRING database, we input the CsRSs protein sequence and searched the protein interaction network with Arabidopsis as a reference. In the plantTFDB database, we input the CsRSs promoter region (2000 bp before ATG) and set the threshold p-value < 1 × 10−7 to predict the transcription factor binding sites and find the upstream regulators. Visualization of transcription factor binding sites was performed using TBtools software [58]. The network diagram of protein interaction and regulation was drawn using the Cytoscape 22.0 software.

4.9. Statistical Analyses

Data analysis was performed using Excel (2020) and SPSS (Version 26, SPSS Inc., Chicago, IL, USA). One-way analysis of variance (ANOVA) was performed, followed by Duncan’s multiple range test at p < 0.05.

5. Conclusions

In this study, nine CsRSs were identified and analyzed. The analysis of phylogenetic relationships, conserved motifs, chromosome locations, gene structures, cis-elements and expression patterns showed that CsRSs are involved in the development process and stress response of tea plants. CsRSs are able to respond to continuous drought stress and cold acclimation. In addition, CsRSs participated in the accumulation of Raf during the withering process of white tea, especially CsRS1 and CsRS8, which had a significant correlation with the accumulation of Raf, and CsDBB might regulate the expression of CsRS8 during the withering process of white tea. This study can provide a reference for exploring the function of CsRSs and studying the quality formation mechanism during the withering process of white tea.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24010202/s1.

Author Contributions

Conceptualization, T.W. and Y.W.; methodology, T.W. and Y.W.; software, T.W., Y.W. and J.Z.; validation, T.W., Y.W., J.Z. and J.K.; formal analysis, T.W. and L.Z.; investigation, T.W., S.Q., J.C., Z.C. and W.S.; resources, W.Z. and W.S.; data curation, T.W., Y.W. and W.Z.; writing—original draft preparation, T.W., Y.W. and W.Z.; writing—review and editing, T.W., Y.W. and W.Z.; visualization, T.W. and Y.W.; supervision, W.Z. and W.S.; project administration, Z.C., W.Z. and W.S.; funding acquisition, W.Z. and W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (2019YFD1001601) and the Excavation, Preservation and Identification of Excellent Tea Plant Germplasm Resources in Fuding (BCY2021K01).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Quanjin Wu (the Open University of Fujian) for her help with the preparation of plant materials in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sonali, S.; Sritama, M.; Papri, B.; Majumder, A.L. Significance of galactinol and raffinose family oligosaccharide synthesis in plants. Front. Plant Sci. 2015, 6, 656. [Google Scholar]
  2. Liu, J.J.; Odegard, W.; LumenM, B.O. Galactinol synthase from kidney bean cotyledon and zucchini leaf. Purification and N-terminal sequences. Plant Physiol. 1995, 109, 505–511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Egert, A.; Keller, F.; Peters, S. Abiotic stress-induced accumulation of raffinose in Arabidopsis leaves is mediated by a single raffinose synthase (RS5, At5g40390). BMC Plant Biol. 2013, 13, 218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lehle, L.; Tanner, W. The Function of myo-Inositol in the Biosynthesis of Raffinose. Eur. J. Biochem. 1973, 38, 103–110. [Google Scholar] [CrossRef] [PubMed]
  5. Leopold, B.L.; LeopoldV, A.C. Changes in soluble carbohydrates during seed storage. Plant Physiol. 1992, 98, 1207–1210. [Google Scholar]
  6. Koster, K.L.; Leopold, A.C. Sugars and desiccation tolerance in seeds. Plant Physiol. 1988, 88, 829–832. [Google Scholar] [CrossRef] [Green Version]
  7. Smythe, B.M. Sucrose crystal growth. III. The relative growth rates of faces and their effect on sucrose crystal shape. Aust. J. Chem. 1967, 20, 1115–1131. [Google Scholar]
  8. Crowe, L.M.; Mouradian, R.; Crowe, J.H.; Jackson, S.A.; Womersley, C. Effects of carbohydrates on membrane stability at low water activities. Biochim. Biophys. Acta BBA Biomembr. 1984, 769, 141–150. [Google Scholar] [CrossRef]
  9. Mccarty, D.R. A core metabolome response of maize leaves subjected to long-duration abiotic stresses. Metabolites 2021, 11, 797. [Google Scholar]
  10. Rohde, P.; Hincha, D.K.; Heyer, A.G. Heterosis in the freezing tolerance of crosses between two Arabidopsis thaliana accessions (Columbia-0 and C24) that show differences in non-acclimated and acclimated freezing tolerance. Plant J. 2004, 38, 790–799. [Google Scholar] [CrossRef]
  11. Schneider, T.; Keller, F. Raffinose in chloroplasts is synthesized in the cytosol and transported across the chloroplast envelope. Plant Cell Physiol. 2009, 50, 2174–2182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ayako, N.; Yukinori, Y.; Shigeru, S. Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008, 147, 1251–1263. [Google Scholar]
  13. Tortuero, F.; Fernández, E.; Rupérez, P.; Moreno, M. Raffinose and lactic acid bacteria influence caecal fermentation and serum cholesterol in rats. Nutr. Res. 1997, 17, 41–49. [Google Scholar] [CrossRef] [Green Version]
  14. Collins, S.L.; Mcmillan, A.; Seney, S.; Veer, C.; Kort, R.; SumarahV, M.W.; Reid, K.; Mark, W.S.; Gregor, R. Promising prebiotic candidate established by evaluation of lactitol, lactulose, raffinose, and oligofructose for maintenance of a lactobacillus-dominated vaginal microbiota. Appl. Environ. Microb. 2017, 84, e2200–e2217. [Google Scholar] [CrossRef] [Green Version]
  15. Kim, H.S.; Cha, E.; Kim, Y.H.; Jeon, Y.H.; Olson, B.H.; Byun, Y.; ParkV, H.D. Raffinose, a plant galactoside, inhibits Pseudomonas aeruginosa biofilm formation via binding to LecA and decreasing cellular cyclic diguanylate levels. Sci. Rep. 2016, 6, 25318. [Google Scholar] [CrossRef]
  16. Zhou, M.L.; Zhang, Q.; Zhou, M.; Sun, Z.M.; Zhu, X.M.; Shao, J.R.; Tang, Y.X.; Wu, Y.M. Genome-wide identification of genes involved in raffinose metabolism in Maize. Glycobiology 2012, 22, 1775–1785. [Google Scholar]
  17. Kannan, U.; Sharma, R.; Gangola, M.P.; Ganeshan, S.; Båga, M.; Chibbar, R.N. Sequential expression of raffinose synthase and stachyose synthase corresponds to successive accumulation of raffinose, stachyose and verbascose in developing seeds of Lens culinaris Medik. J. Plant Physiol. 2021, 265, 153494. [Google Scholar] [CrossRef]
  18. Cui, R.F.; Wang, X.G.; Malik, W.A.; Lu, X.K.; Chen, X.G.; Wang, D.L.; Wang, J.J.; Wang, S.; Chen, C.; Guo, L.X.; et al. Genome-wide identification and expression analysis of Raffinose synthetase family in cotton. BMC Bioinform. 2021, 22, 356. [Google Scholar] [CrossRef]
  19. Zhou, J.; Yang, Y.; Yu, J.; Wang, L.K.; Yu, X.; Ohtani, M.; Kusano, M.; Saito, K.; Demura, T.; Zhuge, Q. Responses of Populus trichocarpa galactinol synthase genes to abiotic stresses. J. Plant Res. 2014, 127, 347–358. [Google Scholar] [CrossRef] [Green Version]
  20. Koning, R.D.; Kiekens, R.; Toili, M.E.M.; Angenon, G. Identification and expression analysis of the genes involved in the raffinose family oligosaccharides pathway of phaseolus vulgaris and glycine max. Plants 2021, 10, 1465. [Google Scholar] [CrossRef]
  21. Li, T.; Zhang, Y.M.; Wang, D.; Liu, Y.; Dirk, L.M.A.; Goodman, J.; Downie, A.B.; Wang, J.M.; Wang, G.Y.; Zhao, T.Y. Regulation of seed vigor by manipulation of raffinose family oligosaccharides in maize and Arabidopsis thaliana. Mol. Plant. 2017, 10, 1540–1555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jing, Y.; Lang, S.; Wang, D.M.; Xue, H.; Wang, X.F. Functional characterization of galactinol synthase and raffinose synthase in desiccation tolerance acquisition in developing Arabidopsis seeds. J. Plant Physiol. 2018, 230, 109–121. [Google Scholar] [CrossRef] [PubMed]
  23. Saito, M.; Yoshida, M. Expression analysis of the gene family associated with raffinose accumulation in rice seedlings under cold stress. J. Plant Physiol. 2011, 168, 2268–2271. [Google Scholar] [CrossRef] [PubMed]
  24. Guo, Q.Q.; Li, X.; Li, N.; Jameson, P.E.; Zhou, W.B. Transcription-associated metabolomic adjustments in maize occur during combined drought and cold stress. Plant Physiol. 2021, 186, 677–695. [Google Scholar] [CrossRef] [PubMed]
  25. Kito, K.; Yamane, K.; Yamamori, T.; Matsuhira, H.; Takabe, T. Isolation, functional characterization and stress responses of raffinose synthase genes in sugar beet. J. Plant Biochem. Biot. 2017, 27, 1–10. [Google Scholar] [CrossRef]
  26. Li, T.; Zhang, Y.; Liu, Y.; Li, X.; Hao, G.; Han, Q.; Dirk, L.; Downie, A.B.; Ruan, Y.L.; Wang, J.M.; et al. Raffinose synthase enhances drought tolerance through raffinose synthesis or galactinol hydrolysis in maize and Arabidopsis plants. J. Biol. Chem. 2020, 295, 8064–8077. [Google Scholar] [CrossRef]
  27. Gu, H.; Lu, M.; Zhang, Z.P.; Xu, J.J.; Cao, W.H.; Miao, M.M. Metabolic process of raffinose family oligosaccharides during cold stress and recovery in cucumber leaves. J. Plant Physiol. 2017, 224–225, 112–120. [Google Scholar] [CrossRef]
  28. Gu, L.; Jiang, T.; Zhang, C.X. Maize HSFA2 and HSBP2 antagonistically modulate raffinose biosynthesis and heat tolerance in Arabidopsis. Plant J. 2019, 100, 128–142. [Google Scholar] [CrossRef]
  29. Wu, X.L.; Kishitani, S.; Ito, Y.; Toriyama, K. Accumulation of raffinose in rice seedlings overexpressing OsWRKY11 in relation to desiccation tolerance. Plant Bio. 2018, 26, 431–434. [Google Scholar] [CrossRef] [Green Version]
  30. Han, Q.H.; Qi, J.L.; Hao, G.L.; Zhang, C.X.; Wang, C.M.; Dirk, L.M.A.; Bruce, D.A.; Zhao, T.Y. ZmDREB1A regulates raffinose synthase controlling raffinose accumulation and plant chilling stress tolerance in maize. Plant Cell Physiol. 2019, 61, 331–341. [Google Scholar] [CrossRef]
  31. Khan, M.; Hu, J.; Dahro, B.; Ming, R.; Zhang, Y.; Wang, Y.; Alhag, A.; Li, C.; Liu, J.H. ERF108 from Poncirus trifoliata (L.) Raf. functions in cold tolerance by modulating raffinose synthesis through transcriptional regulation of PtrRafS. Plant J. 2021, 108, 705–724. [Google Scholar] [CrossRef] [PubMed]
  32. Zou, L.; Shen, S.S.; Wei, Y.M.; Jia, H.Y.; Li, T.H.; Yin, X.C.; Lu, C.Y.; Cui, Q.Q.; He, F.; Deng, W.W.; et al. Evaluation of the effects of solar withering on nonvolatile compounds in white tea through metabolomics and transcriptomics. Food Res. Int. 2022, 162, 112088. [Google Scholar] [CrossRef]
  33. Zhang, L.; Ho, C.T.; Zhou, J.; Santos, J.S.; Armstrong, L.; Granato, D. Chemistry and Biological Activities of Processed Camellia sinensis Teas: A Comprehensive Review. Compr. Rev. Food Sci. Food Saf. 2019, 18, 1474–1495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dias, T.R.; Tomás, G.; Teixeira, N.F.; Alves, M.G.; Silva, B.M. White tea (Camellia Sinensis (L.)): Antioxidant properties and beneficial health effects. Int. J. Food Sci. Nutr. 2013, 2, 19–26. [Google Scholar]
  35. Deng, X.M.; Shang, H.; Chen, J.J.; Wu, J.; Wang, T.; Wang, Y.Q.; Zhu, C.S.; Sun, W.J. Metabolomics combined with proteomics provide a novel interpretation of the changes in flavonoid glycosides during white tea processing. Foods 2022, 11, 1226. [Google Scholar] [CrossRef] [PubMed]
  36. Sharma, E.; Joshi, R.; Gulati, A. L-Theanine: An astounding sui generis integrant in tea. Food Chem. 2018, 242, 601–610. [Google Scholar] [CrossRef] [PubMed]
  37. Zhou, C.Z.; Zhu, C.; Li, X.Z.; Chen, L.; Xie, S.Y.; Chen, G.W.; Zhang, H.; Lai, Z.X.; Lin, Y.L.; Guo, Y.Q. Transcriptome and phytochemical analyses reveal the roles of characteristic metabolites in the taste formation of white tea during the withering process. J. Integr. Agric. 2022, 21, 862–877. [Google Scholar] [CrossRef]
  38. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  39. Xia, E.H.; Li, F.D.; Tong, W.; Li, P.H.; Wu, Q.; Zhao, H.J.; Ge, R.H.; Li, R.P.; Li, Y.Y.; Zhang, Z.Z.; et al. Tea Plant Information Archive: A comprehensive genomics and bioinformatics platform for tea plant. Plant Biotechnol. J. 2019, 17, 1938–1953. [Google Scholar] [CrossRef]
  40. Zhang, Q.; Cai, M.C.; Yu, X.M.; Wang, L.S.; Guo, C.F.; Ming, R.; Zhang, J.S. Transcriptome dynamics of Camellia sinensis in response to continuous salinity and drought stress. Tree Genet. Genomes 2017, 13, 78. [Google Scholar] [CrossRef]
  41. Taji, T.; Ohsumi, C.; Iuchi, S.; Seki, M.; Kasuga, M.; Kobayashi, M.; Yamaguchi, S.K.; Shinozaki, K. Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 2002, 29, 417–426. [Google Scholar] [CrossRef] [PubMed]
  42. Gang, R.; Behmüller, R.; Tenhaken, R. Molecular cloning of AtRS4, a seed specific multifunctional RFO synthase/galactosylhydrolase in Arabidopsis thaliana. Front. Plant Sci. 2015, 6, 789. [Google Scholar]
  43. Li, R.M.; Yuan, S.; He, Y.D.; Fan, J.; Zhou, Y.J.; Qiu, T.T.; Lin, X.J.; Yao, Y.; Liu, J.; Fu, S.P.; et al. Genome-wide identification and expression profiling analysis of the galactinol synthase gene family in Cassava (Manihot esculenta Crantz). Agron 2018, 8, 250. [Google Scholar] [CrossRef] [Green Version]
  44. You, J.; Wang, Y.Y.; Zhang, Y.J.; Dossa, K.; Li, D.H.; Zhou, R.; Wang, L.H.; Zhang, X.R. Genome-wide identification and expression analyses of genes involved in raffinose accumulation in sesame. Sci. Rep. 2018, 8, 4331. [Google Scholar] [CrossRef] [Green Version]
  45. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef] [Green Version]
  46. Hernandez-Garcia, C.M.; Finer, J.J. Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 2014, 217–218, 109–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Sui, X.L.; Meng, F.Z.; Wang, H.Y.; Wei, Y.X.; Li, R.F.; Wang, Z.Y.; Hu, L.P.; Wang, S.H.; Zhang, Z.X. Molecular cloning, characteristics and low temperature response of raffinose synthase gene in Cucumis sativus L. J. Plant Physiol. 2012, 169, 1883–1891. [Google Scholar] [CrossRef]
  48. Yue, C.; Cao, H.L.; Wang, L.; Zhou, Y.H.; Huang, Y.T.; Hao, X.Y.; Wang, Y.C.; Wang, B.; Yang, Y.J.; Wang, X.C. Effects of cold acclimation on sugar metabolism and sugar-related gene expression in tea plant during the winter season. Plant Mol. Biol. 2015, 88, 591–608. [Google Scholar] [CrossRef]
  49. Singh, K.B.; Foley, R.C.; Oñate, S.L. Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol. 2002, 5, 430–436. [Google Scholar] [CrossRef]
  50. Lang, S.R.; Liu, X.X.; Xue, H.; Li, X.; Wang, X.F. Functional characterization of BnHSFA4a as a heat shock transcription factor in controlling the re-establishment of desiccation tolerance in seeds. J. Exp. Bot. 2017, 68, 2361–2375. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, J.H.; Zhao, T.Y. ZmDREB2A regulates ZmGH3.2 and ZmRAFS, shifting metabolism towards seed aging tolerance over seedling growth. Plant J. 2020, 104, 268–282. [Google Scholar]
  52. Wei, W.; Zhang, Y.X.; Han, L.; Guan, Z.Q.; Chai, T.Y. A novel WRKY transcriptional factor from Thlaspi caerulescens negatively regulates the osmotic stress tolerance of transgenic tobacco. Plant Cell Rep. 2008, 27, 795–803. [Google Scholar] [CrossRef] [PubMed]
  53. Khanna, R.; Kronmiller, B.; Maszle, D.R.; Coupland, G.; Holm, M.; Mizuno, T.; Wu, S.H. The Arabidopsis B-box zinc finger family. Plant Cell 2009, 21, 3416–3420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Xia, E.H.; Tong, W.; Hou, Y.; An, Y.L.; Chen, L.B.; Wu, Q.; Liu, Y.L.; Yu, J.; Li, F.D.; Li, R.P.; et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into its genome evolution and adaptation. Mol. Plant 2020, 13, 1013–1026. [Google Scholar] [CrossRef] [PubMed]
  55. Jaina, M.; Sara, C.; Lowri, W.; Matloob, Q.; Gustavoa, S.; Sonnhammer, E.; Tosatto, S.; Lisanna, P.; Shriya, R.; Richardson, L.J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2020, 49, D412–D419. [Google Scholar]
  56. Xu, D.; Lu, Z.C.; Jin, K.M.; Qiu, W.M.; Qiao, G.R.; Han, X.J.; Zhuo, R.Y. Data and text mining SPDE: A multi-functional software for sequence processing and data extraction. Bioinformatics 2021, 37, 3686–3687. [Google Scholar] [CrossRef]
  57. Lu, S.; Wang, J.; Farideh, C.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; Marc, G.; Hurwitz, D.I.; Marchler, G.H.; Song, J.S.; et al. CDD/SPARCLE: The conserved domain database. Nucleic Acids Res. 2019, 48, D265–D268. [Google Scholar] [CrossRef] [Green Version]
  58. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools-an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  59. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. The Proteomics Protocols Handbook-Chapter 52: Protein Identification and Analysis Tools on the ExPASy Server. Proteom. Protoc. Handb. 2005, 571–607. [Google Scholar] [CrossRef]
  60. Xiong, E.H.; Zheng, C.Y.; Wu, X.L.; Wang, W. Protein Subcellular Location: The gap between prediction and experimentation. Plant Mol. Biol. Rep. 2015, 34, 1–10. [Google Scholar] [CrossRef]
  61. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  62. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar]
  63. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Cryst. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  64. Laskowski, R.A.; Thornton, J.M. PDBsum extras: SARS-CoV-2 and AlphaFold models. Protein Sci. 2022, 31, 283–289. [Google Scholar] [CrossRef] [PubMed]
  65. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef] [PubMed]
  66. Ma, S.; Sun, L.L.; Sui, X.L.; Li, Y.X.; Chang, Y.W.; Fan, J.; Zhang, Z.X. Phloem loading in cucumber: Combined symplastic and apoplastic strategies. Plant J. 2019, 98, 391–404. [Google Scholar] [CrossRef] [PubMed]
  67. Livak, K.J.; Schmittgen, T.J.M. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCt method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Ijms 24 00202 g001 550
Figure 1. The phylogenetic tree, distribution and duplication events of CsRs. (A) Phylogenetic tree of RS proteins constructed using the neighbor-joining(NJ) phylogenetic tree by MEGA X. In total, 57 RS protein sequences from different species were used to construct the phylogenetic tree by MEGA X. The seven different groups are indicated by different colors. Am: Alonsoa meridionalis, Cm: Cucumis melo, Cs: Camellia sinensis, Os: Oryza sativa, Ps: Pisum sativum, Pt: Populus trichocarpa, Rc: Ricinus communis, Si: Sesamum indicum, Va: Vigna angularis, Vv: Vitis vinifera, Zm: Zea mays. The red star represents the RSs in Camellis sinensis. (B) Chromosomal localization and gene duplication events in CsRSs. From outside to inside, the first circle is the chromosome number with the scale on top indicating the coordinate position of the chromosome, and the second circle is the chromosome number and density. Red line indicates gene duplication.
Figure 1. The phylogenetic tree, distribution and duplication events of CsRs. (A) Phylogenetic tree of RS proteins constructed using the neighbor-joining(NJ) phylogenetic tree by MEGA X. In total, 57 RS protein sequences from different species were used to construct the phylogenetic tree by MEGA X. The seven different groups are indicated by different colors. Am: Alonsoa meridionalis, Cm: Cucumis melo, Cs: Camellia sinensis, Os: Oryza sativa, Ps: Pisum sativum, Pt: Populus trichocarpa, Rc: Ricinus communis, Si: Sesamum indicum, Va: Vigna angularis, Vv: Vitis vinifera, Zm: Zea mays. The red star represents the RSs in Camellis sinensis. (B) Chromosomal localization and gene duplication events in CsRSs. From outside to inside, the first circle is the chromosome number with the scale on top indicating the coordinate position of the chromosome, and the second circle is the chromosome number and density. Red line indicates gene duplication.
Ijms 24 00202 g001
Ijms 24 00202 g002 550
Figure 2. The molecular structure and cis-acting element analyses of CsRSs. (A) Phylogenetic analysis of RS genes in Camellia sinensis. CsRSs are distributed in the three subgroups on the left. (B) Conserved motif of CsRSs. 10 motifs are predicted by MEME and different colored boxes represent different motifs. The box lengths represent motif lengths. (C) Conserved domain structure of CsRSs. The green box represents the RS domain. (D) Exon-intron structures of CsRSs. The green boxes, yellow boxes and black lines represent the UTRs, exons and introns, respectively. (E) Heat map of cis-acting element of CsRSs. Different color blocks represent different promoter types.
Figure 2. The molecular structure and cis-acting element analyses of CsRSs. (A) Phylogenetic analysis of RS genes in Camellia sinensis. CsRSs are distributed in the three subgroups on the left. (B) Conserved motif of CsRSs. 10 motifs are predicted by MEME and different colored boxes represent different motifs. The box lengths represent motif lengths. (C) Conserved domain structure of CsRSs. The green box represents the RS domain. (D) Exon-intron structures of CsRSs. The green boxes, yellow boxes and black lines represent the UTRs, exons and introns, respectively. (E) Heat map of cis-acting element of CsRSs. Different color blocks represent different promoter types.
Ijms 24 00202 g002
Ijms 24 00202 g003 550
Figure 3. The heat map of CsRSs expression mode. (A) Expression pattern of CsRSs in different parts of the tea plant. (B) Expression pattern of CsRSs during continuous drought stress and cold acclimation. The FPKM of the continuous drought stress and cold acclimation converted by log2. PEG-CK: untreated, PEG-24: PEG treatment for 24 h, PEG-48: PEG treatment for 48 h, PEG-72: PEG treatment for 72 h, Cold-CK: nonacclimated at 25~20 °C, CA1-6 h: fully acclimated at 10 °C for 6 h, CA1-7 d: fully acclimated at 10~4 °C for 7 days, DA-7 d: recovering under 25~20 °C for 7 days. (C) Expression pattern of CsRSs in the withering process of white tea.
Figure 3. The heat map of CsRSs expression mode. (A) Expression pattern of CsRSs in different parts of the tea plant. (B) Expression pattern of CsRSs during continuous drought stress and cold acclimation. The FPKM of the continuous drought stress and cold acclimation converted by log2. PEG-CK: untreated, PEG-24: PEG treatment for 24 h, PEG-48: PEG treatment for 48 h, PEG-72: PEG treatment for 72 h, Cold-CK: nonacclimated at 25~20 °C, CA1-6 h: fully acclimated at 10 °C for 6 h, CA1-7 d: fully acclimated at 10~4 °C for 7 days, DA-7 d: recovering under 25~20 °C for 7 days. (C) Expression pattern of CsRSs in the withering process of white tea.
Ijms 24 00202 g003
Ijms 24 00202 g004 550
Figure 4. Raffinose content and CsDBB expression pattern. (A) Changes of Raf content in white tea during withering. (B) Expression pattern of CsDBB in different parts of the tea plant. (C) Expression pattern of CsDBB in the withering process of white tea. Different letters represent significant differences (p < 0.05).
Figure 4. Raffinose content and CsDBB expression pattern. (A) Changes of Raf content in white tea during withering. (B) Expression pattern of CsDBB in different parts of the tea plant. (C) Expression pattern of CsDBB in the withering process of white tea. Different letters represent significant differences (p < 0.05).
Ijms 24 00202 g004
Ijms 24 00202 g005 550
Figure 5. Transcription factor binding site prediction and protein interaction regulation network of CsRSs. (A) Transcription factor binding sites predicted in the promoters of CsRSs. Protein Interaction Regulatory Network of CsRSs. Boxes of different colors represent different transcription factor families. “+” and “−” represent positive and negative strands, respectively. Dotted lines indicate protein interactions, solid lines indicate that transcription factors regulate CsRSs, and the solid red line indicates that it has been verified by experiments. (B) Protein interactions and regulatory networks of CsRSs. The outermost ring represents RS in Camellia sinensis. The circles represent functional genes and the diamonds represent transcription factors, different color blocks represent different families of transcription factors, the dotted line represents the interaction between genes, the solid line represents the regulation of genes by transcription factors, and the red line represents preliminary validation.
Figure 5. Transcription factor binding site prediction and protein interaction regulation network of CsRSs. (A) Transcription factor binding sites predicted in the promoters of CsRSs. Protein Interaction Regulatory Network of CsRSs. Boxes of different colors represent different transcription factor families. “+” and “−” represent positive and negative strands, respectively. Dotted lines indicate protein interactions, solid lines indicate that transcription factors regulate CsRSs, and the solid red line indicates that it has been verified by experiments. (B) Protein interactions and regulatory networks of CsRSs. The outermost ring represents RS in Camellia sinensis. The circles represent functional genes and the diamonds represent transcription factors, different color blocks represent different families of transcription factors, the dotted line represents the interaction between genes, the solid line represents the regulation of genes by transcription factors, and the red line represents preliminary validation.
Ijms 24 00202 g005
Table
Table 1. Secondary structure of CsRSs protein.
Table 1. Secondary structure of CsRSs protein.
Gene IDAlpha HelixExtended StrandBeta TurnRandom Coil
CsRS128.33%21.73%6.86%43.08%
CsRS232.77%16.85%6.34%44.05%
CsRS327.87%21.53%7%43.59%
CsRS428.92%19.62%6.25%45.20%
CsRS530.68%18.41%5.44%45.47%
CsRS626.99%21.68%7.31%44.02%
CsRS728.36%21.84%8.52%41.28%
CsRS827.89%20.62%6.23%45.27%
CsRS929.73%21.62%7%41.65%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, T.; Wang, Y.; Zhao, J.; Kong, J.; Zhang, L.; Qi, S.; Chen, J.; Chen, Z.; Zeng, W.; Sun, W. Identification, Characterization and Expression Profiling of the RS Gene Family during the Withering Process of White Tea in the Tea Plant (Camellia sinensis) Reveal the Transcriptional Regulation of CsRS8. Int. J. Mol. Sci. 2023, 24, 202. https://doi.org/10.3390/ijms24010202

AMA Style

Wang T, Wang Y, Zhao J, Kong J, Zhang L, Qi S, Chen J, Chen Z, Zeng W, Sun W. Identification, Characterization and Expression Profiling of the RS Gene Family during the Withering Process of White Tea in the Tea Plant (Camellia sinensis) Reveal the Transcriptional Regulation of CsRS8. International Journal of Molecular Sciences. 2023; 24(1):202. https://doi.org/10.3390/ijms24010202

Chicago/Turabian Style

Wang, Tao, Yiqing Wang, Jiamin Zhao, Jiumei Kong, Lingzhi Zhang, Siyu Qi, Jiajia Chen, Zhidan Chen, Wen Zeng, and Weijiang Sun. 2023. "Identification, Characterization and Expression Profiling of the RS Gene Family during the Withering Process of White Tea in the Tea Plant (Camellia sinensis) Reveal the Transcriptional Regulation of CsRS8" International Journal of Molecular Sciences 24, no. 1: 202. https://doi.org/10.3390/ijms24010202

APA Style

Wang, T., Wang, Y., Zhao, J., Kong, J., Zhang, L., Qi, S., Chen, J., Chen, Z., Zeng, W., & Sun, W. (2023). Identification, Characterization and Expression Profiling of the RS Gene Family during the Withering Process of White Tea in the Tea Plant (Camellia sinensis) Reveal the Transcriptional Regulation of CsRS8. International Journal of Molecular Sciences, 24(1), 202. https://doi.org/10.3390/ijms24010202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop