Phytoglobin Expression Alters the Na+/K+ Balance and Antioxidant Responses in Soybean Plants Exposed to Na2SO4

Soybean (Glycine max) is an economically important crop which is very susceptible to salt stress. Tolerance to Na2SO4 stress was evaluated in soybean plants overexpressing or suppressing the phytoglobin GmPgb1. Salt stress depressed several gas exchange parameters, including the photosynthetic rate, caused leaf damage, and reduced the water content and dry weights. Lower expression of respiratory burst oxidase homologs (RBOHB and D), as well as enhanced antioxidant activity, resulting from GmPgb1 overexpression, limited ROS-induced damage in salt-stressed leaf tissue. The leaves also exhibited higher activities of the H2O2-quenching enzymes, catalase (CAT) and ascorbate peroxidase (APX), as well as enhanced levels of ascorbic acid. Relative to WT and GmPgb1-suppressing plants, overexpression of GmPgb1 attenuated the accumulation of foliar Na+ and exhibited a lower Na+/K+ ratio. These changes were attributed to the induction of the Na+ efflux transporter SALT OVERLY SENSITIVE 1 (SOS1) limiting Na+ intake and transport and the inward rectifying K+ channel POTASSIUM TRANSPORTER 1 (AKT1) required for the maintenance of the Na+/K+ balance.


Introduction
Salinity influences plant growth and development and has negative effects on productivity [1]. Occurring virtually under all climatic conditions as a result of natural events and human practices, salinity is more prominent in arid environments where the discrepancy between the plant demand for water and the water provided by precipitation results in the accumulation of minerals in the soil [2]. It is estimated that about 20% of irrigated land is salt-affected [3], and plant survival is dependent on the salinity tolerance of the species. Soybean (Glycine max), a valuable source of proteins and oils, is especially susceptible to high salt concentrations, and several desirable agronomic traits during both vegetative and reproductive stages are severely affected [4]. Among the different stages of development, seedlings are the most vulnerable to salt stress, as demonstrated by substantial growth retardation which compromises the subsequent phases of development and, ultimately, plant productivity [5].
Plants exposed to salinity experience several types of stress including osmotic, ionic, and oxidative stress. Osmotic stress arises from a reduction of the soil's osmotic potential, which impairs the ability of plants to absorb water [6]. Ionic stress occurs as a result of toxic accumulation of ions, especially of Na + , which contributes to mineral imbalance and and drought [38] through reduction of NO, ROS, and oxidative stress. To further our knowledge on the function of Pgb1 as a stress attenuator, we examined the behavior of previously characterized soybean plants overexpressing or downregulating GmPgb1 [37] to Na 2 SO 4 stress. It is hypothesized that during salt stress, overexpression of GmPgb1 attenuates the depression of several morphophysiological parameters through limitation of ROS-induced oxidative damage and Na + accumulation. To address this hypothesis, we measured biomass, water content, leaf damage, and K + and Na + concentrations, as well as photosynthetic gas exchange parameters and the enzymes related to oxidative stress known to be influenced by salinity.

Overexpression of GmPgb1 Augments Tolerance to Salinity Stress
The transcript levels of GmPgb1 during Na 2 SO 4 stress were measured in the WT line and in previously characterized lines [37] overexpressing (GmPgb1(S)17) or downregulating (GmPgb1(RNAi)23) GmPgb1. Salinity gradually augmented the GmPgb1 transcripts in WT roots and leaves, with a marked increment observed in the latter between 6 and 12 h ( Figure 1). The increase in GmPgb1 in leaves trailed the increase in roots by approximately 6 h. The GmPgb1 overexpression resulted in an approximately six-fold increase in root transcripts and an approximately twelve-fold increase in leaf transcripts, likely reflecting the relatively low base levels of transcripts in leaves compared to roots of non-transformed plant materials [38]. The RNAi line did not completely repress the induction of GmPgb1 in roots whereas it did repress the leaf GmPgb1 expression.
waterlogging [37] and drought [38] through reduction of NO, ROS, and oxidative stre To further our knowledge on the function of Pgb1 as a stress attenuator, we examined t behavior of previously characterized soybean plants overexpressing or downregulati GmPgb1 [37] to Na2SO4 stress. It is hypothesized that during salt stress, overexpression GmPgb1 attenuates the depression of several morphophysiological parameters throu limitation of ROS-induced oxidative damage and Na + accumulation. To address this h pothesis, we measured biomass, water content, leaf damage, and K + and Na + concent tions, as well as photosynthetic gas exchange parameters and the enzymes related to o dative stress known to be influenced by salinity.

Overexpression of GmPgb1 Augments Tolerance to Salinity Stress
The transcript levels of GmPgb1 during Na2SO4 stress were measured in the WT li and in previously characterized lines [37] overexpressing (GmPgb1(S)17) or downreg lating (GmPgb1(RNAi)23) GmPgb1. Salinity gradually augmented the GmPgb1 transcrip in WT roots and leaves, with a marked increment observed in the latter between 6 and h ( Figure 1). The increase in GmPgb1 in leaves trailed the increase in roots by appro mately 6 h. The GmPgb1 overexpression resulted in an approximately six-fold increase root transcripts and an approximately twelve-fold increase in leaf transcripts, likely flecting the relatively low base levels of transcripts in leaves compared to roots of no transformed plant materials [38]. The RNAi line did not completely repress the inducti of GmPgb1 in roots whereas it did repress the leaf GmPgb1 expression. Differences in photosynthetic parameters between the lines were clearly detectab after 24 h of Na2SO4 treatment. In leaves of WT plants, salinity stress depressed the ph tosynthetic rate, the transpiration rate, the stomatal conductance, and the water use e ciency by more than 50% (Figure 2). With the exclusion of transpiration, leaves overe pressing GmPgb1 exhibited a higher retention of gas exchange parameters. T GmPgb1(RNAi)23 line had the lowest retention of the photosynthetic rate, the transpi tion, and the stomatal conductance relative to the WT (Figure 2). Salinity also induced t closure of the stomata (Supplemental Figure S1), as evidenced by the reduced stoma conductance. This effect was mitigated in leaves of the GmPgb1(S)17 line. During stre reduction in the photosynthetic rate is often concomitant to the degradation of thylako Differences in photosynthetic parameters between the lines were clearly detectable after 24 h of Na 2 SO 4 treatment. In leaves of WT plants, salinity stress depressed the photosynthetic rate, the transpiration rate, the stomatal conductance, and the water use efficiency by more than 50% (Figure 2). With the exclusion of transpiration, leaves overexpressing GmPgb1 exhibited a higher retention of gas exchange parameters. The GmPgb1(RNAi)23 line had the lowest retention of the photosynthetic rate, the transpiration, and the stomatal conductance relative to the WT ( Figure 2). Salinity also induced the closure of the stomata (Supplemental Figure S1), as evidenced by the reduced stomatal conductance. This effect was mitigated in leaves of the GmPgb1(S)17 line. During stress, reduction in the photosynthetic rate is often concomitant to the degradation of thylakoidassociated protein complexes [39]. The transcript levels of the selected genes known to be sensitive to abiotic stress and involved in the formation of photosystem I (PsaA, PsaB, and LHCA), cytochrome B6F (cylB6F), and photosystem II (PsbA, PsbB, PsbC, and PsbD) complexes were measured at 24 h. Except for LHCA, the relative transcript levels of all the genes, comparable between genotypes under control conditions, decreased in leaves of the WT plants exposed to Na 2 SO 4 (Supplemental Figure S2). This decrease was mitigated by the overexpression of GmPgb1. Relative to the WT, GmPgb1(RNAi)23 leaves exhibited lower levels of PsaA, PsbB, and PsbC (Supplemental Figure S2). associated protein complexes [39]. The transcript levels of the selected genes known to be sensitive to abiotic stress and involved in the formation of photosystem I (PsaA, PsaB, and LHCA), cytochrome B6F (cylB6F), and photosystem II (PsbA, PsbB, PsbC, and PsbD) complexes were measured at 24 h. Except for LHCA, the relative transcript levels of all the genes, comparable between genotypes under control conditions, decreased in leaves of the WT plants exposed to Na2SO4 (Supplemental Figure S2). This decrease was mitigated by the overexpression of GmPgb1. Relative to the WT, GmPgb1(RNAi)23 leaves exhibited lower levels of PsaA, PsbB, and PsbC (Supplemental Figure S2).  Morphological differences between the lines were visible after seven days of stress. Overexpression of GmPgb1 limited the Na 2 SO 4 depression of the shoot water content (WC) and the root and shoot dry weight, which were more severe in the WT and GmPgb1(RNAi)23 lines ( Figure 3). Leaf injury, not noticeable in leaves of the GmPgb1(S)17 line exposed to salinity for seven days, was apparent in the other two lines, and especially in leaves of the GmPgb1(RNAi)23 plants (Figure 3). nt. J. Mol. Sci. 2022, 23, x FOR PEER REVIEW Morphological differences between the lines were visible after seven Overexpression of GmPgb1 limited the Na2SO4 depression of the shoot water and the root and shoot dry weight, which were more severe in GmPgb1(RNAi)23 lines (Figure 3). Leaf injury, not noticeable in l GmPgb1(S)17 line exposed to salinity for seven days, was apparent in the o and especially in leaves of the GmPgb1(RNAi)23 plants (Figure 3).

Na + and K + Accumulation during Salinity Stress Is Influenced by GmPgb1
Depression of the photosynthetic rate and leaf injury during salinity s related to the toxicity of the Na + [13,14] accumulating in leaves and disturb ostasis of K + needed in several photosynthetic functions [16]. Salt stress au and leaf Na + in all the genotypes, but the increment in leaves was reduced expressing GmPgb1 ( Figure 4A). The level of K + was generally reduced by N organs, and the reduction was abolished when GmPgb1 was overexpresse Because of these changes, the foliar Na + /K + ratio after 24 h of salt stress wa leaves overexpressing GmPgb1 relative to those of WT and GmPgb1(RNAi) ure 4C) To further understand the mechanisms for the differential accumulati

Na + and K + Accumulation during Salinity Stress Is Influenced by GmPgb1
Depression of the photosynthetic rate and leaf injury during salinity stress are often related to the toxicity of the Na + [13,14] accumulating in leaves and disturbing the homeostasis of K + needed in several photosynthetic functions [16]. Salt stress augmented root and leaf Na + in all the genotypes, but the increment in leaves was reduced in plants overexpressing GmPgb1 ( Figure 4A). The level of K + was generally reduced by Na 2 SO 4 in both organs, and the reduction was abolished when GmPgb1 was overexpressed ( Figure 4B). Because of these changes, the foliar Na + /K + ratio after 24 h of salt stress was lower in the leaves overexpressing GmPgb1 relative to those of WT and GmPgb1(RNAi)23 plants ( Figure 4C) significant increase also occurred in WT leaves. Salinity also elevated the expression of AKT1 in the roots of GmPgb1(S)17 plants.  To further understand the mechanisms for the differential accumulation of Na + and K + , we measured the expression of two key genes involved in Na + and K + transport: the Na + efflux transporter SALT-OVERLY SENSITIVE 1 (SOS1) limiting the Na + intake in roots and its long-distance movement [16] and the inward-rectifying K + channel AKT1 required for the maintenance of the Na + /K + ratio [40]. The expression of SOS1, which remained unchanged in GmPgb1(RNAi)23 roots and leaves following imposition of Na 2 SO 4 stress, increased in both organs of the plants overexpressing GmPgb1 ( Figure 4D). A small significant increase also occurred in WT leaves. Salinity also elevated the expression of AKT1 in the roots of GmPgb1(S)17 plants.

Na 2 SO 4 Accumulation of Foliar H 2 O 2 and Lipid Peroxidation Are Mitigated by the Overexpression of GmPgb1
A consequence of Na + toxicity is the overproduction of ROS causing severe cellular damage including peroxidation of lipid membranes, which can be estimated by the level of MDA [41]. In both the WT and GmPgb1(RNAi)23 leaves, exposure to salinity for 24 h elevated foliar H 2 O 2 and MDA ( Figure 5A). Salinity increased both H 2 O 2 and MDA, whereas in the stressed plants overexpressing GmPgb1, the increase in H 2 O 2 was limited, and the level of MDA was not different from the control ( Figure 5).

Na2SO4 Accumulation of Foliar H2O2 and Lipid Peroxidation Are Mitigated by the Overexpression of GmPgb1
A consequence of Na + toxicity is the overproduction of ROS causing severe cellular damage including peroxidation of lipid membranes, which can be estimated by the level of MDA [41]. In both the WT and GmPgb1(RNAi)23 leaves, exposure to salinity for 24 h elevated foliar H2O2 and MDA ( Figure 5A). Salinity increased both H2O2 and MDA, whereas in the stressed plants overexpressing GmPgb1 , the increase in H2O2 was limited, and the level of MDA was not different from the control ( Figure 5).
A similar pattern was observed for the expression of two respiratory burst oxidase homologs (RBOHB and D), components of the ROS generator NADPH oxidase enzyme [42] and indicators of ROS production in soybean [37]. The increase in foliar expression of both genes, more apparent after 6 h of treatment, was in fact attenuated by the overexpression of GmPgb1 ( Figure 5B). No differences between the genotypes were observed for RBOHG.  A similar pattern was observed for the expression of two respiratory burst oxidase homologs (RBOHB and D), components of the ROS generator NADPH oxidase enzyme [42] and indicators of ROS production in soybean [37]. The increase in foliar expression of both genes, more apparent after 6 h of treatment, was in fact attenuated by the overexpression of GmPgb1 ( Figure 5B). No differences between the genotypes were observed for RBOHG.

Antioxidant Responses in Leaves of Salt-Stressed Plants Are Induced by GmPgb1
Limitation of oxidative damage is a tolerance trait in several forms of abiotic stress, including salinity [43]. The cellular level of ROS is determined not only by their production, but also by the antioxidant system composed by an array of enzymes. These include superoxide dismutase (SOD), scavenging superoxide; catalases (CAT) and ascorbate peroxidase (APX) participating in the detoxification of H 2 O 2 ; and enzymes required for the turnover of ascorbic acid, dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR), and glutathione reductase (GR).
Expression analyses of the genes encoding these enzymes and responsive to abiotic stress [37] were performed. The foliar expression of SOD1 and CAT1 increased after 6 h of Na 2 SO 4 stress ( Figure 6A). This increase was more pronounced in the GmPgb1(S)17 plants and lower in the the GmPgb1(RNAi)23 line relative to the WT. A similar profile was also observed for CAT2 despite an earlier (0-6 h) increase in expression. Salinity induced SOD3 in the three genotypes, and after 24 h, the expression of this gene was greatest in the GmPgb1(S)17 plants, intermediate in the WT plants, and lowest in those suppressing GmPgb1 ( Figure 6A).
Transcript levels of the enzymes participating in glutathione and ascorbate metabolism were also measured. Salinity increased the foliar expression of APX1, especially in the GmPgb1(S)17 line ( Figure 6B). No major differences between the lines were observed for MDHAR1 and 4 despite an increase in expression in the WT leaves after 12 h. The expression of DHAR1 increased between 0 and 6 h independently of the genotype; at 24 h, a higher expression level was observed in leaves, upregulating GmPgb1. Of the two GR genes, the expression of GR1 did not exhibit pronounced fluctuations during the Na 2 SO 4 treatment, while a gradual suppression of GR2 was recorded in the GmPgb1(S)17 line. This contrasted with the WT and the GmPgb1(RNAi)23 line which showed an increase in GR2 transcripts after 6 h ( Figure 6B).
To further interpret the changes in transcript levels, the activities of the antioxidant enzymes were measured in leaves of plants grown for 24 h under control or Na 2 SO 4 conditions. The activities of many enzymes, similar under control conditions among the three genotypes, were induced by salt stress (Figure 7). Relative to the WT and GmPgb1(S)23 lines, the foliar activities of APX, DHAR, and CAT were higher in the GmPgb1(S)17 line under salinity stress. The three lines exhibited similar activity values for MDHAR and GR, while the activity of SOD was slightly reduced in the GmPgb1(RNAi)23 line when compared to the values of the other two lines exposed to stress (Figure 7).

GmPgb1 Expression in Commercial Soybean Cultivars
The increased salinity tolerance in plants overexpressing GmPgb1 prompted an analysis of the natural germplasm to determine if a relationship exists between GmPgb1 expression and salt tolerance. Based on the combined analyses of several morphological parameters used on 18 commercial cultivars (Supplemental Figure S3), we measured the level of GmPgb1 in three cultivars (1-3) exhibiting the highest tolerance to Na 2 SO 4 and three cultivars (16)(17)(18)

GmPgb1 Expression in Commercial Soybean Cultivars
The increased salinity tolerance in plants overexpressing GmPgb1 prompted an analysis of the natural germplasm to determine if a relationship exists between GmPgb1 expression and salt tolerance. Based on the combined analyses of several morphological parameters used on 18 commercial cultivars (Supplemental Figure S3), we measured the level of GmPgb1 in three cultivars (1-3) exhibiting the highest tolerance to Na2SO4 and three cultivars (16)(17)(18) exhibiting the highest susceptibility ( Figure 8). The expression of GmPgb1 increased during the first 12 h of stress in roots and especially shoots of tolerant lines while it remained low in the susceptible cultivars ( Figure 8). . Expression level of GmPgb1 in roots and leaves of commercial soybean cultivars exposed to Na2SO4. Cultivars were denoted as tolerant (1, 2, 3) or susceptible (16,17,18) to Na2SO4 stress. Values are the means ± SE of three biological replicates, each consisting of at least three plants.

Discussion
Overexpression of GmPgb1 alleviated salinity stress in soybean by mitigating the saltinduced damage and the reduction of several morphophysiological parameters. In both WT and GmPgb1-downregulating plants, exposure to Na2SO4 caused a reduction in bio- Figure 8. Expression level of GmPgb1 in roots and leaves of commercial soybean cultivars exposed to Na 2 SO 4 . Cultivars were denoted as tolerant (1,2,3) or susceptible (16,17,18) to Na 2 SO 4 stress. Values are the means ± SE of three biological replicates, each consisting of at least three plants.

Discussion
Overexpression of GmPgb1 alleviated salinity stress in soybean by mitigating the saltinduced damage and the reduction of several morphophysiological parameters. In both WT and GmPgb1-downregulating plants, exposure to Na 2 SO 4 caused a reduction in biomass and a depression of the photosynthetic rate, the stomatal conductance, the transpiration, and the water use efficiency (Figure 2). These effects, typical of salt-treated soybean plants [5,6], may be the consequence of the combinatorial effects of salinity resulting from the osmotic stress reducing the plant's hydraulic conductivity and, most importantly, ionic stress caused by the toxic accumulation of ions, especially of Na + [6,7,22,44]. Noxious levels of Na + depress photosynthesis [13] by disrupting the function of the photosystem II reaction center and the oxygen-evolving complex [3]. In soybean, this can be further aggravated by the concomitant reduction in foliar K + (Figure 4) attributable to the Na + interference with its translocation [17,45]. Many photosynthetic reactions, in fact, require K + [16].
Salt-treated WT plants also experience oxidative stress, as evidenced by the increase in H 2 O 2 and MDA, a marker of lipid peroxidation ( Figure 5A), and these effects can be linked to excessive levels of Na + [46,47]. Consistent with this observation is the transcriptional induction of RBOHs ( Figure 5B), components of the NADPH oxidase multicomplex system [41] and markers of ROS production [37], and the activation of several ROS-detoxifying enzymes, including APX, CAT, and SOD ( Figure 5A). Thus, the responses observed in Na 2 SO 4 in the WT leaves are typical in many respects to those observed in many species [46], including soybean [5].
Overexpression of GmPgb1 ameliorates plant response to salinity through the retention of photosynthetic gas exchange parameters and mitigates foliar tissue damage (Figures 2 and 3). These effects are associated to the ability of these plants to limit accumulation of Na + (and retain a low Na + /K + ratio) and mitigate oxidative stress. A plausible reason for the reduced Na + levels in tissue overexpressing GmPgb1 is the increase in SOS1 ( Figure 4D), a member of the SOS pathway contributing to salinity tolerance [48]. Strategically located in epidermal, cortical, and pericycle cells, SOS1 limits the Na + intake in root cells, restricts the radial movement of the ion across cells layers, thereby reducing its long-distance transport through the xylem, and excludes Na + in leaf cells [16]. The induction of GmSOS1 in both roots and leaves of the Na 2 SO 4 -treated GmPgb1(S)17 plants, consistent with the higher expression levels observed in species tolerant to salt stress [49] and the acquired tolerance of transgenic plants overexpressing SOS1 [20], suggests limited Na + intake in root cells and transport to the aerial organs.
The increased expression of GmATK1 in the roots of the salt stress-exposed GmPgb1(S)17 plants ( Figure 4D) is also an indication of increased tolerance [45,50]. As an inward-rectifying K + channel, ATK1 is required for the maintenance of the Na + /K + balance by ensuring sustained K + intake, especially in high salt environments. The crucial function of this transporter during salinity stress was demonstrated by Wang et al. [45] documenting enhanced tolerance and retention of the Na + /K + balance in Arabidopsis plants overexpressing the soybean GmATK1.
Besides restricting Na + from leaves, overexpression of GmPgb1 alleviates salt damage by relieving oxidative stress, specifically reducing the level of ROS and ROS-induced damage. This protective role, documented under developmental and stress events [26,29,30], is linked to the function of the protein to scavenge NO [34], as also observed for GmPgb1 (Supplemental Figure S4). Here, we show that expression of GmPgb1 during salt stress attenuates the increase in H 2 O 2 and MDA, possibly by limiting the Na 2 SO 4 induction of RBOHB and RBOHD, markers of ROS production in soybean [37]. The mitigation of ROS production in GmPgb1-overexpressing leaves is a contributing factor for a decreased leaf injury and a higher retention of the photosynthetic rate. Additionally, overexpression of GmPgb1 also attenuates the Na 2 SO 4 depression of the transcripts involved in the formation of photosystem I (PsaA and B), cytochrome B6F (cylB6F), and photosystem II (PsbB, PsbC, PsbD). Retention of these transcripts during conditions of stress has been associated to abiotic stress tolerance [37].
Reduced oxidative stress in leaves overexpressing GmPgb1 can be attributed to an improved antioxidant system. Relative to the WT, overexpression of GmPgb1 is associated to higher levels of the CAT1, 2 and APX1 transcripts, as well as the CAT and APX activity. This behavior might be accountable for the lower level of H 2 O 2 measured in the GmPgb1(S)17 leaves during salt stress. Elevated activity of CAT and APX are desirable tolerance traits for salinity stress [22]. While the direct mechanisms through which GmPgb1 modulates the expression and activity of CAT and APX1 are unknown, regulation of the APX1 activity might be the result of post-translational modifications mediated by NO. Nitration of APX by peroxynitrate, formed in environments enriched in NO, inhibits the enzyme [51], an observation consistent with the higher APX activity (Figure 7) in GmPgb1 leaves which are depleted in NO (Supplemental Figure S4).
Sustained activity of APX and removal of H 2 O 2 during conditions of stress is dependent upon the availability of the substrate ascorbate produced de novo and regenerated through the ascorbate-glutathione cycle [43]. Our work suggests that during salt stress, overexpression of GmPgb1 might result in higher levels of ascorbate required for the operation of APX. First, the expression of the three ascorbate biosynthetic genes GDP-L-galactose phosphorylase GmGPP and L-galactose-dehydrogenase (GmGDH1 and 2), good indicators of ascorbate production in soybean [37,52], is induced in leaves of Na 2 SO 4 -treated plants overexpressing GmPgb1 (Supplemental Figure S5). Second, the same leaves exhibit a higher activity of DHAR, which is required for ascorbate recycling through the glutathionedependent dehydroascorbate reduction [23]. The relevance of this enzyme in replenishing the ascorbate pool during salinity stress was documented [53]. The slightly lower ascorbate level in the salt-stressed GmPgb1(S)17 leaves (Supplemental Figure S6) is difficult to interpret, but it might be a consequence of higher utilization by APX. While we suggest that removal of ROS by ascorbate is mediated by APX, the operation of APX-independent mechanisms, which also involve ascorbate [54], cannot be discounted. Activation of the antioxidant components and similar Na + /K + metabolic adaptations to those described here were associated to the ability of thymol to increase salinity tolerance in tobacco plants [15]. In the same study, it was observed that thymol elevates the level of NO, an observation which is, however, not consistent with our study.
The relevance of high expression of GmPgb1 as a mechanism to cope with salinity stress is not limited to transgenic material. Our study shows that this gene is preferentially induced in commercial cultivars exhibiting tolerance to Na 2 SO 4 , relative to those susceptible to the same stress. This observation, consistent with a previous work linking induction of Pgbs to the ability of soybean plants to cope with hypoxic stress [37], could be exploited as a tool to screen plant material at early stages of development for salt stress tolerance.
In conclusion, this study demonstrates that overexpression of the soybean GmPgb1 enhances tolerance to salt stress. The Na 2 SO 4 depression of gas exchange parameters, water use efficiency, leaf damage, and reduction in water content are attenuated by the overexpression of GmPgb1. These beneficial effects can be attributed to the ability of plants overexpressing GmPgb1 to limit the accumulation of Na + in leaf tissue (thus maintaining a lower Na + /K + ratio) and mitigate oxidative stress. This is achieved by reducing production of ROS as well as by activating antioxidant responses centered around the activities of CAT and APX.
Collectively, these data suggest GmPgb1 as an effective modulator of soybean response to salt stress.

Plant Material and Growth Conditions
Transformed soybean (Glycine max) plants with altered levels of the Phytoglobin1 gene GmPgb1, described by Mira et al. [37], were selected for this study. Soybean cultivars grown in Western Canada were provided by Dennis Lange (Manitoba Agriculture, Winnipeg, Canada). These cultivars were the same as the ones utilized previously [55].
Seed sterilization was carried out with chlorine gas for 20 h [56]. The sterilized seeds were germinated on moistened filter paper for about one week. Uniform seedlings characterized by unfolded unifoliate leaves (VC stage of development) were selected and grown further hydroponically in a solution containing half-strength Hoagland [57]. Acclimation to the hydroponic conditions was conducted for 7 days under a 16 h photoperiod of 22 • C (light) and 18 • C (dark), with 334 µmol m −2 s −1 light intensity. As soon as the plants developed their first trifoliate leaf (V1 stage of development) they were transferred to a half-strength Hoagland solution containing 100 mM Na 2 SO 4 . The level of Na 2 SO 4 and the timing of the measurements were determined empirically to ensure consistent manifestation of symptoms. Plants in the control group were given identical treatments as the experimental plants, except for Na 2 SO 4 . The plants were harvested at different timepoints for morphophysiological parameter assessments, as indicated below. Unless otherwise specified, all the experiments carried out used at least five biological replicates, each with a minimum of five plants.

Determination of Morphophysiological Parameters
Fresh and oven-dry weights of shoot or root tissues were measured after 7 days of treatment. The extent of leaf injury, also measured at 7 days, was assessed using a visual scale adopted in a previous work [55], based on an index (0-100) reflecting the percentage of foliar damage.
Gas exchange parameters (photosynthetic rate, stomatal conductivity, water use efficiency, and transpiration rate) were measured on the first trifoliate leaves after 24 h of treatment using a gas exchange infrared gas analyzer (IRGA; LI-6400, LI-COR, Inc., Lincoln, NE, USA), as reported by Youssef et al. [55]. The measurements were conducted between 9 am and 3 pm. Gas exchange parameters were measured under photosynthetically active radiation of 400 µmol m −2 s −1 , atmospheric CO 2 of 400 µmol mol −1 , relative humidity of 50%, and temperature of 22 • C. The measurements were performed on three biological replicates, each consisting of at least five plants.

Measurements of K + and Na +
Leaf and root tissues (2 g) were harvested at 24 h under control or Na 2 SO 4 treatment and oven-dried for 5 days. Analysis of K + and Na + was performed using Agvise (Agvise Laboratories Inc., Northwood, ND, USA; www.agvise.com, 29 March 2022). Extraction of both elements was performed using 1 M ammonium acetate (pH 7), and determination was carried out by means of inductively coupled plasma atomic emission spectroscopy (ICP-AES).

Stomata Measurements
Analysis of stomata was performed after 24 of stress [55]. Briefly, a thin layer of translucent nail polish was applied to a leaf's abaxial side and left to dry for 5 min. A transparent tape was applied on the dry nail polish and carefully peeled off before being placed on a clean microscope slide for microscopic examination. Stomata were categorized into three groups based on their length-to-width ratio (L/W): fully open (L/W = 1-2), partially open (L/W = 2.1-2.5), or closed (L/W 2.6) (Supplemental Figure S1). The analyses were conducted using three biological replicates each, with at least 100 stomata.

Measurements of Transcript Levels
Transcript levels were measured on the first trifoliate leaves at 0, 6, 12, and 24 h of stress. Tissue was harvested and stored at −80 • C. Total RNA was extracted with a TRI Reagent solution (Invitrogen, Waltham, MA, USA) following the manufacturer's instructions. The RNA samples were treated with DNase (Thermo Fisher Scientific, Waltham, MA, USA), and a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Waltham, MA, USA) was used for cDNA synthesis. The relative levels of transcripts were measured by means of quantitative PCR using the primers listed in Supplemental Table S1 and analyzed with the 2 −∆∆CT method [58] using helicase (Gm04g07180.4) as the reference gene.

Antioxidant Enzyme Activity and Quantification of Ascorbic Acid and Glutathione
The enzyme activity of APX, DHAR, MDHAR, and GR were measured as reported by Zhang and Kirkham [59]. Trifoliate leaf tissue (100-150 mg fresh weight), collected at 24 h under control or Na 2 SO 4 treatment, was immediately homogenized on ice with 50 mM phosphate buffer (pH 7.8) containing 0.2 mM EDTA and 1% (w/v) polyvinylpyrrolidone. The homogenate was then centrifuged at 18,000× g at 4 • C for 15 min. The supernatant was collected for the analysis of enzyme activities [59] using a Bio-Tek Synergy H4 Microplate Reader.
The APX activity was estimated by measuring the H 2 O 2 -dependent oxidation of ascorbic acid following the decrease in absorbance at 290 nm. The reaction mixture contained 0.5 mM ASC, 0.1 mM H 2 O 2 , 0.1 mM EDTA, and 50 mM sodium phosphate buffer (pH 7.0).
The DHAR activity was measured by following the GSH-dependent production of ASC at 265 nm. The reaction mixture contained 2.5 mM GSH, 0.1 mM EDTA, 0.2 mM dehydroascorbate (DHA), and 50 mM sodium phosphate buffer (pH 7.0).
The MDHAR activity was measured by observing the change in absorbance at 340 nm. The reaction mixture contained 0.1 mM NADH, 2.5 mM ASC, 50 mM sodium phosphate buffer (pH 7.6), and four units of ASC oxidase.
The GR activity was determined by observing the NADPH-dependent oxidation of GSSG following the decrease in absorbance at 340 nm. The reaction mixture contained 1 mM EDTA, 1 mM GSSG, 0.2 mM NADPH, and 0.1 M sodium phosphate buffer (pH 7.8).
For CAT and SOD, leaf tissue (100-150 mg fresh weight) was homogenized on icecold 50 mM sodium phosphate buffer (pH 7.0) containing 0.2 mM EDTA and 1% (w/v) polyvinylpyrrolidone. The homogenate was centrifuged at 15,000× g at 4 • C for 20 min. The supernatant was then collected for the enzyme activity analysis [59].
For the analysis of ascorbic acid (ASC) and dehydroascorbate (DHA), about 100-150 mg (fresh weight) of leaf tissue was first powdered in liquid nitrogen and then homogenized with 2 mL of ice-cold 5% metaphosphoric acid. The homogenate was centrifuged at 21,000× g at 4 • C for 15 min. The supernatant was collected and the absorbance (525 nm) was read [59].

Lipid Peroxidation and H 2 O 2 Measurements
Lipid peroxidation was estimated by measuring the amount of malondialdehyde (MDA) [60]. The leaves (about 0.5 g) collected at 24 h under control or Na 2 SO 4 treatment were homogenized in 10 mL of 0.1% (w/v) trichloroacetic acid (TCA), centrifuged at 15,000× g for 10 min, and 1 mL of the supernatant was added to 4 mL thiobarbituric acid (TBA) in 20% (w/v) TCA. The solution was heated at 95 • C for 30 min centrifuged at 10,000× g for 10 min and used to measure the absorbance at 532 nm. A standard curve was generated to determine the content of MDA in the tissue.
Hydrogen peroxide was measured in trifoliate leaves using the ferrous oxidation protocol with xylenol orange [61]. The level of H 2 O 2 was determined using a standard curve.

Nitric Oxide (NO) Localization
The accumulation of intracellular NO in the leaf tissue was monitored using the NO-specific fluorescent probe DAF-FM DA (MilliporeSigma, Burlington, MA, USA), as described previously [62,63]. The leaves, collected at 24 h under control or Na 2 SO 4 treatment, were incubated for 1 h at 37 • C with gentle agitation in the dark in a buffer (50 mM Tris and 50 mM KCl, pH 7.2) containing 1% (v/v) Triton X-100 and 50 µM DAF-FM DA. The leaves were rinsed twice in the same buffer (50 mM Tris and 50 mM KCl, pH 7.2). The fluorescent signal was captured with a GFP light cube in a EVOS M5000 Imaging System, with excitation and emission peaks at 470 and 510 nm. Image J was used to determine the DAF-FM DA fluorescent pixel intensity in leaf segments.

Statistical Analysis
The data were analysed by means of one-way analysis of variance (ANOVA) using the SPSS program (version 19.0). Treatment means were compared using Duncan's test (α = 0.05) to differentiate the significance of the difference between various parameters.