A Heat Shock Transcription Factor TrHSFB2a of White Clover Negatively Regulates Drought, Heat and Salt Stress Tolerance in Transgenic Arabidopsis

Heat shock transcription factors (HSF) are divided into classes A, B and C. Class A transcription factors are generally recognized as transcriptional activators, while functional characterization of class B and C heat shock transcription factors have not been fully developed in most plant species. We isolated and characterized a novel HSF transcription factor gene, TrHSFB2a (a class B HSF) gene, from the drought stress-sensitive forage crop species, white clover (Trifolium repens). TrHSFB2a was highly homologous to MtHSFB2b, CarHSFB2a, AtHSFB2b and AtHSFB2a. The expression of TrHSFB2a was strongly induced by drought (PEG6000 15% w/v), high temperature (35 °C) and salt stresses (200 mM L−1 NaCl) in white clover, while subcellular localization analysis showed that it is a nuclear protein. Overexpression of the white clover gene TrHSFB2a in Arabidopsis significantly reduced fresh and dry weight, relative water contents (RWC), maximum photosynthesis efficiency (Fv/Fm) and performance index on the absorption basis (PIABS), while it promoted leaf senescence, relative electrical conductivity (REC) and the contents of malondialdehyde (MDA) compared to a wild type under drought, heat and salt stress conditions of Arabidopsis plants. The silencing of its native homolog (AtHSFB2a) by RNA interference in Arabidopsis thaliana showed opposite trends by significantly increasing fresh and dry weights, RWC, maximum photosynthesis efficiency (Fv/Fm) and performance index on the absorption basis (PIABS) and reducing REC and MDA contents under drought, heat and salt stress conditions compared to wild type Arabidopsis plants. These phenotypic and physiological indicators suggested that the TrHSFB2a of white clover functions as a negative regulator of heat, salt and drought tolerance. The bioinformatics analysis showed that TrHSFB2a contained the core B3 repression domain (BRD) that has been reported as a repressor activator domain in other plant species that might repress the activation of the heat shock-inducible genes required in the stress tolerance process in plants. The present study explores one of the potential causes of drought and heat sensitivity in white clover that can be overcome to some extent by silencing the TrHSFB2a gene in white clover.


Introduction
White clover, belonging to the genus Trifolium, is an excellent worldwide cultivated forage due to its high-yield characteristics, better forage quality and strong biological nitrogen fixation ability. It is an integral part of the grass industry. The development of white clover can promote grain conversion to feed and encourage grass and animal products in competitive markets. However, it is vulnerable to drought and heat stresses silenced by RNA interference (RNAi) to determine how HSFB2a protein responds under drought, heat and salt stress conditions. The current study results would improve the understanding of the molecular mechanism of the responses to extreme environmental conditions in white clover and explore the regulatory role of HSFs during ecological stress responses in plants.

TrHSFB2a Was Predicated as One of the Class B HSF Proteins
The sequence of the isolated TrHSFB2a gene comprises an open reading frame (ORF) of 924 bp (Figure 1a), which constitutes a protein of 307 amino acids (Figure 1c). Interprotein scan results showed that this protein belonged to the heat shock transcription family and was predicated as the heat stress transcription factor B-2a (HSFB2-a). It contained a heat shock factor-type, DNA-binding domain (HSF-DNA-bd) of 93 amino acids , commonly associated with DNA-binding transcription factor activity, sequence-specific DNA binding and regulation of DNA-dependent cellular transcription. The N-terminal region contained a helix-turn-helix motif (H2-T-H3) in the center of the DNA binding domain (DND) (Figure 1c). Based on the homology of amino acids, TrHSFB2a showed an overall 80.46%, 70.61% and 73.74% homology with TpHSFB2b, MtHSFB2b and CarHSFB2a, respectively ( Figure 1b). Additionally, amplification and sequencing of TrHSFB2a from white clover genomes showed it has one intron of 107bp, similar to other HSF genes of plants. Based on conserved domains of HSF genes, these are divided into Class A, B and C members. Phylogenetic tree constructed using 27 HSF proteins showed TrHSFB2a clustered with class B HSF proteins of closely related species and Arabidopsis, indicating that it is one of class B HSF proteins (Figure 1d).    . TrHSFB2  206  MtHSFB2b  221  AtHSFB2a  240  AtHSFB2b  202  AtHSFB1 Consensus l TrHSFB2  241  MtHSFB2b  253  AtHSFB2a  280  AtHSFB2b 238 AtHSFB1 Consensus TrHSFB2  280  MtHSFB2b  293  AtHSFB2a  320  AtHSFB2b 278 TrHSFB2  302  MtHSFB2b  299  AtHSFB2a  342  AtHSFB2b 284 AtHSFB1 Consensus Bioinformatics tools were used to understand physiochemical and structural characteristics of TrHSFB2a and results are presented in Supplementary Figure S1. Using pro-param predication [31], it was found that the molecular formula of TrHSFB2a was C1499H2383N425O475S11, molecular weight was 34,311.63, theoretical isoelectric point was 8.7, instability index was 49.92, and it was an unstable protein. DeepTMHMM [32] predicted that it was a globular protein and had no transmembrane structure (Supplementary Figure S1a). Hydrophobicity analysis using ProtScale revealed that it was a hy-drophobic protein (Supplementary Figure S1b), with an average hydropathicity −0.566. The secondary structure of the encoded protein predicted by SOPMA [33] showed that the secondary structure was composed of 33. 88% α helix, 2.93% β-turn, 55.70% random coil and 7.49% extended strand (Supplementary Figure S1d). The phosphorylation sites of protein sequence determined by Netphos 3.1 server [34] showed TrHSFB2a protein had 58 phosphorylation sites, including 32 serine sites, 19 threonine sites and 7 tyrosine sites (Supplementary Figure S1c). Protein signal peptide analysis using Signalp 4.1 server [35] showed that the protein had no signal peptide (Supplementary Figure S1e). The spatial protein structure predicted by SWISS-Mod. server [36] revealed that the three-dimensional structure had a monomer structure (Supplementary Figure S1f).

TrHSFB2a Is a Nuclear Protein
TrHSFB2a was predicated as transcription factors, and transcription factors generally function in the nucleus. To confirm this, we initially used the WoLF PSORT program and predicted its location in the nucleus, then ligated the coding region of TrHSFB2a to 3 end of GFP gene in pSUPER1300 containing 35S promoter, while GPP gene alone with promoter served as control. The subcellular location of TrHSFB2a protein was determined by transiently expressing TrHSFB2a::GFP fusion protein in the leaves of Nicotiana benthamiana. As expected, the TrHSFB2a::GFP protein is localized in the nucleus ( Figure 2).

Relative Expression of TrHSFB2a in White Clover under Different Conditions
A previous transcriptome study of our research group showed that TrHSFB2a is strongly expressed in white clover under drought stress [28]. After isolating TrHSFB2a from white clover, its expression pattern was studied. The four-week-old white clover plants were subjected to drought (PEG6000 15% w/v), heat (35 °C) and salt (NaCl 200 mM) stresses. Total RNA of whole plants (root + shoot) was isolated after 0 h, 3 h, 6 h, 12 h and 24 h and used to construct cDNA for determining relative gene expression levels. A relative expression of TrHSFB2a was determined under normal conditions (0 h) in white clover, indicating that this gene is also expressed under normal conditions, at least in one-

Relative Expression of TrHSFB2a in White Clover under Different Conditions
A previous transcriptome study of our research group showed that TrHSFB2a is strongly expressed in white clover under drought stress [28]. After isolating TrHSFB2a from white clover, its expression pattern was studied. The four-week-old white clover plants were subjected to drought (PEG6000 15% w/v), heat (35 • C) and salt (NaCl 200 mM) stresses. Total RNA of whole plants (root + shoot) was isolated after 0 h, 3 h, 6 h, 12 h and 24 h and used to construct cDNA for determining relative gene expression levels. A relative expression of TrHSFB2a was determined under normal conditions (0 h) in white clover, indicating that this gene is also expressed under normal conditions, at least in one-day-old white clover plants. Under polyethylene glycol (PEG6000 15% w/v) and salt stress of 3 h, its increased expression was observed, which kept increasing for 12 h under salt stress and, at least, up to 24 h under drought stress. However, the expression pattern of TrHSFB2a under heat stress was not the same. In heat stress, a strong expression of TrHSFB2a was detected after 3 h, and then it reduced continuously but remained relatively higher than the highest levels observed under salt and drought stresses ( Figure 3). These results indicated that TrHSFB2a was involved in abiotic stress responses in white clover and more strongly induced under heat stress.

Relative Expression of TrHSFB2a in White Clover under Different Conditions
A previous transcriptome study of our research group showed that TrHSFB2a is strongly expressed in white clover under drought stress [28]. After isolating TrHSFB2a from white clover, its expression pattern was studied. The four-week-old white clover plants were subjected to drought (PEG6000 15% w/v), heat (35 °C) and salt (NaCl 200 mM) stresses. Total RNA of whole plants (root + shoot) was isolated after 0 h, 3 h, 6 h, 12 h and 24 h and used to construct cDNA for determining relative gene expression levels. A relative expression of TrHSFB2a was determined under normal conditions (0 h) in white clover, indicating that this gene is also expressed under normal conditions, at least in oneday-old white clover plants. Under polyethylene glycol (PEG6000 15% w/v) and salt stress of 3 h, its increased expression was observed, which kept increasing for 12 h under salt stress and, at least, up to 24 h under drought stress. However, the expression pattern of TrHSFB2a under heat stress was not the same. In heat stress, a strong expression of TrHSFB2a was detected after 3 h, and then it reduced continuously but remained relatively higher than the highest levels observed under salt and drought stresses ( Figure  3). These results indicated that TrHSFB2a was involved in abiotic stress responses in white clover and more strongly induced under heat stress.

TrHSFB2a Negatively Regulates Drought Tolerance in Transgenic Overexpression Arabidopsis
The detection of TrHSFB2a expression under abiotic stress observing its role in stress responses in detail. Therefore, the overexpression vector of TrHSFB2a was constructed and transferred into Arabidopsis thaliana, and thereby, constructed the TrHSFB2a overexpression Arabidopsis. After developing homozygous T3 lines, the relative expression of TrHSFB2a was analyzed in transgenic Arabidopsis (Figure 4g), and two lines highly expressing TrHSFB2a were selected for a comparative analysis along with the wild type and T-DNA lines of Arabidopsis gene AtHSFB2a, purchased from Arabidopsis share platform "https://www.arashare.cn/index/ (accessed 18 December 2020)", in which TrHSFB2a has been silenced. Under normal conditions, there was no significant difference in plant root length among groups after 7 d of transplanting in 1 /2 MS medium (Figure 4d,e,h). However, the soil-grown overexpression Arabidopsis were significantly lower in fresh and dry weight compared to wild-type Arabidopsis, and T-DNA lines of Arabidopsis were significantly higher in fresh and dry weights compared to wild type (Figure 4i,j), suggesting that TrHSFB2a produced negative effects on early plant growth and development, which was also obvious in phenotypic observations (Figure 4f), suggesting that overexpression or silencing of TrHSFB2a at least do not cause detrimental effects on growth and developments of Arabidopsis. Two overexpression lines (OE3 and OE4), two T-DNA lines (denoted as T77 and T95) and the wild type (Wt) were subjected to different drought treatments. Initially, 5-day-old seedlings of Arabidopsis (wild type and transgenic) were transferred to ½ MS plates containing 200 or 300 mM mannitol. After two weeks, the overexpression lines (OE-3 and OE-4) under 200 mM mannitol showed a significantly higher leaf-yellowing rate compared to the wild type, while most of the plants of overexpression lines under 300 mM mannitol dried; the wild type although had curled leaves, all were green (Figure 5a,b). Moreover, root length was significantly higher in T- Initially, 5-day-old seedlings of Arabidopsis (wild type and transgenic) were transferred to 1 /2 MS plates containing 200 or 300 mM mannitol. After two weeks, the overexpression lines (OE-3 and OE-4) under 200 mM mannitol showed a significantly higher leaf-yellowing rate compared to the wild type, while most of the plants of overexpression lines under 300 mM mannitol dried; the wild type although had curled leaves, all were green (Figure 5a,b). Moreover, root length was significantly higher in T-DNA lines and was significantly lower in overexpression lines under 200 mM mannitol stress, indicating that TrHSFB2a possibly negatively regulates drought stress (Figure 5c-e). Three-week-old wildtype and transgenic Arabidopsis plants were also subjected to drought stress in nutrient soil by withholding water. Plants were well watered, and then, irrigation was stopped. After 12 days of water withholding, all the leaves of overexpression lines OE3 and OE4 wilted 7 of 19 and dried, while wild-type plants showed less wilting compared to the overexpression lines, while T-DNA lines did not exhibit leaf wilting at this stage ( Figure 6a). Three days after rewatering, more than 70% of wild-type plants recovered, while only about 30% of overexpression plants could survive. These observations indicated that the overexpression of TrHSFB2a significantly reduced drought tolerance, while suppressing its expression significantly increased drought tolerance in transgenic Arabidopsis at the seedling stage. Fresh weight, dry weight and photochemical efficiency determine how plants grow under stressed conditions. Relative water contents (RWC) of the plant showed their water status in terms of physiological consequences under water-deficient conditions. The cell membrane permeability and lipid peroxidation are the general indicators of damage in plants caused by different stresses and are measured as REC of electrolyte leakage and malondialdehyde (MDA) contents, respectively. Lipid peroxidation is generally associated with the oxidation of lipids of the cell membrane; thus, its quantification determines the extent of cell membrane damage. In the current study, the phytochemical efficiency of photosystem II and performance index was determined after 10 days of water withholding. Results showed that overexpression lines had significantly reduced photochemical efficiency of Photosystem II and the performance index compared to wild type, while T-DNA lines had significantly higher photochemical efficiency of photosystem II and performance index than that of wild type. Similarly, fresh weight, dry weight and RWC were significantly lower in overexpression lines compared to wild type and were significantly higher in T-DNA lines compared to wild type ( Figure 6). In contrast, REC and MDA contents were significantly higher in overexpressed lines than that of wild type and significantly lower in T-DNA lines than that of wild type ( Figure 6). These results, along with phenotypic data, demonstrated that TrHSFB2a negatively regulates drought stress, and RNA interference in T-DNA lines improved the drought tolerance in Arabidopsis.

TrHSFB2a Negatively Regulates the Heat Tolerance in Arabidopsis
Heat shock proteins and transcription factors are generally involved in regulating heat stress responses. One-month-old soil-grown plants were transferred to an incubator adjusted to 35 °C temperature for two weeks in order to validate the role of TrHSFB2a in Arabidopsis under high-temperature stress. Photochemical efficiency, REC and MDA contents were determined after one week, while plants' dying rate and dry weight were measured after two weeks when phenotypic differences appeared significantly. As shown in Figure 7, around 60% of plants of overexpression lines completely dried and could not recover when transferred in normal conditions. The rate of drying was significantly higher in overexpression lines than that of wild type. In comparison, T-DNA lines showed a

TrHSFB2a Negatively Regulates the Heat Tolerance in Arabidopsis
Heat shock proteins and transcription factors are generally involved in regulating heat stress responses. One-month-old soil-grown plants were transferred to an incubator adjusted to 35 • C temperature for two weeks in order to validate the role of TrHSFB2a in Arabidopsis under high-temperature stress. Photochemical efficiency, REC and MDA contents were determined after one week, while plants' dying rate and dry weight were measured after two weeks when phenotypic differences appeared significantly. As shown in Figure 7, around 60% of plants of overexpression lines completely dried and could not recover when transferred in normal conditions. The rate of drying was significantly higher in overexpression lines than that of wild type. In comparison, T-DNA lines showed a significantly lower drying rate compared to wild type (Figure 7a). Photochemical efficiency and the plant's dry weight were significantly lower for overexpression lines compared to wild type and significantly higher for T-DNA lines. In contrast, overexpression lines showed significantly increased REC % and MDA contents than that of wild type, and T-DNA lines had significantly reduced REC % and MDA contents compared to wild types, indicating that RNA interference lines had lower damage at the cellular level compared to wild type, while overexpression lines showed opposite trends. Collectively, phenotypic and physiological data demonstrated that TrHSFB2a negatively regulates heat tolerance in Arabidopsis at the seedling stage. significantly lower drying rate compared to wild type (Figure 7a). Photochemical efficiency and the plant's dry weight were significantly lower for overexpression lines compared to wild type and significantly higher for T-DNA lines. In contrast, overexpression lines showed significantly increased REC % and MDA contents than that of wild type, and T-DNA lines had significantly reduced REC % and MDA contents compared to wild types, indicating that RNA interference lines had lower damage at the cellular level compared to wild type, while overexpression lines showed opposite trends. Collectively, phenotypic and physiological data demonstrated that TrHSFB2a negatively regulates heat tolerance in Arabidopsis at the seedling stage.

TrHSFB2a Negatively Regulates Salt Tolerance in Arabidopsis
To check how the TrHSFB2a gene functions under salt stress conditions, we observed the overexpression and RNA interference T-DNA lines of Arabidopsis in 1 /2 MS media and in soil. At five days old, seedlings were shifted to 1/2MS+Agar plates, supplemented with 100 mM or 150 mM NaCl, and phenotypic data were recorded at 10 d of transplanting. The results showed that TrHSFB2a had a significantly lower rate of survival (around 20%) compared to wild type (47.33%), and RNA interference T-DNA lines showed a significantly higher survival rate (more than 65%) than that of wild type (Figure 8), suggesting that overexpression of TrHSFB2a significantly reduced tolerance of salt, while RNA interference of TrHSFB2a can enhance salt tolerance of transformed Arabidopsis at the seedling stage. Root length was significantly higher in T-DNA lines compared to wild type under 100 mM and 150 mM salt stress, while one overexpression line (OE-3) showed significantly reduced roots length than that of wild type, and one (OE-3) did not show a significant difference with the wild type under both salt stress treatments. However, the total fresh weight was significantly higher in T-DNA lines and significantly lower in overexpression lines compared to wild type, respectively ( Figure 8). Simultaneously, 10 d old wild and transgenic Arabidopsis lines grown on 1/2MS+Agar media were transferred into the soil, grown additionally for three weeks in nutrient soil, then irrigated with 100 mM, 200 mM and 300 mM NaCl gradients water, each twice after every two days (Figure 9a,b). RWC, EC, photochemical efficiency and performance index were measured after two times of 100 mM and 2 times 200 mM saline water irrigation, on the 8th d of treatment. Results showed that RWC, photochemical efficiency and performance index on an absorption basis were significantly lower in overexpression lines and were significantly higher in T-DNA lines compared to wild type, respectively. In contrast, EC % had significantly increased in overexpression lines and was significantly lower in T-DNA lines compared to wild type (Figure 9c-f). The data of physiological analysis, together with phenotypic observations, indicated that overexpression of TrHSFB2a reduced salt tolerance, while RNA interference for this gene can improve the salt tolerance of Arabidopsis at seedling and developing growth stages.

Discussion
Plant heat shock transcription factors (HSFs) play important roles in regulating biotic and abiotic stress responses, especially heat stress regulation [5,6]. HSFs can be divided into three classes (A, B and C) based on their domain characteristics, which have been reported to function differently in different stress responses [5]. The HSFs of class A have been generally reported as the positive regulator of stress tolerance in different plant species and possess an exclusive C terminal activation domain [16]. The class C HSFs have also been shown to play important functions in mediating the effects of heat, salt and osmotic stress in plants [37][38][39]. The class B HSFs including HSFB1 and HSFB2, show repressive activity and are reported as negative regulators of stress tolerance in Cicer arietinium [40], Oryza sativa [41], Glycine max [42] and Arabidopsis thialiana [26,27]. In the current study, we reported and characterized an HSF of white clover HSFB2a, which belongs to the Class B HSFs family. Current study data suggested that the TrHSFB2a might negatively regulate heat, drought and salt tolerance in transgenic Arabidopsis.

Abiotic Stress Treatments Induce TrHSFB2a Expression
The expression pattern of TrHSFB2a was investigated under different treatments of abiotic stresses, including high temperature, PEG and salt stress. The expression of TrHSFB2a was strongly induced under drought stress, which is in accordance with our previous study of transcriptome analysis of white clover under drought stress [28]. Moreover, TrHSFB2a is strongly upregulated under heat stress like other HSFs members,

Discussion
Plant heat shock transcription factors (HSFs) play important roles in regulating biotic and abiotic stress responses, especially heat stress regulation [5,6]. HSFs can be divided into three classes (A, B and C) based on their domain characteristics, which have been reported to function differently in different stress responses [5]. The HSFs of class A have been generally reported as the positive regulator of stress tolerance in different plant species and possess an exclusive C terminal activation domain [16]. The class C HSFs have also been shown to play important functions in mediating the effects of heat, salt and osmotic stress in plants [37][38][39]. The class B HSFs including HSFB1 and HSFB2, show repressive activity and are reported as negative regulators of stress tolerance in Cicer arietinium [40], Oryza sativa [41], Glycine max [42] and Arabidopsis thialiana [26,27]. In the current study, we reported and characterized an HSF of white clover HSFB2a, which belongs to the Class B HSFs family. Current study data suggested that the TrHSFB2a might negatively regulate heat, drought and salt tolerance in transgenic Arabidopsis.

Abiotic Stress Treatments Induce TrHSFB2a Expression
The expression pattern of TrHSFB2a was investigated under different treatments of abiotic stresses, including high temperature, PEG and salt stress. The expression of TrHSFB2a was strongly induced under drought stress, which is in accordance with our previous study of transcriptome analysis of white clover under drought stress [28]. Moreover, TrHSFB2a is strongly upregulated under heat stress like other HSFs members, as well as strongly expressed in response to NaCl stress. These observations are well in agreement with previous studies that HSFB2 candidates expressed under high temperature, high salinity and drought conditions [23,40,41]. Like stress-related cis-acting elements, including HSE, DRE, ABRE, MYCRS and MYBRS, present in the promoter region of the HSFB2 genes in Arabidopsis and other plant species [5,43,44], they might also be involved in inducing response and the expression of TrHSFB2a in white clover under drought, heat and high salinity stresses, but it needs to be further verified.

TrHSFB2a Negatively Regulates Drought, Heat and Salt Stress in Arabidopsis
The overexpression of TrHSFB2a in Arabidopsis, although it enhanced the transcriptome abundance and RNA interference of AtHSFB2a in Arabidopsis, reduced the accumulation of AtHSFB2a transcripts. There were no significant differences in plant stature, root length or fresh and dry weights under normal conditions. It indicated that TrHSFB2a does not have significant effects on plant development under normal conditions, as similarly observed in the effects of HSBs of class B on other plant species [40,41]. However, under drought, salt and heat stress, the overexpression of TrHSFB2a significantly reduced plant growth compared to wild type, while RNA interference of AtHSFB2a significantly improved salt, heat and drought tolerance both in 1 /2 MS media and in the soil. By increasing the stress levels (NaCl 150 mM and Mannitol 300 mM) or by extending stress duration, the overexpression of TrHSFB2a led to significantly higher drying and dying rates compared to wild type, while RNA interference of AtHSFB2a exhibited opposite trends. These observations agreeing has been observed in the transgenic lines developed by overexpression and RNA interference homolog of this OsHSFB2b in rice [41] and overexpression of CarHSFB2b of Chickpea in Arabidopsis [40].
RWC, photosynthetic activity, REC percentage and MDA contents are important indicators to assess the stress tolerance potential differences in plants. Under drought and heat stress treatment, the overexpression lines of TrHSFB2a showed higher REC % and MDA contents than that of wild type, and in T-DNA lines of AtHSFB2a, the accumulation of MDA contents and REC % was lower compared to wild type. As electrolyte leakage is an indirect indicator of cell membrane damage due to stress conditions [45,46], the higher REC in TrHSFB2a overexpression lines showed higher cell membrane damage, which had occurred there, while there was comparatively less cell membrane damage in T-DNA lines. One of the end products of lipid peroxidation due to external stresses is MDA, which is a result of membrane damage caused by free radicals [47,48], was higher in overexpression lines and lower in RNA interference lines compared to wild types, representing higher lipid damage in overexpression lines and a lower in T-DNA lines. Abiotic stresses on plants mainly reduce the photosynthetic efficiency of the plants under stress as a result of the negative effects of the abiotic stresses on the biosynthesis of photosynthetic pigments, photosystem performance, CO 2 fixation and gases exchanges, and on carbohydrate metabolism, etc. The photosynthetic activity and photosystem performance index significantly reduced overexpression lines and significantly improved in T-DNA lines compared to wild-type user drought, salt and heat stress, indicating better adaptability of T-DNA line and reduced fitness by overexpressing TrHSFB2a gene in Arabidopsis. Moreover, lower RWC, as observed in overexpression lines, can reduce or even completely stop photosynthetic activity in a stressed condition. Leaf wilting, drying, RWC, photosynthetic system performance, REC and MDA-related to plant stress that has been already used to assess the function of HSFs under drought, salt and heat stress in transgenic plants [41,49,50]. These findings provide evidence that TrHSFB2a negatively regulates stress tolerance at least in Arabidopsis.
The HSFB class has been reported to act as a transcriptional repressor in plant species [23,26,41], including Arabidopsis [5,26]. In Arabidopsis, the single mutant HSFB2B and double mutant HSFB1 HSFB2 induced the strong expression of stress defense-related genes, and thereby, significantly enhanced the resistance of Arabidopsis against pathogens [27]. These double mutants also improve heat tolerance in Arabidopsis, which might be because HSFB1 and HSFB2b repress the expression of heat shock (HS) inducible genes required in the stress tolerance process [26]. Likewise, it is speculated that TrHSFB2a suppresses the ex-pression of HS-inducible genes in white clover and as well as in overexpression Arabidopsis lines, while when it was mutated, or its expression decreased by RNA interference, it could not suppress the expression of HS-inducible genes, and full functionality of HS-inducible gens in stress condition in the absence of AtHSFB2 protein might have improved the stress tolerance in T-DNA lines. The differences observed in the tolerance level of different stress between overexpression lines and T-DNA lines might be attributed due to the differential expression of HS-inducible genes because TrHSFB2a may repress the expression of such genes due to the existence of the core B3 repression domain (BRD) domain (Figure 1c). This domain is conserved in plant species and speculated as a core domain that repress the activation of downstream HSPs genes in different plant species [41], as observed in the case of HSFB1 in Arabidopsis, but it needs to be tested.
In conclusion, the current study demonstrated that TrHSFB2a is expressed under various abiotic stresses, such as drought, heat and salt stress, and it negatively regulates drought, heat and salt stress, possibly by repressing the activation of HS-inducible genes under stressed conditions. Moreover, diminishing or reducing the expression of TrHSFB2a may have promising utility for enhancing the drought and heat tolerance of white clover.

Plant Growth Conditions
A total of 0.5 g of healthy and uniform-sized white clover (Trifolium repens L. cv. Landino) seeds were surface-sterilized with NaClO (0.5% (w/v), followed by being rinsed in autoclaved ddH 2 O 6 times and germinated in a plastic pot (24 cm length, 18 cm width and 9 cm deep) filled with sterilized moisturized quartz sand in a controlled growth chamber for seven days to maintain germination. Seven-day-old seedlings were irrigated with Hoagland's nutrient solution for an additional three weeks (21 days) until the second leaves fully expanded, changing the Hoagland's nutrient solution after every 2 days. The growth chamber was maintained at a 12 h photoperiod, with day/night temperatures of 23/19 • C, relative humidity of 75% and 250-µmol m −2 ·s −1 photosynthetic photon flux density. Fourweek-old white clover plants were subjected to drought stress by 15% (w/v) polyethylene glycol (PEG) 6000, which was dissolved in Hoagland's solution and samples were collected after 0 h, 1.5 h, 3 h, 6 h, 12 h and 24 h to isolate the total RNA for the subsequent cDNA library constructions. The total RNA was isolated from white clover plants (root + shoot) using an RNAprep Pure Plant Kit (TIANGEN Biotech (Beijing) Co., Ltd., Beijing, China) and was then reverse-transcribed to cDNA by an iScript TM cDNA Synthesis Kit (Bio-Rad Laboratories, (Shanghai) Co., Ltd., Shanghai, China) according to the manufacturer's manual instructions. Arabidopsis thaliana ecotype Col-0 was used for genetic transformation, and thus, grown as wild type. The Arabidopsis genotypes grown in the current stud consisted of wild-type (Col-0) and transgenic Arabidopsis, including overexpressing the TrHSFB2 and RNA interference T-DNA lines of AtHSFB2a for the comparative study. The seeds of Arabidopsis were surface sterilized with alcohol (75% w/v) and (0.1% w/v) NaClO, and subsequently, propagated on 1 /2 MS medium containing 3% sucrose and 0.7% agar, followed by being placed in the dark at 4 • C for 2-3 days for vernalization. Then, for germination, the seeds were transferred to the growth chamber in the set condition of 21 • C temperature with 65% relative humidity, 16 h photoperiod and 150-µmol m −2 ·s −1 photoactive radiation. After 5 days, the equal-sized seedlings were transferred to square plates containing 1/2MS medium with sucrose (3%) and agar (0.7%) and arranged at 90 degrees in the same growth chamber for observing roots and shoot growth differences under normal and stressed conditions between transgenic and wild type Arabidopsis plants. Moreover, after 10 days, well-grown uniform-sized wild and transgenic Arabidopsis plants were transferred into plastic pots (10 cm length, 10 cm width and 7 cm deep) filled with nutrient soil consisting of peat moss, vermiculite and perlite (3:1:1) and placed in the growth chamber for completing rest growth periods or being subjected to different stress treatments after about 25 d of transplanting.

TrHSFB2 Gene Isolation from White Clover and Bioinformatics Analysis
Initially, the homologous sequences of TrHSFB2a were obtained through the BLAST program of NCBI "http://www.ncbi.nlm.nih.gov/Blast (accessed on 12 April 2020)" [51] by searching and identifying conserved sequences region of TrHSFB2a with closely related species and then aligning these sequences with partially published white clover genomes on the NCBI database https://www.ncbi.nlm.nih.gov/genome/13404 (accessed on 12 April 2020) [52] to obtain gene sequence information for the designing primers. The premier primer-5 software was used to design primers for TrHSFB2a amplification, and thereby, obtained a full-length CDS sequence of TrHSFB2a from the cDNA library by using the primer pair TrHSFB2-F CTCGCGAACCTTCTAGAACTCTCA and TrHSFB2-R TCCCTAATCCATCTAACATCAGGTGTCA through the touch-down PCR method and Phanta Max Master Mix (see the procedure detail in Supplementary Methods and Supplementary Table S1). The full-length amplified PCR fragments were excised and cloned into the pMD19 simple vector for sequencing; after confirming the TrHSFB2 sequence from three sequenced clones through a homology search on the NCBI database, an open reading frame of 924bp was cloned into the pBI21 overexpression vector (detail below). The molecular weight and pI values of primary TrHSFB2 protein sequences were determined using the ExPASy software http://web.expasy.org/protparam (accessed on 1 March 2022) [53]. Clustal Omega "https://www.ebi.ac.uk/Tools/msa/clustalo/ (accessed on 1 March 2022)" [54] and PSORT "https://www.genscript.com/psort.html (accessed on 10 December 2020)" [55] were used for multiple sequence alignment and subcellular protein localization, respectively. The motif and conserved domains were identified using "https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 1 March 2022)" [56], and open reading frames (ORFs) were predicted using the NCBI program https://www.ncbi.nlm.nih.gov/orffinder (accessed on 10 May 2020) [57]. The construction of the phylogenetic tree was carried out through the MEGA 11 bioinformatic tool.

Plasmid Construction, Genetic Transformation and Developing Homozygous Overexpression Lines of TrHSFB2a
pSUPER1300-GFP and pBI121 vectors were used for subcellular localization and for constructing the overexpression of Arabidopsis, respectively. The complete coding region (without stop codon) of TrHSFB2a was amplified from a previously constructed pMD19 simple vector by the PCR reaction using Prime STAR ® Max DNA-Polymerase (Takara-Biology Technology (Beijing) Co., Beijing, China) by following manual instructions, and inserted into the restriction sites of XbaI and BamHI of the linearized pBI121 (containing CaMV 35S promoter) and XbaI and KpnI site of the linearized pSuper1300-GFP using the EasyGeno Assembly Cloning kit (TIANGEN Biotech (Beijing) Co., Ltd., Beijing, China) to generate TrHSFB2::GFP in-frame fused protein according to manual's instruction. The primers for amplifying the coding region and flanking vector sequences for both constructs are enlisted in Supplementary Table S1. The constructed pBI121 vector (CaMV 35S::TrHSFB2a) was then mobilized into the floral tissues of the Arabidopsis thaliana genotype "Col-0" by Agrobacterium-mediated transformation through the floral dip method [58]. The seeds of transformed plants were harvested, surface sterilized and sown on the selective media of 1 /2 MS media with Kanamycine-50. Surviving plants were transferred into the soil at the four-leaf stage, confirmed by PCR for the presence of TrHSFB2a, and thereby developed ten independent transgenic overexpression lines. The expression level of TrHSFB2a in transgenic lines was identified by qRT-PCR and selected overexpression lines OE3 and OE4 and OE6 for further experiments based on relative expression levels. Two T-DNA lines (N675207 and N677095) of the Arabidopsis gene AtHSFB2A (Gene ID: At5g62020) were purchased from Arabidopsis mutants stock share platform AraShare "https://www.arashare.cn (accessed on 18 December 2020)" for the comparative functional analysis of TrHAFB2a in Arabidopsis.

Subcellular Localization
The fusion constructs TrHSFB2::GFP and control GFP vector were inserted into the Agrobacterium strain EHA105 competent cells using a standard protocol. Briefly, the transformed Agrobacterium culture was shaken for 2-3 h in YEB (Yeast Extract Beef) solution without an antibiotic at 28 • C, then centrifuged at 4000 rpm for one min (minute) and plated 100 µL on YEB medium containing rifamycin-20 and kanamycin-50, and incubated at 28 • C for 72 h followed by positive colonies selection using colony PCR. Cell culture was shaken at 28 • C until the OD-600 value reached 0.8, then the 30 mL cell culture was centrifuged at 4000 rpm for 10 min at 4 • C, dissolved in autoclaved water containing 5% sucrose and 0.02% silwet and, finally, introduced into tobacco leaves. Then the subcellular location of TrHSFB2a protein was determined by transiently expressing the TrHSFB2a::GFP fusion protein in the leaves of Nicotiana benthamiana. Tobacco plants were placed in a growth chamber in the dark, set at a day/night temperature of 23/19 • C, with 75% relative humidity for 24 h, followed by observing the GFP localization using a fluorescent microscope (Olympus).

Relative Expression of TrHSFB2a
To determine the relative expression level of TrHSFB2a in white clover under drought, heat and salt stress conditions, the four-week-old white clover plants were subjected to drought or salt stress by irrigating with Hoagland's solution containing PEG6000 15% w/v or 200 mM L −1 , respectively, in a controlled growth chamber. High-temperature stress was induced by placing the four-week-old white clover plants in a growth chamber adjusted with a continuous temperature of 35 • C. The processing times were 0 h, 1.5 h, 3 h, 6 h, 12 h and 24 h. The plant samples were collected, immediately frozen in liquid nitrogen, and stored at −80 • C for the subsequent process of total RNA extraction. Plant tissues (leaf, stem and roots) were ground in liquid nitrogen and 100 mg of tissue powder was used to extract the total RNA using HiPure Universal RNA Mini Kit (Magen, Beijing) by following the manual's instructions. After confirming the quantity and quality on the NanDrop spectrophotometer and the RNA integrity on 1.5 agarose gel, cDNA was synthesized by reverse transcription of the first strand of RNA using MonScript™ RTIII All-in-One Mix with dsDNase (Monad Biotech Co., Ltd., Beijing, China) following kit protocols. Later, the relative gene expression was determined using 2X M5 HiPer SYBR Premix EsTaq (with Tli RNaseH) (Mei5 Biotechnology, Co., Ltd., Beijing, China) and the CFX96 Real-Time PCR detection system (Bio-RAD) according to manual instructions. The qRT-PCR primers are provided in Supplementary Table S1. The TrB-Actin gene was used as an internal control for white clover. The amplification conditions were as follows: initial denaturation at 95 • C for 5 min, then, 40 cycles at 95 • C for 10 sec and 60 • C for 30 s. The relative difference in the expression of the genes was calculated by the 2 −∆∆Ct method [59].
For determining TrHSFB2a expressions in transgenic Arabidopsis, the seeds of 10 overexpression lines were sown on 1 /2 MS media supplemented with 3% sucrose and 0.7% agar, followed by vernalization for 2 days. After germination, seedlings were grown for an additional 10 days in the same 1 /2 MS media plates placed in a growth chamber set at 21 • C with 65% relative humidity, 16 h photoperiod, and 150 µmol m −2 ·s −1 photoactive radiation. Then, the extraction of the total RNA, reverse transcription and qRT-PCR analysis was the same as stated above for white clover.

Drought, Salt and Heat Stress in Arabidopsis
The drought, salt and heat stress treatments were according to Jia et al. (2021) [60]. Briefly, the seeds of homozygous T3 lines of transgenic Arabidopsis (overexpression and T-DNA lines) were used for studying stress responses. For drought and salt stress treatments in 1 /2 MS medium, 5 d old wild and transgenic lines were transferred into 1 /2 MS medium (15 plants of each genotype with three technical repeats) supplemented with 200 mML −1 and 300 mML −1 of mannitol or 100 mML −1 and 150 mML −1 of sodium chloride (NaCl). The phenotypes of plants under stress were observed daily for up to 26 days. Photographs of the plants were taken using a digital camera. The yellowing rate (senescence) was calculated by counting the respective plant at the end of the experiment from each treatment. Onemonth-old wild and transgenic lines were also subjected to drought, heat and salt treatment in nutrient soil (peat moss, vermiculite, perlite, 3:1:1) in controlled growth chamber conditions, as stated above. Drought stress treatment was induced by withholding the irrigation for 22 days. Salt treatment was applied by irrigation with 100 mML −1 , 200 mML −1 and 300 mML −1 of NaCl every two times after every two consecutive days. For inducing heat stress, the one-month-old soil-grown WT and transgenic plants were transferred to another growth chamber adjusted with a continuous temperature of 35 • C for two weeks.

Determination of Relative Water Contents (RWC)
RWC was determined by following Barrs and Weatherley [61]. A total of 0.3 g of leaf blade samples were collected between 9:00-10:00 am, wrapped well in ordinary absorbent paper and placed into centrifuge tubes of 50 mL, followed by filling the tubes with water, covered with a lid and placing in a protected place for 24 h. When the leaves had absorbed water at full saturation level, we took out the leaf samples, wiped off the surface water, weighed the saturated fresh weight, and placed the samples in a blast oven set at 105 • C temperature for 45 min. Finally, we dried the samples at 75 • C to reach a constant weight and measured dry weights. Three biological repeats were used for all treatments. RWC of the leaves were calculated by employing the following formula: RWC % = [(FW − DW)/(TW − DW)] × 100 FW, TW and DW represent fresh weight, saturated fresh weight after drenching the leaves in water for 12 h and dry weight, respectively.

Determination of Relative Electrical Conductivity (EC)
Relative electrical conductivity was determined by following Blum and Ebercon [62]. Briefly, 0.1 g of leaf blade samples of the different materials in the study were placed into test tubes filled with 15 mL of deionized water. The initial conductivity (C initial ) was determined by the conductivity meter (YSI Model 32) after 24 h, then we boiled the samples at 100 centigrade for 15 min, followed by cooling down at room temperature and measuring their final conductivity (C max ). Three biological repeats were used for all treatments. The following formula was used to calculate EL: Relative electrical conductivity EL = C initial /C max × 100 where C initial is initial conductivity and C max is final conductivity.

Determination of Chlorophyll Fluorescence Parameters
Chlorophyll fluorescence parameters were measured according to Jia et al. (2021) [60]. Briefly, for determining chlorophyll fluorescence parameters, we placed the soil-grown WT and transgenic lines in the dark for 30 min, and measured the maximum quantum yield of PSII photochemistry (Fv/Fm) and PI ABS by a pulse-amplitude modulation portable chlorophyll fluorometer (PAM-2500), using 15 replicates for each treatment.

Determination of Malondialdehyde (MDA) Content
The MDA contents of the leaves of different materials in the study were determined by following Heath and Packer (1968) method. Briefly, 0.1 g leaf tissues were immediately frozen in liquid nitrogen, followed by thoroughly grounded in ice within 2 mL of 50 mM pre-cooled phosphate-buffered saline (PBS) solution (pH 7.8) and centrifuged at 12,000× g for up to 30 min at 4 • C. The supernatant was taken and used to extract and determine MDA content by adding 1 mL of reaction solution comprising of trichloroacetic acid (20% w/v) and thiobarbituric acid (0.5% w/v) of a 0.5 mL solution of crude enzyme. Then samples were heated at 95 degrees Celsius for up to 30 min in a water bath, followed by quickly cooling to room temperature in an ice bath, with continuous gentle shaking to avoid bubble formations, and then centrifuged at 10,000× g for 10 min, removing the air bubbles, if any, with the pipette. Finally, we took out the supernatant and measured absorbance values at the wavelength of 532 nm and 600 nm. Three biological repeats were used for all treatments. To calculate MDA contents values, we subtracted the absorbance value of 600 nm wavelength from the absorbance value obtained at 532 nm wavelength by using an extinction coefficient of 155 mM −1 em −1 [63].
MDA (nmole g −1 DW) = (A532 − A600) × V × 1000/155 × W where A532 is the absorbance at 532 nm, A600 is the absorbance at 600 nm, V is the extraction volume and W is the dry weight of the leaf.

Statistical Analysis
Statistical analyses were performed using IBM SPSS and GraphPad Prism 8.3.0 software, and graphs were made using GraphPad Prism 8.3.0 software.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare that they have no competing interest.