Next Article in Journal
Serum and Hepatic Autofluorescence as a Real-Time Diagnostic Tool for Early Cholestasis Assessment
Previous Article in Journal
Grape Seed Proanthocyanidins Improve White Adipose Tissue Expansion during Diet-Induced Obesity Development in Rats
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 9
10.3390/ijms19092633
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

Genome-Wide Identification, Classification and Expression Analysis of the HSP Gene Superfamily in Tea Plant (Camellia sinensis)

by
Jiangfei Chen
,
Tong Gao
,
Siqing Wan
,
Yongheng Zhang
,
Jiankun Yang
,
Youben Yu
and
Weidong Wang
*
College of Horticulture, Northwest A&F University, Yangling 712100, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(9), 2633; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092633
Submission received: 7 August 2018 / Revised: 28 August 2018 / Accepted: 3 September 2018 / Published: 5 September 2018
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Heat shock proteins (HSPs) function as molecular chaperones. These proteins are encoded by a multigene family whose members play crucial roles in plant growth, development and stress response. However, little is known about the HSP gene superfamily in tea plant. In this study, a total of 47 CsHSP genes were identified, including 7 CsHSP90, 18 CsHSP70, and 22 CssHSP genes. Phylogenetic and composition analyses showed that CsHSP proteins in the same subfamily have similar gene structures and conserved motifs, but significant differences exist in the different subfamilies. In addition, expression analysis revealed that almost all CsHSP genes were specifically expressed in one or more tissues, and significantly induced under heat and drought stress, implying that CsHSP genes play important roles in tea plant growth, development, and response to heat and drought stress. Furthermore, a potential interaction network dominated by CsHSPs, including HSP70/HSP90 organizing protein (HOP) and heat shock transcription factor (HSF), is closely related to the abovementioned processes. These results increase our understanding of CsHSP genes and their roles in tea plant, and thus, this study could contribute to the cloning and functional analysis of CsHSP genes and their encoded proteins in the future.

Graphical Abstract

1. Introduction

Plants in the growth stage often have difficulty avoiding threats from various environmental stresses, such as drought, cold, heat, salinity, and heavy metal stress [1,2]. Generally, plants are unable to change their location to escape the abovementioned stresses, but have instead developed various defense mechanisms, including activating the synthesis of heat shock proteins (HSPs), to resist them [3,4]. HSPs are ubiquitous proteins found in plant cells, which were originally detected in plant response to heat stress but are now known to be induced by various biotic and abiotic stresses [5]. HSPs, as molecular chaperones involved in the cellular processes of protein folding and re-folding, assembly, translocation, and degradation, help to maintain the stability of proteins and membranes under stress conditions [6,7,8]. Plant HSPs are divided into five conserved families according to their molecular weights, namely, the HSP100, HSP90, HSP70, HSP60, and small HSP (sHSP) families [9,10]. Among them, the members of the HSP90, HSP70, and sHSP families are the most abundant, and the structures and biological functions of these proteins have received increasing attention [10,11]. Previous studies have confirmed a very high degree of conservation in the structures of HSPs in these three families. HSP90 consists of an N-terminal ATP-binding domain, a C-terminal dimerization domain, and a middle domain that contains an amphipathic loop [12], and HSP70 also contains three distinct domains: an N-terminal ATPase domain, a substrate binding domain, and a highly variable C-terminal domain, and is considered to be the most highly conserved HSP [13]. Notably, HSP110 proteins, a class of macromolecular HSPs, are also included in the HSP70 family, due to their high sequence and structural homology to Hsp70 proteins [14]. In comparison, sHSPs share only a conserved carboxyl-terminal domain of approximately 90 amino acid residues called the a-crystallin domain (ACD), which is characterized by a compact β-strand structure and contains two conserved regions (CRs): CRI with β2, β3, β4, and β5; and CRII with β7, β8, and β9, and a β6 loop [15].
At present, the functional mechanisms of HSP90, HSP70, and sHSP in plant stress response have been widely reported. For instance, the expression of Arabidopsis HSP90 genes is strongly induced under heat, cold, salinity, and heavy metal stress [16,17]; the overexpression of alfalfa MsHSP70 and horsegram MuHSP70 improved the tolerance of transgenic plants to multiple stresses, such as heat, cold, drought, salinity, and oxidative stress [18,19]. In addition, existing studies have confirmed that HSP70 and HSP90 proteins can interact with each other, depending on co-chaperone molecules, such as HOP (HSP70/HSP90 organizing protein), ROF (FK506-binding protein), and HIP (HSP70 interacting protein), to form a complex to jointly respond to external stress [20,21,22]. For example, AtHOP3 interacted with HSP90 and HSP70 to regulate the response of Arabidopsis to different biotic and abiotic stresses [23]. Similarly, the roles of sHSP in plant response to stress have been generally confirmed; for example, rice OsHSP26 could significantly enhance the tolerance of tall fescue to oxidative and heat stresses by protecting photosystem II (PSII) to maintain photosynthesis [24], and a substantial and durable increase in thermotolerance was achieved through sHSP overexpression in hybrid poplar [25]. Furthermore, HSP activity is widely believed to be controlled and regulated by specific types of transcription factors called heat shock transcription factors (HSFs), which not only regulate the expression of HSP genes, but also participate in the protein interactions among HSPs [26,27]. To date, the HSP gene families have been investigated in several plant species, such as Arabidopsis [17], rice [28], soybean [29], populus [30], and grapes [31], and the functional mechanisms of HSPs in plant stress response have gradually become a hot topic of study.
Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops in the world, its young leaves are processed to prepare a popular non-alcoholic beverage known as “tea”, which has positive effects such as delaying ageing and preventing cancer and cardiovascular diseases [32,33]. Tea plant is often affected by harsh environments, and heat and drought stress are important factors that significantly constrain the yield and quality of tea products [34]. Therefore, clarifying the mechanisms of tea plant response to heat and drought stress and identifying resistance genes are of great significance. Recently, four sHSP genes in tea plant were reported to play important roles in the process of tea plant response to heat stress [35,36]. However, other studies of HSPs in tea plant are extremely rare. The genome of tea plant has been sequenced and published, rendering the identification and characterization of the CsHSP gene superfamily easy and reliable [37]. In this study, 47 CsHSP genes were identified, and a comprehensive analysis was performed, including sequence characteristics, phylogenetic relationships, gene structures, conserved motifs, and functional interaction network analysis. In addition, the expression patterns of the CsHSP genes were detected in different tissues of tea plants and under heat and drought stress. The results of this study reveal the molecular characteristics of the CsHSP gene superfamily and provide a theoretical basis for future studies of the biological functions of CsHSPs under abiotic stresses.

2. Results

2.1. Identification and Characterization of the CsHSP Gene Superfamily in Tea Plant

A total of 47 CsHSP genes were isolated from the tea plant genome and transcriptome database, and their corresponding amino acid sequences were confirmed using the NCBI BLAST program (Version 2.4.0) (Table S1). Among them, 7 (14.89%), 18 (38.29%) and 22 (46.82%) members were determined to belong to the HSP90, HSP70, and sHSP families, respectively. The results of physiological and biochemical properties analysis showed that the lengths of the CsHSP90 proteins ranged from 643 (CsHSP90-7) to 823 (CsHSP90-2) amino acids, the molecular weights were between 73.68 kDa (CsHSP90-7) and 94.32 kDa (CsHSP90-2), and the isoelectric point (pI) values ranged from 4.84 (CsHSP90-2) to 5.20 (CsHSP90-1); the lengths of the CsHSP70 proteins ranged from 407 (CsHSP70-17) to 892 (CsHSP70-3) amino acids, the molecular weights were between 45.45 kDa (CsHSP70-17) and 98.88 kDa (CsHSP70-3), and the pI values ranged from 5.03 (CsHSP70-7) to 9.04 (CsHSP70-17); the lengths of the CssHSPs ranged from 136 (CssHSP-1) to 255 (CssHSP-5) amino acids, the molecular weights were between 15.45 kDa (CssHSP-1) and 28.07 kDa (CssHSP-5), and the pI values ranged from 4.83 (CssHSP-15) to 9.81 (CssHSP-11). More detailed information, including the grand average of hydropathicity (GRAVY), instability index, aliphatic index and subcellular localization, is listed in Table 1.
Twenty other species in which HSP superfamily genes have been identified were comparatively analyzed to further evaluate the family composition of the CsHSPs, and the results showed that the distribution of the CsHSP members maintained a similar proportion to those of other higher plants, namely, the HSP90 family had significantly fewer members than the HSP70 and sHSP families did. However, the above trend is not evident in lower plants, which have very few sHSP members. In addition, the total number of HSPs in lower plants is significantly less than that in higher plants. In animals, almost equal numbers of genes were found in the three HSP families, showing a distribution clearly different from those of higher plants (Figure 1).

2.2. Phylogenetic Analysis of CsHSPs

By clustering the HSP90 proteins of Arabidopsis, rice, tomato, soybean, and tea plant, we found that the CsHSP90 proteins were divided into 3 groups. Group I included 4 members (CsHSP90-4, -5, -6, -7), while group II and group III contained 1 member (CsHSP90-2) and 2 members (CsHSP90-1 and CsHSP90-3), respectively (Figure 2A). Similarly, the CsHSP70 proteins were classified into 5 subfamilies, including 9 in the cytoplasm (CsHSP70-1, -8, -9, -10, -12, -13, -14, -16, and -17), 2 in the ER (CsHSP70-7 and CsHSP70-11), 2 in the mitochondria (CsHSP70-6 and CsHSP70-15), 2 in the plastid (CsHSP70-4 and CsHSP70-5), and 3 in the HSP110/SSE (CsHSP70-2, -3 and-18) subfamily (Figure 2B). According to the classification of sHSPs in Arabidopsis [38], the 22 CssHSPs were classified into only 9 protein subfamilies, including CI (9), CII (2), CV (1), CVI (4), CVII (1), ER (1), MI (1), P (3), and Px (1), while no CssHSPs existed in the CIII, CIV, and MII subfamilies (Figure 2C). Interestingly, the P subfamily clustered into two branches, and was closely related to the MI and MII subfamilies, suggesting that they have similar origins and evolutionary relationships.

2.3. Gene Structures of the CsHSP Gene Superfamily

To further analyze the evolutionary relationships of the CsHSP gene superfamily, exon–intron diagrams of the CsHSP90, CsHSP70, and CssHSP genes were generated according to their genome and coding sequences. Group II and III CsHSP90 genes are characterized by a large number of introns (14–19), while group I CsHSP90 genes contain only a few introns (2–3), except for CsHSP90-4 that has 11 introns (Figure 3A). In the CsHSP70 genes, there are 0–3 introns in the cytoplasm subfamily genes, 3–5 introns in the mitochondria subfamily genes, and at least 6 introns in the ER, P, and HSP110/SSE subfamily genes (Figure 3B). In comparison, most CssHSP genes have no introns or one intron; only CssHSP-5, in the P subfamily, contains 4 introns (Figure 3C).

2.4. Conserved Protein Motif Analysis of CsHSPs

To investigate the protein sequence features of the CsHSPs, 10 motifs in the CsHSP90, CsHSP70, and CssHSP families were individually predicted by the MEME tool, and the regular expression levels of the conserved motifs are listed in Table S3. In the CsHSP90 family, all members have a similar motif composition, although small differences between members in different group are also present (Figure 4A). In the CsHSP70 and CssHSP families, members belonging to the same subfamily contain similar motifs, while the type and number of motifs in the different subfamilies of CsHSP70s and CssHSPs differ significantly (Figure 4B,C). In addition, protein motifs 1, 2, and 3 in the CssHSPs are highly conserved and distributed across almost all members; together, they form the ACD protein domain of the sHSP (Figure S1), which is the structural basis for the biological function of sHSPs.

2.5. Protein Interaction Network of CsHSPs

In this study, the functional and physical protein interactions of the CsHSPs were validated by constructing an Arabidopsis association model using STRING software (version 10.5, https://string-db.org/). The results showed that 47 CsHSPs associated with 26 known Arabidopsis HSPs participated in the interaction network (Table S4 and Figure 5). As expected, the network showed general and complex interactions between HSP70 and HSP90 family proteins; for example, HSP70 (corresponding to CsHSP70-10, -12, -16, and -17), HSP70T-2 (corresponding to CsHSP70-5), HSC70-1 (corresponding to CsHSP70-1, -8, -9, -13, and -14) and MTHSC70-2 (corresponding to CsHSP70-6 and -15) closely interact with HSP81.4 (corresponding to CsHSP90-4 and -5), SHD (corresponding to CsHSP90-2), and HSP90.1 (corresponding to CsHSP90-6 and -7). Simultaneously, HOP3, an HSP70/HSP90 organizing protein, was closely connected to many HSP70 and HSP90 family proteins, implying that it plays important roles in their functions. In addition, 3 HSP110/SSE subfamily proteins (corresponding to CsHSP70-2, -3, and -18) showed very positive interactions with other HSP70 family proteins, especially the members of the cytoplasm subfamily. Furthermore, an extensive interaction was observed between the transcription factor HSFA2 and the HSP70 and HSP90 proteins. By contrast, some sHSPs (corresponding to CssHSP-2, -3, -6, 10, -11, -12, -13, -17, -20, and -21) were associated with only HSP70T-2 proteins, and others were independent.

2.6. Expression Patterns of CsHSP Genes in Different Tissues of Tea Plant

To elucidate the tissue-specific expression patterns of CsHSP genes, qRT-PCR was used to determine the expression levels of 47 CsHSP genes in four tissues, including the root, stem, leaf, and flower of tea plant (Figure 6 and Table S5). All CsHSP90 genes were highly expressed in the leaf, and CsHSP90-2, -4, and -5 were also highly expressed in the stem (Figure 6A). Most CsHSP70 genes are differentially expressed in different tissues, although they have no obvious regularity (Figure 6B). In addition, we found that 7 CssHSP genes had higher expression levels in the stem, 10 CssHSP genes had higher expression levels in the leaf, and 5 CssHSP genes had higher expression levels in the flower, although some of these genes are also significantly expressed in other tissues (Figure 6C). In general, almost all CsHSP genes were specifically expressed in one or more tissues, implying that these genes play different roles in the growth and development of tea plant.

2.7. Expression Profiles of CsHSP Genes in Response to Heat and Drought Stress

As shown in Figure 7 and Table S6, the expression of all CsHSP90 genes was significantly upregulated under heat stress, although the expression levels of certain genes briefly and slightly decreased at individual time points (Figure 7A). Among the CsHSP70 genes, 12 members showed highly upregulated expression throughout the heat treatment process, 5 members were upregulated in the early stage and downregulated in the late stage, and the expression of the CsHSP70-4 gene was significantly suppressed (Figure 7B). Similarly, the expression of 18 CssHSP genes was significantly upregulated and maintained at a very high level, while the remaining 4 CssHSP genes, including CssHSP-1, -2, -11, and -12, were significantly downregulated, although the expression levels of certain genes briefly increased at 2 h or 24 h (Figure 7C). By contrast, the expression trends of the CsHSP90, CsHSP70, and CssHSP genes under drought stress were complex and diverse, although most of them also showed upregulated expression (Figure 8 and Table S7). Specifically, most CsHSP90 genes were downregulated at 2 and 4 h, and significantly upregulated at other time points (Figure 8A). Among the CsHSP70 genes, the expression levels of 5 members showed a trend of increasing first and then decreasing, and other members were upregulated overall, although the change trends of certain genes were weak at 2 to 24 h (Figure 8B). All CssHSP genes, except CssHSP-12, were upregulated under drought stress, and the expression of most of them also showed a trend of increasing first and then decreasing (Figure 8C). In addition, the expression trends of most CsHSP genes were identical under both drought and heat stress, but the opposite was observed in a few genes, including CsHSP90-2, CsHSP90-4, CsHSP70-4, CssHSP-2, and CssHSP-11 (Figure 7 and Figure 8). These results reveal that almost all CsHSP genes are involved in the responses of tea plant to heat and drought stress, and the response mechanisms of different CsHSP genes to different stresses are complex and diverse.

3. Discussion

HSPs are ubiquitous molecular chaperones that play vital roles during plant growth and development, and protect the plant cellular machinery under stress conditions [5,6,10]. The biological function of HSPs has been widely studied in many plants, such as Arabidopsis [14], rice [28], and wheat [39]. However, only a few efforts have been made to elucidate HSPs in tea plant. In this study, a comprehensive genome-wide analysis of the CsHSP gene superfamily in tea plant was carried out, and the results will provide a powerful theoretical foundation for future functional studies.
In the current study, 47 CsHSP genes were identified in tea plant, including 7 CsHSP90, 18 CsHSP70, and 22 sHSP genes, demonstrating a similar proportion to the distribution of HSP members in other higher plants. Notably, the number of HSP genes is much greater in higher plants than in lower plants, which may be because higher plants generally require more genes for biological processes [40]. However, HSP gene numbers do not correspond with evolutionary level in higher plants, implying diversity in plant HSP gene evolution [41]. It has been reported that the phylogenetic relationships of HSPs are closely correlated with their subcellular localizations, and subfamilies are usually named for their protein localization in this group [9,15]. Our phylogenetic analysis of the CsHSP90, CsHSP70, and CssHSP families showed results highly consistent with the predictions of protein localization, which further confirms the above conclusions. Notably, three CsHSP70 proteins (CsHSP70-2, -3, and -18) were clustered into the HSP110/SSE subfamily, which generally act as nucleotide exchange factors for cytoplasmic HSP70 proteins, and participate in HSP70-mediated protein folding [42], suggesting that CsHSP70-2, -3, and -18 may have similar biological functions. This hypothesis has been confirmed to some extent by the results of the protein interaction network prediction, which showed positive interactions of CsHSP70-2, -3, and -18 with cytoplasmic CsHSP70 proteins. In addition, the 22 CssHSPs were grouped into only 9 subfamilies, fewer than the 12 subfamilies in Arabidopsis and similar to the sHSP classifications in pepper [43], tomato [44], and rice [26]. Intriguingly, the CssHSPs in the P and M (MI and MII) subfamilies were closely related to each other, which was consistent with the sHSPs of the M subfamily having evolved later from the P subfamily [39,43,45].
Previous studies have shown that gene organization plays an important role in the evolution of multiple gene families and in response to stress conditions [46]. In our investigation, we found that the most closely related CsHSP genes in the same subfamily share similar gene structures, indicating that gene structures shed light on phylogenetic relationships. In addition, cytoplasmic CsHSP genes generally have few introns, whereas organelle-specific CsHSPs have more introns, similar to the HSP genes in other plants [16,31,43], suggesting that intron pattern is closely related to gene function. Additionally, there are indications that genes that must be rapidly activated in response to stress tend to evolve a decreased intron density [43]. Our results showed that most CssHSP genes have no introns or just one intron, which contributes to their transcriptional regulation under stress conditions; this feature provides a possible explanation for how almost all of these CssHSP genes are rapidly and strongly induced under heat and drought stress [46]. Furthermore, our conserved motif analysis showed that CsHSPs in the same subfamily contain similar protein motifs but differ significantly among the different subfamilies, especially the members of the CsHSP70 and CssHSP families; this result supports our abovementioned results from the phylogenetic analysis of CsHSPs. Notably, 3 highly conserved protein motifs were detected in the majority of CssHSPs; these motifs form the ACD domain that is the structural basis for the biological function of sHSPs, which is consistent with earlier findings in pepper and switchgrass [43,45].
Numerous studies have shown significant histological specificity in the transcription of plant HSP genes [13,47]. For example, 11 tomato sHSP genes were shown to have tissue- and development-specific expression in the leaf, root, and hypocotyl [44]; the transcription of AtHSP90-1 could be detected only in roots; and several AtHSP90 genes are highly expressed in the root, but have very low or no expression in other organs [48]. Similarly, most PtHSP90 genes in populus are mainly expressed in the stem [12], while the OsHSP90 genes of rice are specifically expressed in different tissues [49]. In the present study, we also found that almost all CsHSP genes were specifically expressed in one or more tissues, especially the leaf and stem of tea plant, similar to the results from recent reports in rice [50], tomato [44], and switchgrass [49], suggesting that the CsHSP genes are extensively involved in the growth and development of tea plant, and may play different roles in different tissues. The name of the HSP family is well known to be derived from their rapid and efficient expression under heat stress, which has been widely confirmed in various plant species [9,11,13]. As expected, all CsHSP90 genes and most CsHSP70 and CssHSP genes were significantly induced under heat stress, and the expression levels of the CssHSP genes were maintained at an extremely high level, similar to the results of earlier investigations in other plants [44,51,52], suggesting that CsHSPs are critical in tea plant response to heat stress. On the other hand, more recent studies have indicated that HSPs also play important roles in plant responses to other biotic and abiotic stresses in addition to heat stress, especially drought stress [4,30]. For example, five switchgrass sHSP genes were significantly upregulated under severe drought stress [45], and some HSP70 and HSP90 genes are also induced by drought stress in rice and populus [30,49,51]. Our results showed that most genes in the CsHSP90, CsHSP70, and CssHSP families also showed upregulated expression under drought stress, although their expression trends were complex and diverse, implying that they actively participated in the drought stress response of tea plant. In addition, we found differences in the response mechanisms of different CsHSP genes to different stresses, implying a diversity of CsHSP functions in the stress response mechanisms of tea plant [45].
Previous investigations have confirmed that HSPs interact with each other to regulate plant growth, development, and stress response [13,53]. HSP90 and HSP70 interact, depending on the action of their co-chaperone molecules, to form a complex to jointly respond to external stress [20,21,22]. In the present study, we also found that the interaction between CsHSP90 and CsHSP70 family proteins is ubiquitous, and may play important roles in tea plant responses to heat and drought stress. Simultaneously, we speculated that HOPs participate in the above process as the co-chaperones connecting the CsHSP90 and CsHSP70 proteins, similar to the recent reports that HOP interacting with HSP90 and HSP70 regulated stress responses in Arabidopsis and orchardgrass [54,55]. In addition, our data confirm that HSP110/SSE subfamily proteins positively interact with cytoplasm subfamily HSP70 proteins, as previously demonstrated in both plants and yeast [56,57]. In addition, some CssHSPs were also found to interact with a few CsHSP70s, similar to previous reports in pea and tobacco [58,59]; thus, they may cooperate to act in the stress response of tea plant. Furthermore, numerous studies have confirmed that the expression of HSP genes is regulated by HSF transcription factors, such as Arabidopsis HSPA2, and populus HSFs transcriptionally regulate the expression of different HSP genes participating in plant development and response to abiotic stresses [60,61]. HSFs have also been confirmed to be involved in the protein interactions of HSPs in plants, which are vital in plant response to multiple stresses [60,61]. Our results showed that most CsHSP90 and CsHSP70 members were closely related to HSFA2, implying that CsHSPs are also regulated by HSFs; this regulation may include the transcript and protein levels. Overall, these results suggest that the interaction network dominated by CsHSP, including HOP and HSF, is ubiquitous and plays important roles in tea plant growth, development and stress response.
In conclusion, this study provides the first comprehensive and systematic analysis of the HSP gene superfamily in tea plant. In total, 47 CsHSP genes from the HSP90, HSP70, and sHSP families were identified, and bioinformatics and expression profile analyses of these genes were performed to determine their potential functions in the growth, development, and stress responses of tea plant. The results indicate that CsHSP genes are actively involved in regulating tea plant growth, development and responses to heat and drought stress, and these processes are closely related to the interaction network involving CsHSP, HOP, and HSF. These results provide an important foundation for further functional studies investigating the CsHSPs and contribute to illuminating the mechanisms of stress tolerance in tea plant.

4. Materials and Methods

4.1. Plant Materials, Growth Conditions, and Stress Treatments

Cutting seedlings of 2-year-old tea plant (C. sinensis cv. “Longjingchangye”) were grown in a chamber at Northwest A&F University (Yangling, China) under a 12 h light (300 μmol·m−2·s−1)/12 h dark cycle at 23 ± 2 °C ambient temperature and 70 ± 5% relative humidity. Then, the tea plants were exposed to heat and drought stress, respectively. For heat treatment, the tea plants grown in normal conditions were transferred to an artificial climate chamber maintained at 40 °C for 48 h. For drought treatment, the root of tea plants together with the medium was immersed completely in the solution containing 20% (w/v) polyethylene glycol (PEG) 6000 for 48 h. All stress treatments were completed under other consistent growth conditions. The first or second tender leaves of the treated tea plants were respectively collected at 0, 1, 2, 4, 8, 12, 24, and 48 h, immediately frozen in liquid nitrogen, and then stored at −80 °C for further analysis.

4.2. Identification of the CsHSP Gene Superfamily in Tea Plant

The amino acid sequences of the HSPs of Arabidopsis thaliana and other plants were downloaded from the TAIR (http://www.arabidopsis.org/) and NCBI databases; the detailed information is shown in Table S2. All the candidate CsHSP genes of tea plant were retrieved from the tea plant genome [37,62] and transcriptome (NCBI SRA: SRP128078) database. The amino acid sequences of the AtHSPs were selected as query sequences to identify the CsHSPs. The online program SMART (http://smart.embl-heidelberg.de/) and the NCBI Conserved Domain Database (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/Structure/cdd/wrpsb.cgi) were used to survey the conserved domains of the candidate proteins.

4.3. Sequence Analysis and Phylogenetic Tree Construction

The exon–intron structures of the CsHSP genes were displayed through Gene Structure Display Server 2.0 (http://gsds.cbi.pku.edu.cn). The sequences of the CsHSPs were analyzed with ExPASy ProtParam (http://www.expasy.org/tools/protparam.html) to obtain the number of amino acids, molecular weight, theoretical isoelectric point (pI), and instability index. The WoLF PSORT program (https://wolfpsort.hgc.jp/) was used to predict the subcellular localization of CsHSPs. Multiple sequence alignment of CsHSPs with other plants’ HSP proteins (Table S2) were performed using DNAMAN 6.0. Unrooted neighbor-joining phylogenetic trees were constructed using MEGA 7.0 and EvolView (http://www.evolgenius.info/evolview).

4.4. Conserved Motif Analysis and Protein Interaction Network Prediction of CsHSPs

The MEME program (version 4.10.0, http://meme-suite.org/) was used to identify the conserved protein motifs of all CsHSPs, with the following parameters: number of repetitions any, maximum number of motifs 10, and optimum motif widths from 6 to 200 amino acid residues. The functional interacting network models of CsHSPs were integrated in STRING (https://string-db.org/) with the confidence parameter set at a 0.40 threshold.

4.5. RNA Isolation and Expression Analysis of CsHSP Genes

Total RNA from tea plant leaves and other tissues was extracted using the Plant RNA Kit (Omega, Norcross, GA, USA), and the concentration and integrity of the RNA samples were checked by a NanoDrop ND 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and 1.2% agarose gel electrophoresis. Subsequently, 1 μg of total RNA was reverse transcribed to first-strand cDNA by the 5× All-In-One RT MasterMix Kit (ABM, Richmond, BC, Canada) according to the manufacturer’s protocol, and the cDNA was diluted 10-fold for standby. Quantitative real-time PCR (qRT-PCR) was performed using SYBR® Premix Ex Taq™ II (TaKaRa, Dalian, China) on an iQ5 real-time PCR platform (Bio-Rad, Hercules, CA, USA). Briefly, each reaction was performed in a total volume of 20 µL containing 10 µL SYBR® Premix Ex Taq II, 1 μL diluted cDNA template, 0.8 μL each primer, and 7.4 μL ddH2O using the following PCR program: 95 °C for 3 min, followed by 40 cycles of 95 °C for 30 s and 60 °C for 1 min 30 s. Melting curves were obtained to verify the amplification specificity through a stepwise heating of the amplicon from 60 to 95 °C. Primer pairs (Table S8) were designed by OligoAnalyzer 3.1 (https://sg.idtdna.com/calc/analyzer) and tested by NCBI Primer BLAST. The CsPTB gene was used as an internal control gene [63]. Three independent biological replicates were performed, and the relative expression levels of the CsHSP genes were calculated with the 2−∆∆Ct method [64]. The heatmap was created using MultiExperiment Viewer (MeV).

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/19/9/2633/s1.

Author Contributions

W.W. conceived and designed the experiments; J.C., T.G. and S.W. performed the experiments; J.C., T.G., Y.Z. and J.Y. analyzed the data; J.C. and W.W. wrote the paper; W.W. and Y.Y. participated in the data analysis and the paper revision.

Funding

This research was funded by the China Postdoctoral Science Foundation (2016M602873), the Fundamental Research Funds for the Central Universities (2452017074), the earmarked fund for Modern Agro-industry Technology Research System (CARS-19), the special fund for University-Supported Extension Model (XTG2018-04).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mittler, R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 2006, 11, 15–19. [Google Scholar] [CrossRef] [PubMed]
  2. Mittler, R.; Finka, A.; Goloubinoff, P. How do plants feel the heat? Trends Biochem. Sci. 2012, 37, 118–125. [Google Scholar] [CrossRef] [PubMed]
  3. Montero-Barrientos, M.; Hermosa, R.; Nicolás, C.; Cardoza, R.E.; Gutiérrez, S.; Monte, E. Overexpression of a Trichoderma HSP70 gene increases fungal resistance to heat and other abiotic stresses. Fungal Genet. Biol. 2008, 45, 1506–1513. [Google Scholar] [CrossRef] [PubMed]
  4. Balfagóna, D.; Zandalinasa, S.I.; Baliñob, P.; Muriachb, M.; Gómez-Cadenas, A. Involvement of ascorbate peroxidase and heat shock proteins on citrus tolerance to combined conditions of drought and high temperatures. Plant Physiol. Biochem. 2018, 127, 194–199. [Google Scholar] [CrossRef] [PubMed]
  5. Park, C.; Seo, Y. Heat Shock Proteins: A Review of the molecular chaperones for plant immunity. Plant Pathol. J. 2015, 31, 323–333. [Google Scholar] [CrossRef] [PubMed]
  6. Vierling, E. The Roles of heat shock proteins in plants. Annu. Rev. Plant Physiol. 1991, 42, 579–620. [Google Scholar] [CrossRef]
  7. Sung, D.Y.; Kaplan, F.; Lee, K.J.; Guy, C.L. Acquired tolerance to temperature extremes. Trends Plant Sci. 2003, 8, 179–187. [Google Scholar] [CrossRef] [Green Version]
  8. Charng, Y.; Liu, H.; Liu, N.; Hsu, F.; Ko, S. Arabidopsis Hsa32, a novel heat shock protein, is essential for acquired thermotolerance during long recovery after acclimation. Plant Physiol. 2014, 140, 1297–1305. [Google Scholar] [CrossRef] [PubMed]
  9. Sarkar, N.K.; Kim, Y.K.; Grover, A. Rice sHsp genes: Genomic organization and expression profiling under stress and development. BMC Genom. 2009, 10, 393. [Google Scholar] [CrossRef] [PubMed]
  10. Waters, E.R. The evolution, function, structure, and expression of the plant sHSPs. J. Exp. Bot. 2013, 64, 391–403. [Google Scholar] [CrossRef] [PubMed]
  11. Timperio, A.M.; Egidi, M.G.; Zolla, L. Proteomics applied on plant abiotic stresses: Role of heat shock proteins (HSP). J. Proteom. 2008, 71, 391–411. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, J.; Li, J.; Liu, B.; Zhang, L.; Chen, J.; Lu, M. Genome-wide analysis of the populus Hsp90 gene family reveals differential expression patterns, localization, and heat stress responses. BMC Genom. 2013, 14, 532. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, W.; Vinocur, B.; Shoseyov, O.; Altman, A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9, 244–252. [Google Scholar] [CrossRef] [PubMed]
  14. Lin, B.; Wang, J.; Liu, H.; Chen, R.; Meyer, Y.; Barakat, A.; Delseny, M. Genomic analysis of the Hsp70 superfamily in Arabidopsis thaliana. Cell Stress Chaperon. 2001, 6, 201–208. [Google Scholar] [CrossRef]
  15. Scharf, K.D.; Siddique, M.; Vierling, E. The expanding family of Arabidopsis thaliana small heat stress proteins and a new family of proteins containing a-crystallin domains (acd proteins). Cell Stress Chaperon. 2001, 6, 225–237. [Google Scholar] [CrossRef]
  16. Krishna, P.; Gloor, G. The HSP90 family of proteins in Arabidopsis thaliana. Cell Stress Chaperon. 2001, 6, 238–246. [Google Scholar] [CrossRef]
  17. Dimitra Milioni, D.; Hatzopoulos, P. Genomic organization of hsp90 gene family in Arabidopsis. Plant Mol. Biol. 1997, 35, 955–961. [Google Scholar] [CrossRef]
  18. Li, Z.; Long, R.; Zhang, T.; Wang, Z.; Zhang, F.; Yang, Q.; Kang, J.; Sun, Y. Molecular cloning and functional analysis of the drought tolerance gene MsHSP70 from alfalfa (Medicago sativa L.). J. Plant Res. 2017, 130, 387–396. [Google Scholar] [CrossRef] [PubMed]
  19. Masand, S.; Yadav, S.K. Overexpression of MuHSP70 gene from Macrotyloma uniflorum confers multiple abiotic stress tolerance in transgenic Arabidopsis thaliana. Mol. Biol. Rep. 2016, 43, 1–12. [Google Scholar] [CrossRef] [PubMed]
  20. Aviezer-Hagai, K.; Skovorodnikova, J.; Galigniana, M.; Farchi-Pisanty, O.; Maayan, E.; Bocovza, S.; Efrat, Y.; von Koskull-Döring, P.; Ohad, N.; Breiman, A. Arabidopsis immunophilins ROF1 (AtFKBP62) and ROF2 (AtFKBP65) exhibit tissue specificity, are heat-stress induced, and bind HSP90. Plant Mol. Biol. 2007, 63, 237–255. [Google Scholar] [CrossRef] [PubMed]
  21. Meiri, D.; Breiman, A. Arabidopsis ROF1 (FKBP62) modulates thermotolerance by interacting with HSP90.1 and affecting the accumulation of HsfA2-regulated sHSPs. Plant J. 2009, 59, 387–399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Li, Z.; Hartl, F.; Bracher, A. Structure and function of Hip, an attenuator of the Hsp70 chaperone cycle. Nat. Struct. Mol. Biol. 2013, 20, 929–935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Fernández-Bautista, N.; Fernández-Calvino, L.; Muñoz, A.; Castellano, M.M. HOP3, a member of the HOP family in Arabidopsis, interacts with BiP and plays a major role in the ER stress response. Plant Cell Environ. 2017, 40, 1341. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, K.H.; Alam, I.; Kim, Y.G.; Sharmin, S.A.; Lee, K.W.; Lee, S.H.; Lee, B.H. Overexpression of a chloroplast-localized small heat shock protein OsHSP26 confers enhanced tolerance against oxidative and heat stresses in tall fescue. Biotechnol. Lett. 2012, 34, 371–377. [Google Scholar] [CrossRef] [PubMed]
  25. Merino, I.; Contreras, A.; Jing, Z.; Gallardo, F.; Francisco, M.; Cánovas, F.M.; Gómez, L. Plantation forestry under global warming: Hybrid poplars with improved thermotolerance provide new insights on the in vivo function of small heat shock protein chaperones. Plant Physiol. 2014, 164, 978–991. [Google Scholar] [CrossRef] [PubMed]
  26. Zhu, Y.; Yan, H.; Wang, Y.; Feng, L.; Chen, Z.; Xiang, Y. Genome duplication and evolution of heat shock transcription factor (HSF) gene family in four model angiosperms. J. Plant Growth Regul. 2016, 35, 903–920. [Google Scholar] [CrossRef]
  27. Jacob, P.; Hirt, H.; Bendahmane, A. The heat-shock protein/chaperone network and multiple stress resistance. Plant Biotechnol. J. 2017, 15, 405–414. [Google Scholar] [CrossRef] [PubMed]
  28. Ouyang, Y.; Chen, J.; Xie, W.; Wang, L.; Zhang, Q. Comprehensive sequence and expression profile analysis of Hsp20 gene family in rice. Plant Mol. Biol. 2009, 70, 341–357. [Google Scholar] [CrossRef] [PubMed]
  29. Lopes-Caitar, V.S.; de Carvalho, M.C.; Darben, M.; Kuwahara, M.K.; Nepomuceno, A.L.; Dias, W.P.; Abdelnoor, R.V.; Marcelino-Guimarães, F.C. Genome-wide analysis of the Hsp20 gene family in soybean: Comprehensive sequence, genomic organization and expression profile analysis under abiotic and biotic stresses. BMC Genom. 2013, 14, 577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Yer, E.N.; Baloglu, M.C.; Ziplar, U.T.; Ayan, S.; Unver, T. Drought-responsive Hsp70 gene analysis in populus at genome-wide level. Plant Mol. Biol. Rep. 2016, 34, 483–500. [Google Scholar] [CrossRef]
  31. Banilas, G.; Korkas, E.; Engezos, V.; Nisiotou, A.A.; Hatzopouos, P. Genome-wide analysis of the heat shock protein 90 gene family in grapevine (Vitis vinifera L.). Aust. J. Grape Wine Res. 2012, 18, 29–38. [Google Scholar] [CrossRef]
  32. Song, D.; Jung, Y.; Chay, K.; Chung, M.A.; Lee, K.H.; Yang, S.Y.; Shin, B.A.; Ahn, B.W. Effect of drinking green tea on age-associated accumulation of maillard-type fluorescence and carbonyl groups in rat aortic and skin collagen. Arch. Biochem. Biophys. 2002, 397, 424–429. [Google Scholar] [CrossRef] [PubMed]
  33. Bordonia, A.; Hrelia, S.; Angelonia, C.; Giordanob, E.; Guarnierib, C.; Caldarerab, C.M.; Biagia, P.L. Green tea protection of hypoxia/reoxygenation injury in cultured cardiac cells. J. Nutr. Biochem. 2002, 13, 103–111. [Google Scholar] [CrossRef]
  34. Wang, W.; Xin, H.; Wang, M.; Ma, Q.; Wang, L.; Kaleri, N.A.; Wang, Y.; Li, X. Transcriptomic analysis of Camellia sinensis reveals the molecular mechanisms of drought-stress-induced decreases in leaf quality. Front. Plant Sci. 2016, 7, 385. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, M.; Zou, Z.; Li, Q.; Xin, H.; Zhu, X.; Chen, X.; Li, X. Heterologous expression of three Camellia sinensis small heat shock protein genes confers temperature stress tolerance in yeast and Arabidopsis thaliana. Plant Cell Rep. 2017, 36, 1125–1135. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, M.; Zou, Z.; Li, Q.; Sun, K.; Chen, X.; Li, X. The CsHSP17.2 molecular chaperone is essential for thermotolerance in Camellia sinensis. Sci. Rep. 2016, 7, 1237. [Google Scholar] [CrossRef] [PubMed]
  37. Wei, C.; Yang, H.; Wang, S.; Zhao, J.; Liu, C.; Gao, L.; Xia, E.; Lu, Y.; Tai, Y.; She, G.; et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl. Acad. Sci. USA 2018, 115, E4151–E4158. [Google Scholar] [PubMed]
  38. Siddique, M.; Gernhard, S.; von Koskull-Döring, P.; Vierling, E.; Scharf, K.D. The plant sHSP superfamily: Five new members in Arabidopsis thaliana with unexpected properties. Cell Stress Chaperon. 2008, 13, 183–197. [Google Scholar] [CrossRef] [PubMed]
  39. Muthusamy, S.D.; Dalala, M.; Chinnusamyc, V.; Bansal, K.C. Genome-wide identification and analysis of biotic and abiotic stress regulation of small heat shock protein (HSP20) family genes in bread wheat. J. Plant Physiol. 2017, 211, 100–113. [Google Scholar] [CrossRef] [PubMed]
  40. Wu, Z.; Li, X.; Liu, Z.; Li, H.; Wang, Y.; Zhuang, J. Transcriptome-wide identification of Camellia sinensis WRKY transcription factors in response to temperature stress. Mol. Genet. Genom. 2016, 291, 255–269. [Google Scholar] [CrossRef] [PubMed]
  41. Wang, X.; Wang, C.; Ban, F.; Zhu, D.; Liu, S.; Wang, X. Genome-wide identification and characterization of HSP gene superfamily in whitefly (Bemisia tabaci) and expression profiling analysis under temperature stress. Insect Sci. 2017. [Google Scholar] [CrossRef] [PubMed]
  42. Dragovic, Z.; Broadley, S.A.; Shomura, Y.; Bracher, A.; Hartl, F.U. Molecular chaperones of the Hsp110 family act as nucleotide exchange factors of Hsp70s. EMBO J. 2006, 25, 2519–2528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Guo, M.; Liu, J.; Lu, J.; Zhai, Y.; Wang, H.; Gong, Z.; Wang, S.; Lu, M. Genome-wide analysis of the CaHsp20 gene family in pepper: Comprehensive sequence and expression profile analysis under heat stress. Front. Plant Sci. 2015, 6, 806. [Google Scholar] [CrossRef] [PubMed]
  44. Yu, J.; Cheng, Y.; Feng, K.; Ruan, M.; Ye, Q.; Wang, R.; Li, Z.; Zhou, G.; Yao, Z.; Yang, Y.; et al. Genome-wide identification and expression profiling of tomato Hsp20 gene family in response to biotic and abiotic stresses. Front. Plant Sci. 2016, 7, 1215. [Google Scholar] [CrossRef] [PubMed]
  45. Yan, H.; Zhang, A.; Chen, J.; He, X.; Xu, B.; Xie, G.; Miao, Z.; Zhang, X.; Huang, L. Genome-wide analysis of the PvHsp20 family in switchgrass: Motif, genomic organization, and identification of stress or developmental-related Hsp20s. Front. Plant Sci. 2017, 8, 1024. [Google Scholar] [CrossRef] [PubMed]
  46. Xua, G.; Guo, C.; Shan, H.; Kong, H. Divergence of duplicate genes in exon-intron structure. Proc. Natl. Acad. Sci. USA 2012, 109, 1187–1192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Wehmeyer, N.; Vierling, E. The expression of small heat shock proteins in seeds responds to discrete developmental signals and suggests a general protective role in desiccation tolerance. Plant Physiol. 2000, 122, 1099–1108. [Google Scholar] [CrossRef] [PubMed]
  48. Yabe, N.; Takahashi, T.; Komeda, T. Analysis of tissue-specific expression of Arabidopsis thaliana HSP90-family gene HSP81. Plant Cell Physiol. 1994, 35, 1207. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, H.; Li, L.; Ye, T.; Chen, R.; Gao, X.; Xu, Z. Molecular characterization, expression pattern and function analysis of the OsHSP90 family in rice. Biotechnol. Biotechnol. Equip. 2016, 30, 1–8. [Google Scholar] [CrossRef] [Green Version]
  50. Mertz-Henning, L.M.; Pegoraro, C.; Maia, L.C.; Venske, E.; Rombaldi, C.V.; de Oliveira, C.A. Expression profile of rice Hsp genes under anoxic stress. Genet. Mol. Res. 2016, 15. [Google Scholar] [CrossRef] [PubMed]
  51. Sarkar, N.K.; Kundnani, P.; Grover, A. Functional analysis of Hsp70 superfamily proteins of rice (Oryza sativa). Cell Stress Chaperon. 2013, 18, 427–437. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, J.; Liu, B.; Li, J.; Zhang, L.; Wang, Y.; Zheng, H.; Lu, M.; Chen, J. Hsf and Hsp gene families in populus: Genome-wide identification, organization and correlated expression during development and in stress responses. BMC Genom. 2015, 16, 181. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, L.; Zhao, H.; Dong, Q.; Zhang, Y.; Wang, Y.; Li, H.; Xing, G.; Li, Q.; Dong, Y. Genome-wide analysis and expression profiling under heat and drought treatments of HSP70 gene family in soybean (Glycine max L.). Front. Plant Sci. 2015, 6, 773. [Google Scholar] [CrossRef] [PubMed]
  54. Fernández-Bautista, N.; Fernández-Calvino, L.; Muñoz, A.; Toribio, R.; Mock, H.P.; Castellano, M.M. HOP family plays a major role in long term acquired thermotolerance in Arabidopsis. Plant Cell Environ. 2018, 41, 1852–1869. [Google Scholar] [CrossRef] [PubMed]
  55. Cha, J.Y.; Baek, D.; Lee, K.W.; Cheong, M.S.; Hong, S.; Son, D. Identification of heat shock-inducible transcripts and Hop (Hsp-organizing protein)-interacting proteins in orchardgrass (Dactylis glomerata L.). Acta Physiol. Plant. 2018, 40, 132. [Google Scholar] [CrossRef]
  56. Makhnevych, T.; Houry, W.A. The control of spindle length by Hsp70 and Hsp110 molecular chaperones. FEBS Lett. 2013, 587, 1067–1072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yam, A.Y.W.; Albanèse, V.; Lin, H.T.J.; Frydman, J. HSP110 cooperates with different cytosolic HSP70 systems in a pathway for de novo folding. J. Biol. Chem. 2005, 280, 41252. [Google Scholar] [CrossRef] [PubMed]
  58. Lee, G.J.; Vierling, E. A small heat shock protein cooperates with heat shock protein 70 systems to reactivate a heat-denatured protein. Plant Physiol. 2002, 122, 189–197. [Google Scholar] [CrossRef]
  59. Mogk, A.; Deuerling, E.; Vorderwülbecke, S.; Vierling, E.; Bukau, B. Small heat shock proteins, ClpB and the DnaK system form a functional triade in reversing protein aggregation. Mol. Microbiol. 2003, 50, 585–595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Liu, H.; Charng, Y. Common and distinct functions of Arabidopsis cass A1 and A2 heat shock factors in diverse abiotic stress responses and development. Plant Physiol. 2013, 163, 276–290. [Google Scholar] [CrossRef] [PubMed]
  61. Hahn, A.; Bublak, D.; Schleiff, E.; Scharf, K.D. Crosstalk between HSP90 and HSP70 chaperones and heat stress transcription factors in tomato. Plant Cell 2011, 23, 741–755. [Google Scholar] [CrossRef] [PubMed]
  62. Xia, E.; Zhang, H.; Sheng, J.; Li, K.; Zhang, Q.; Kim, C.; Zhang, Y.; Liu, Y.; Zhu, T.; Li, W.; et al. The tea tree genome provides insights into tea flavor and independent evolution ofcaffeine biosynthesis. Mol. Plant. 2017, 6, 866–877. [Google Scholar] [CrossRef] [PubMed]
  63. Hao, X.; Horvath, D.P.; Chao, W.; Yang, Y.; Wang, X.; Xiao, B. Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plant (Camellia sinensis (L.) O. Kuntze). Int. J. Mol. Sci. 2014, 15, 22155–22172. [Google Scholar] [CrossRef] [PubMed]
  64. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−△△Ct method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02633 g001 550
Figure 1. Summary of the HSP (heat shock protein) gene superfamily among 21 species.
Figure 1. Summary of the HSP (heat shock protein) gene superfamily among 21 species.
Ijms 19 02633 g001
Ijms 19 02633 g002 550
Figure 2. Phylogenetic tree of CsHSP90 (A), CsHSP70 (B), and CssHSP (C). Phylogenetic tree of HSP proteins in tea plant and other plant species were generated by MEGA 7 using neighbor-joining. Cs: Camellia sinensis, At: Arabidopsis thaliana, Os: Oryza sativa, Gm: Glycine max, and Sl: Solanum lycopersicum; CI–CVII: cytoplasm I–VII, ER: endoplasmic reticulum, MI–MII: mitochondria I–II, P: plastid, and Px: peroxide. The protein ID of all species is listed in Table S2.
Figure 2. Phylogenetic tree of CsHSP90 (A), CsHSP70 (B), and CssHSP (C). Phylogenetic tree of HSP proteins in tea plant and other plant species were generated by MEGA 7 using neighbor-joining. Cs: Camellia sinensis, At: Arabidopsis thaliana, Os: Oryza sativa, Gm: Glycine max, and Sl: Solanum lycopersicum; CI–CVII: cytoplasm I–VII, ER: endoplasmic reticulum, MI–MII: mitochondria I–II, P: plastid, and Px: peroxide. The protein ID of all species is listed in Table S2.
Ijms 19 02633 g002
Ijms 19 02633 g003 550
Figure 3. Exon–intron structure analyses of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes. The yellow sections represent exons, and the grey parts indicate introns. CI, CII, CV, and CVI: cytoplasm I, II, V and VI; ER: endoplasmic reticulum, MI: mitochondria I; P: plastid; and Px: peroxide.
Figure 3. Exon–intron structure analyses of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes. The yellow sections represent exons, and the grey parts indicate introns. CI, CII, CV, and CVI: cytoplasm I, II, V and VI; ER: endoplasmic reticulum, MI: mitochondria I; P: plastid; and Px: peroxide.
Ijms 19 02633 g003
Ijms 19 02633 g004 550
Figure 4. Distribution of conserved motifs in CsHSP90 (A), CsHSP70 (B), and CssHSP (C). Putative motifs are represented by a number in a colored box. CI, CII, CV, and CVI: cytoplasm I, II, V and VI; ER: endoplasmic reticulum, MI: mitochondria I; P: plastid, and Px: peroxide.
Figure 4. Distribution of conserved motifs in CsHSP90 (A), CsHSP70 (B), and CssHSP (C). Putative motifs are represented by a number in a colored box. CI, CII, CV, and CVI: cytoplasm I, II, V and VI; ER: endoplasmic reticulum, MI: mitochondria I; P: plastid, and Px: peroxide.
Ijms 19 02633 g004
Ijms 19 02633 g005 550
Figure 5. Putative interaction network of CsHSPs in tea plant. Homologous proteins in tea plant and Arabidopsis are shown in red and black, respectively.
Figure 5. Putative interaction network of CsHSPs in tea plant. Homologous proteins in tea plant and Arabidopsis are shown in red and black, respectively.
Ijms 19 02633 g005
Ijms 19 02633 g006 550
Figure 6. Expression patterns of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in different tissues of tea plant. R indicates root; S, stem; L, leaf; F, flower. The heatmap was generated by MeV software (TIGR, Rockville, MD, USA) using the CsHSP genes’ expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Figure 6. Expression patterns of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in different tissues of tea plant. R indicates root; S, stem; L, leaf; F, flower. The heatmap was generated by MeV software (TIGR, Rockville, MD, USA) using the CsHSP genes’ expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Ijms 19 02633 g006
Ijms 19 02633 g007 550
Figure 7. Expression profiles of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in tea plant under heat stress (40 °C). The heatmap was generated by MeV software using the CsHSP genes expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Figure 7. Expression profiles of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in tea plant under heat stress (40 °C). The heatmap was generated by MeV software using the CsHSP genes expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Ijms 19 02633 g007
Ijms 19 02633 g008 550
Figure 8. Expression profiles of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in tea plant under drought stress (20% polyethylene glycol (PEG) 6000). The heatmap was generated by MeV software using the CsHSP genes expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Figure 8. Expression profiles of CsHSP90 (A), CsHSP70 (B), and CssHSP (C) genes in tea plant under drought stress (20% polyethylene glycol (PEG) 6000). The heatmap was generated by MeV software using the CsHSP genes expression data, and normalized log2 transformed values were used with hierarchical clustering. The red and green colors indicate higher or lower transcript abundances than those of the relevant control, respectively.
Ijms 19 02633 g008
Table
Table 1. Summary information of physiological and biochemical properties of the CsHSP proteins.
Table 1. Summary information of physiological and biochemical properties of the CsHSP proteins.
NameUnigene IDAnimo AcidsMW (kDa)pIGRAVYInstability IndexAliphatic IndexSubcellular Localization
CsHSP90-1TEA014732.179089.905.20−0.5834.3777.62mito: 7, chlo: 4, nucl: 2, cysk: 1
CsHSP90-2TEA027510.182394.334.84−0.7537.0179.62E.R.: 14
CsHSP90-3TEA017408.181092.315.01−0.6043.7780.26chlo: 11.5, chlo_mito: 7.5, mito: 2.5
CsHSP90-4TEA027510.170180.445.11−0.6035.0884.95cyto: 7, nucl: 2, plas: 2, mito: 1, vacu: 1, golg: 1
CsHSP90-5TEA008112.169979.994.99−0.5636.6884.23cyto: 8, plas: 2, E.R.: 2, mito: 1, golg: 1
CsHSP90-6TEA007238.170781.145.04−0.5942.7582.26nucl: 6, cyto: 5, chlo: 2, cysk: 1
CsHSP90-7TEA027790.164373.685.11−0.5039.3385.99cyto: 9, nucl: 3, plas: 1, golg: 1
CsHSP70-1TEA005970.164871.035.16−0.4234.5681.87cyto: 7, cysk: 6, chlo: 1
CsHSP70-2TEA029156.184693.335.21−0.4244.2080.18cyto: 8, chlo: 2, nucl: 2, cysk: 2
CsHSP70-3TEA023352.189298.895.22−0.4641.3383.55E.R.: 14
CsHSP70-4TEA003665.170975.515.38−0.2933.7085.51chlo: 14
CsHSP70-5TEA016684.156961.815.540.0945.27100.39cyto: 11, chlo: 2, nucl: 1
CsHSP70-6TEA022542.167973.096.02−0.3341.8286.07mito: 12, chlo: 2
CsHSP70-7TEA023884.166973.735.03−0.4531.4086.43E.R.: 14
CsHSP70-8CL5657.Contig2_All▲65171.195.20−0.4534.6680.06cyto: 12, chlo: 1, nucl: 1
CsHSP70-9TEA025820.167775.125.54−0.2734.5687.70plas: 11, nucl: 1, vacu: 1, E.R.: 1
CsHSP70-10TEA024518.174482.425.79−0.4332.9682.34cyto: 8, cysk: 5, chlo: 1
CsHSP70-11TEA023111.166974.025.14−0.4231.8687.43plas: 11.5, golg_plas: 7, golg: 1.5, E.R.: 1
CsHSP70-12CSA029398.1◆51957.326.15−0.3139.7688.86cysk: 11, chlo: 2, cyto: 1
CsHSP70-13TEA014652.161467.305.43−0.3731.3685.94chlo: 9, nucl: 1, cyto: 1, extr: 1, vacu: 1, pero: 1
CsHSP70-14TEA024522.157463.627.99−0.4337.8185.45cyto: 11, chlo: 2, nucl: 1
CsHSP70-15TEA029776.162769.615.070.0136.23106.67cyto: 9, chlo: 2, nucl: 2, cysk: 1
CsHSP70-16TEA029256.144949.875.32−0.5240.9979.73extr: 6, cyto: 3, chlo: 2, vacu: 1, E.R.: 1, mito_plas: 1
CsHSP70-17TEA024511.140745.469.04−0.3832.4187.00chlo: 11, mito: 1, extr: 1, E.R.: 1
CsHSP70-18TEA017380.176785.855.60−0.4143.6483.42nucl: 6, chlo: 3, cyto: 2, cysk: 2, vacu: 1
CssHSP-1TEA023567.113615.465.60−0.1858.9188.16cyto: 10, chlo: 1, nucl: 1, extr: 1, golg: 1
CssHSP-2TEA025531.115517.425.13−0.5752.6080.52cyto: 13, chlo: 1
CssHSP-3TEA030322.119922.455.99−0.5652.9288.14mito: 11, chlo: 3
CssHSP-4TEA028966.124327.408.29−0.7143.7862.55chlo: 13, nucl: 1
CssHSP-5TEA031630.125528.088.76−0.4360.8269.25chlo: 10, nucl: 2, extr: 2
CssHSP-6TEA004366.115918.255.68−0.7062.2171.70cyto: 14
CssHSP-7TEA031692.115317.475.53−0.6154.3873.20cyto: 14
CssHSP-8TEA033542.122325.195.83−0.5346.8981.26chlo: 10, mito: 2, vacu: 2
CssHSP-9TEA015431.115017.255.56−0.7144.1068.20cyto: 13, extr: 1
CssHSP-10CL9235.Contig2_All▲14416.086.76−0.4941.0983.12pero: 8, nucl: 2, cyto: 2, chlo: 1, golg: 1
CssHSP-11CSA020588.1◆15918.349.81−0.5645.2283.96cyto: 6, mito: 5, nucl: 3
CssHSP-12CSA025700.1◆19621.665.30−0.7152.6768.67cyto: 8, nucl: 6
CssHSP-13TEA032865.114917.058.83−0.7239.9869.93cyto: 14
CssHSP-14TEA015431.114917.115.38−0.6544.2768.66cyto: 7, nucl: 3, extr: 2, chlo: 1, plas: 1
CssHSP-15TEA010603.119822.774.83−0.5846.1469.34chlo: 5, nucl: 5, cyto: 3, plas: 1
CssHSP-16TEA021291.119522.155.29−0.5943.4370.36nucl: 5, cyto: 5, chlo: 3, plas: 1
CssHSP-17TEA033850.116418.085.78−0.3739.3083.84cyto: 8, nucl: 5, chlo: 1
CssHSP-18TEA005909.120522.919.52−0.4132.1384.73cyto: 3, golg: 3, plas: 2.5, E.R.plas: 2.5, mito: 2, E.R.: 1.5, chlo: 1, vacu: 1
CssHSP-19TEA005935.119920.958.93−0.2847.3476.58E.R.: 3.5, golg: 3, cyto: 2.5, E.R.plas: 2.5, mito: 2, cyto_pero: 2, chlo: 1, vacu:1
CssHSP-20TEA004364.124927.729.49−0.4348.2874.90mito: 9, nucl: 3, chlo: 2
CssHSP-21TEA028346.114016.256.77−0.8842.5574.43cyto: 5, extr: 4, nucl: 2, golg: 2, plas: 1
CssHSP-22TEA017741.118821.786.47−0.4949.1392.82chlo: 7, mito: 3, cyto: 2, extr: 1, E.R.: 1
Note: Sequenced IDs marked with triangle are from transcriptome, with diamond are from C. sinensis var. assamica (CSA) genome, others are from C. sinensis var. sinensis (CSS) genome. chlo: chloroplast, cyto: cytoplasm, ER: endoplasmic reticulum, extr: extracellular, golg: golgi apparatus, mito: mitochondria, nucl: nucleus, pero: peroxide, plas: plasma membrane, vacu: vacuole.

Share and Cite

MDPI and ACS Style

Chen, J.; Gao, T.; Wan, S.; Zhang, Y.; Yang, J.; Yu, Y.; Wang, W. Genome-Wide Identification, Classification and Expression Analysis of the HSP Gene Superfamily in Tea Plant (Camellia sinensis). Int. J. Mol. Sci. 2018, 19, 2633. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092633

AMA Style

Chen J, Gao T, Wan S, Zhang Y, Yang J, Yu Y, Wang W. Genome-Wide Identification, Classification and Expression Analysis of the HSP Gene Superfamily in Tea Plant (Camellia sinensis). International Journal of Molecular Sciences. 2018; 19(9):2633. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092633

Chicago/Turabian Style

Chen, Jiangfei, Tong Gao, Siqing Wan, Yongheng Zhang, Jiankun Yang, Youben Yu, and Weidong Wang. 2018. "Genome-Wide Identification, Classification and Expression Analysis of the HSP Gene Superfamily in Tea Plant (Camellia sinensis)" International Journal of Molecular Sciences 19, no. 9: 2633. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092633

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