Next Article in Journal
Matrix Metalloproteinase-3 is Key Effector of TNF-α-Induced Collagen Degradation in Skin
Next Article in Special Issue
Comparative Genome-wide Analysis and Expression Profiling of Histone Acetyltransferase (HAT) Gene Family in Response to Hormonal Applications, Metal and Abiotic Stresses in Cotton
Previous Article in Journal
Hydrogen Sulfide in Bone Tissue Regeneration and Repair: State of the Art and New Perspectives
Previous Article in Special Issue
Genome-Wide Identification and Characterization of Cucumber BPC Transcription Factors and Their Responses to Abiotic Stresses and Exogenous Phytohormones
 
 
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 Hsf Gene Family in Carnation (Dianthus caryophyllus)

1
College of Landscape Architecture and Forestry, Qingdao Agricultural University, Qingdao 266000, China
2
State Key Laboratory of Tree Genetics and Breeding, Key Laboratory of Tree Breeding and Cultivation of the National Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China
3
Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2019, 20(20), 5233; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20205233
Submission received: 21 September 2019 / Revised: 16 October 2019 / Accepted: 18 October 2019 / Published: 22 October 2019
(This article belongs to the Special Issue Mapping Abiotic Stress-Tolerance Genes in Plants)

Abstract

:
Heat shock transcription factors (Hsfs) are a class of important transcription factors (TFs) which play crucial roles in the protection of plants from damages caused by various abiotic stresses. The present study aimed to characterize the Hsf genes in carnation (Dianthus caryophyllus), which is one of the four largest cut flowers worldwide. In this study, a total of 17 non-redundant Hsf genes were identified from the D. caryophyllus genome. Specifically, the gene structure and motifs of each DcaHsf were comprehensively analyzed. Phylogenetic analysis of the DcaHsf family distinctly separated nine class A, seven class B, and one class C Hsf genes. Additionally, promoter analysis indicated that the DcaHsf promoters included various cis-acting elements that were related to stress, hormones, as well as development processes. In addition, cis-elements, such as STRE, MYB, and ABRE binding sites, were identified in the promoters of most DcaHsf genes. According to qRT-PCR data, the expression of DcaHsfs varied in eight tissues and six flowering stages and among different DcaHsfs, even in the same class. Moreover, DcaHsf-A1, A2a, A9a, B2a, B3a revealed their putative involvement in the early flowering stages. The time-course expression profile of DcaHsf during stress responses illustrated that all the DcaHsfs were heat- and drought-responsive, and almost all DcaHsfs were down-regulated by cold, salt, and abscisic acid (ABA) stress. Meanwhile, DcaHsf-A3, A7, A9a, A9b, B3a were primarily up-regulated at an early stage in response to salicylic acid (SA). This study provides an overview of the Hsf gene family in D. caryophyllus and a basis for the breeding of stress-resistant carnation.

1. Introduction

Plant growth and production are affected by abiotic stresses such as heat, cold, drought, and salinity [1,2,3]. Unlike animals, plants are sessile organisms. Consequently, to cope with environmental stresses, plants have evolved a series of defense or signaling mechanisms. Furthermore, each process involves different types of transcription factors (TFs). These include heat shock transcription factors (Hsfs), such as WRKY, MYB, AP2/ERF, and NAC, which regulate the expression of thousands of genes under various stress conditions [4,5,6]. In plants, the Hsf family is one of the most important TF families in plants involved in resistance to heat [7] and other abiotic stresses or chemical stressors, such as abscisic acid (ABA) and salicylic acid (SA) [8,9]. Hsfs regulate the expression of Heat shock proteins (Hsps) as well as other stress-responsive proteins, such as reactive oxygen species (ROS)-scavenging enzymes [ascorbate peroxidase (APX) and catalase (CAT)] [4]. Besides their roles in stress responses, Hsfs are also involved in plant growth and development [10,11,12].
Similar to many other TFs, Hsfs are a part of an evolutionarily conserved gene family. Hsf genes are composed of several structurally and functionally conserved domains, including DNA-binding domains (DBD), N-terminal adjacent bipartite oligomerization domains (HR-A/B), nuclear localization signals (NLS), nuclear export signals (NES), C-terminal activator peptide proteins (AHA), and repressor domains (RD) [13]. Among these conserved domains, DBD is characterized by a central helix–turn–helix motif and is responsible for binding to the heat shock elements (HSEs) of the target genes [7]. Notably, the HSEs are palindromic binding motifs (5’-AGAAnnTTCT-3’) conserved in the promoters of heat stress-inducible genes [7,14]. According to the flexible linker of variable lengths (about 15–80 amino acids) and HR-A/B regions, plant Hsfs can be divided into at least three types, i.e., class A (subclasses A1, A2, A3, A4, A5, A6, A7, A8, and A9), class B (subclasses B1, B2, B3, and B4), and class C (subclasses C1 and C2) [15,16,17].
The size of the Hsf gene family varies significantly in different plant species. For instance, there are 22 Hsf members in the model plant Arabidopsis thaliana [18], 25 members in Oryza sativa [18], 16 members in Medicago truncatula [19], 26 members in Glycine max [20], 25 members in Zea mays [21], 25 members in Malus domestica [22], 21 members in Cucumis sativus [23], 28 members in Populus trichocarpa [24], 40 members in Gossypium hirsutum [25], and 56 members in members in Triticum aestivum [26]. To date, the largest Hsf gene family has been identified in Brassica napus, with 64 Hsfs [27].
Carnation (Dianthus caryophyllus L.) is a major floricultural crop and one of the four largest cut flowers [28,29]. Until now, more than 300 Dianthus species have been identified worldwide. Carnations are cultivated widely for their attractive characteristics such as flower color, flower size, fragrance, and flower longevity. However, the vegetative and reproductive growth of carnations are severely impaired in heat stress conditions, resulting in flower wilting and quality decline [30]. The completion of the draft genome sequence of D. caryophyllus L. has greatly facilitated the identification of Hsfs at the whole-genome level and it is extremely important to study the heat-resistant mechanism of carnation [31]. To our knowledge, there are no reports on the identification and functional analysis of carnation Hsfs to date. In this study, we aimed to comprehensively study the structural and expression profiles of the Hsf gene family in D. caryophyllus. A total of 17 putative genes were identified and characterized as members of the Hsf gene family from D. caryophyllus. Additionally, we performed bioinformatic analyses of phylogenetic relationships, conserved domains, motifs, and other. Furthermore, the expression level of these genes in various tissues and in response to abiotic stresses were compared. Our results will be a reference and provide valuable information for the functional analysis of the Hsf genes in D. caryophyllus.

2. Results

2.1. Identification of DcaHsfs in Carnation

The amino acid sequences of putative Hsf proteins were examined using the conserved Hsf domain (PF00447) from the carnation database (DB, http://carnation.kazusa.or.jp). Additionally, searches using the BLASTP program resulted in the identification of putative Hsf gene candidates. In total, 17 proteins were retrieved as DcaHsfs in D. caryophyllus.
The physical and chemical properties of the 17 DcaHsfs were analyzed (Table 1). The DcaHsfs ranged from 133 amino acids (aa; DcaHsf-C1, incomplete) to 495 aa (DcaHsf-A5) in length. The predicted isoelectric points (pI) varied from 4.74 (DcaHsf-A1) to 8.88 (DcaHsf-B3a and DcaHsf-B3b), and the molecular weight (MW) varied from 15.89 kDa (DcaHsf-C1) to 54.52 kDa (DcaHsf-A5). The instability index, i.e., the stability of the protein in a test tube, indicated that all DcaHsfs were unstable, except for DcaHsf-A9b and DcaHsf-B1. The GRAVY value reflects the hydropathicity of a protein; the low GRAVY values (<0) of DcaHsfs suggest that all DcaHsfs are hydrophilic. The total number of negatively charged residues (Asp + Glu, n.c.r.) and the total number of positively charged residues (p.c.r.) of class A were all greater than those of class B and C. These differences might be caused by differences in the amino acid composition of the non-conserved region. We determined the scaffold locations of DcaHsfs on the basis of the information from the Carnation genomic database. We mapped 17 DcaHsfs to 17 scaffolds, and these genes were distributed evenly in the Carnation genome (Figure 1).

2.2. Phylogenetic and Sequence Conservation Analysis of DcaHsfs

To explore the phylogenetic relationship of Hsfs in D. caryophyllus and other species, the amino acid sequences of Hsfs from A. thaliana, O. sativa, and P. trichocarpa were used, together with those of DcaHsfs, as a means to construct a phylogenetic tree. In this study, 21 Hsf proteins from A. thaliana [18], 25 from O. sativa [18], 31 from P. trichocarpa [32], and 17 from D. caryophyllus were utilized for the phylogenetic analysis. A total of 94 Hsf proteins from the four species were clearly divided into three classes (class A, B, and C) with well-supported bootstrap values (Figure 2).
In D. caryophyllus, 9 DcaHsfs out of 17 proteins belonged to class A, making it the largest subclass, followed by 7 DcaHsfs belonging to class B. The number of class B Hsfs in D. caryophyllus (7) was greater than that in Arabidopsis (5). The class C Hsf was present as a single copy in D. caryophyllus, Arabidopsis, and P. trichocarpa, whereas four copies of class C Hsfs were discovered in O. sativa. However, none of the DcaHsfs belonged to subclasses A6 and A8. Sequence conservation among DcaHsfs was also supported by their identity at the amino acid level. Detailed information on the identity of AtHsfs, OsHsfs, and PtHsfs amino acid sequences is illustrated in Table S1.

2.3. Structural and Motif Analysis of DcaHsfs

The structural diversity of the DcaHsf family was analyzed in terms of the exon/intron arrangement of the coding sequences. The number of introns in DcaHsfs ranged from one to three. The detailed gene structure of DcaHsfs is pictured in Figure 3a. Three introns were identified in DcaHsfs-A7, whereas all the other DcaHsfs had only one intron. Most closely related DcaHsfs in the same class or subfamily shared a similar gene structure in terms of intron number and intron and exon length (Figure 3a).
To investigate the protein sequence features of DcaHsfs, 20 different motifs were identified in DcaHsfs, with lengths ranging from 10 to 50 aa. (Figure 3b, Table 2). All members showed similar motif composition, but small differences between different groups were also found (Figure 3b). The conserved motifs in Hsf genes indicated that all DcaHsfs contained motif 1, motif 2, motif 3, and motif 4, except for DcaHsf-C1. Additionally, some motifs were only discovered in a certain subfamily of DcaHsfs. For instance, motif 6, motif 8, motif 9, and motif 14 were present in the B2 subfamily, whereas motif 10 and motif 12 were present in the B3 subfamily. Specifically, the phylogenetic analysis showed that the same clusters of DcaHsfs shared similar conserved domain composition. This indicates that Hsf genes are evolutionarily well conserved or possess similar regulatory functions in D. caryophyllus (Figure 3b). Additionally, motifs 1–3 represent the Hsf DBD domains (~100 aa). The DBD domain contains three α-helices and a four-stranded antiparallel β-sheet (α1-β1-β2-α2-α3-β3-β4) (Figure S1).
The 20 motifs consist of six different domains, including DBD, HR-A/B, NLS, NES, RD, and AHA domains. Among these, the highly structured DBD domain is the most conserved section in the DcaHsf family (Table 3). In addition to DBD, HR-A/B is critical for Hsf–Hsf interactions in the formation of a trimer [7], HR-A/B is also present in all DcaHsfs, and class A Hsfs have longer HR-A/B regions compared with class B and class C Hsfs (Figure S1, Table 3) Meanwhile, the other four conserved domains were only identified in specific DcaHsf members. The majority of class A DcaHsfs contained an NLS sequence rich in basic amino acid residues (K/R), except DcaHsf-A1/A3, whereas two or three NLS domains were located in seven DcaHsfs (A2a, A2b, A7, A9a, A9b, B3a, and B3b). NLS domains were not identified in DcaHsf-C1 and in some class B proteins (DcaHsf-B1, B2a, B2b, B2c). NES motifs were found in nine DcaHsfs. Also, five Class B Hsfs, except DcaHsf-B2a and DcaHsf-B2b, contained an RD in the C-terminus, characterized by the tetrapeptide LFGV. Transcription activator AHA motifs were located in class A DcaHsfs (Figure S1, Table 3). Sequence conservation among DcaHsfs was also supported by their identity at the amino acid level (0.023–0.83, Table S2). Four pairs of DcaHsfs (A2a–A2b, A9a–A9b, B2a–B2b, and B3a–B3b) exhibited high sequence identity (Table S2).

2.4. Cis-Acting Element Analysis in the Promoters of DcaHsfs

To predict the biological function of DcaHsfs, 1500-bp upstream sequences from the translation start sites of DcaHsfs were analyzed through the PlantCARE database. The promoter of each DcaHsfs consists of several cis-acting elements, such as phytohormone-, abiotic stress-, and developmental process-related elements. As illustrated in Figure 4, the MYB element, ARE element (essential for anaerobic induction), and STRE element (activated by heat shock, osmotic stress, low pH, and nutrient starvation) were identified in the promoters of 15, 12, and 12 DcaHsf genes, respectively. The promoters of 11 DcaHsfs contained the ABA-responsive element (ABRE), methyl jasmonate (MeJA)-responsive element (CGTCA-motif), and TGACG motif involved in MeJA responsiveness (Figure 4). Also, the ethylene-responsive element (ERE), cis-acting element involved in salicylic acid responsiveness (TCA-element), stress-inducible element (TCA), wounding and pathogen responsiveness elements (W-box) were all found in 10, 10, 8, and 8 DcaHsfs, respectively (Figure 4). In total, 17 DcaHsf promoters contained 30 MYB, 30 STRE, 23 ARE, 21 ABRE, and 17 CGTCA-motif elements (Figure 4). These findings demonstrate that DcaHsfs might be associated with various transcriptional regulations involving development, hormones, and stress responses.

2.5. The Expression Pattern of DcaHsfs in Different Tissues and Flower Development

To elucidate the tissue-specific expression patterns of DcaHsfs, qRT-PCR was utilized to determine the expression levels of 17 DcaHsfs in 8 carnation tissues [root (R), stem (S), calyx (CA), young leaf (YL), mature leaf (ML), stigma (ST), ovary (OV), and flower (F)] and at 6 flowering stages (FS1, FS2, FS3, FS4, FS5, and FS6) (Figure 5, Table S4). Interestingly, the expression levels differed in different tissues and flowering stages, and the expression patterns of different members of DcaHsfs also differed, even for the same class. Among the different tissues, DcaHsf-A1, A2a, 2b, A7, A9b, B2a, 2c, B4 were up-regulated in S, CA, YL, and ML. Meanwhile, all DcaHsfs were down-regulated in ST, OV, and F (Figure 5). Twelve DcaHsfs were more highly expressed in CA, and 10 out of 17 DcaHsfs had higher expression levels in ML (Figure 5). Some genes demonstrated tissue-specific expression patterns. For instance, DcaHsf-A9a, B1, B3a, B3b, C1 were up-regulated in R, and DcaHsf-A3, A4, A5 were expressed at high levels in CA.
During the six flowering stages of carnations, all DcaHsfs showed relatively high expression levels at FS1. Additionally, DcaHsfA1, A2a, A9a, B2a, B3a were up-regulated at FS2 (Figure 5), implying that these genes may be involved in the early development of carnation flowers. In contrast, DcaHsf-A5 and DcaHsf-B2b exhibited a high expression level at FS6 (Figure 5).

2.6. DcaHsfs Response to Various Stresses

To determine the potential roles of the DcaHsfs in plant responses to various environmental stresses, qRT-PCR was conducted on the 17 DcaHsfs using the leaves of carnations exposed to heat, cold, drought, salt, ABA, and SA. The results illustrated that almost all DcaHsfs revealed three types of expression patterns under different stress conditions: (1) the expression of all genes was up-regulated; (2) the expression of all genes was down-regulated; and (3) some genes were expressed at higher levels in the early stage of stress, while others were up-regulated in the later stage of stress (Figure 6). Regarding the first category (1), all DcaHsfs were up-regulated after leaf exposure to heat and polyethylene glycol (PEG) treatments (Figure 6). For the second category of genes (2), almost all DcaHsfs displayed a decrease in their expression levels under cold, salt, or ABA stresses (except for individual DcaHsfs) (Figure 6). For example, four DcaHsfs (DcaHsf-A2a, A5, B2b, C1) were slightly induced at different time points at 4 °C., while the transcription levels of the remaining 13 genes were down-regulated at the tested time points (Figure 6). DcaHsf-A5 demonstrated higher transcript accumulation compared to the other genes at 12 h under 200 mM NaCl treatment. Meanwhile, DcaHsf-A3 was slightly up-regulated at 12 h under ABA treatment (Figure 6). Finally, the genes in category (3) and the expression of DcaHsf-A3, A7, A9a, A9b, B3a were primarily up-regulated at the earlier stage of SA treatment, whereas other DcaHsfs were strongly up-regulated after 12 h of SA treatment (Figure 6). These findings indicate that DcaHsf genes might play crucial roles in different stress response pathways.

3. Discussion

3.1. Characterization of the Carnation Hsf Genes Family

Hsfs exist extensively in all plant species and act as the key regulatory components involved in various abiotic stresses to protect the plant cellular machinery under stress conditions [4,13,26]. In this study, a comprehensive genome-wide analysis of the DcaHsf family in carnations was carried out for the first time. A total of 17 DcaHsf genes were identified from the Carnation genome database [31]. The size of the carnation Hsf gene family is smaller compared with that of three other plant species, i.e., A. thaliana, O. sativa, and P. trichocarpa. Meanwhile, all four species have a similar subfamily distribution, which indicates that parallel evolutionary events of Hsf genes occurred in dicots and monocots. Additionally, the subclasses A6 and A8 are absent in carnation, and the diversification of Hsf members could provide some clues about the biological function of the corresponding Hsf counterparts in carnation. This suggests that gene loss and gene duplication events occurred at different stages of the evolutionary process, resulting in Hsf diversity [18] (Table S1).

3.2. Cis-Element Analysis in the Promoters of DcaHsfs

The number and form of cis-elements in promoter regions might play an essential function in the regulation of gene expression related to metabolic pathways [33]. The results illustrate that abiotic stress-related cis-elements, including MYB, STRE, ARE, ABRE, CGTCA-motif element, ERE, TCA-element, and W-box, are major regulatory elements in DcaHsfs promoters activated by heat shock or other abiotic stress. The presence of these stress-related elements is related to the expression response of DcaHsfs to heat, drought, ABA, and SA treatments. STRE is a marker element for plant Hsf proteins, which has been located in the promoters of the 17 DcaHsfs (Figure 4). In our study, a large number of STRE elements were identified in the promoter of 12 DcaHsf genes, which coincides with their expression (Figure 6). These findings suggest that STRE plays a vital role in transcriptional regulation under heat conditions in carnation. DcaHsf subclass A promoters contained MYB binding sites which participate in drought, low temperature, salt, ABA, and gibberellic acid (GA) stress responses [34]. We found that 15 DcaHsfs included 30 MYB binding sites in their promoter regions. However, the presence of MYB elements seems to be correlated with the positive regulation of DcaHsfs during drought and the negative regulation of DcaHsfs in response to salt and ABA treatments (Figure 6, Table S3). Other abiotic stress-related cis-elements, including the CGTCA-motif, TGACG-motif, ERE element, and TCA-element, were also major regulatory elements identified in DcaHsfs. Furthermore, the presence of these stress-related elements appears to be correlated with MeJA, SA, and stress responsiveness, suggesting their potential roles in the response to pathogen infections. Consequently, DcaHsfs could be taken as candidate genes to understand the responses to drought and other biotic stresses.

3.3. Structural Analysis of DcaHsfs

The detailed knowledge of A. thaliana, O. sativa, and P. trichocarpa Hsf functional domains enabled us to analyze similar domains in D. caryophyllus Hsf gene family. It has been reported that the number of introns both regulate gene expression and participate in gene evolution [35]. Analysis of Hsf gene structure revealed that 16 of 17 DcaHsfs have one intron in their DBD domain (Figure 3a), which is an evolutionarily conserved intron [36]. However, DcaHsf-A7 contains three introns (Figure 3a), which might affect its expression under stress conditions. All 17 DcaHsfs proteins contain the necessary DBD domain and specific protein domains (HR-A/B, NLS, NES, RD, and AHA) (Table 3, Figure 3b), which provide the structural basis for their conserved function [22]. The Hsf DBD domain of approximately 100 amino acid residues is highly conserved in different organisms, from plants to animals [7]. However, the DBD of DcaHsf-C1 contains only 44 aa and is shorter than the other DcaHsfs, lacking the full α1-helix, β1-sheet, β1-sheet, and α2-helix. Notably, this might be caused by the current genome assembly. It is interesting that AHA, an essential domain for the activator function in the HsfA class [7], was not found in several members of class A DcaHsfs (A1, A3, A9a, and A9b) (Table 3). The members of Hsfs lacking AHA domains might contribute differently to the activator function or bind to other HsfAs to form hetero-oligomers [18].

3.4. DcaHsfs Involvement in Carnation Development Processes

The expression patterns of DcaHsfs in seven different organs or tissues uncovered that DcaHsfs have different expression profiles in carnation. This suggests that they may participate in various developmental processes or regulatory pathways. In this study, nearly all the DcaHsfs were found to display high transcription levels in ML and at FS1 (Figure 5). Within the potato HsfA1 group, StHsf002 is highly expressed in flowers, petals, and sepals, whereas StHsf003 is highly expressed in roots, flowers, carpels, and sepals [37]. In our study, DcaHsf-A1 demonstrated up-regulation in S, CA, YL, and ML (Figure 5). Phyllostachys edulis PheHsfA2a-2 is predicted to play an important role in flower and shoot development [38] and Cicer arietinum CarHsfA2 is up-regulated in shoot, root, and flower [12]. Their orthologs in carnations, DcaHsf-A2a, A2b, were constitutively expressed in S, CA, YL, and ML at relatively high levels (Figure 5). These findings indicate that the members of the Hsf-A1 and A2 sub-families are conserved and involved in the development of vegetative organs.
HsfA5 has been reported to play a vital role in stress tolerance during anther/pollen development as well as in other stages of plant reproduction in tomato and Arabidopsis [22,39]. In this study, DcaHsf-A5 was highly expressed in CA and OV (Figure 5). This implies that the function of DcaHsf-A5 might be conserved for regulating reproductive organ development and the growth of carnation. Salix suchowensis SsuHsf-A9 is specifically expressed in the female catkin [32]. In Populus female catkin development, PtHsf-A9 displays relatively high transcription levels [40]. Our results indicate that carnation DcaHsf-A9b was up-regulated in S and CA and is possibly widely involved in the development of both vegetative and reproductive tissues.
For Class B Hsfs, Chickpea CarHsfB2c is highly expressed in the late flowering stages, while CarHsfB2a is expressed in root, flower, pod wall, and grain. CarHsfB4b is specifically expressed in flower and grain [12]. In carnations, DcaHsfs-B2a, B2c, B4 are highly expressed in S, CA, YL, and ML. However, DcaHsf-B2a/B3a and DcaHsf B2b are highly expressed in FS2 and FS6, respectively (Figure 5). This indicates that members of the Class B DcaHsfs might be widely involved in the development of both vegetative and reproductive organs and tissues. The expression patterns of Hsf-C1 genes were diverse in different tissues. For example, the transcripts of Vitis pseudoreticulata VpHsfC1a remain at relatively lower levels (even undetectable) in roots, stems, leaves, and tendrils [41]. Similarly, SaHsfC1a is expressed at low levels in all the tested tissues in Sedum alfredii [42]. Carnation DcaHsf-C1 was down-regulated in almost all tested tissues (except for R) (Figure 5), which is consistent with the expression pattern of SaHsfC1b [42]. Specifically, it may be attributed to the fact that DcaHsf-C1 acts as a negative regulator in the development of organs.

3.5. DcaHsfs are Involved in Carnation Stress Response

The genome-wide expression profile analyses indicated that the majority of the Hsf genes are involved in heat, cold, drought, and salt stress responses [22,26]. Under heat or other stress conditions, plant Hsfs display diversity in patterns of expression [14]. In our study, all 17 DcaHsfs were found to be induced by a high temperature of 42 °C (Figure 6), which is in agreement with a previous study [43]. All DcaHsfs accumulated during drought treatment (Figure 6), and a previous study revealed that approximately 90% of sesame Hsfs are drought-responsive [44]. DcaHsf-A2a, A2b revealed to be strongly induced under heat stress conditions. This indicates that HsfA2 is a dominant regulator during the heat stress response in carnation, which is consistent with the studies of Arabidopsis, tomato, apple, Populus euphratica, and Phyllostachys edulis [11,22,38,39,45]. HsfA3 has been identified as an important player in the responses to heat, high salinity, and drought stresses in Solanum lycopersicum [46], whereas a similar function for HsfA3 is not detected in tomato [14]. In this study, DcaHsfA3 was also up-regulated in response to four analyzed abiotic stresses (heat, drought, ABA, and SA) (Figure 6). Group A4 Hsfs are involved in controlling reactive oxygen species homeostasis in plants, and group A5 Hsfs act as specific repressors of HsfA4 [47,48]. Fragaria vesca FvHsfA4a, A5a were both distinctly up-regulated in response to abiotic stresses such as cold, drought, and salt and hormone treatments (ABA, Eth, MeJA, and SA) [49]. Our data are highly similar, indicating DcaHsf-A4 accumulation during heat, drought, and SA treatment, as well as DcaHsf-A5 upregulation in response to cold, heat, drought, and salt and SA treatments (Figure 6). Arabidopsis Hsf-A9a is associated with ABA-mediated stress signaling and drought resistance [50]. Similarly, DcaHsfA9a was also induced in response to salt, ABA, and SA (Figure 6). Compared to Class A Hsfs, the members in Class B and C still have not been well studied. Arabidopsis AtHsfB1a and F. vesca FvHsfB1a were highly induced and accumulated in response to SA treatment [45,49,51]. Additionally, additional evidence demonstrated that AtHsfB1a, B2b are crucial components in primed defense gene activation and pathogen-induced acquired immune response [51]. Similarly, in this study, DcaHsf-B1, B2a, B2b accumulated at high levels at the later stage of SA treatment (Figure 6) Therefore, it is reasonable to speculate that DcaHsf-B1, B2a, B2b play a crucial role in the acquired immune response to pathogens. Additionally, DcaHsf-C1 acts as a positive regulator of heat shock proteins under heat stress conditions or PEG stress. The expression of rice OsHsfC1b was induced by salt, mannitol, and ABA [52]. In V. pseudoreticulata, VpHsfC1a was up-regulated in response to ABA treatment but significantly down-regulated during both MeJA and Eth treatments [41]. However, the expression of DcaHsf-C1 was not up-regulated by cold, ABA, or salt stress (Figure 6). We can speculate that DcaHsf-C1 might be involved in ABA-independent pathways in carnations. However, gene expression is a complex biological process, and more thorough studies are required to decipher the regulatory mechanisms.

4. Materials and Methods

4.1. Identification and Characterization of Hsf Genes in D. caryophyllus

The protein and nucleotide sequences of D. caryophyllus were downloaded from the carnation DB (http://carnation.kazusa.or.jp). The conserved domain of Hsf DBD (Pfam: PF00447) was submitted as a query in a BLASTP search of the D. caryophyllus proteome. The SMART 7 software (http://smart.embl-heidelberg.de/) was used to identify integrated DBD domain and (HR-A/B) domain in the putative Hsfs. Candidate proteins without integrated DBD domain and HR-A/B domain were removed. The ExPaSy-Protparam tool (https://www.expasy.org/tools/ProtParam.html) was used to analyze the physical properties of the predicted Hsf proteins.

4.2. Phylogenetic Analysis

Multiple sequence alignments of full-length Hsf proteins from D. caryophyllus and other three model species, i.e., A. thaliana, O. sativa, and P. trichocarpa, were performed using Clustal W2 (https://www.ebi.ac.uk/Tools/msa/clustalw2/, Dublin, Ireland). An unrooted neighbor-joining (NJ) phylogenetic tree was constructed using MEGA7.0 (Philadelphia, PA, U.S.A.) with 1000 bootstrap replicates. Distinctive names for each of the Hsfs identified in D. caryophyllus were given according to the classification of Hsfs in classes A, B, and C, referred to as DcaHsf genes.

4.3. Structural and Motif Analyses of DcaHsf Genes

The gene structures including exons and introns were displayed using Gene Structure Display Server (GSDS, http://gsds.cbi.pku.edu.cn/index.php, Beijing, China). The conserved motifs of DcaHsfs were defined by Multiple Em for Motif Elicitation (MEME, http://meme-suite.org/, U.S.A.) using the following parameters: number of repetitions = any, maximum number of motifs = 20, minimum width ≥10, maximum width ≤200, and only motifs with an E-value < 0.01 were retained for further analysis. NLS domains were predicted using cNLS Mapper software (http://nls-mapper.iab.keio.ac.jp/cgi-bin/NLS_Mapper_form.cgi, Tsuruoka, Japan). NES domains in the DcaHsfs were predicted with the NetNES 1.1 server software (http://www.cbs.dtudk/services/NetNES/, Lyngby Denmark).

4.4. Cis-acting Element Analysis of DcaHsfs

The 1500-bp sequence upstream from the initiation codon of each DcaHsf gene was obtained from the D. caryophyllus genome database. These sequences were used to identify cis-acting regulatory elements with the online program PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, Ghent, Belgium).

4.5. Plant Materials, Growth Conditions, and Stress Treatments

Tissue culture seedlings of carnation were grown in a chamber at Qingdao Agriculture University (Qingdao, China) under a 12 h light (300 μmol·m−2·s−1)/12 h dark cycle at 23–25 °C ambient temperature and 70% relative humidity. Various tissues, including the root (R), shoot (S), calyx (CA), young leaf (YL), mature leaf (ML), stigma (ST), ovary (OV), and flower petals (F), and six flowering stages (FS1, FS2, FS3, FS4, FS5, and FS6) were collected from the carnation seedlings. For abiotic stress and hormone treatments, the seedlings were treated at 42 °C (for heat stress), with 20% (w/v) polyethylene glycol (PEG) 6000 (for drought stress), 200 mM NaCl (for salt stress), 100 μM ABA, or 100 μM SA. The first or second tender leaves of the seedlings were collected at 0, 1, 6, and 12 h, immediately frozen in liquid nitrogen, and then stored at −80 °C for further analysis. Three biological replicates were performed for each sample.

4.6. RNA Isolation and Expression Analysis of DcaHsf Genes

Total RNA from carnation leaves was extracted using the Plant RNA Kit (Omega, Norcross, GA, USA) according to the instructions. Subsequently, 500 ng of total RNA was reverse-transcribed to first-strand cDNA by the PrimeScrip RT reagent Kit with gDNA Eraser (TaKaRa, Dalian, China) according to the manufacturer’s protocol, and the cDNA was diluted 10-fold for quantitative real-time PCR (qRT-PCR). qRT-PCR was performed using 2 μL of cDNA in a 20 μL reaction volume with SYBR® Premix Ex Taq™ II (TaKaRa, Dalian, China) on a StepOnePlus Real-Time PCR System (ABI, USA), using the following PCR program: 95 °C for 3 min, followed by 40 cycles at 95 °C for 30 s and at 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 were designed by Primer Premier 5.0 (Table S4). The GAPDH gene was used as an internal control gene. Three independent biological replicates were performed, and the relative expression levels of the DcaHsf genes were calculated with the 2−∆∆Ct method [53].

5. Conclusions

In this study, 17 DcaHsf genes were identified in the carnation genome for the first time. Comprehensive analyses of these genes, including phylogeny, genes structure, conserved motifs, and expression profiles in various tissues and under abiotic stresses were performed. Structural characteristics and comparisons with A. thaliana, O. sativa, and P. trichocarpa assisted in classifying these genes into three major classes (A, B, and C), with members of class A being the most abundant. The DcaHsf members were expressed in at least one tissue among root, stem, calyx, young leaf, mature leaf, stigma, ovary, and flower. In addition, DcaHsfA1, A2a, A9a, B2a, B3a revealed their putative involvement in the early flowering stages. The results of qRT-PCR revealed that all DcaHsfs responded to heat and drought, and many DcaHsfs were also regulated by cold, salt, and osmotic stress, as well as by the phytohormones ABA and SA. Our research suggests that DcaHsf A2a, A2b may be used as candidate genes for the breeding of heat-resistant carnation. Meanwhile, DcaHsf-B2a, B3a and DcaHsfA5, B2b could be considered as probable candidate genes for promoting early blooming and prolonging florescence in carnations.

Supplementary Materials

Author Contributions

J.Z. conceived and designed the research. W.L., X.-L.W., and J.-Y.Y. performed the experiments. K.-L.W. contributed with valuable discussions. W.L. wrote the original draft. J.Z. revised the manuscript. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ABAAbscisic acid
ABREABA-responsive element
APXAscorbate peroxidase
CACalyx
CATCatalase
DBDDNA-binding domains
EREEthylene-responsive element
FFlower
FSFlowering stages
HSEHeat shock element
HsfHeat shock transcription factor
MEMEMultiple Em for Motif Elicitation
MLMature leaf
MWMolecular weight
NESNuclear export signals
NLSNuclear localization signals
OVOvary
PEGPolyethylene glycol
qRT-PCRQuantitative real-time PCR
RRoot
RDRepressor domains
ROSReactive oxygen species
SStem
SASalicylic acid
STStigma
TFsTranscription factors
YLYoung leaf

References

  1. Hu, H.; Xiong, L. Genetic engineering and breeding of drought-resistant crops. Annu. Rev. Plant Biol. 2014, 65, 715–741. [Google Scholar] [CrossRef] [PubMed]
  2. Pereira, A. Plant abiotic stress challenges from the changing environment. Front. Plant Sci. 2016, 7, 1123. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, J.K. Abiotic stress signaling and responses in plants. Cell 2016, 167, 313–324. [Google Scholar] [CrossRef] [PubMed]
  4. Ohama, N.; Sato, H.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Transcriptional regulatory network of plant heat stress response. Trends Plant Sci. 2017, 22, 53–65. [Google Scholar] [CrossRef]
  5. Rushton, P.J.; Somssich, I.E.; Ringler, P.; Shen, Q.J. WRKY transcription factors. Trends Plant Sci. 2010, 15, 247–258. [Google Scholar] [CrossRef]
  6. Puranik, S.; Sahu, P.P.; Srivastava, P.S.; Prasad, M. NAC proteins: Regulation and role in stress tolerance. Trends Plant Sci. 2012, 17, 369–381. [Google Scholar] [CrossRef]
  7. Scharf, K.D.; Berberich, T.; Ebersberger, I.; Nover, L. The plant heat stress transcription factor (Hsf) family: Structure, function and evolution. Biochim. Biophys. Acta Gene Regul. Mech. 2012, 1819, 104–119. [Google Scholar] [CrossRef]
  8. Kotak, S.; Larkindale, J.; Lee, U.; Koskull-Döring, P.V.; Vierling, E.; Scharf, K.D. Complexity of the heat stress response in plants. Curr. Opin. Plant Biol. 2007, 10, 310–316. [Google Scholar] [CrossRef]
  9. Swindell, W.R.; Huebner, M.; Weber, A.P. Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genom. 2007, 8, 125. [Google Scholar] [CrossRef]
  10. Wei, Y.X.; Hu, W.; Xia, F.Y.; Zeng, H.Q.; Li, X.L.; Yan, Y.; He, C.Z.; Shi, H.T. Heat shock transcription factors in banana: Genome-wide characterization and expression profile analysis during development and stress response. Sci. Rep. 2016, 6, 36864. [Google Scholar] [CrossRef]
  11. Zhang, J.; Jia, H.X.; Li, J.B.; Li, Y.; Lu, M.Z.; Hu, J.J. Molecular evolution and expression divergence of the Populus euphratica Hsf genes provide insight into the stress acclimation of desert poplar. Sci. Rep. 2016, 6, 30050. [Google Scholar] [CrossRef] [PubMed]
  12. Chidambaranathan, P.; Jagannadham, P.T.K.; Satheesh, V.; Kohli, D.; Basavarajappa, S.H.; Chellapilla, B.; Kumar, J.; Jain, P.K.; Srinivasan, R. Genome-wide analysis identifies chickpea (Cicer arietinum) heat stress transcription factors (Hsfs) responsive to heat stress at the pod development stage. J. Plant Res. 2017, 131, 525–542. [Google Scholar] [CrossRef] [PubMed]
  13. Guo, M.; Liu, J.H.; Ma, X.; Luo, D.X.; Gong, Z.H.; Lu, M.H. The plant Heat Stress Transcription Factors (HSFs): Structure, regulation, and gunction in response to abiotic stresses. Front. Plant Sci. 2016, 7, 114. [Google Scholar] [CrossRef] [PubMed]
  14. Koskull-DöRing, P.V.; Scharf, K.D.; Nover, L. The diversity of plant heat stress transcription factors. Trends Plant Sci. 2007, 12, 452–457. [Google Scholar] [CrossRef]
  15. Singh, A.; Mittal, D.; Lavania, D.; Agarwal, M.; Mishra, R.C.; Grover, A. OsHsfA2c and OsHsfB4b are involved in the transcriptional regulation of cytoplasmic OsClpB (Hsp100) gene in rice (Oryza sativa L.). Cell Stress Chaperones 2012, 17, 243–254. [Google Scholar] [CrossRef]
  16. Cheng, Q.; Zhou, Y.; Liu, Z.; Zhang, L.; Song, G.; Guo, Z.; Wang, W.; Qu, X.; Zhu, Y.; Yang, D. An alternatively spliced heat shock transcription factor, OsHSFA2dI, functions in the heat stress-induced unfolded protein response in rice. Plant Biol. 2015, 17, 419–429. [Google Scholar] [CrossRef]
  17. Giesguth, M.; Sahm, A.; Simon, S.; Dietz, K.J. Redox-dependent translocation of the heat shock transcription factor AtHSFA8 from the cytosol to the nucleus in Arabidopsis thaliana. FEBS Lett. 2015, 589, 718–725. [Google Scholar] [CrossRef]
  18. Guo, J.K.; Wu, J.; Ji, Q.; Wang, C.; Luo, L.; Yuan, Y.; Wang, Y.H.; Wang, J. Genome-wide analysis of heat shock transcription factor families in rice and Arabidopsis. J. Genet. Genom. 2008, 35, 105–118. [Google Scholar] [CrossRef]
  19. Wang, F.M.; Dong, Q.; Jiang, H.Y.; Zhu, S.W.; Chen, B.J.; Xiang, Y. Genome-wide analysis of the heat shock transcription factors in Populus trichocarpa and Medicago truncatula. Mol. Biol. Rep. 2012, 39, 1877–1886. [Google Scholar] [CrossRef]
  20. Chung, E.; Kim, K.M.; Lee, J.H. Genome-wide analysis and molecular characterization of heat shock transcription factor family in Glycine max. J. Genet. Genom. 2013, 40, 127–135. [Google Scholar] [CrossRef]
  21. Lin, Y.X.; Jiang, H.Y.; Chu, Z.X.; Tang, X.L.; Zhu, S.W.; Cheng, B.J. Genome-wide identification, classification and analysis of heat shock transcription factor family in maize. BMC Genom. 2011, 12, 76. [Google Scholar] [CrossRef] [PubMed]
  22. Giorno, F.; Guerriero, G.; Baric, S.; Mariani, C. Heat shock transcriptional factors in Malus domestica: Identification, classification and expression analysis. BMC Genom. 2012, 13, 639. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, S.J.; Zhang, P.; Jing, Z.; Shi, J.L. Genome-wide identification and analysis of heat shock transcription factor family in cucumber (Cucumis sativus L.). Plant Omics 2013, 6, 449–455. [Google Scholar]
  24. Zhang, J.; Liu, B.B.; Li, J.B.; Zhang, L.; Wang, Y.; Zheng, H.Q.; Lu, M.Z.; 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]
  25. Wang, J.; Sun, N.; Deng, T.; Zhang, L.D.; Zuo, K.J. Genome-wide cloning, identification, classification and functional analysis of cotton heat shock transcription factors in cotton (Gossypium hirsutum). BMC Genom. 2014, 15, 961. [Google Scholar] [CrossRef]
  26. Xue, G.P.; Sadat, S.; Drenth, J.; Mcintyre, C.L. The heat shock factor family from Triticum aestivum in response to heat and other major abiotic stresses and their role in regulation of heat shock protein genes. J. Exp. Bot. 2014, 65, 539–557. [Google Scholar] [CrossRef]
  27. Zhu, X.Y.; Huang, C.Q.; Zhang, L.; Liu, H.F.; Yu, J.H.; Hu, Z.Y.; Hua, W. Systematic analysis of Hsf family genes in the Brassica napus genome reveals novel responses to heat, drought and high CO2 stresses. Front. Plant Sci. 2017, 8, 1174. [Google Scholar] [CrossRef]
  28. Boxriker, M.; Boehm, R.; Krezdorn, N.; Rotter, B.; Piepho, H.P. Comparative transcriptome analysis of vase life and carnation type in Dianthus caryophyllus L. Sci. Hortic. 2017, 217, 61–72. [Google Scholar] [CrossRef]
  29. Onozaki, T.; Ikeda, H.; Yamaguchi, T. Genetic improvement of vase life of carnation flowers by crossing and selection. Sci. Hortic. 2001, 87, 107–120. [Google Scholar] [CrossRef]
  30. Muneer, S.; Soundararajan, P.; Jeong, B.R. Proteomic and antioxidant analysis elucidates the underlying mechanism of tolerance to hyperhydricity stress in in vitro shoot cultures of Dianthus caryophyllus. J. Plant Growth Regul. 2016, 35, 667–679. [Google Scholar] [CrossRef]
  31. Yagi, M.; Kosugi, S.; Hirakawa, H.; Ohmiya, A.; Tanase, K.; Harada, T.; Kishimoto, K.; Nakayama, M.; Ichimura, K.; Onozaki, T.; et al. Sequence analysis of the genome of carnation (Dianthus caryophyllus L.). DNA Res. 2013, 21, 231–241. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, J.; Li, Y.; Jia, H.X.; Li, J.B.; Huang, J.; Lu, M.Z.; Hu, J.J. The heat shock factor gene family in Salix suchowensis: A genome-wide survey and expression profiling during development and abiotic stresses. Front. Plant Sci. 2015, 6, 748. [Google Scholar] [CrossRef] [PubMed]
  33. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; de Peer, Y.V.; 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]
  34. He, Y.; Li, W.; Lv, J.; Jia, Y.B.; Wang, M.C.; Xia, G.M. Ectopic expression of a wheat MYB transcription factor gene, TaMYB73, improves salinity stress tolerance in Arabidopsis thaliana. J. Exp. Bot. 2012, 63, 1511–1522. [Google Scholar] [CrossRef] [PubMed]
  35. Rose, A.B. Intron-mediated regulation of gene expression. Curr. Top. Microbiol. Immunol. 2008, 326, 277–290. [Google Scholar] [PubMed]
  36. Pelham, H.R.; Bienz, M. A synthetic heat-shock promoter element confers heat-inducibility on the herpes simplex virus thymidine kinase gene. EMBO J. 1981, 1, 1473–1477. [Google Scholar] [CrossRef]
  37. Tang, R.M.; Zhu, W.J.; Song, X.Y.; Lin, X.Y.; Cai, J.H.; Man, W.; Yang, Q. Genome-wide identification and function analyses of heat shock transcription factors in potato. Front. Plant Sci. 2016, 7, 490. [Google Scholar] [CrossRef]
  38. Xie, L.H.; Li, X.Y.; Hou, D.; Cheng, Z.C.; Liu, J.; Li, J.; Mu, S.H.; Gao, J. Genome-wide analysis and expression profiling of the heat shock factor gene family in Phyllostachys edulis during development and in response to abiotic stresses. Forests 2019, 10, 100. [Google Scholar] [CrossRef]
  39. Fragkostefanakis, S.; Röth, S.; Schleiff, E.; Scharf, K.D. Prospects of engineering thermotolerance in crops through modulation of heat stress transcription factor and heat shock protein networks. Plant Cell Environ. 2015, 38, 1881–1895. [Google Scholar] [CrossRef]
  40. Liu, B.; Hu, J.; Zhang, J. Evolutionary divergence of duplicated Hsf genes in Populus. Cells 2019, 8, 438. [Google Scholar] [CrossRef]
  41. Hu, Y.; Han, Y.T.; Zhang, K.; Zhao, F.L.; Li, Y.J.; Zheng, Y.; Wang, Y.J.; Wen, Y.Q. Identification and expression analysis of heat shock transcription factors in the wild Chinese grapevine (Vitis pseudoreticulata). Plant Physiol. Biochem. 2016, 99, 1–10. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, S.S.; Jiang, J.; Han, X.J.; Zhang, Y.X.; Zhuo, R.Y. Identification, expression analysis of the Hsf family, and characterization of class A4 in Sedum Alfredii Hance under cadmium stress. Int. J. Mol. Sci. 2018, 19, 1216. [Google Scholar] [CrossRef] [PubMed]
  43. Mittal, D.; Chakrabarti, S.; Sarkar, A.; Singh, A.; Grover, A. Heat shock factor gene family in rice: Genomic organization and transcript expression profiling in response to high temperature, low temperature and oxidative stresses. Plant Physiol. Biochem. 2009, 47, 785–795. [Google Scholar] [CrossRef] [PubMed]
  44. Dossa, K.; Diouf, D.; Cisse, N. Genome-wide investigation of Hsf genes in sesame reveals their segmental duplication expansion and their active role in drought stress response. Front. Plant Sci. 2016, 7, 1522. [Google Scholar] [CrossRef]
  45. Ikeda, M.; Mitsuda, N.; Ohme-Takagi, M. Arabidopsis HsfB1 and HsfB2b act as repressors of the expression of heat-inducible Hsfs but positively regulate the acquired thermotolerance. Plant Physiol. 2011, 157, 1243–1254. [Google Scholar] [CrossRef]
  46. Li, Z.J.; Zhang, L.L.; Wang, A.X.; Xu, X.Y.; Li, J.F. Ectopic overexpression of SlHsfA3, a heat stress transcription factor from tomato, confers increased thermotolerance and salt hypersensitivity in germination in transgenic Arabidopsis. PLoS ONE 2013, 8, e54880. [Google Scholar] [CrossRef]
  47. Yamanouchi, U.; Yano, M.; Lin, H.X.; Ashikari, M.; Yamada, K. A rice spotted leaf gene, Spl7, encodes a heat stress transcription factor protein. Proc. Natl. Acad. Sci. USA 2002, 99, 7530–7535. [Google Scholar] [CrossRef]
  48. Baniwal, S.K.; Chan, K.Y.; Scharf, K.D.; Nover, L. Role of heat stress transcription factor HsfA5 as specific repressor of HsfA4. J. Biol. Chem. 2007, 282, 3605–3613. [Google Scholar] [CrossRef]
  49. Hu, Y.; Han, Y.T.; Wei, W.; Li, Y.J.; Zhang, K.; Gao, Y.R.; Zhao, F.L.; Feng, J.Y. Identification, isolation, and expression analysis of heat shock transcription factors in the diploid woodland strawberry Fragaria vesca. Front. Plant Sci. 2015, 6, 736. [Google Scholar] [CrossRef]
  50. Schramm, F.; Ganguli, A.; Kiehlmann, E.; Englich, G.; Walch, D.; Koskull-Döring, P.V. The heat stress transcription factor HsfA2 serves as a regulatory amplifier of a subset of genes in the heat stress response in Arabidopsis. Plant Mol. Biol. 2006, 60, 759–772. [Google Scholar] [CrossRef]
  51. Pick, T.; Jaskiewicz, M.; Peterhänsel, C.; Conrath, U. Heat shock factor HsfB1 primes gene transcription and systemic acquired resistance in Arabidopsis. Plant Physiol. 2012, 159, 52–55. [Google Scholar] [CrossRef] [PubMed]
  52. Schmidt, R.; Schippers, J.H.M.; Welker, A.; Mieulet, D.; Guiderdoni, E.; Mueller-Roeber, B. Transcription factor OsHsfC1b regulates salt tolerance and development in Oryza sativa ssp. japonica. AoB Plants 2012, 2012, pls011. [Google Scholar] [CrossRef] [PubMed]
  53. Livak, K.J.; Schmittgenb, 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]
Figure 1. Scaffold locations of DcaHsfs. Bars represent the scaffolds, DcaHsfs are marked by redlines.
Figure 1. Scaffold locations of DcaHsfs. Bars represent the scaffolds, DcaHsfs are marked by redlines.
Ijms 20 05233 g001
Figure 2. The phylogenetic tree of Hsf proteins. The phylogenetic tree of Hsf proteins in carnation and other plant species was generated by MEGA 7 using the neighbor-joining method. Dca, D. caryophyllus; At, Arabidopsis thaliana; Os, Oryza sativa and Pt, Populus trichocarpa.
Figure 2. The phylogenetic tree of Hsf proteins. The phylogenetic tree of Hsf proteins in carnation and other plant species was generated by MEGA 7 using the neighbor-joining method. Dca, D. caryophyllus; At, Arabidopsis thaliana; Os, Oryza sativa and Pt, Populus trichocarpa.
Ijms 20 05233 g002
Figure 3. Intron and exon structure (a) and amino acid motifs (b) of members of the DcaHsf family. (a) Boxes and lines represent exons and introns, respectively. (b) A total of 20 conserved motifs were identified using Multiple Em for Motif Elicitation (MEME).
Figure 3. Intron and exon structure (a) and amino acid motifs (b) of members of the DcaHsf family. (a) Boxes and lines represent exons and introns, respectively. (b) A total of 20 conserved motifs were identified using Multiple Em for Motif Elicitation (MEME).
Ijms 20 05233 g003
Figure 4. Cis-regulatory elements in the promoter region of DcaHsfs. (a) Number of each cis-acting element in the promoter region (1.5 kb upstream of the translation initiation site) of DcaHsfs. Statistics of the total number of DcaHsfs including the corresponding cis-acting elements (red dot) and the total number of cis-acting elements in the DcaHsf gene family (gray box). (b) Frequency of the cis-acting elements in each gene. Based on the functional annotation, the cis-acting elements were classified into three major classes: stress-, hormone-, and development-related cis-acting elements.
Figure 4. Cis-regulatory elements in the promoter region of DcaHsfs. (a) Number of each cis-acting element in the promoter region (1.5 kb upstream of the translation initiation site) of DcaHsfs. Statistics of the total number of DcaHsfs including the corresponding cis-acting elements (red dot) and the total number of cis-acting elements in the DcaHsf gene family (gray box). (b) Frequency of the cis-acting elements in each gene. Based on the functional annotation, the cis-acting elements were classified into three major classes: stress-, hormone-, and development-related cis-acting elements.
Ijms 20 05233 g004
Figure 5. The expression levels of DcaHsfs in different tissues and flowering stages. The different colors correspond to log2-transformed fold change, green indicates down-regulation, and red represents up-regulation.
Figure 5. The expression levels of DcaHsfs in different tissues and flowering stages. The different colors correspond to log2-transformed fold change, green indicates down-regulation, and red represents up-regulation.
Ijms 20 05233 g005
Figure 6. Expression levels of DcaHsfs under various abiotic stresses, determined by qRT-PCR. The different colors correspond to log2-transformed fold change, green indicates down-regulation, and red represents up-regulation.
Figure 6. Expression levels of DcaHsfs under various abiotic stresses, determined by qRT-PCR. The different colors correspond to log2-transformed fold change, green indicates down-regulation, and red represents up-regulation.
Ijms 20 05233 g006
Table 1. Summary information of Dianthus caryophyllus heat shock transcription factors (DcaHsfs) in carnation. Notes: I.I., instability index; Stability, (U: unstable protein, S: stable protein); A.I., aliphatic index; n.c.r., total number of negatively charged residues (Asp + Glu); p.c.r., total number of positively charged residues (Arg + Lys); GRAVY, grand average of hydropathicity; pI, isoelectric point; MW, molecular weight.
Table 1. Summary information of Dianthus caryophyllus heat shock transcription factors (DcaHsfs) in carnation. Notes: I.I., instability index; Stability, (U: unstable protein, S: stable protein); A.I., aliphatic index; n.c.r., total number of negatively charged residues (Asp + Glu); p.c.r., total number of positively charged residues (Arg + Lys); GRAVY, grand average of hydropathicity; pI, isoelectric point; MW, molecular weight.
Protein NameGene IDSubfamilySizeI.I.StabilityA.I.n.c.r. (%)p.c.r.(%)GRAVYpIMW (kDa)
DcaHsf-A1Dca57201.1A148855.97U66.726947−0.6904.7454.49
DcaHsf-A2aDca14360.1A238059.24U83.555945−0.5035.1242.96
DcaHsf-A2bDca52568.1A235958.97U80.035743−0.5485.0540.55
DcaHsf-A3Dca41810.1A324464.58U67.336750−0.5664.9851.75
DcaHsf-A4Dca23163.1A439047.75U71.955947−0.8595.744.99
DcaHsf-A5Dca19769.1A548948.53U72.176752−0.7455.4554.52
DcaHsf-A7Dca4574.1A742548.11U67.366361−0.7716.7248.91
DcaHsf-A9aDca9629.1A940147.92U69.736746−0.691546.12
DcaHsf-A9bDca41703.1A933136.63S69.434544−0.8056.2338.21
DcaHsf-B1Dca60410.1B127634.47S68.773839−0.8047.6131.0
DcaHsf-B2aDca22545.1B233754.72U61.643529−0.721636.48
DcaHsf-B2bDca48996.1B233760.89U60.983932−0.709636.48
DcaHsf-B2cDca54105.1B231861.35U68.573332−0.5925.9133.80
DcaHsf-B3aDca44175.1B324448.80U70.333137−0.7638.8828.38
DcaHsf-B3bDca48010.1B345748.40U68.963137−0.7848.8828.33
DcaHsf-B4Dca39623.1B428754.07U72.303032−0.8748.5233.72
DcaHsf-C1 (incomplete)Dca24054.1C113358.77U76.242122−0.8818.0115.89
Table 2. Analysis and distribution of conserved motifs in carnation DcaHsfs.
Table 2. Analysis and distribution of conserved motifs in carnation DcaHsfs.
MotifE-ValueWidthBest Possible Match
Motif12.70 × 10−42641VWDPAEFARDLLPRYFKHNNFSSFVRQLNTYGFRKVVPDRW
Motif25.20 × 10−22529PFLTKTYDMVDDPSTDDIVSWSEDGTSFV
Motif31.10 × 10−18729EFANEGFLRGQKHLLKNIKRRKTTTAHSQ
Motif44.20 × 10−11541GLEGENERLRRENEVLMSELVKLKQQQQNTFSLLQAMESRL
Motif54.80 × 10−5434QSTEWKQKQMMTFLAKAMQNPTFVQQLVQKKDER
Motif61.60 × 10−2249PPQQQPSTAAPTNSSDEQVISNSNSPPLAIPSVIMHRHHHHHHLYHNNN
Motif71.20 × 10−1941KSVKAIRVSMKRRLTSTLSAPNLNDVVEPELVRSMAVSSDN
Motif82.60 × 10−1650CGGGGGGSPMIFGVSIGGKRGREGGDDGGGEVVGGGEGLGATEVHDDHMH
Motif92.60 × 10−1350ATLDNGTDGDIKEQKVDDSMPPEIDTNVGDVSQTSWEELLWAEDEEGFRQ
Motif107.70 × 10−1029MRELSIKGLFDDHDDDDECGIIMRRKMTK
Motif111.50 × 10−913PKPMEGLNEMNPP
Motif121.80 × 10−541DSDGDDGNNKNRPKLFGVRLDLQDESERKRRKKLALDYTRT
Motif135.20 × 10−821NNNNNNNVVITRKNNENEMNN
Motif141.40 × 10−429PVENVVPESGNWGEDVEDLIEQLGFLGPM
Motif151.80 × 10−421MTAVLVTVSDLVSSSTTSSSS
Motif162.10 × 10−429MSPPPSPPAEEKPEKLTAVVVGGGGGETQ
Motif171.70 × 10−321AAPSRVNDAVWTQLLTLPRGS
Motif182.90 × 10−310GFRKVDPDKW
Motif191.30 × 10−223TSTTCTCTPLSTESPQLGLQLSP
Motif201.00 × 10−219YWYDFDGEDEVELEERVPC
Table 3. Functional domains of DcaHsfs.
Table 3. Functional domains of DcaHsfs.
Gene Name DBDHR-A/BNLSNESRDAHA
DcaHsf-A131–124162–182/201–212N.D.(403) LN.D.N.D.
DcaHsf-A2a40–154183–201/222–233(147–156) KTIKRRRNVT
(258–267) AGMKRRLTST
N.D.N.D.(329–338) QTSWEELLWA
DcaHsf-A2b40–133162–180/201–212(126–137) LLKTIKRRRNVT
(237–246) AGMKRRLTST
(232–237) LDITHLN.D.(308–317) QTSWEELLWA
DcaHsf-A337–148181–199/220–230N.D.(319) LN.D.N.D.
DcaHsf-A411–104139–157/178–189(204–213) HDRKRRFSRPN.D.N.D.(325–334) DVFWEQFLTE
DcaHsf-A520–113138–156/177–187(206–217) LSAYNKKRRLPPN.D.ND(438–447) DLFWEQFLTE
DcaHsf-A740–133164–182/203–213(126–140) LLKNIKRRKNPSQTF
(237–246) LSKKRRRPIE
N.D.ND(322–331) DDFWEDLLNE
DcaHsf-A9a87–180202–220/241–251(173–184) LLKSIKRKRHGS
(275–285) RVSKKRRLAST
(204) LDQEALKVEINDN.D.
DcaHsf-A9b95–188210–228/249–259(181–192) LLKSIKRKRHGS
(283–293) RVSKKRRLAST
(217–221) LKVEINDN.D.
DcaHsf-B16–99125–131N.D.(155–157) LEL(220-226) KLFGVWLN.D.
DcaHsf-B2a32–125151–169/190–200N.D.N.D.NDN.D.
DcaHsf-B2b32–125151–169/190–200N.D.N.D.NDN.D.
DcaHsf-B2c26–119145–153/172–178N.D.(197–202) KENMSL(284–290) KLFGVSIN.D.
DcaHsf-B3a29–122147–152(5–30) SIKGLFDDHDDDDECGIIMRRKMTKP
(178–187) NAMKRKCQEL
(207–235) KNRPKLFGVRLDLQDESERKRRKKLALDY
(222–238) LDLQDESERKRRKKLAL(216–222) KLFGVRLN.D.
DcaHsf-B3b29–122147–152(5–30) SIKGLFDDHDDDDECGIIMRRKMTKP
(178–187) NAMKRKCQEL
(207–235) KNRPKLFGVRLDLQDESERKRRKKLALDY
N.D.(216–222) KLFGVRLN.D.
DcaHsf-B412–105131–149/163–166(275–283) HSKKRLHLAN.D.(268–274) RLFGVPLN.D.
DcaHsf-C1 (not full)1–4473–89/98–108N.D.(66) LN.D.N.D.
N.D., not detected.

Share and Cite

MDPI and ACS Style

Li, W.; Wan, X.-L.; Yu, J.-Y.; Wang, K.-L.; Zhang, J. Genome-Wide Identification, Classification, and Expression Analysis of the Hsf Gene Family in Carnation (Dianthus caryophyllus). Int. J. Mol. Sci. 2019, 20, 5233. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20205233

AMA Style

Li W, Wan X-L, Yu J-Y, Wang K-L, Zhang J. Genome-Wide Identification, Classification, and Expression Analysis of the Hsf Gene Family in Carnation (Dianthus caryophyllus). International Journal of Molecular Sciences. 2019; 20(20):5233. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20205233

Chicago/Turabian Style

Li, Wei, Xue-Li Wan, Jia-Yu Yu, Kui-Ling Wang, and Jin Zhang. 2019. "Genome-Wide Identification, Classification, and Expression Analysis of the Hsf Gene Family in Carnation (Dianthus caryophyllus)" International Journal of Molecular Sciences 20, no. 20: 5233. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20205233

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