1. Introduction
Heat shock factors (HSFs) are well known for their functions as the main regulator of heat stress response [1]. Plant HSFs are encoded by a multi-gene family, and their protein structure is highly conserved and generally contains five core domains: DNA binding domain (DBD), oligomerization domain (OD), nuclear localization signal (NLS), nuclear export signal (NES), and C-terminal short activator peptide motif (AHA motif) [2]. Three evolutionary-conserved classes of HSFs (classes A, B, and C) are divided primarily based on the structural features of OD [3]. DBD is located close to the N-terminus of all HSFs, and consists of three α-helices and four anti-parallel β-sheets. The hydrophobic core of DBD is a helix-turn-helix motif (H2-T-H3), the functional region for transcriptional activation to recognize and bind highly conserved heat stress elements (HSE: 5′-nGAAnnTTCn-3′) [3,4]. Though binding to the HSE motif, plant HSFs can regulate the transcription of potential downstream genes, including, but not limited to, classical heat shock protein genes (HSPs) [5,6,7]. In addition, Arabidopsis HsfA1a also has a high affinity to stress responsive element (STRE: 5′-AGGGG-3′) and modified HSE motifs (gap-type and TTC-rich-type) in the promoters of target genes [5,7].
Plant HSFs play essential roles in both heat-stress defenses and other abiotic stress tolerance in an HSPs-dependent manner [8,9]. HSFs have also shown their versatile roles in biotic stress, such as powdery mildew infection [10]. Besides their roles in stress responses, HSFs are also involved in plant growth and development [11,12]. Plenty of research has revealed that even under normal growth conditions, plant HSFs participate in organ development, secondary mechanisms, and other life activities. For example, MtHSFA9 from Medicago truncatula regulates the depth of dormancy seeds and the biosynthesis of flavonoid, carotenoid, and other secondary metabolites [13]. HSFs have also been found to perform an important function in flower development in lettuce (Lactuca sativa), Asteraceae [6]; moso bamboo (Phyllostachys edulis), Poaceae [11]; carnation (Dianthus caryophyllus), Caryophyllaceae [12]; and wintersweet (Chimonanthus praecox), Calycanthaceae [14]. What is more, in Narcissus tazetta, 10 members from HSF transcription factor family are highly expressed in corona at full-opening stage, when the emission of floral scent (mainly including phenethyl acetate and benzyl acetate) reaches a maximum. The transcriptome-based WGCNA analysis indicates that those HSFs are possibly involved in the regulation of biosynthesis and the emission of benzenoid/phenylpropanoid metabolites in Narcissus tazetta flowers [15]. Up until now, the regulation of HSFs regarding floral scent in other species has not been clear.
Petunia (Petunia hybrida) is a kind of ornamental plant widely used as a garden and potted flower. Some petunia cultivars, such as “Mitchell”, emit a strong fragrance at night [16], which is dominated by benzenoid/phenylpropanoid compounds derived from the amino acid phenylalanine [17,18]. As high-quality whole genome sequences of petunia wild parents (P. axillaris and P. inflata) have been published [19], it was easily found that HSF-binding elements HSE and/or STRE exist in the promoters of most benzenoid biosynthetic genes from their genomes (SNG,
PAL catalyzes the first committed step in the phenylpropanoid pathway. Considering the contribution of PAL to the regulation of the carbon flux from the amino acid phenylalanine to downstream products of the secondary metabolism, PAL genes were predicted to be the key targets of many transcription factors, such as ODORANT1 (ODO1), EMISSION OF BENZENOIDS I (EOBI), and EOBII, which are positive regulators of benzenoid/phenylpropanoid biosynthesis and emissions [16,20,21]. PAL1 (GenBank accession number AY705976) and PAL2 (CO805160), two gene members of petunia PALs, have been proven to be target genes of positive transcription factors ODORANT1 (ODO1) and EMISSION OF BENZENOIDS I (EOBI) [16,21]. Additionally, PAL2 can be regulated by EOBII, but PAL1 cannot [20].
The present study was carried out to prove the hypothesis of whether HSF genes have a potential function in the regulation of benzenoid/phenylpropanoid biosynthesis in petunia. First, we systematically identified and characterized the petunia HSF genes. Based on transcriptomic data and qRT-PCR analysis, PhHSF19 was found as a potential effector involved in petunia benzenoid/phenylpropanoid biosynthesis. Furthermore, PhHSF19 was proven to positively regulate the expression of PhPAL2. Our results will provide valuable additional information for the functional analysis of HSF genes in petunia.
2. Results
2.1. Identification of HSF Gene Family Members in Petunia
Through domain verification, 34 HSF gene family members were obtained in petunia. The length of the 34 HSF proteins varied from 69 (Peaxi162Scf00205g00820) to 512 (Peaxi162Scf00045g00131) amino acids, while the molecular weight ranged from 8.02 kD to 56.56 kD (Table 1). The theoretical isoelectric point was between 4.32 to 9.42. Furthermore, except one HSF protein (Peaxi162Scf00205g00820), the total average hydrophobic index of other HSF proteins was negative, suggesting that most HSF proteins belonged to hydrophilic proteins. The aliphatic index is the relative content of aliphatic side chains, including alanine, valine, isoleucine, and leucine. The aliphatic index of HSF proteins was from 54.71 (Peaxi162Scf00450g00422) to 94.49 (Peaxi162Scf00205g00820). All HSFs were predicted to be located in the nucleus.
2.2. Gene Structure, Phylogeny Analysis, and Conserved Motifs of Petunia HSFs
The intron and exon structures were analyzed using GSDS 2.0 to reveal the gene structural features of petunia HSFs (Figure 1). Out of 34 petunia HSFs, 28 genes contained two exons and one intron. The other six genes had three to five exons. Peaxi162Scf00045g00131 had the longest sequence (6430 bp), while Peaxi162Scf00205g00820 had the shortest (477 bp).
Based on well-established Arabidopsis HSF family classification [22], petunia HSFs were also divided into three major subfamilies with 21 HSF proteins from A. thaliana, HSFA (red), HSFB (blue), and HSFC (green), which is shown in Figure 2. The HSFA subfamily is the largest with 23 petunia HSFs, while there are only two petunia HSFs in the HSFC subfamily. Among the HSFA subfamily, eight HSFs belong to class A9, followed by four HSFs in classes A1. However, none of the petunia HSFs belongs to class A7.
Fifteen conserved motifs in petunia HSFs were identified and analyzed using the online software MEME (Figure 3). Motif 1, motif 2, and motif 3 in the N-terminals contained 37, 41, and 29 amino acids, respectively. Motifs 1–3 constituted the DBD domain (~100aa) and were widely distributed and highly conserved in most petunia HSFs. In addition, 20 out of 34 HSF proteins had motif 13, referring to the NLS motif. In addition, motif 4 was observed in the HSFs from subfamily A and subfamily C, but was not included in subfamily B. Motif 9 was only present in subfamily B. Motif 5 was found in HSFs from subfamily A and subfamily B, rather than subfamily C. AHA motif (motif 10) and NES motif (motif 11) only existed in subfamily A.
2.3. Cis-Acting Elements Analysis of Petunia HSFs
The PlantCARE online tool was used to analyze the promoter cis-acting elements of 34 petunia HSFs, and the core cis-acting elements in each promoter are shown in Figure 4. It was predicted that promoters of petunia HSFs contained a variety of hormone-response cis-acting elements, such as ABA-response element ABRE, methyl jasmonate response element TGACG, salicylic acid response element TCA, and ethylene response element ERE. Several well-known stress-related elements, including HSE, STRE, TC-rich, LTR-motif, ARE, and MBS, were found in the promoters of 34 petunia HSFs. In addition, promoters of petunia HSFs also possessed light stress response elements (Box 4, G-box, and Spl), wounding responsive element WUN, and circadian element. Among these mentioned elements, 28 out of 34 petunia HSF promoters contained STRE motifs, the highest number of cis-acting elements. There were six STRE motifs in the promoter of Peaxi162Scf00560g00011, followed by Peaxi162Scf00170g01112 and Peaxi162Scf00316g00518 with five STRE motifs. Except for the STRE motif, the other three cis-acting elements, ARE, Box 4, and HSE, had a higher frequency in the promoters of 34 HSFs. There were four HSE motifs in the promoters of Peaxi162Scf00004g00431, Peaxi162Scf00327g00413, and Peaxi162Scf00461g00026.
2.4. Selection of Potential Key Petunia HSF Genes Based on RNA Sequencing and qRT-PCR Analysis
Based on the reported RNA sequencing data of the day 0 (bud) or day 2 (post-anthesis) corolla tissue of P. hybrida “Mitchell” [23], 33 assembled sequences could be mapped to the HSF genes of Petunia axillaris, and were named PhHSF01-PhHSF33 (Table S1). Among them, six PhHSFs with a fold-change of day 2 vs. day 0 ≥ 2 and the expression level (RPKM) at day 2 ≥ 1000, including PhHSF08, PhHSF09, PhHSF13, PhHSF19, PhHSF20, and PhHSF25, were predicted to have potential roles in regulating benzenoid/phenylpropanoid biosynthesis (Table S1). The expression levels of six potential key PhHSF genes in different tissues and during the flower opening period were detected by qRT-PCR. The results showed that PhHSF19 is highly expressed in petal tubes, followed by stems and petal limbs (Figure 5a). However, the expression level of the other five PhHSF genes was much lower in the petal limbs—the main tissues emitting volatile benzenoid/phenylpropanoid in petunia [24]. During flower development, the expression level of PhHSF19 in D1 increased 19.4-fold compared to that in D0 (Figure 5b). It gradually decreased from D1 to D4, but dramatically increased in D5, and then sharply decreased. Taking the expression patterns of PhHSF19 in different tissues and at different flower development into consideration, PhHSF19 probably had a positive effect on volatile benzenoid/phenylpropanoid biosynthesis and was chosen for further biological function analysis.
2.5. Y1H and Dual Luciferase Assays
When all kinds of benzenoids/phenylpropanoids accumulated to high levels during early flower opening period [25], the expression level of PhPAL2 dramatically increased, while that of PhPAL1 changed slightly (Figure 6a), revealing that the PhPAL2 gene would play a more important role in volatile benzenoid/phenylpropanoid biosynthesis, compared with PhPAL1. The expression pattern of PhHSF19 was positively related with PhPAL2 (R2 = 0.8423), rather than PhPAL1 (Figure 6b). There were three STRE elements in the promoter of PhPAL2 located at nucleotides −1951, −750, and −172 from the initial ATG codon (Figure 6c). We first performed the Y1H assay to address whether PhHSF19 could regulate the expression of PhPAL2. Only co-transformed yeasts of pGADT7-Rec2-PhHSF19 and pAbAi-PhPAL2-Pro could survive and grow normally on the selective medium of SD/-Leu/100 ng mL−1 AbA, but the control did not survive (Figure 6d), revealing that PhHSF19 can interact with PhPAL2 promoter in the yeast system. In addition, a dual-luciferase assay was employed to further verify the regulation of PhHSF19 on the PhPAL2 promoter. As shown in Figure 6e, a stronger fluorescent signal was detected in the tobacco leaf co-transformed with pCNHP-PhHSF19 plus pGreenII-PhPAL2-Pro-LUC compared to the control, and the LUC/REN ratio of tobacco co-transformed with the effector containing PhHSF19 and the reporter containing promoter of PhPAL2 were drastically enhanced with an approximate 2.8-fold induction.
3. Discussion
HSFs act as essential regulators in various abiotic stresses [1]. There are 21 HSF gene members in Arabidopsis [22], 25 members in rice [22], 31 members in Populus trichocarpa [26], and 17 members in Dianthus caryophyllus [12]. In this study, a total of 34 HSFs were identified in the petunia genome, which is greater than the number of HSF gene members in the four above-mentioned plant species. Similar to other species, most petunia HSFs belong to the HSFA subfamily. There are eight petunia HSF members in class A9, while only one or two members of the same class are found in Arabidopsis, P. trichocarpa, and D. caryophyllus [12,22,26]. Additionally, class A7 is absent in petunia. These results suggest that gene loss and gene duplication of HSFs occur in the evolution process of petunia.
Besides motifs 1–3 (DBD domain) and motif 13 (NLS motif), which are in almost all petunia HSFs, some conserved motifs are present in specific HSF subfamilies. AHA motif (motif 10), which functions as a transcription activator [1], only exists in subfamily A, suggesting that HSF members from subfamily B and subfamily C might probably lack transcriptional activation function. Similar to the AHA motif, the NES motif (motif 11) is also only found in the HSFs from subfamily A. OD, called the HR-A/B region, can be used to distinguish subfamily B from subfamilies A and C [1], as this domain in petunia is longer (as motif 4) in subfamilies A and C, but is less (as motif 9) in subfamily B, which is more similar to other species [27,28,29].
In addition to their roles in stress responses [30,31], HSFs are also involved in plant development [6,32] and secondary mechanism [5,15]. Several HSFs in Narcissus tazetta are predicted to have an effect on flower scent emission [15]. HSF-binding elements HSE and/or STRE are present in the promoters of many benzenoid biosynthetic genes, which advances the hypothesis that HSFs probably play a vital part in the regulation of benzenoid/phenylpropanoid biosynthesis and emissions in petunia. A previous study revealed that flowers at one day or two days post-anthesis emitted the greatest amount of 12 benzenoids/phenylpropanoids, such as methyl benzoate, benzaldehyde, benzyl benzoate, phenylacetaldehyde, and isoeugenol [25]. PhHSF19, with an up-regulated expression during this period based on RNA sequencing and qRT-PCR analysis (Table S1 and Figure 5), was identified as a potential key HSF gene in benzenoid/phenylpropanoid biosynthesis for further research.
PAL is involved in the production of precursors which contribute to benzenoid/phenylpropanoid biosynthesis [17]. Three STRE elements, the important cis-acting regulatory elements that HSF can bind to, are present in the promoter of the PhPAL2 gene (Figure 6c). What is more, PhPAL2 was chosen as a candidate target gene of PhHSF19 because of their similar expression pattern. The expression level of PhPAL2 is significantly higher at the early flower opening stage than that at the bud stage (Figure 6a), and the expression pattern of PhHSF19 is positively related with PhPAL2, rather than PhPAL1 (Figure 6b). Both PhPAL1 and PhPAL2 from the PAL gene family are supposed to have a similar and redundant function, but the expression pattern of these two genes is not the same, suggesting that PhPAL1 and PhPAL2 exhibit a different response to the regulation of different transcription factors. In this study, PhPAL2 was the candidate target gene of PhHSF19, while PhPAL1 was not. Similarly, PhPAL2 can be regulated by EOBII, but PhPAL1 cannot [20].
The phylogeny analysis showed that PhHSF19 belonging to class A3 without AHA motifs does not perform as a transcription activator [1]. However, a previous study found that tryptophane residues in the C-terminal domains of HSFs from class A3 provide an additive contribution to the activator function [33]. The important tryptophane residues are also present in the C-terminal domains of PhHSF19, and we speculated that PhHSF19 should possess a transcriptional activation function. Then, the Y1H and dual luciferase assays not only prove this speculation, but also confirmed the interaction between PhHSF19 and PhPAL2 in our study. Based on the dual luciferase assay, about a 2.8-fold induction in the LUC/REN ratio verified that PhHSF19 could positively regulate the expression of PhPAL2. Therefore, an up-regulated expression of PhPAL2 under the regulation of PhHSF19 would enhance the transformation of phenylalanine to trans-cinnamic acid [17], which probably provided more metabolic substrates for further benzenoid biosynthesis and made the flower fragrance more intense. What is more, LsHSFB2A-1 can bind to the LsMADS55 promoter to regulate floral organ specification and flower development in lettuce [6]. The expression levels of DcaHsf-A5 and DcaHsf-B2b are closely related with flower development and are the highest at the full-opening stage (FS6) [12]. This suggests that PhHSF19 might have a similar function in the flower opening process, which would also affect flower scent to a certain extent. More research needs be carried out to prove this assumption.
In summary, 34 HSF gene family members were identified in petunia. The expression analysis suggested that PhHSF19 possesses potential roles in regulating benzenoid/phenylpropanoid biosynthesis. In addition, PhHSF19 could positively regulate the expression of PhPAL2 based on the yeast one hybrid and dual luciferase assays.
4. Materials and Methods
4.1. Plant Materials
Petunia hybrida “Mitchell” plants were cultured in an environmental chamber in Zhejiang Agriculture and Forestry University, Hangzhou, China. The chamber was programmed for 16 h light/8 h dark with a temperature of 24 °C, with the light period starting at 07:00 a.m. Petal limbs were sampled at 16:00 p.m. from three individual flowers before anthesis (day 0, D0) and from the first to seventh day of anthesis (D1–D7). In addition, different tissues, including leaf, peduncle, anther, ovary, petal tube, petal limb, stigma, and sepal of three individual plants at 16:00 p.m. were collected on the second day of anthesis. All of the samples were immediately frozen in liquid nitrogen and stored at −80 °C.
4.2. Identification HSF Genes in Petunia
The Hidden Markov Model (HMM) of the HSF-type DBD domain (PF00447) was obtained from the Pram database (
4.3. Gene Structure and Phylogenetic Analysis of Petunia HSFs
The exon/intron structures of PhHSFs were displayed by the Gene Structure Displayer Server v2.0 (GSDS,
4.4. Conserved Motifs and Cis-Acting Elements Analysis of Petunia HSFs
Multiple Expectation maximization for Motif Elicitation (MEME,
4.5. Selection of Potential Key PhHSF Genes Based on RNA Sequencing
RNA sequencing data for the day 0 (bud) or day 2 (post-anthesis) corolla tissue of P. hybrida “Mitchell” were obtained from the NBCI Gene Expression Omnibus (GSE70948) [23]. RNA sequencing analysis was conducted using the same method as in a previous study [34]. The assembled sequences were assigned to Petunia axillaris (v1.6.2) genomes [19,35] with Local BLAST to identify PhHSF genes, and the corresponding gene expressions were acquired with the reads per kilobases per million reads (RPKM) statistics. Twenty PhHSF genes were screened under a threshold of identify of ≥95% (Table S1). Then, six differentially expressed PhHSF genes were further determined with a fold-change of day 2 vs. day 0 ≥ 2 and the expression level (RPKM) at day 2 ≥ 1000.
4.6. Total RNA Extraction, cDNA Reverse Transcription and qRT-PCR Analysis
The total RNA was extracted from different samples according to the manufacturer’s instruction for the FastPure® Plant Total RNA Isolation Kit (Polysaccharides and Polyphenolics-rich) (Vazyme, Nanjing, China). The first-strand cDNA was synthesized from 500 ng total RNA, 2 μL 5× PrimerScript RT Master Mix, and RNase Free ddH2O to 10 μL reaction mixture using the PrimeScript™ RT Master Mix (Takara, Dalian, China). The relative expression of the related genes (six differentially expressed PhHSF genes and two PhPALs) was determined by qRT-PCR using ChamQ SYBR qPCR Master Mix (Vazyme, Nanjing, China) on LightCycler® 480 II (Roche Diagnostics, Mannheim, Germany). The 10 μL reaction condition was configured according to the following instruction: 95 °C for 30 s, 40 cycles of amplification (95 °C for 10 s, 60 °C for 30 s, and 60 °C for 30 s), and melting curve (95 °C for 5 s, 60 °C for 1 min, and then 95 °C). Primer sequences are shown in Table S2. The relative expression of the target genes was calculated with 2−∆CT using the reference gene PhEF1-α.
4.7. Y1H Assay
The yeast one hybrid assay was used to identify the transcriptional regulators of PhHSF19. The promoter sequence of PhPAL2 (PhPAL2-Pro) and was amplified and inserted into the pAbAi vector, and the CDS of PhHSF19 was cloned into the pGADT7-Rec2 vector with the primers listed in Table S2. pAbAi-PhPAL2-Pro was transformed into Y1H Gold according to the protocol of Yeastmaker™ Yeast Transformation System 2 User Manual, and the lowest inhibitory concentration of aureobasidin A (AbA) for the bait strain was tested on SD/-Ura media. The prey was transformed into yeast cells containing recombinant pAbAi vector and then incubated on SD/-Leu medium at 30 °C for 3 days. The interactions between pGADT7-Rec2-PhHSF19 and pAbAi-PhPAL2-Pro were detected on an SD/-Leu medium with 100 ng mL−1 AbA at 30 °C for 3–5 days.
4.8. Dual Luciferase Assay
A dual-luciferase assay was used to further test the interaction between PhHSF19 and the promoter of the PhPAL2 in vivo. Full-length sequence of regulatory genes and target promoters were cloned into pCNHP and pGreenII 0800-LUC vectors, respectively, using the primers listed in Table S2. Then, the recombinant vectors were transformed into Agrobacterium strain GV3101 (pSoup). The transformed strains were cultured in 28 °C and then resuspended in an infiltration medium (10 mM MgCl2 and 200 μM acetosyringone) to an OD 600 of 1.6. The cell suspension was mixed at a ratio of 1:1 (v/v) and injected into Nicotiana benthamiana leaves with needleless syringes. After 3 days, the transformed tobacco leaves, which had been sprayed with D-luciferin sodium salt (Coolaber, Beijing, China), were observed under Fuison-FX7 multi-imaging apparatus (Vilber, Collégien, France).
Conceptualization, J.F.; data curation, S.H., J.Q., H.Q., Z.W. and H.C.; funding acquisition, J.F.; methodology, S.H., J.Q., H.Q. and C.Z.; supervision, J.F. and C.Z.; writing—original draft, S.H., J.Q. and C.Z.; writing—review and editing, J.F. and C.Z. All authors have read and agreed to the published version of the manuscript.
This research was funded by the National Natural Science Foundation of China (grant no. 32002077).
Not applicable.
Not applicable.
Not applicable.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Gene structure analysis of 34 petunia HSFs. The exons and introns are indicated by the yellow boxes and black lines, respectively.
Enlarge this image.Figure 2. Phylogeny analysis of petunia HSFs with AtHSFs from Arabidopsis thaliana. HSFs were divided into three major subfamilies HSFA (red), HSFB (blue), and HSFC (green).
Enlarge this image.Figure 3. Conserved motif analysis of petunia HSFs. Proteins are organized according to the groups in Figure 2. Fifteen motifs are found in the protein sequences.
Enlarge this image.Figure 4. Cis-acting elements analysis of promoter sequences of petunia HSFs. The position of cis-acting elements in the 2kb upstream region is represented by different colored bars.
Enlarge this image.Figure 5. qRT-PCR analysis of six candidate PhHSF genes. All transcript levels were determined by qRT-PCR and were obtained relative to the reference gene (PhEF1-α). (a) Tissue-specific expression of PhHSFs in leaf, peduncle, anther, ovary, petal tube, petal limb, stigma, and sepal collected at 1500 h of day 1 post-anthesis (D1). (b) Stage-specific expression of PhHSFs from mature buds (D0) to day 7 post-anthesis (D7) at 1500 h. Representative histograms from multiple biological replications using the 2−∆CT method are shown (mean ± SE, n = 3).
Enlarge this image.Figure 6. PhHSF19 could bind to the PhPAL2 promoter. (a) Relative expression level of PhPAL1 and PhPAL2 during flower opening. (b) Correlation analysis between the relative expression of PhHSF19 and PhPAL1/PhPAL2. (c) Location of the STRE cis-element in the PhPAL2 promoter. (d) Yeast one-hybrid assays of the PhHSF19 and PhPAL2 promoters. (e) Dual luciferase assay between the PhHSF19 and PhPAL2 promoter, including the schematic diagram of the double-reporter and effector plasmids, qualitative results, and quantitative results. T-test was used for statistical analyses of quantitative results compared with corresponding control (**, p < 0.01).
Characteristics of petunia heat shock transcription factors, i.e., protein length (aa), molecular weight (MW), theoretical isoelectric point (pI), instability index, hydrophilic index, aliphatic index, and predicted subcellular localization.
| Loci | Protein Length |
Molecular Weight |
Theoretical pI | Instability Index | Hydrophilic Index | Aliphatic Index | Subcellular Localization |
|---|---|---|---|---|---|---|---|
| Peaxi162Scf00002g02718 | 250 | 28.83 | 5.77 | 55.15 | −0.868 | 64.40 | nucleus |
| Peaxi162Scf00004g00431 | 253 | 29.36 | 6.56 | 41.27 | −0.745 | 84.31 | nucleus |
| Peaxi162Scf00010g00269 | 403 | 45.91 | 6.34 | 60.36 | −0.734 | 70.89 | nucleus |
| Peaxi162Scf00036g00820 | 314 | 35.82 | 6.12 | 47.69 | −0.805 | 72.07 | nucleus |
| Peaxi162Scf00045g00131 | 512 | 56.56 | 4.71 | 68.62 | −0.636 | 68.57 | nucleus |
| Peaxi162Scf00112g00104 | 400 | 45.40 | 5.48 | 47.29 | −0.649 | 79.20 | nucleus |
| Peaxi162Scf00112g00814 | 302 | 34.25 | 5.99 | 42.22 | −0.728 | 75.56 | nucleus |
| Peaxi162Scf00130g00326 | 249 | 28.91 | 6.68 | 54.38 | −0.865 | 65.30 | nucleus |
| Peaxi162Scf00131g00415 | 359 | 41.74 | 5.00 | 51.96 | −0.835 | 64.32 | nucleus |
| Peaxi162Scf00145g00167 | 224 | 25.98 | 9.18 | 41.12 | −0.728 | 74.06 | nucleus |
| Peaxi162Scf00154g01021 | 405 | 46.04 | 5.39 | 54.92 | −0.820 | 71.23 | nucleus |
| Peaxi162Scf00170g01112 | 500 | 55.63 | 4.75 | 67.61 | −0.647 | 66.08 | nucleus |
| Peaxi162Scf00205g00818 | 390 | 45.15 | 4.63 | 54.33 | −0.673 | 73.46 | nucleus |
| Peaxi162Scf00205g00820 | 69 | 8.02 | 4.32 | 42.50 | 0.333 | 94.49 | nucleus |
| Peaxi162Scf00221g00015 | 360 | 41.84 | 6.64 | 57.17 | −0.759 | 70.92 | nucleus |
| Peaxi162Scf00252g00226 | 314 | 35.81 | 5.35 | 48.60 | −0.803 | 71.43 | nucleus |
| Peaxi162Scf00259g00710 | 369 | 42.10 | 5.42 | 51.46 | −0.751 | 65.77 | nucleus |
| Peaxi162Scf00316g00518 | 378 | 44.14 | 5.08 | 49.36 | −0.883 | 67.57 | nucleus |
| Peaxi162Scf00327g00413 | 399 | 45.58 | 5.19 | 53.32 | −0.406 | 87.39 | nucleus |
| Peaxi162Scf00433g00411 | 461 | 51.58 | 4.90 | 54.28 | −0.526 | 66.18 | nucleus |
| Peaxi162Scf00450g00422 | 221 | 24.71 | 5.85 | 41.11 | −0.765 | 54.71 | nucleus |
| Peaxi162Scf00461g00026 | 199 | 23.31 | 9.42 | 49.99 | −0.735 | 69.10 | nucleus |
| Peaxi162Scf00498g00438 | 357 | 41.59 | 5.08 | 55.61 | −0.777 | 70.73 | nucleus |
| Peaxi162Scf00503g00003 | 348 | 39.75 | 7.75 | 52.24 | −0.664 | 64.91 | nucleus |
| Peaxi162Scf00516g00645 | 503 | 55.67 | 5.67 | 47.46 | −0.604 | 74.77 | nucleus |
| Peaxi162Scf00546g00118 | 319 | 36.05 | 5.35 | 69.94 | −0.600 | 63.86 | nucleus |
| Peaxi162Scf00560g00011 | 263 | 30.25 | 6.62 | 47.95 | −0.651 | 64.52 | nucleus |
| Peaxi162Scf00619g00912 | 385 | 45.37 | 5.43 | 42.70 | −0.907 | 60.08 | nucleus |
| Peaxi162Scf00652g00018 | 443 | 49.37 | 5.39 | 57.72 | −0.794 | 65.80 | nucleus |
| Peaxi162Scf00714g00331 | 305 | 33.23 | 5.24 | 56.40 | −0.341 | 80.89 | nucleus |
| Peaxi162Scf00815g00357 | 469 | 52.79 | 5.40 | 49.77 | −0.550 | 72.69 | nucleus |
| Peaxi162Scf00915g00029 | 323 | 35.67 | 4.94 | 53.00 | −0.667 | 68.85 | nucleus |
| Peaxi162Scf01009g00222 | 407 | 46.63 | 5.20 | 48.51 | −0.754 | 70.42 | nucleus |
| Peaxi162Scf16358g00003 | 224 | 25.90 | 9.19 | 41.46 | −0.729 | 74.06 | nucleus |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Scharf, K.D.; Berberich, T.; Ebersberger, I.; Nover, L. The plant heat stress transcription factor (Hsf) family: Structure, function and evolution. Biochim. Biophys. Acta; 2012; 1819, pp. 104-119. [DOI: https://dx.doi.org/10.1016/j.bbagrm.2011.10.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22033015]
2. 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, pp. 785-795. [DOI: https://dx.doi.org/10.1016/j.plaphy.2009.05.003]
3. Nover, N.; Bharti, K.; DRing, P.; Mishra, S.K.; Scharf, G. Arabidopsis and the heat stress transcription factor world: How many heat stress transcription factors do we need?. Cell Stress Chaperones; 2001; 6, pp. 177-189. [DOI: https://dx.doi.org/10.1379/1466-1268(2001)006<0177:AATHST>2.0.CO;2]
4. Pelham, H.R.B. A regulatory upstream promoter element in the Drosophila Hsp 70 heat-shock gene. Cell; 1982; 30, pp. 517-528. [DOI: https://dx.doi.org/10.1016/0092-8674(82)90249-5]
5. Li, S.; Liu, S.; Lin, X.; Grierson, D.; Yin, X.; Chen, K. Citrus heat shock transcription factor CitHsfA7-mediated citric acid degradation in response to heat stress. Plant Cell Environ.; 2021; 45, pp. 95-104. [DOI: https://dx.doi.org/10.1111/pce.14207] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34705284]
6. Ning, K.; Han, Y.; Chen, Z.; Luo, C.; Wang, S.; Zhang, W.; Li, L.; Zhang, X.; Fan, S.; Wang, Q. Genome-wide analysis of MADS-box family genes during flower development in lettuce. Plant Cell Environ.; 2019; 42, pp. 1868-1881. [DOI: https://dx.doi.org/10.1111/pce.13523] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30680748]
7. Guo, L.; Chen, S.; Liu, K.; Liu, Y.; Ni, L.; Zhang, K.; Zhang, L. Isolation of heat shock factor HsfA1a-binding sites in vivo revealed variations of heat shock elements in Arabidopsis thaliana. Plant Cell Physiol.; 2008; 9, pp. 1306-1315. [DOI: https://dx.doi.org/10.1093/pcp/pcn105]
8. Jiang, Y.; Zheng, Q.; Chen, L.; Liang, Y.; Wu, J. Ectopic overexpression of maize heat shock transcription factor gene ZmHsf04 confers increased thermo and salt-stress tolerance in transgenic Arabidopsis. Acta Physiol. Plant.; 2017; 40, pp. 9-20. [DOI: https://dx.doi.org/10.1007/s11738-017-2587-2]
9. Zhu, B.; Ye, C.; Lü, H.; Chen, X.; Chai, G.; Chen, J.; Wang, C. Identification and characterization of a novel heat shock transcription factor gene, GmHsfA1, in soybeans (Glycine max). J. Plant Res.; 2006; 119, pp. 247-256. [DOI: https://dx.doi.org/10.1007/s10265-006-0267-1]
10. 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. [DOI: https://dx.doi.org/10.3389/fpls.2015.00736] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26442049]
11. Xie, L.; Li, X.; Hou, D.; Cheng, Z.; 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. [DOI: https://dx.doi.org/10.3390/f10020100]
12. 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. [DOI: https://dx.doi.org/10.3390/ijms20205233]
13. Zinsmeister, J.; Berriri, S.; Basso, D.; Vu, B.L.; Buitink, J. The seed-specific heat shock factor A9 regulates the depth of dormancy in Medicago truncatula seeds via ABA signaling. Plant Cell Environ.; 2020; 43, pp. 2508-2522. [DOI: https://dx.doi.org/10.1111/pce.13853]
14. Shang, J.; Tian, J.; Cheng, H.; Yan, Q.; Chen, L. The chromosome-level wintersweet (Chimonanthus praecox) genome provides insights into floral scent biosynthesis and flowering in winter. Genome Biol.; 2020; 21, pp. 200-227. [DOI: https://dx.doi.org/10.1186/s13059-020-02088-y]
15. Yang, J.W.; Ren, Y.J.; Zhang, D.Y.; Chen, X.W.; Huang, J.Z.; Xu, Y. Transcriptome-Based WGCNA Analysis Reveals Regulated Metabolite Fluxes between Floral Color and Scent in Narcissus tazetta Flower. Int. J. Mol. Sci.; 2021; 15, 8249. [DOI: https://dx.doi.org/10.3390/ijms22158249] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34361014]
16. Verdonk, J.C.; Haring, M.; van Tunen, A.J.; Schuurink, R.C. ODORANT1 Regulates Fragrance Biosynthesis in Petunia Flowers. Plant Cell; 2005; 17, pp. 1612-1624. [DOI: https://dx.doi.org/10.1105/tpc.104.028837] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15805488]
17. Boatright, J.; Negre, F.; Chen, X.; Kish, C.M. Understanding in Vivo Benzenoid Metabolism in Petunia Petal Tissue. Plant Physiol.; 2004; 135, pp. 1993-2011. [DOI: https://dx.doi.org/10.1104/pp.104.045468] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15286288]
18. Verdonk, J.C.; Vos, C.; Verhoeven, H.A.; Haring, M.A.; Schuurink, R.C. Regulation of floral scent production in Petunia revealed by targeted metabolomics. Phytochemistry; 2003; 62, pp. 997-1008. [DOI: https://dx.doi.org/10.1016/S0031-9422(02)00707-0]
19. Aureliano, B.; Michel, M.; Avichai, A.; Manzhu, B.; Laure, B.; Cornelius, S.B.; Mattijs, B. Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat. Plants; 2016; 2, 16074.
20. Spitzer-Rimon, B.; Marhevka, E.; Barkai, O.; Marton, I.; Edelbaum, O.; Masci, T.; Prathapani, N.K.; Shklarman, E.; Vainstein, O.A. EOBII, a Gene Encoding a Flower-Specific Regulator of Phenylpropanoid Volatiles’ Biosynthesis in Petunia. Plant Cell; 2010; 22, pp. 1961-1976. [DOI: https://dx.doi.org/10.1105/tpc.109.067280] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20543029]
21. Spitzer-Rimon, B.; Farhi, M.; Albo, B.; Cna’Ani, A.; Zvi, M.M.B. The R2R3-MYB–Like Regulatory Factor EOBI, Acting Downstream of EOBII, Regulates Scent Production by Activating ODO1 and Structural Scent-Related Genes in Petunia. Plant Cell; 2012; 24, pp. 5089-5105. [DOI: https://dx.doi.org/10.1105/tpc.112.105247] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23275577]
22. Guo, J.; Jian, W.; Qian, J.; Chao, W.; Lei, L.; Yi, Y.; Wang, Y. Genome-wide analysis of heat shock transcription factor families in rice and Arabidopsis. J. Genet. Genom.; 2008; 035, pp. 105-118. [DOI: https://dx.doi.org/10.1016/S1673-8527(08)60016-8]
23. Widhalm, J.R.; Gutensohn, M.; Yoo, H.; Adebesin, F.; Qian, Y.; Guo, L.; Jaini, R.; Lynch, J.H.; Mccoy, R.M.; Shreve, J.T. Identification of a plastidial phenylalanine exporter that influences flux distribution through the phenylalanine biosynthetic network. Nat. Commun.; 2015; 6, 8142. [DOI: https://dx.doi.org/10.1038/ncomms9142] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26356302]
24. Underwood, B.A.; Tieman, D.M.; Shibuya, K.; Dexter, R.J.; Loucas, H.M.; Simkin, A.J.; Clark, D.G. Ethylene-Regulated Floral Volatile Synthesis in Petunia Corollas. Plant Physiol.; 2005; 138, pp. 255-266. [DOI: https://dx.doi.org/10.1104/pp.104.051144] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15849311]
25. Colquhoun, T.A.; Verdonk, J.C.; Schimmel, B.C.J.; Tieman, D.M.; Underwood, B.A.; Clark, D.G. Petunia floral volatile benzenoid/phenylpropanoid genes are regulated in a similar manner. Phytochemistry; 2010; 71, pp. 158-167. [DOI: https://dx.doi.org/10.1016/j.phytochem.2009.09.036] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19889429]
26. Wang, F.; Dong, Q.; Jiang, H.; Zhu, S.; Chen, B.; Xiang, Y. Genome-wide analysis of the heat shock transcription factors in Populus trichocarpa and Medicago truncatula. Mol. Biol. Rep.; 2012; 39, pp. 1877-1886. [DOI: https://dx.doi.org/10.1007/s11033-011-0933-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21625849]
27. Liu, M.; Huang, Q.; Sun, W.; Ma, Z.; Huang, L.; Wu, Q.; Tang, Z.; Bu, T.; Li, C.; Chen, H. Genome-wide investigation of the heat shock transcription factor (Hsf) gene family in Tartary buckwheat (Fagopyrum tataricum). BMC Genom.; 2019; 20, pp. 871-887. [DOI: https://dx.doi.org/10.1186/s12864-019-6205-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31730445]
28. Lohani, N.; Golicz, A.A.; Singh, M.B.; Bhalla, P.L. Genome-wide analysis of the Hsf gene family in Brassica oleracea and a comparative analysis of the Hsf gene family in B. oleracea, B. rapa and B. napus. Funct. Integr. Genom.; 2019; 18, pp. 649-665. [DOI: https://dx.doi.org/10.1007/s10142-018-0649-1]
29. Guo, M.; Lu, J.P.; Zhai, Y.F.; Chai, W.G.; Gong, Z.H.; Lu, M.H. Genome-wide analysis, expression profile of heat shock factor gene family (CaHsfs) and characterisation of CaHsfA2 in pepper (Capsicum annuum L.). BMC Plant Biol.; 2015; 15, 151. [DOI: https://dx.doi.org/10.1186/s12870-015-0512-7]
30. Sun, T.; Wang, C.; Liu, R.; Zhang, Y.; Wang, Y.; Wang, L. Thhsfa1 Confers Sa1t Stress Tolerance through Modulation of Reactive Oxygen Species Scavenging by Directly Regulating ThWRKY4. Int. J. Mol. Sci.; 2021; 22, 5048. [DOI: https://dx.doi.org/10.3390/ijms22095048]
31. Rao, S.; Das, J.R.; Mathur, S. Exploring the master regulator heat stress transcription factor HSFA1a-mediated transcriptional cascade of HSFs in the heat stress response of tomato. J. Plant Biochem. Biot.; 2021; 30, pp. 878-888. [DOI: https://dx.doi.org/10.1007/s13562-021-00696-8]
32. Mishra, S.K.; Poonia, A.; Chaudhary, B.; Baranwal, V.K.; Arora, D.; Kumar, R.; Chauhan, H. Genome-wide identification, phylogeny and expression analysis of HSF gene family in barley during abiotic stress response and reproductive development. Plant Gene; 2020; 23, 100231. [DOI: https://dx.doi.org/10.1016/j.plgene.2020.100231]
33. Bharti, K.; Schmidt, E.; Lyck, R.; Heerklotz, D.; Scharf, K.D. Isolation and characterization of HsfA3, a new heat stress transcription factor of Lycopersicon peruvianum. Plant J.; 2010; 22, pp. 355-365. [DOI: https://dx.doi.org/10.1046/j.1365-313x.2000.00746.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10849352]
34. Patrick, R.M.; Huang, X.Q.; Natalia, D.; Li, Y. Dynamic histone acetylation in floral volatile synthesis and emission in petunia flowers. J. Exp. Bot.; 2021; 10, pp. 3704-3722. [DOI: https://dx.doi.org/10.1093/jxb/erab072] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33606881]
35. Noe, F.P.; Naama, M.; Edwards, J.D.; Surya, S.; Tecle, I.Y.; Strickler, S.R.; Aureliano, B.; Thomas, F.Y.; Anuradha, P.; Hartmut, F. The Sol Genomics Network (SGN)—From genotype to phenotype to breeding. Nucleic Acids Res.; 2015; 10, pp. 1036-1041.
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
