GmWRKY45 Enhances Tolerance to Phosphate Starvation and Salt Stress, and Changes Fertility in Transgenic Arabidopsis

WRKY Transcription Factor Responses and Tolerance to Abiotic Stresses in Plants

Plants are subjected to abiotic stresses throughout their developmental period. Abiotic stresses include drought, salt, heat, cold, heavy metals, nutritional elements, and oxidative stresses. Improving plant responses to various environmental stresses is critical for plant survival and perpetuation. WRKY transcription factors have special structures (WRKY structural domains), which enable the WRKY transcription factors to have different transcriptional regulatory functions. WRKY transcription factors can not only regulate abiotic stress responses and plant growth and development by regulating phytohormone signalling pathways but also promote or suppress the expression of downstream genes by binding to the W-box [TGACCA/TGACCT] in the promoters of their target genes. In addition, WRKY transcription factors not only interact with other families of transcription factors to regulate plant defence responses to abiotic stresses but also self-regulate by recognising and binding to W-boxes in their own target genes to regulate their defence responses to abiotic stresses. However, in recent years, research reviews on the regulatory roles of WRKY transcription factors in higher plants have been scarce and shallow. In this review, we focus on the structure and classification of WRKY transcription factors, as well as the identification of their downstream target genes and molecular mechanisms involved in the response to abiotic stresses, which can improve the tolerance ability of plants under abiotic stress, and we also look forward to their future research directions, with a view of providing theoretical support for the genetic improvement of crop abiotic stress tolerance.

Recent advances in unraveling the mystery of combined nutrient stress in plants

Efficiently regulating growth to adapt to varying resource availability is crucial for organisms, including plants. In particular, the acquisition of essential nutrients is vital for plant development, as a shortage of just one nutrient can significantly decrease crop yield. However, plants constantly experience fluctuations in the presence of multiple essential mineral nutrients, leading to combined nutrient stress conditions. Unfortunately, our understanding of how plants perceive and respond to these multiple stresses remains limited. Unlocking this mystery could provide valuable insights and help enhance plant nutrition strategies. This review focuses specifically on the regulation of phosphorous homeostasis in plants, with a primary emphasis on recent studies that have shed light on the intricate interactions between phosphorous and other essential elements, such as nitrogen, iron, and zinc, as well as non-essential elements like aluminum and sodium. By summarizing and consolidating these findings, this review aims to contribute to a better understanding of how plants respond to and cope with combined nutrient stress.

Unregulated GmAGL82 due to Phosphorus Deficiency Positively Regulates Root Nodule Growth in Soybean

Nitrogen fixation, occurring through the symbiotic relationship between legumes and rhizobia in root nodules, is crucial in sustainable agriculture. Nodulation and soybean production are influenced by low levels of phosphorus stress. In this study, we discovered a MADS transcription factor, GmAGL82, which is preferentially expressed in nodules and displays significantly increased expression under conditions of phosphate (Pi) deficiency. The overexpression of GmAGL82 in composite transgenic plants resulted in an increased number of nodules, higher fresh weight, and enhanced soluble Pi concentration, which subsequently increased the nitrogen content, phosphorus content, and overall growth of soybean plants. Additionally, transcriptome analysis revealed that the overexpression of GmAGL82 significantly upregulated the expression of genes associated with nodule growth, such as GmENOD100, GmHSP17.1, GmHSP17.9, GmSPX5, and GmPIN9d. Based on these findings, we concluded that GmAGL82 likely participates in the phosphorus signaling pathway and positively regulates nodulation in soybeans. The findings of this research may lay the theoretical groundwork for further studies and candidate gene resources for the genetic improvement of nutrient-efficient soybean varieties in acidic soils.

Functional Characterisation of the Transcription Factor GsWRKY23 Gene from Glycine soja in Overexpressed Soybean Composite Plants and Arabidopsis under Salt Stress

WRKY proteins are a superfamily of transcription factors (TFs) that play multiple roles in plants' growth, development, and environmental stress response. In this study, a novel WRKY gene called GsWRKY23 that is specifically upregulated in salt-tolerant Glycine soja accession BB52 seedlings was identified by transcriptomic analysis under salt stress. How the physiological functions and mechanisms of the GsWRKY23 gene affect salt tolerance was investigated using transformations of soybean hairy roots and Arabidopsis, including wild-type (WT) and atwrky23-mutant plants. The results showed that GsWRKY23 in the roots, stems, and leaves of BB52, along with its promoter in the cotyledons and root tips of GsWRKY23pro::GUS Arabidopsis seedlings, displayed enhanced induction under salt stress. GsWRKY23 localises to the nucleus and shows transcriptional activation ability in yeast cells. Compared to GsWRKY23-RNAi wild soybean hairy-root composite plants under salt stress, obvious improvements, such as superior growth appearance, plant height and fresh weight (FW), and leaf chlorophyll and relative water content (RWC), were displayed by GsWRKY23-overexpressing (OE) composite plants. Moreover, their relative electrolytic leakage (REL) values and malondialdehyde (MDA) contents in the roots and leaves declined significantly. Most of the contents of Na+ and Cl- in the roots, stems, and leaves of GsWRKY23-OE plants decreased significantly, while the content of K+ in the roots increased, and the content of NO3- displayed no obvious change. Ultimately, the Na+/K+ ratios of roots, stems, and leaves, along with the Cl-/NO3- ratios of roots and stems, decreased significantly. In the transgenic WT-GsWRKY23 and atwrky23-GsWRKY23 Arabidopsis seedlings, the salt-induced reduction in seed germination rate and seedling growth was markedly ameliorated; plant FW, leaf chlorophyll content, and RWC increased, and the REL value and MDA content in shoots decreased significantly. In addition, the accumulation of Na+ and Cl- decreased, and the K+ and NO3- levels increased markedly to maintain lower Na+/K+ and Cl-/NO3- ratios in the roots and shoots. Taken together, these results highlight the role of GsWRKY23 in regulating ionic homeostasis in NaCl-stressed overexpressed soybean composite plants and Arabidopsis seedlings to maintain lower Na+/K+ and Cl-/NO3- ratios in the roots and shoots, thus conferring improved salt tolerance.

Laser capture microdissection-based spatiotemporal transcriptomes uncover regulatory networks during seed abortion in seedless Ponkan (Citrus reticulata)

Seed abortion is an important process in the formation of seedless characteristics in citrus fruits. However, the molecular regulatory mechanism underlying citrus seed abortion is poorly understood. Laser capture microdissection-based RNA-seq combined with Pacbio-seq was used to profile seed development in the Ponkan cultivars 'Huagan No. 4' (seedless Ponkan) (Citrus reticulata) and 'E'gan No. 1' (seeded Ponkan) (C. reticulata) in two types of seed tissue across three developmental stages. Through comparative transcriptome and dynamic phytohormone analyses, plant hormone signal, cell division and nutrient metabolism-related processes were revealed to play critical roles in the seed abortion of 'Huagan No. 4'. Moreover, several genes may play indispensable roles in seed abortion of 'Huagan No. 4', such as CrWRKY74, CrWRKY48 and CrMYB3R4. Overexpression of CrWRKY74 in Arabidopsis resulted in severe seed abortion. By analyzing the downstream regulatory network, we further determined that CrWRKY74 participated in seed abortion regulation by inducing abnormal programmed cell death. Of particular importance is that a preliminary model was proposed to depict the regulatory networks underlying seed abortion in citrus. The results of this study provide novel insights into the molecular mechanism across citrus seed development, and reveal the master role of CrWRKY74 in seed abortion of 'Huagan No. 4'.

Twenty years of mining salt tolerance genes in soybean

Current combined challenges of rising food demand, climate change and farmland degradation exert enormous pressure on agricultural production. Worldwide soil salinization, in particular, necessitates the development of salt-tolerant crops. Soybean, being a globally important produce, has its genetic resources increasingly examined to facilitate crop improvement based on functional genomics. In response to the multifaceted physiological challenge that salt stress imposes, soybean has evolved an array of defences against salinity. These include maintaining cell homeostasis by ion transportation, osmoregulation, and restoring oxidative balance. Other adaptations include cell wall alterations, transcriptomic reprogramming, and efficient signal transduction for detecting and responding to salt stress. Here, we reviewed functionally verified genes that underly different salt tolerance mechanisms employed by soybean in the past two decades, and discussed the strategy in selecting salt tolerance genes for crop improvement. Future studies could adopt an integrated multi-omic approach in characterizing soybean salt tolerance adaptations and put our existing knowledge into practice via omic-assisted breeding and gene editing. This review serves as a guide and inspiration for crop developers in enhancing soybean tolerance against abiotic stresses, thereby fulfilling the role of science in solving real-life problems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01383-3.

Improving phosphorus acquisition efficiency through modification of root growth responses to phosphate starvation in legumes

Phosphorus (P) is one of the essential macronutrients for plant growth and development, and it is an integral part of the major organic components, including nucleic acids, proteins and phospholipids. Although total P is abundant in most soils, a large amount of P is not easily absorbed by plants. Inorganic phosphate (Pi) is the plant-available P, which is generally immobile and of low availability in soils. Hence, Pi starvation is a major constraint limiting plant growth and productivity. Enhancing plant P efficiency can be achieved by improving P acquisition efficiency (PAE) through modification of morpho-physiological and biochemical alteration in root traits that enable greater acquisition of external Pi from soils. Major advances have been made to dissect the mechanisms underlying plant adaptation to P deficiency, especially for legumes, which are considered important dietary sources for humans and livestock. This review aims to describe how legume root growth responds to Pi starvation, such as changes in the growth of primary root, lateral roots, root hairs and cluster roots. In particular, it summarizes the various strategies of legumes to confront P deficiency by regulating root traits that contribute towards improving PAE. Within these complex responses, a large number of Pi starvation-induced (PSI) genes and regulators involved in the developmental and biochemical alteration of root traits are highlighted. The involvement of key functional genes and regulators in remodeling root traits provides new opportunities for developing legume varieties with maximum PAE needed for regenerative agriculture.

Genome-wide annotation and expression analysis of WRKY and bHLH transcriptional factor families reveal their involvement under cadmium stress in tomato (Solanum lycopersicum L.)

The WRKY and bHLH transcription factors have been implicated in the regulation of gene expression during various physiological processes in plants, especially in plant stress responses. However, little information about the heavy metal-responsive SlWRKY and SlbHLH in tomato (Solanum lycopersicum) is available. We performed a genome-wide investigation for these two TF families in S. lycopersicum and determined their role in cadmium (Cd) stress tolerance. Furthermore, ortholog analysis with the Arabidopsis genome led to classifying WRKY and bHLH ortholog genes into nine and 11 clusters, respectively. The comparative phylogenetic analysis revealed duplication events and gene loss in Arabidopsis and S. lycopersicum, which occurred during evolution both before and after the last common ancestor of the two species. Orthologous relationships are also supported by additional evidence, such as gene structure, conserved motif compositions, and protein-protein interaction networks for the majority of genes, suggesting their similar functions. A comprehensive transcriptomics analysis revealed that both WRKY and bHLH genes were differentially expressed in response to cadmium stress as compared with control plants. A gene ontology analysis revealed that most WRKYs and bHLHs are DNA-binding essential proteins that regulate gene expression positively and negatively. Analyses of interaction networks revealed that both WRKYs and bHLHs mediate networks implicated in several stress-signaling pathways. The findings of this work may help us to comprehend the intricate transcriptional control of WRKY and bHLH genes and identify potential stress-responsive genes relevant to tomato genetic improvement. Moreover, identifying heavy metal stress-responsive WRKY and bHLH genes in S. lycopersicum will provide fundamental insights for developing new heavy metal stress-tolerant varieties of tomato crops.

GmWRKY21, a Soybean WRKY Transcription Factor Gene, Enhances the Tolerance to Aluminum Stress in Arabidopsis thaliana

The WRKY transcription factors (TFs) are one of the largest families of TFs in plants and play multiple roles in plant growth and development and stress response. In this study, GmWRKY21 encoding a WRKY transcription factor was functionally characterized in Arabidopsis and soybean. The GmWRKY21 protein containing a highly conserved WRKY domain and a C₂H₂ zinc-finger structure is located in the nucleus and has the characteristics of transcriptional activation ability. The GmWRKY21 gene presented a constitutive expression pattern rich in the roots, leaves, and flowers of soybean with over 6-fold of relative expression levels and could be substantially induced by aluminum stress. As compared to the control, overexpression of GmWRKY21 in Arabidopsis increased the root growth of seedlings in transgenic lines under the AlCl₃ concentrations of 25, 50, and 100 μM with higher proline and lower MDA accumulation. The results of quantitative real-time polymerase chain reaction (qRT-PCR) showed that the marker genes relative to aluminum stress including ALMT, ALS3, MATE, and STOP1 were induced in GmWRKY21 transgenic plants under AlCl₃ treatment. The stress-related genes, such as KIN1, COR15A, COR15B, COR47, GLOS3, and RD29A, were also upregulated in GmWRKY21 transgenic Arabidopsis under aluminum stress. Similarly, stress-related genes, such as GmCOR47, GmDREB2A, GmMYB84, GmKIN1, GmGST1, and GmLEA, were upregulated in hair roots of GmWRKY21 transgenic plants. In summary, these results suggested that the GmWRKY21 transcription factor may promote the tolerance to aluminum stress mediated by the pathways regulating the expression of the acidic aluminum stress-responsive genes and abiotic stress-responsive genes.

Maize Transcription Factor ZmARF4 Confers Phosphorus Tolerance by Promoting Root Morphological Development

Plant growth and development are closely related to phosphate (Pi) and auxin. However, data regarding auxin response factors (ARFs) and their response to phosphate in maize are limited. Here, we isolated ZmARF4 in maize and dissected its biological function response to Pi stress. Overexpression of ZmARF4 in Arabidopsis confers tolerance of Pi deficiency with better root morphology than wild-type. Overexpressed ZmARF4 can partially restore the absence of lateral roots in mutant arf7 arf19. The ZmARF4 overexpression promoted Pi remobilization and up-regulated AtRNS1, under Pi limitation while it down-regulated the expression of the anthocyanin biosynthesis genes AtDFR and AtANS. A continuous detection revealed higher activity of promoter in the Pi-tolerant maize P178 line than in the sensitive 9782 line under low-Pi conditions. Meanwhile, GUS activity was specifically detected in new leaves and the stele of roots in transgenic offspring. ZmARF4 was localized to the nucleus and cytoplasm of the mesophyll protoplast and interacted with ZmILL4 and ZmChc5, which mediate lateral root initiation and defense response, respectively. ZmARF4 overexpression also conferred salinity and osmotic stress tolerance in Arabidopsis. Overall, our findings suggest that ZmARF4, a pleiotropic gene, modulates multiple stress signaling pathways, and thus, could be a candidate gene for engineering plants with multiple stress adaptation.

GmWRKY46, a WRKY transcription factor, negatively regulates phosphorus tolerance primarily through modifying root morphology in soybean

Phosphorus (P) deficiency affects soybean growth and development, resulting in significant reduction of yields. However, the regulatory mechanism of P deficiency tolerance in soybean is still largely unclear. WRKY transcription factors are a family of regulators involved in a variety of abiotic stresses in plants while rarely reported in P deficiency. Here, we demonstrated that a soybean GmWRKY46 gene, belonging to group III of WRKY TF family, was involved in the regulation of P deficiency tolerance in soybean. The expression of GmWRKY46 in low P sensitive soybean varieties was significantly higher than that in tolerant soybean varieties. It was primarily expressed in roots and strongly induced by P deprivation. GmWRKY46 was localized in the nucleus. Compared with the control expressing the empty vector, overexpression of GmWRKY46 in soybean hairy roots exhibited more sensitive phenotypes to low P stress, while the RNA interfered GmWRKY46 significantly enhanced P deficiency tolerance by increasing the proliferation, elongation and P absorption efficiency of hairy roots. Expression patterns of a number of P-responsive genes (GmPht1;1, GmPht1;4, GmPTF1, GmACP1, GmPAP21 and GmExpansin-A7) were altered in both overexpression and gene silenced plants. The results provided a novel insight into how soybean responds to low P stress and new gene that may be used to improve soybean low P tolerance through gene editing approach.

Transcription Factor GmWRKY46 Enhanced Phosphate Starvation Tolerance and Root Development in Transgenic Plants

Phosphorus (P) is one of the essential macronutrients, whose deficiency limits the growth and development of plants. In this study, we investigated the possible role of GmWRKY46 in the phosphate (Pi) starvation stress tolerance of soybean. GmWRKY46 belonged to the group III subfamily of the WRKY transcription factor family, which was localized in the nucleus and had transcriptional activator activity. GmWRKY46 could be strongly induced by Pi starvation, especially in soybean roots. Overexpression of GmWRKY46 significantly enhanced tolerance to Pi starvation and lateral root development in transgenic Arabidopsis. RNA-seq analysis showed that overexpression of GmWRKY46 led to change in many genes related to energy metabolisms, stress responses, and plant hormone signal transduction in transgenic Arabidopsis. Among these differential expression genes, we found that overexpression of AtAED1 alone could enhance the tolerance of transgenic Arabidopsis to Pi starvation. Y1H and ChIP-qPCR analyses showed that GmWRKY46 could directly bind to the W-box motif of the AtAED1 promoter in vitro and in vivo. Furthermore, results from intact soybean composite plants with GmWRKY46 overexpression showed that GmWRKY46 was involved in hairy roots development and subsequently affected plant growth and Pi uptake. These results provide a basis for the molecular genetic breeding of soybean tolerant to Pi starvation.

Gm6PGDH1, a Cytosolic 6-Phosphogluconate Dehydrogenase, Enhanced Tolerance to Phosphate Starvation by Improving Root System Development and Modifying the Antioxidant System in Soybean

Phosphorus plays an important role in plant growth and development, and is an important limiting factor for crop yield. Although previous studies have shown that 6-phosphogluconate dehydrogenase (6PGDH) plays an important role in plant resistance to adversity, its response to low phosphorus (P) stress remains unknown. In this study, we reported the cloning and characterization of a cytosolic 6PGDH gene, Gm6PGDH1, which enhanced the tolerance to phosphate (Pi) starvation by improving root system development and modifying the antioxidant system in transgenic plants. Gm6PGDH1 was highly expressed in the root at full bloom stage, and strongly induced by Pi starvation. The results from intact soybean composite plant and soybean plant, both containing a Gm6PGDH1-overexpressing construct, showed that Gm6PGDH1 was involved in root system development, and subsequently affected P uptake under Pi-deficient conditions. Meanwhile, the accumulation of reactive oxygen species (ROS) in the root tip of transgenic soybean was reduced, and the activity of ROS-scavenging enzymes was enhanced compared with those of the wild type under Pi-deficient conditions. Interestingly, we found that the overexpression of Gm6PGDH1 weakened the response of several other important Pi-answer genes to Pi starvation, such as some purple acid phosphatases (PAPs) and redox-related genes. In addition, the results from a virus-induced gene silencing (VIGS) indicated that Gm6PGDH1 might have functional redundancy in soybean, and the results from a heterogeneous transformation system showed that overexpressing Gm6PGDH1 also enhanced tolerance to Pi starvation in transgenic Arabidopsis. Together, these results suggested the great potential of Gm6PGDH1 in crop breeding for low Pi tolerance.

Genome-wide identification and comparative analysis of the WRKY gene family in aquatic plants and their response to abiotic stresses in giant duckweed (Spirodela polyrhiza)

WRKY is one of the largest transcription factor families across higher plant species and is involved in important biological processes and plant responses to various biotic/abiotic stresses. However, only a few investigations on WRKYs have been conducted in aquatic plants. This study first systematically analyzed the gene structure, protein properties, and phylogenetic relationship of 693 WRKYs in nine aquatic and two wetland plants at the genome-wide level. The pattern of WRKY groups in two aquatic ferns provided new evidence for the origin and evolution of WRKY genes. ARE cis-regulatory elements show an unusual high frequency in the promoter region of WRKY genes, indicating the adaptation to the aquatic habitat in aquatic plants. The WRKY gene family experienced a series of gene loss events in aquatic plants, especially group III. Further studies were conducted on the interaction network of SpWRKYs, their target genes, and non-coding RNAs. The expression profile of SpWRKYs under phosphate starvation, cold, and submergence conditions revealed that most SpWRKYs are involved in the response to abiotic stresses. Our investigations lay the foundation for further study on the mechanism of WRKYs responding to abiotic stresses in aquatic plants.

Genome-wide identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

Genome-wide identification and characterization of GRAS genes in soybean (Glycine max)

BACKGROUND

GRAS proteins are crucial transcription factors, which are plant-specific and participate in various plant biological processes. Thanks to the rapid progress of the whole genome sequencing technologies, the GRAS gene families in different plants have been broadly explored and studied. However, comprehensive research on the soybean (Glycine max) GRAS gene family is relatively lagging.

RESULTS

In this study, 117 Glycine max GRAS genes (GmGRAS) were identified. Further phylogenetic analyses showed that the GmGRAS genes could be categorized into nine gene subfamilies: DELLA, HAM, LAS, LISCL, PAT1, SCL3, SCL4/7, SCR and SHR. Gene structure analyses turned out that the GmGRAS genes lacked introns and were relatively conserved. Conserved domains and motif patterns of the GmGRAS members in the same subfamily or clade exhibited similarities. Notably, the expansion of the GmGRAS gene family was driven both by gene tandem and segmental duplication events. Whereas, segmental duplications took the major role in generating new GmGRAS genes. Moreover, the synteny and evolutionary constraints analyses of the GRAS proteins among soybean and distinct species (two monocots and four dicots) provided more detailed evidence for GmGRAS gene evolution. Cis-element analyses indicated that the GmGRAS genes may be responsive to diverse environmental stresses and regulate distinct biological processes. Besides, the expression patterns of the GmGRAS genes were varied in various tissues, during saline and dehydration stresses and during seed germination processes.

CONCLUSIONS

We conducted a systematic investigation of the GRAS genes in soybean, which may be valuable in paving the way for future GmGRAS gene studies and soybean breeding.

Comparative Transcriptomics Analysis and Functional Study Reveal Important Role of High-Temperature Stress Response Gene GmHSFA2 During Flower Bud Development of CMS-Based F1 in Soybean

High-temperature (HT) is one of the most important environmental factors that negatively impact the yield of some soybean cytoplasmic male sterility (CMS)-based hybrid (F₁) combinations. The response of soybean to HT, especially at the male organ development stage, is poorly understood. To investigate the molecular mechanisms of the response from soybean CMS-based F₁ male organ to HT, a detailed transcriptomics analysis was performed during flower bud development of soybean HT-tolerant and HT-sensitive CMS-based F₁ combinations (NF₁ and YF₁) under normal-temperature and HT conditions. Obvious HT damage was observed by subjecting YF₁ with HT, such as indehiscent anthers and decreased pollen fertility, whereas the male fertility of NF₁ was normal. In total, 8,784 differentially expressed genes (DEGs) were found to respond to HT stress, which were mainly associated with anther/pollen wall development, carbohydrate metabolism and sugar transport, and auxin signaling. The quantitative real-time PCR (qRT-PCR) analysis and substance content detection also revealed that HT caused male fertility defects in YF₁ by altering pectin metabolism, auxin, and sugar signaling pathways. Most importantly, the sugar signaling-PIF-auxin signaling pathway may underlie the instability of male fertility in YF₁ under HT. Furthermore, HT induced the expression of heat shock factor (HSF) and heat shock protein (HSP) gene families. Overexpression of GmHSFA2 in Arabidopsis can promote the expression of HT protective genes (such as HSP20) by binding to the HSE motifs in their promoters, so as to improve the HT tolerance during flowering. Our results indicated that GmHSFA2 acted as a positive regulator, conferring HT tolerance improvement in soybean CMS-based F₁. GmHSFA2 may be directly involved in the activation of male fertility protection mechanism in the soybean CMS-based F₁ under HT stress.