Genome-wide characterization of the HpHsf gene family and expression analysis under heat stress in Herpetospermum pedunculosum

Introduction: Heat shock transcription factors (Hsfs) are key regulators of plant responses to heat stress and other environmental challenges. However, the Hsf gene family of Herpetospermum pedunculosum, an alpine medicinal plant valued for its hepatoprotective properties, remains poorly characterized. Investigating the characteristics of the Hsf genes in H. pedunculosum will enhance our understanding of its thermotolerance mechanisms and provide a theoretical basis for improving its environmental adaptability.

Methods: In this study, we conducted a genome-wide identification and characterization of the Hsf gene family in H. pedunculosum. The study included analyses of protein physicochemical properties, chromosome locations, phylogenetic classification, conserved motifs, gene structures, collinearity, and cis-acting promoter elements. qRT-PCR was performed to assess the transcription levels of HpHsf genes in H. pedunculosum leaves under heat stress over a time course. Three representative HpHsfs were selected for subcellular localization analysis.

Results: Here, a total of 21 HpHsf genes were identified. Phylogenetic analysis classified the HpHsfs into three main groups (A, B, and C), consistent with classifications in other plant species. Promoter analysis revealed enrichment of stress, hormone signaling, and development-related elements. Furthermore, expression profiling under heat stress revealed significant upregulation of several HpHsfs, suggesting their involvement in heat stress response. Subcellular localization assays of three representative proteins (HpHsf3, HpHsf5, and HpHsf17) in Nicotiana benthamiana epidermal cells confirmed their nuclear localization, supporting their function as transcription factors.

Discussion: These findings provide new insights into the HpHsf gene family and lay a foundation for future functional studies on heat tolerance in alpine plants.

1 Introduction

Global climate change poses a significant threat to agricultural productivity by altering environmental conditions, highlighting the urgent need to elucidate plant responses to abiotic stressors, particularly heat stress (Sato et al., 2024). Elevated temperatures compromise plant health and performance, leading to stunted growth, reduced yields, and disrupted physiological functions. Alpine medicinal plants, adapted to cooler environments, are especially vulnerable to rising temperatures that impair their development and reproductive capacity (Zhu et al., 2024a). Understanding heat stress responses mechanisms in these high-altitude species is essential for formulating adaptive strategies to enhance their resilience and ensure sustainable cultivation in a warming climate.

Heat shock transcription factors (Hsfs) constitute a key family of transcription factors that serve as the central regulators of the heat shock response (Scharf et al., 2012). This response represents a highly conserved protective mechanism ubiquitously employed by organisms to mitigate various stresses, particularly thermal stress (Ohama et al., 2017; Wang and Xu, 2025). Structurally, Hsfs are characterized by several conserved functional domains, including a DNA-binding domain (DBD), an oligomerization domain (OD), a nuclear localization signal (NLS), a nuclear export signal (NES), a repressor domain (RD), and an aromatic, hydrophobic, and acidic (AHA) motif (Wang and Xu, 2025). Based on these structural features, Hsfs are generally classified into three major groups: A, B, and C. Class A Hsfs possess a C-terminal AHA motif that enables direct transcriptional activation of downstream target genes (Von Koskull-Döring et al., 2007). In contrast, Class B Hsfs contain a repressor domain with a highly conserved LFGV tetrapeptide motif (Ikeda and Ohme-Takagi, 2009). Class C Hsfs are generally shorter and structurally simpler than those in Classes A and B. Interestingly, the expansion of Class C Hsfs in monocots suggests potential functional specialization among members within this group (Andrási et al., 2020).

Plant Hsfs maintain a dynamic equilibrium, switching between active and inactive states in response to environmental conditions and adaptive mechanisms (Bakery et al., 2024). According to the activation cycle model, Hsfs are initially maintained in an inactive state within cytoplasmic complexes bound to HSP70 and HSP90. Under heat stress, misfolded proteins compete for and sequester HSPs, leading to Hsf release. Once liberated, Hsfs oligomerize, translocate into the nucleus, and activate the transcription of heat-responsive genes, including those encoding HSPs and additional Hsfs (Scharf et al., 2012; Bakery et al., 2024). In Arabidopsis thaliana, HsfA1 is a master regulator of high-temperature signaling, coordinating heat shock and warm-temperature responses. Specifically, HsfA1d links these pathways by stabilizing phytochrome-interacting factor 4 (PIF4), which promotes thermomorphogenesis, while simultaneously activating HsfA2 to reinforce thermotolerance (Raturi and Zinta, 2024; Li et al., 2024; Tan et al., 2023). Nevertheless, HsfBs are occasionally regarded as repressors of the heat shock response, capable of modulating the activity of class A members (Andrási et al., 2020; Wang and Xu, 2025). Loss-of-function mutations in AtHsfB1 and AtHsfB2b confer enhanced heat tolerance in Arabidopsis, accompanied by elevated expression of AtHsfA2 and AtHsfA7a under heat stress (Pick et al., 2012; Ikeda et al., 2011). In tomato, SlHsfb1 acts as a co-activator of SlHsfA1a, promoting the accumulation of heat shock proteins and enhancing thermotolerance. In addition, SlHsfb1 functions as a transcriptional repressor of other SlHsfs, including SlHsfA1b and SlHsfA2, thereby maintaining a balance between growth and stress responses under heat stress conditions (Fragkostefanakis et al., 2019). Furthermore, class C Hsfs contribute to the regulation of plant heat tolerance. For example, in lilies, high-temperature stress suppresses LlHsfC2 homodimerization while enhancing its heterodimerization with HsfAs, thereby enabling sustained co-activation of heat tolerance mechanisms (Wu et al., 2024).

Herpetospermum pedunculosum, commonly known as ‘Bo-Leng-Gua’, is an annual species of Cucurbitaceae family, endemic to the Tibetan Plateau. It has garnered significant attention due to its high pharmacological value, particularly in traditional Tibetan medicine (Zhao et al., 2019). The seeds of H. pedunculosum are notably rich in lignan compounds (Zhu et al., 2024b), which exhibit diverse medicinal properties, including hepatoprotective, anti-cholestatic, and anti-inflammatory effects (Fang et al., 2007; Wei et al., 2021; Li et al., 2023a). These attributes underscore its importance within traditional Tibetan medicine. Ecologically, H. pedunculosum is well adapted to the extreme conditions of high-altitude habitats, typically occurring between 2,300 and 3,500 m, where it withstands such as strong winds, low temperatures, drought, and hypoxic stress (Zhao et al., 2019; Chen et al., 2025). However, with ongoing climate change and increasing global temperatures, its long-term survival and medicinal utility may be at risk. Thus, elucidating the mechanisms underlying its thermotolerance is crucial not only for understanding its ecological resilience but also for informing future conservation, cultivation, and sustainable utilization strategies.

Hsfs play a central role in plant responses to high-temperature stress by regulating the expression of heat shock proteins and other protective genes (Scharf et al., 2012; Kan et al., 2023). Investigating the Hsf gene family in H. pedunculosum offers valuable insights into its thermotolerance mechanisms and adaptive evolution, while providing a molecular basis for developing stress-resilient cultivars. In this study, we comprehensively characterize the HpHsf gene family in H. pedunculosum and analyze their expression patterns under heat stress conditions, thereby advancing the broader understanding of thermotolerance mechanisms in plants.

2 Materials and methods

2.1 Identification of Hsf genes from H. pedunculosum

The complete genome sequence of H. pedunculosum was obtained from Figshare (https://doi.org/10.6084/m9.figshare.21626153.v2) (Yang et al., 2023). Hidden Markov Model (HMM)-based screening was employed to identify Hsf genes. The Hsf domain profile (PF00447) was retrieved from the Pfam database, and putative Hsfs were identified using HMMER v3.0 with a significance threshold of E-value < 1 × 10^-5 (Potter et al., 2018; Mistry et al., 2020). Candidate Hsf protein sequences were subsequently validated using the NCBI Batch CD-Search tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi).

2.2 Sequence analyses of HpHsfs

Protein physicochemical properties and gene structures were analyzed using TBtools (Chen et al., 2023). Subcellular localization was predicted using WOLF PSORT (https://www.genscript.com/wolf-psort.html). Conserved motifs were identified with MEME Suite v5.5.8 (https://meme-suite.org/meme/) (Bailey et al., 2009), while conserved domains within protein sequences were predicted using the NCBI Batch Web CD-Search tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi). Key amino acid motifs were identified from a ClustalW-based multiple sequence alignment of all HpHsf proteins, with reference to the established characteristics of Hsf signature motifs (Scharf et al., 2012).

2.3 Phylogenetic analysis of the HpHsf gene family

For phylogenetic analysis and subfamily classification, full-length Hsf amino acid sequences from Cucumis sativus, Arabidopsis thaliana, and Oryza sativa were retrieved using a consistent approach. Multiple sequence alignment was performed with ClustalW, and phylogenetic trees were constructed in MEGA 12 using the neighbor-joining method with 1000 bootstrap replicates and default settings (Substitution type: Amino acid; Model/Method: Jones-Taylor-Thornton (JTT) model; Rates among Sites: Uniform Rates; Gaps/Missing Data Treatment: Patial deletion, 50%) (Tamura et al., 2021). The resulting subfamily clustering was consistent with the previously established classification of the C. sativus Hsf family (Chen et al., 2021). Phylogenetic trees were visualized using the ChiPlot online tool (https://www.chiplot.online/tvbot.html) (Xie et al., 2023).

2.4 Analysis of cis-regulatory elements of HpHsfs

Using the Gtf/Gff3 Sequences Extractor in TBtools, 2,000 bp promoter regions upstream of the ATG start codon of all H. pedunculosum Hsf genes were extracted from the genome. Cis-acting regulatory elements within these sequences were identified using the PlantCARE web server (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (Lescot et al., 2002), and the results were visualized with TBtools’ Simple BioSequence Viewer.

2.5 Chromosome distribution and gene collinearity analysis of HpHsf gene family

Based on the H. pedunculosum genome annotation, we retrieved chromosomal localization information for HpHsf genes, generated a chromosomal distribution map with TBtools, and designated gene names according to their positional order along the chromosomes. Intraspecific collinearity analysis of Hsf genes in H. pedunculosum was performed using MCScanX integrated into TBtools, and the results were visualized with the Advanced Circos module. Genome data for A. thaliana were obtained from The Arabidopsis Information Resource (TAIR; http://www.arabidopsis.org) (Lamesch et al., 2011), while those for C. sativus were retrieved from the Cucurbit Genomics Database (http://cucurbitgenomics.org/) (Li et al., 2019). Interspecific collinearity analysis between H. pedunculosum and the two reference species was also conducted using MCScanX in TBtools.

2.6 Plant materials and heat stress

Seeds of H. pedunculosum were sourced from the cultivation base of Tibet Rhodiola Pharmaceutical Holding Company in Lhasa (29°79′ N, 94°09′ E, Nyingchi, Xizang, China). These seeds were morphologically authenticated by Prof. Qi Zhao from Chengdu University and subsequently stored at the Engineering Research Center of Sichuan-Xizang Traditional Medicinal Plants, Chengdu, China. Seedlings were hydroponically grown in Hoagland solution under controlled growth chamber conditions (24°C day/22°C night, 16-h photoperiod) for 40 days. Uniform seedlings were subjected to heat stress at 42°C (Xu et al., 2006; Chen et al., 2021), and leaf samples were collected at 0 (control), 1, 3, 6, and 12 hours post-treatment. All samples were immediately frozen in liquid nitrogen and stored at -80°C until RNA extraction.

2.7 RNA extraction, qRT-PCR, and RT-PCR analysis

Leaf samples of H. pedunculosum were homogenized using an ALLSHENG Bioprep-24R homogenizer (Hangzhou, China). Subsequently, total RNA was extracted from the homogenized tissue using a plant total RNA extraction kit (Foregene Biotechnology, Chengdu, China) in accordance with the manufacturer’s instructions. cDNA was synthesized using the HiFiScript gDNA Removal RT MasterMix (Cwbio, Jiangsu, China). The expression levels of HpHsfs were quantified using quantitative real-time PCR (qRT-PCR) using the Cwbio SYBR Master Mix kit (Cwbio, Jiangsu, China) on a Bio-Rad CFX96 Real-time PCR system, with the H. pedunculosum Actin gene serving as the internal control. The primer sequences are provided in Supplementary Table 1. Statistical analysis was performed using one-way analysis of variance (ANOVA) in SPSS 26.0. Data obtained from three or four independent technical replicates are presented as the mean ± standard deviation (SD). The results were verified in two independent biological experiments. Significant differences (p < 0.05) were determined using Tukey’s multiple comparisons test and are indicated by distinct superscript letters. The raw data are provided in Supplementary Table 3.

For RT-PCR, first-strand cDNA synthesis was performed with the BioRT All-in-One RT Master Mix for qPCR (Bioer Technology, Hangzhou, China). We amplified the PCR products of HpHsfs splice variants with Green Taq Mix (Vazyme, Nanjing, China) and resolved them on a 2% agarose gel. The corresponding primer sequences can be found in Supplementary Table 1. All assays were repeated at least twice and yielded concordant results.

2.8 Subcellular localization analysis

The open reading frames (ORFs) of HpHsf3, HpHsf5, and HpHsf17 without stop codons were amplified using primers containing BamH I and Kpn I sites (Supplementary Table 1) and subcloned in-frame into the pCAMBIA1300-35S-YFP vector. The RFP-NLS_SV40 fusion protein was utilized to label the nucleus (Huang et al., 2014), while the empty vector served as a negative control. All constructs were individually transformed into the Agrobacterium tumefaciens strain GV3101. Following the protocol described by (Li et al., 2023b), Agrobacterium strains harboring the plasmids were infiltrated into Nicotiana benthamiana leaves, and fluorescence signals were observed under a laser confocal microscope two days post-transient expression (Nikon Ti2-E AXR+NSPARC, Tokyo, Japan).

3 Results

3.1 Identification and characterization of the physical properties of Hsf gene family members in H. pedunculosum

A total of 21 HpHsf genes were identified in the H. pedunculosum genome and designated HpHsf1 through HpHsf21 according to their chromosomal starting positions (Table 1, Figure 1). These genes were unevenly distributed across eight of the 10 chromosomes. Analysis of their physicochemical properties revealed predicted ORF lengths ranging from 549 bp (HpHsf6) to 1566 bp (HpHsf11), corresponding to protein lengths of 182–521 amino acids. The molecular weights of the encoded proteins varied from 21.43 kDa to 56.99 kDa, while their theoretical isoelectric points (pI) ranged from 4.68 (HpHsf1) to 9.23 (HpHsf8). All HpHsf proteins exhibited negative grand average of hydropathicity (GRAVY) values, suggesting hydrophilic characteristics. Subcellular localization analysis further predicted that all HpHsf proteins are localized in the nucleus.

Table 1

Table 1. The detailed information of HpHsf members.

Figure 1

Figure 1. Chromosomal distribution of HpHsfs in H. pedunculosum. The ruler on the left represents the length of the chromosome. Regions of higher gene density are depicted in blue, contrasting with lower-density regions rendered in orange.

3.2 Chromosomal localization and phylogenetic analysis of the HpHsf gene family

Chromosomal distribution analysis revealed that the 21 HpHsf genes are unevenly distributed across chromosomes 1 and 3 to 9. Among these, chromosome 8 exhibited the highest number, with five HpHsf genes, whereas chromosome 3 harbored only a single gene (Figure 1). We performed a phylogenetic analysis of the Hsf gene family in H. pedunculosum using amino acid sequences from 21 HpHsf proteins, 21 AtHsf proteins (A. thaliana), 25 OsHsf proteins (O. sativa), and 23 CsHsf proteins (C. sativus) (Figure 2; Supplementary Table 2). Based on this analysis, Hsf proteins were classified into three subfamilies: HsfA, HsfB, and HsfC (Shamshad et al., 2023; Guo et al., 2008; Chen et al., 2021). Class A was further divided into nine subtypes (A1-A9) comprising 11 members, while Class B was segregated into five subtypes (B1-B5) containing eight members. Class C was separated into two subtypes (C1-C2), of which only C1 was represented, with a single member. The phylogenetic tree showed that multiple HpHsfs clustered closely with CsHsfs. In addition, members of the AtHsf, HpHsf, and CsHsf subfamilies were absent from subtype C2.

Figure 2

Figure 2. Phylogenetic analysis of Hsf gene family members in H. pedunculosum (Hp), A.thaliana (At), O. sativa (Os), and Cucumis sativus (Cs). Distinct colors represent different Hsf subfamilies.

3.3 Gene duplication and synteny analysis of HpHsf genes

Gene duplication is a key evolutionary mechanism in plants, contributing to genetic diversity and driving the development of novel traits such as stress resistance and floral complexity (Panchy et al., 2016). Intraspecific collinearity analysis identified eight segmental duplication pairs among the HpHsfs genes, including HpHsf2/HpHsf16, HpHsf3/HpHsf9, HpHsf3/HpHsf12, HpHsf7/HpHsf14, HpHsf9/HpHsf12, HpHsf14/HpHsf15, and HpHsf18/HpHsf20 (Figure 3A). Notably, HpHsf3, HpHsf9, and HpHsf12 exhibited linear correlations, while HpHsf14 correlated with HpHsf7 and HpHsf15. These findings suggest that gene duplication events have contributed to the expansion of the HpHsf family and likely played a significant role in its evolutionary diversification.

Figure 3

Figure 3. Collinearity analysis of H pedunculosum. (A) Collinearity among HpHsf members in H pedunculosum. Red lines highlight duplicated HpHsf gene pairs, whereas the gray line represents all duplicate gene pairs. (B) Comparative synteny analysis of Hsf genes between H pedunculosum, A thaliana, andC sativus. Background gray lines indicate collinear blocks between H pedunculosum and other plant genomes, with Hsf-related gene pairs highlighted in red.

To further investigate the evolutionary patterns of the HpHsf gene family, synteny analysis was conducted between H. pedunculosum and two representative species, A. thaliana and C. sativus (Figure 3B). The analysis identified 22 collinear gene pairs with A. thaliana and 32 with C. sativus, suggesting a closer phylogenetic relationship between H. pedunculosum and cucumber than with Arabidopsis. In addition, 15 HpHsf genes showed collinearity with both A. thaliana and C. sativus, indicating that these genes represent conserved orthologs inherited from a common ancestor. Notably, HpHsf9 and HpHsf12 were associated with two to three collinear gene pairs in both species, implying potential gene duplication events and a shared evolutionary origin for these loci.

3.4 Analysis of the HpHsf gene structures, conserved motifs, and conserved domains

We conducted a comprehensive characterization of the HpHsf family by analyzing phylogenetic relationships, motif patterns, gene structures, and conserved domains (Table 1; Figure 4). A total of ten conserved motifs were identified, of which motifs 2, 1, and 4 were present in all members (Figure 4B). Motifs 6, 7, and 8 were exclusively detected in Class A Hsfs, whereas Motifs 5 and 10 were specific to Class B, with the exception of HpHsf15. In contrast, Class B members lacked motif 3, which was consistently observed in both Class A and Class C. Further characterization of canonical Hsf motifs revealed distinct architectural features among the three classes (Table 1). Class A proteins were defined by a canonical suite of motifs (DBD, OD, NLS, AHA, and NES), with the exceptions of HpHsf1 and HpHsf19, which lacked the NLS and NES, respectively. The RD motif was ubiquitous in Class B members, except for HpHsf6 and HpHsf13. Class C members exhibited the simplest architecture, typically possessing only the DBD and OD, occasionally accompanied by an NLS. Taken together, these results demonstrate that the Hsf family in H. pedunculosum conforms to the fundamental structural framework typical of plant HSF proteins. Moreover, the distinct motif signatures observed among the three subfamilies reflect clear structural and functional divergence, supporting their phylogenetic classification.

Figure 4

Figure 4. Phylogenetic relationship, motif distribution, gene structure and conserved domain of HpHsf proteins. (A) Phylogenetic relationships of HpHsf proteins, with different colors representing distinct subgroups. (B) Motif composition of HpHsf members, with each of the 10 motifs depicted by uniquely colored boxes. (C) Gene structures of HpHsf gene family. (D) Conserved domain structures of HpHsf proteins.

Gene structure analysis revealed that HpHsf17 was the only member containing three exons, whereas the other HpHsf genes consisted of two exons (Figure 4C). Furthermore, conserved domain analysis showed that all Class A members contain the HSF domain, while Class B members display a broader array of domains, including HSF, HSF_DNA-bind, HSF1 superfamily, and BRLZ. This pattern implies that Class A domains have remained more conserved throughout evolution, whereas Class B members exhibit greater domain diversification (Figure 4D).

3.5 Analysis of cis−acting elements in the promoter of HpHsf genes

To investigate the potential regulatory mechanisms and functional roles of HpHsf genes, the 2,000 bp upstream sequences from the translation start site were analyzed using the PlantCARE database to identify putative cis-acting regulatory elements. The promoter analysis revealed that the HpHsf gene family contains various cis-acting elements associated with stress response (low-temperature responsiveness, defense and stress responsiveness, anoxic specific inducibility, and anaerobic induction), hormonal signaling pathways (auxin responsiveness, MeJA-responsiveness, gibberellin-responsiveness, abscisic acid responsiveness, zein metabolism regulation, and salicylic acid responsiveness), and developmental processes (light responsiveness, meristem expression, and endosperm expression) (Figure 5). Notably, all 21 HpHsf gene promoters contained multiple light-responsive elements, while the distribution of stress- and development-related elements varied among members. These findings suggest that the HpHsf gene family plays a crucial role in light-response regulation and possesses diverse regulatory potential in other biological processes.

Figure 5

Figure 5. Analysis of cis-acting elements in the HpHsf promoters. Different colored squares represent distinct cis-acting elements within the HpHsf promoter regions.

3.6 Expression patterns of HpHsfs under heat stress and subcellular localization of three HpHsfs

Extensive evidence indicates that Hsfs play key roles in plant responses to high temperatures (Kan et al., 2023; Ding et al., 2020). Here, we characterized the expression profiles of all 21 HpHsfs under heat stress at different time points (0, 1, 3, 6, and 12 h) using qRT-PCR (Figure 6). HpHsf2, HpHsf3, HpHsf5, and HpHsf17 exhibited transient upregulation, with transcript levels peaking at 1 h and declining by 3 h, indicating a rapid transcriptional response. Among these, HpHsf3 and HpHsf17 belong to the Class A subfamily, while HpHsf2 and HpHsf5 are categorized into the Class B subfamily. Notably, the expression of HpHsf3 and HpHsf17 showed signs of reactivation at 12 h after heat stress. These results imply that these genes may function as critical regulators in the heat stress response of H. pedunculosum. To further validate their functional roles, we examined the subcellular localization of HpHsf3, HpHsf5, and HpHsf17 (Figure 7). The YFP-tagged HpHsf proteins were exclusively localized in the nucleus, colocalizing with a nuclear marker. In contrast, the free YFP control was distributed throughout the entire cell. These results confirm the nuclear localization of these proteins, which is consistent with the predicted localization patterns.

Figure 6

Figure 6. The relative expression levels of HpHsfs in response to heat stress. Transcript levels of HpHsfs were analyzed by qRT-PCR in H. pedunculosum leaves subjected to heat stress (42 °C) at different time points (0 h, 1 h, 3 h, 6 h, and 12 h). Samples collected at 0 h served as the normalization control. Data are mean ± SD from three or four independent experiments. Significance was determined via one-way ANOVA with Tukey’s multiple comparisons test. Different letters indicate significant differences at p < 0.05. All experiments were performed with two biological replicates.

Figure 7

Figure 7. Subcellular localization of HpHsf3, HpHsf5, and HpHsf17 in N. benthamiana leaves. The transient expression of p35S::HpHsf3-YFP, p35S::HpHsf5-YFP, and p35S::HpHsf17-YFP, was observed in N. benthamiana leaves, using the YFP empty vector as a control. The RFP-NLS_SV40 fusion protein served as a nuclear marker. Scale bars, 50 μm.

To investigate alternative splicing of HpHsf genes under heat stress, we analyzed their transcript variants by RT-PCR (Supplementary Figure 1). Multiple amplicons were detected for HpHsf1, HpHsf3, HpHsf5, HpHsf6, HpHsf9, HpHsf10, and HpHsf19, suggesting that these genes potentially undergo alternative splicing in H. pedunculosum. However, no heat stress-induced splicing variants were observed. This indicates that alternative splicing is unlikely to be a primary regulatory mechanism for these HpHsf genes in response to heat stress.

4 Discussion

Given the pivotal role of the Hsf gene family in mediating plant thermotolerance, its composition and functions have been extensively characterized in several model plants and major crops, including Arabidopsis (Guo et al., 2008), rice (Shamshad et al., 2023), wheat (Zhou et al., 2019), and maize (Lin et al., 2011). Systematic identification and analysis of Hsf genes have been reported within the cucurbit family in species such as cucumber (Chen et al., 2021) and pumpkin (Shen and Yuan, 2020). However, in the medicinal cucurbit H. pedunculosum, the Hsf gene family remains largely unexplored, and its genomic organization and expression patterns under heat stress are still unclear.

In this study, 21 members of the HpHsf family were identified in H. pedunculosum and designated HpHsf1 to HpHsf21 according to their chromosomal locations (Table 1). Compared with cucumber (23 genes) and pumpkin (36 genes), H. pedunculosum harbors a smaller number of Hsf genes. This difference is likely due to variations in whole-genome duplication events, followed by lineage-specific gene loss or retention (Jiang et al., 2013). We further characterized the physicochemical properties of HpHsf (Table 1). The isoelectric point (pI) of HpHsf proteins ranged from 4.68 to 9.23, a distribution comparable to that reported for CsHsf proteins from cucumber (Chen et al., 2021) and CmHsf proteins from pumpkin (Shen and Yuan, 2020). This conserved pI range suggests that the charge variability may be a common feature of the Hsf proteins across cucurbitaceous plants. Furthermore, all HpHsf proteins exhibited negative GRAVY values, suggesting their hydrophilic nature and likely solubility in the aqueous cellular environment. This observation aligns with the water-soluble properties of Hsf proteins reported in other plant species (Liao et al., 2023; Li et al., 2025). All HpHsfs proteins are predicted to localize in the nucleus, consistent with the canonical subcellular distribution of transcription factors.

Phylogenetic analyses have shown that the plant Hsf gene family is evolutionarily conserved and can be classified into three major groups: A, B, and C (Guo et al., 2016; Scharf et al., 2012). In cucumber, the CsHsf family contains 23 genes, categorized into classes A (12), B (9), and C (2), with each class comprising several phylogenetically distinct subclades (Chen et al., 2021). In pumpkin, the classification is more intricate: 36 CmHsf genes are divided into three subfamilies (I, II, and III) according to phylogenetic analysis, which correspond to the conventional groups A, B, and C based on conserved structural characteristics. Specifically, subfamily II aligns with group A, subfamily III corresponds to group B, and subfamily I includes genes from groups A and C (Shen and Yuan, 2020). Herein, phylogenetic analysis of Hsf members from H. pedunculosum, A. thaliana, O. sativa, and C. sativus revealed a conserved clustering pattern, with genes grouped into three major classes (A, B, and C) and 15 subgroups (A1-A9, B1-B5, and C1) (Figure 2). This classification is consistent with previous reports in H. pedunculosum and cucumber (Chen et al., 2021). The results suggest that orthologous genes with similar motif composition may exhibit functional redundancy, while paralogous genes within the same lineage may retain overlapping functions. Notably, no Class C2 subfamily members were identified in H. pedunculosum. This indicates that C2-type Hsfs are absent in this dicot species, a finding consistent with phylogenetic patterns where C2-type Hsfs are typically restricted to monocots (Andrási et al., 2020). Collectively, these results reinforce the link between the functional diversification of Hsf subfamilies and the evolutionary history of plant taxa, offering valuable insights for future research on the molecular mechanisms underlying stress resistance in H. pedunculosum.

Eukaryotic genomes vary in the extent to which genes are retained on homologous chromosomes (synteny) and in the preservation of their linear order (collinearity) over evolutionary timescales (Tang et al., 2008). Gene duplication promotes genetic diversification within species, while collinearity analyses provide insights into the evolutionary trajectories and functional divergence of genes. This study identified eight segmentally duplicated pairs among the HpHsf genes in H. pedunculosum, encompassing 10 genes distributed across eight chromosomes (Figure 3A). This number is slightly lower than that reported in Cucurbita moschata (nine pairs) and C. sativus (13 pairs) (Shen and Yuan, 2020; Chen et al., 2021). These results suggest that segmental duplication has likely contributed significantly to the expansion of the Hsf gene family in H. pedunculosum and potentially across the Cucurbitaceae family. Moreover, collinearity analysis revealed 22 and 32 homologous gene pairs between HpHsf genes and those of A. thaliana and C. sativus, respectively (Figure 3B). These conserved syntenic relationships highlight the presence of shared ancestral genes and provide valuable clues for inferring the functional roles of HpHsf genes.

Plant Heat Shock Factors (HSFs) have a modular structure, typically consisting of several functional domains (Guo et al., 2016). Here, our detailed characterization of canonical Hsf motifs reveals a distinct modular architecture that underscores the functional divergence among the three phylogenetic classes (Table 1). The comprehensive suite of motifs in Class A proteins, including the critical AHA activation domain, is consistent with their established role as potent transcriptional activators of the heat stress response (Döring et al., 2000). Notably, HpHsf1 and HpHsf19 lack a canonical NLS and NES, respectively. Given that a defective NES causes nuclear retention in SlHsfA2 (Heerklotz et al., 2001), the absence of these signals in HpHsf1 and HpHsf19 suggests their subcellular localization may be regulated by alternative mechanisms. In Class B, the RD repressor motif (-LFGV-) is conserved at the C-terminus but is notably absent in HpHsf6 and HpHsf13. Phylogenetic analysis placed HpHsf6 in the HsfB5 subfamily (Figure 2). This subfamily includes members like StHsfB5 from potato, which lack the RD motif and act as coactivators in heat stress responses (Zhu et al., 2023). Given this and the close evolutionary relationship between HpHsf6 and HpHsf13 (Figure 3), both are predicted to function as transcriptional activators regulating stress-responsive genes. However, this prediction requires experimental validation.

Analysis of the motif composition, domain architecture, and gene structure revealed that members within the same subfamily exhibit highly conserved characteristics, indicative of conserved functional roles. The distribution of motifs and conserved domains among HpHsf members follows distinct subfamily-specific patterns. For example, Class A members exclusively contain motifs 6, 7, or 8 and consistently harbor the HSF_DNA-bind domain. In contrast, Class B members are characterized by the presence of motifs 5 or 10, with the exception of HpHsf15 (Figure 4). Similar subfamily-specific patterns have also been reported in pumpkin (Shen and Yuan, 2020), maize (Lin et al., 2011), and tea plant (Li et al., 2025). Furthermore, analysis of cis-acting elements in the promoter regions provided insights into the potential regulatory mechanisms of HpHsf genes (Hernandez-Garcia and Finer, 2014). Hsfs are integral components of complex signaling networks that orchestrate plant responses to diverse abiotic stresses, particularly heat (Andrási et al., 2020). They also act as modulators of transcriptional dynamics, fine-tuning stress responses while maintaining a balance between adaptation and normal developmental processes (Bakery et al., 2024). Consistent with these roles, our analysis revealed that HpHsf promoters are enriched in cis-acting elements associated not only with stress responses, but also with hormone signaling and developmental processes (Figure 5). Research has shown that red and blue light regulate Hsf gene expression in species such as medicinal plant Astragalus mongholicus and Cannabis sativa L., with corresponding effects on cannabinoids accumulation. In addition, light-responsive cis-elements have been identified in the promoters of these Hsf genes (Wang et al., 2024; Qian et al., 2023). Similarly, promoters of all Hphsfs in H. pedunculosum also contain such light-responsive elements, suggesting that their expression may likewise be regulated by light. It is noteworthy that alpine species such as H. pedunculosum are routinely exposed to concurrent low-temperature and high-light stress. Our cis-regulatory element analysis revealed that the HpHsf genes not only harbor light-responsive elements but also contain low-temperature response elements in seven members (33%). In contrast, while six (26%) of the cucumber CsHsfs possess low-temperature responsive elements, no notable distribution of light-responsive elements has been reported in their promoters (Chen et al., 2021). The prevalence of both element types in H. pedunculosum suggests a potential coordinated regulatory mechanism for adapting to the complex high-altitude environment. However, these findings are currently based on in silico predictions, and the functional role of these cis-elements requires experimental validation. These findings underscore the multifaceted regulatory roles of HpHsf genes in integrating environmental cues with endogenous developmental programs.

High temperatures represent a significant threat to the survival of the alpine plant H. pedunculosum. To elucidate the roles of HpHsf genes under heat stress, we examined the transcriptional expression profiles of all family members following exposure to 42 °C stress at multiple time points. Two Class A genes (HpHsf3 and HpHsf17) and two Class B genes (HpHsf2 and HpHsf5) were rapidly upregulated within 1 hour of treatment, followed by a noticeable decline by 3 hours (Figure 6). Notably, the two Class A members showed further induction at 12 h, suggesting that multiple Hsf genes contribute to thermoregulation, with Class A members potentially mediating multi-phase regulatory functions. Interestingly, homologs of HpHsf3 and HpHsf17 in cucumber (CsHsf23 and CsHsf7, respectively) also demonstrated rapid responses to high-temperature stress (Chen et al., 2021), implying that these members may have conserved regulatory functions across cucurbitaceous species. Furthermore, the nuclear localization of HpHsf3, HpHsf5, and HpHsf17 provides cytological evidence supporting their roles as transcription factors (Figure 7).

5 Conclusion

In summary, this study identified 21 Hsf genes in the H. pedunculosum genome through comprehensive genomic analysis. These genes were systematically characterized by their protein physicochemical properties, chromosomal locations, phylogenetic relationships, conserved motifs, exon-intron structures, protein domains, and cis-regulatory elements. Moreover, we investigated the expression patterns of all HpHsf genes under heat stress and confirmed the subcellular localization of three representative HpHsf proteins via transient expression in tobacco epidermal cells. Collectively, these findings offer valuable insights into the HpHsf gene family and establish a robust foundation for future functional studies to elucidate their roles in heat stress response.

Data availability statement

The raw genomic dataset of H. pedunculosum presented in the study are deposited in the Figshare repository, accession number: https://doi.org/10.6084/m9.figshare.21626153.v2. The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

ZZ: Conceptualization, Methodology, Writing – original draft. HL: Investigation, Writing – original draft. MS: Conceptualization, Methodology, Writing – review & editing. YY: Methodology, Resources, Writing – review & editing. XQ: Resources, Writing – review & editing. QZ: Resources, Supervision, Writing – review & editing. RL: Supervision, Writing – review & editing. YT: Conceptualization, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was funded by the Natural Science Foundation of Sichuan (Grant No.: 2025ZNSFSC1125), the Open Project Program (Grant No.: SKL-KF202402) of State Key Laboratory of Crop Gene Exploration and Utilization in Southwest China, and the Sichuan Province “three districts” talent support project.

Conflict of interest

Author XQ was employed by the company Sichuan Liangyuan Food Co., Ltd.

The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1701203/full#supplementary-material

References

Andrási, N., Pettkó-Szandtner, A., and Szabados, L. (2020). Diversity of plant heat shock factors: regulation, interactions, and functions. J. Exp. Bot. 72, 1558–1575. doi: 10.1093/jxb/eraa576

PubMed Abstract | Crossref Full Text | Google Scholar

Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). Meme suite: tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208. doi: 10.1093/nar/gkp335

PubMed Abstract | Crossref Full Text | Google Scholar

Bakery, A., Vraggalas, S., Shalha, B., Chauhan, H., Benhamed, M., and Fragkostefanakis, S. (2024). Heat stress transcription factors as the central molecular rheostat to optimize plant survival and recovery from heat stress. New Phytol. 244, 51–64. doi: 10.1111/nph.20017

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, D., Sun, M., Yang, Y., Tan, B., Ren, D., Tao, Y., et al. (2025). Comprehensive genome-wide analysis and functional characterization of the dir gene family in herpetospermum pedunculosum: insights from HpDIR16 and HpDIR17. Plant Physiol. Biochem. 226, 110074. doi: 10.1016/j.plaphy.2025.110074

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, X., Wang, Z., Tang, R., Wang, L., Chen, C., and Ren, Z. (2021). Genome-wide identification and expression analysis of Hsf and Hsp gene families in cucumber (Cucumis sativus L.). Plant Growth Regul. 95, 223–239. doi: 10.1007/s10725-021-00739-z

Crossref Full Text | Google Scholar

Chen, C., Wu, Y., Li, J., Wang, X., Zeng, Z., Xu, J., et al. (2023). Tbtools-ii: A “One for all, all for one. Bioinf. Platform For Biol. Big-Data Mining. Mol. Plant 16, 1733–1742. doi: 10.1016/j.molp.2023.09.010

PubMed Abstract | Crossref Full Text | Google Scholar

Ding, Y., Shi, Y., and Yang, S. (2020). Molecular regulation of plant responses to environmental temperatures. Mol. Plant 13, 544–564. doi: 10.1016/j.molp.2020.02.004

PubMed Abstract | Crossref Full Text | Google Scholar

Döring, P., Treuter, E., Kistner, C., Lyck, R., Chen, A., and Nover, L. (2000). The role of AHA motifs in the activator function of tomato heat stress transcription factors HsfA1 and HsfA2. Plant Cell 12, 265–278. doi: 10.1105/tpc.12.2.265

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, Q. M., Zhang, H., Cao, Y., and Wang, C. (2007). Anti-inflammatory and free radical scavenging activities of ethanol extracts of three seeds used as “Bolengguazi. J. Ethnopharmacol. 114, 61–65. doi: 10.1016/j.jep.2007.07.024

PubMed Abstract | Crossref Full Text | Google Scholar

Fragkostefanakis, S., Simm, S., El-Shershaby, A., Hu, Y., Bublak, D., Mesihovic, A., et al. (2019). The repressor and co-activator HsfB1 regulates the major heat stress transcription factors in tomato. Plant Cell Environ. 42, 874–890. doi: 10.1111/pce.13434

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, M., Liu, J. H., Ma, X., Luo, D. X., Gong, Z. H., and Lu, M. H. (2016). The plant heat stress transcription factors (Hsfs): structure, regulation, and function in response to abiotic stresses. Front. Plant Sci. 7, 114. doi: 10.3389/fpls.2016.00114

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, J., Wu, J., Ji, Q., Wang, C., Luo, L., Yuan, Y., et al. (2008). Genome-wide analysis of heat shock transcription factor families in rice and Arabidopsis. J. Genet. Genomics 35, 105–118. doi: 10.1016/S1673-8527(08)60016-8

PubMed Abstract | Crossref Full Text | Google Scholar

Heerklotz, D., Döring, P., Bonzelius, F., Winkelhaus, S., and Nover, L. (2001). The balance of nuclear import and export determines the intracellular distribution and function of tomato heat stress transcription factor HsfA2. Mol. Cell. Biol. 21, 1759–1768. doi: 10.1128/MCB.21.5.1759-1768.2001

PubMed Abstract | Crossref Full Text | Google Scholar

Hernandez-Garcia, C. M. and Finer, J. J. (2014). Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 217-218, 109–119. doi: 10.1016/j.plantsci.2013.12.007

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, Y.-Y., Shi, Y., Lei, Y., Li, Y., Fan, J., Xu, Y.-J., et al. (2014). Functional identification of multiple nucleocytoplasmic trafficking signals in the broad-spectrum resistance protein RPW8.2. Planta 239, 455–468. doi: 10.1007/s00425-013-1994-x

PubMed Abstract | Crossref Full Text | Google Scholar

Ikeda, M., Mitsuda, N., and Ohme-Takagi, M. (2011). Arabidopsis HsfB1 and HsfB2b act as repressors of the expression of heat-inducible hsfs but positively regulate the acquired thermotolerance. Plant Physiol. 157, 1243–1254. doi: 10.1104/pp.111.179036

PubMed Abstract | Crossref Full Text | Google Scholar

Ikeda, M. and Ohme-Takagi, M. (2009). A novel group of transcriptional repressors in Arabidopsis. Plant Cell Physiol. 50, 970–975. doi: 10.1093/pcp/pcp048

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, W. K., Liu, Y. L., Xia, E. H., and Gao, L. Z. (2013). Prevalent role of gene features in determining evolutionary fates of whole-genome duplication duplicated genes in flowering plants. Plant Physiol. 161, 1844–1861. doi: 10.1104/pp.112.200147

PubMed Abstract | Crossref Full Text | Google Scholar

Kan, Y., Mu, X.-R., Gao, J., Lin, H.-X., and Lin, Y. (2023). The molecular basis of heat stress responses in plants. Mol. Plant 16, 1612–1634. doi: 10.1016/j.molp.2023.09.013

PubMed Abstract | Crossref Full Text | Google Scholar

Lamesch, P., Berardini, T. Z., Li, D., Swarbreck, D., Wilks, C., Sasidharan, R., et al. (2011). The Arabidopsis Information Resource (Tair): improved gene annotation and new tools. Nucleic Acids Res. 40, D1202–D1210. doi: 10.1093/nar/gkr1090

PubMed Abstract | Crossref Full Text | Google Scholar

Lescot, M., Déhais, P., Thijs, G., Marchal, K., Moreau, Y., Van De Peer, Y., et al. (2002). PlantCARE, A database of plant cis-acting regulatory elements and A portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi: 10.1093/nar/30.1.325

PubMed Abstract | Crossref Full Text | Google Scholar

Li, B., Jiang, S., Gao, L., Wang, W., Luo, H., Dong, Y., et al. (2024). Heat shock factor A1s are required for phytochrome-interacting factor 4-mediated thermomorphogenesis in arabidopsis. J. Integr. Plant Biol. 66, 20–35. doi: 10.1111/jipb.13579

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Q., Li, H., Huang, W., Xu, Y., Zhou, Q., Wang, S., et al. (2019). A chromosome-scale genome assembly of cucumber (Cucumis sativus L.). Gigascience 8, Giz072. doi: 10.1093/gigascience/giz072

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Lu, Q., Peng, M., Liao, J., Zhang, B., Yang, D., et al. (2023a). Water extract from Herpetospermum pedunculosum attenuates oxidative stress and ferroptosis induced by acetaminophen via regulating Nrf2 and NF-κb pathways. J. Ethnopharmacol. 305, 116069. doi: 10.1016/j.jep.2022.116069

PubMed Abstract | Crossref Full Text | Google Scholar

Li, G., Shi, X., Lin, Q., Lv, M., Chen, J., Wen, Y., et al. (2025). Genome-wide identification and expression analysis of heat shock transcription factors in camellia sinensis under abiotic stress. Plants 14, 697. doi: 10.3390/plants14050697

PubMed Abstract | Crossref Full Text | Google Scholar

Li, R., Xiao, M., Li, J., Zhao, Q., Wang, M., and Zhu, Z. (2023b). Transcriptome analysis of cyp450 family members in Fritillaria cirrhosa D. Don and profiling of key cyp450s related to isosteroidal alkaloid biosynthesis. Genes 14, 219. doi: 10.3390/genes14010219

PubMed Abstract | Crossref Full Text | Google Scholar

Liao, W., Tang, X., Li, J., Zheng, Q., Wang, T., Cheng, S., et al. (2023). Genome wide investigation of hsf gene family in phoebe bournei: identification, evolution, and expression after abiotic stresses. J. For. Res. 35, 11. doi: 10.1007/s11676-023-01661-y

Crossref Full Text | Google Scholar

Lin, Y.-X., Jiang, H.-Y., Chu, Z.-X., Tang, X.-L., Zhu, S.-W., and Cheng, B.-J. (2011). Genome-wide identification, classification and analysis of heat shock transcription factor family in maize. BMC Genomics 12, 76. doi: 10.1186/1471-2164-12-76

PubMed Abstract | Crossref Full Text | Google Scholar

Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Gustavo, A., Sonnhammer, E. L. L., et al. (2020). Pfam: the protein families database in 2021. Nucleic Acids Res. 49, D412–D419. doi: 10.1093/nar/gkaa913

PubMed Abstract | Crossref Full Text | Google Scholar

Ohama, N., Sato, H., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2017). Transcriptional regulatory network of plant heat stress response. Trends Plant Sci. 22, 53–65. doi: 10.1016/j.tplants.2016.08.015

PubMed Abstract | Crossref Full Text | Google Scholar

Panchy, N., Lehti-Shiu, M., and Shiu, S.-H. (2016). Evolution of gene duplication in plants. Plant Physiol. 171, 2294–2316. doi: 10.1104/pp.16.00523

PubMed Abstract | Crossref Full Text | Google Scholar

Pick, T., Jaskiewicz, M., Peterhansel, C., and Conrath, U. (2012). Heat shock factor hsfb1 primes gene transcription and systemic acquired resistance in arabidopsis. Plant Physiol. 159, 52–55. doi: 10.1104/pp.111.191841

PubMed Abstract | Crossref Full Text | Google Scholar

Potter, S. C., Luciani, A., Eddy, S. R., Park, Y., Lopez, R., and Finn, R. D. (2018). Hmmer web server: 2018 update. Nucleic Acids Res. 46, W200–W204. doi: 10.1093/nar/gky448

PubMed Abstract | Crossref Full Text | Google Scholar

Qian, G., Meng, X., Wang, S., Mi, Y., Qin, Z., Liu, T., et al. (2023). Genome-wide identification of hsf gene family and their expression analysis in vegetative tissue of young seedlings of hemp under different light treatments. Ind. Crops Prod. 204, 117375. doi: 10.1016/j.indcrop.2023.117375

Crossref Full Text | Google Scholar

Raturi, V. and Zinta, G. (2024). Hsfa1 heat shock factors integrate warm temperature and heat signals in plants. Trends Plant Sci. 29, 1165–1167. doi: 10.1016/j.tplants.2024.07.002

PubMed Abstract | Crossref Full Text | Google Scholar

Sato, H., Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2024). Complex plant responses to drought and heat stress under climate change. Plant J. 117, 1873–1892. doi: 10.1111/tpj.16612

PubMed Abstract | Crossref Full Text | Google Scholar

Scharf, K.-D., Berberich, T., Ebersberger, I., and Nover, L. (2012). The plant heat stress transcription factor (Hsf) family: structure, function and evolution. Bba-Gene Regul. Mech. 1819, 104–119. doi: 10.1016/j.bbagrm.2011.10.002

PubMed Abstract | Crossref Full Text | Google Scholar

Shamshad, A., Rashid, M., and Zaman, Q. U. (2023). In-silico analysis of heat shock transcription factor (OsHSF) gene family in rice (Oryza sativa L.). BMC Plant Biol. 23, 395. doi: 10.1186/s12870-023-04399-1

PubMed Abstract | Crossref Full Text | Google Scholar

Shen, C. and Yuan, J. (2020). Genome-wide characterization and expression analysis of the heat shock transcription factor family in pumpkin (Cucurbita moschata). BMC Plant Biol. 20, 471. doi: 10.1186/s12870-020-02683-y

PubMed Abstract | Crossref Full Text | Google Scholar

Tamura, K., Stecher, G., and Kumar, S. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38, 3022–3027. doi: 10.1093/molbev/msab120

PubMed Abstract | Crossref Full Text | Google Scholar

Tan, W., Chen, J., Yue, X., Chai, S., Liu, W., Li, C., et al. (2023). The heat response regulators hsfa1s promote arabidopsis thermomorphogenesis via stabilizing pif4 during the day. Sci. Adv. 9, Eadh1738. doi: 10.1126/sciadv.adh1738

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, H., Bowers, J. E., Wang, X., Ming, R., Alam, M., and Paterson, A. H. (2008). Synteny and collinearity in plant genomes. Science 320, 486–488. doi: 10.1126/science.1153917

PubMed Abstract | Crossref Full Text | Google Scholar

Von Koskull-Döring, P., Scharf, K.-D., and Nover, L. (2007). The diversity of plant heat stress transcription factors. Trends Plant Sci. 12, 452–457. doi: 10.1016/j.tplants.2007.08.014

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Z., Wang, P., He, J., Kong, L., Zhang, W., Liu, W., et al. (2024). Genome-wide analysis of the HSF gene family reveals its role in astragalus mongholicus under different light conditions. Biology 13, 280. doi: 10.3390/biology13040280

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y.-X. and Xu, J.-H. (2025). The heat shock transcription factors regulate response mechanisms to abiotic stresses in plants. Crop Design 4, 100109. doi: 10.1016/j.cropd.2025.100109

Crossref Full Text | Google Scholar

Wei, X., Ma, Y., Dong, Z., Wang, G., Lan, X., Liao, Z., et al. (2021). Dehydrodiconiferyl alcohol, A lignan from Herpetospermum pedunculosum, alleviates cholestasis by activating pathways associated with the farnesoid X receptor. Phytomedicine 80, 153378. doi: 10.1016/j.phymed.2020.153378

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, Z., Li, T., Ding, L. P., Wang, C. P., Teng, R. D., Xu, S. J., et al. (2024). Lily LlHSFC2 coordinates with HSFAs to balance heat stress response and improve thermotolerance. New Phytol. 241, 2124–2142. doi: 10.1111/nph.19507

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, J., Chen, Y., Cai, G., Cai, R., Hu, Z., and Wang, H. (2023). Tree visualization by one tab le (Tvbot): A web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res. 51, W587–W592. doi: 10.1093/nar/gkad359

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, S., Li, J., Zhang, X., Wei, H., and Cui, L. (2006). Effects of heat acclimation pretreatment on changes of membrane lipid peroxidation, antioxidant metabolites, and ultrastructure of chloroplasts in two cool-season turfgrass species under heat stress. Environ. Exp. Bot. 56, 274–285. doi: 10.1016/j.envexpbot.2005.03.002

Crossref Full Text | Google Scholar

Yang, Y., Zhang, B., Bao, Y., Huang, P., Li, J., Li, R., et al. (2023). Chromosome-level genome assembly of Herpetospermum pedunculosum (Cucurbitaceae). Genome Biol. Evol. 15, Evad005. doi: 10.1093/gbe/evad005

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Xu, F., Liu, J., Guan, F., Quan, H., and Meng, F. (2019). The adaptation strategies of Herpetospermum pedunculosum (Ser.) baill at altitude gradient of the tibetan plateau by physiological and metabolomic methods. BMC Genomics 20, 451. doi: 10.1186/s12864-019-5778-y

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, M., Zheng, S., Liu, R., Lu, J., Lu, L., Zhang, C., et al. (2019). Genome-wide identification, phylogenetic and expression analysis of the heat shock transcription factor family in bread wheat (Triticum aestivum L.). BMC Genomics 20, 505. doi: 10.1186/s12864-019-5876-x

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, Z., Bao, Y., Yang, Y., Zhao, Q., and Li, R. (2024a). Research progress on heat stress response mechanism and control measures in medicinal plants. Int. J. Mol. Sci. 25, 8600. doi: 10.3390/ijms25168600

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, Z., Chen, D., Sun, M., Xiao, M., Huang, P., Ren, D., et al. (2024b). Integrative analysis of the metabolome and transcriptome provides insights into the mechanisms of lignan biosynthesis in Herpetospermum pedunculosum (Cucurbitaceae). BMC Genomics 25, 421. doi: 10.1186/s12864-024-10306-1

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, W., Xue, C., Chen, M., and Yang, Q. (2023). StHsfB5 promotes heat resistance by directly regulating the expression of hsp genes in potato. Int. J. Mol. Sci. 24, 16528. doi: 10.3390/ijms242216528

PubMed Abstract | Crossref Full Text | Google Scholar