The Manifestation of the Dual ROS-Processing and Redox Signaling Roles of Glutathione Peroxidase-Like Enzymes in Development of Arabidopsis Seedlings

1. Introduction

Adverse environmental conditions may increase the production of reactive oxygen species (ROS), such as superoxide radicals (O₂^•-), hydrogen peroxide (H₂O₂), hydroxyl radicals (OH^•) and singlet oxygen (¹O₂). The elevated level of ROS may induce detrimental oxidation of macromolecules including lipids, proteins and nucleic acids thus they may cause damage or irreversible effects on development of tissues and organs [1,2]. However, a lot of evidence support that change in endogenous oxidant levels can fulfil signalling functions regulating normal plant growth and responses to stresses [3,4,5,6]. ROS play a central role in integrating external and intracellular signals and, in addition to stress responses, are important regulators of plant growth and development [7,8]. Although the role of ROS in signal transduction is widely accepted, the precise signaling mechanisms leading to specific responses remain poorly understood.

To keep ROS levels tightly controlled and minimize ROS-derived impairments, different non-enzymatic antioxidants (such as ascorbate, glutathione, carotenoids, tocopherols) and ROS-processing enzymes have evolved in aerobic organisms [2]. The antioxidant enzymes include superoxide dismutase (SOD), catalase (CAT) and different peroxidases, such as ascorbate peroxidase (APX), guaiacol peroxidase (POX) and glutathione peroxidase. The mammalian glutathione peroxidases (GPXs) are pivotal ROS-eliminating antioxidant enzymes that can reduce H₂O₂ and organic peroxide substrates to water or corresponding alcohols using reduced glutathione (GSH) as electron donor [9,10,11]. Beside participating in the maintenance of membrane integrity they can oxidase cysteine-containing proteins involved in signaling, such as phosphatases, kinases, and transcription factors and thus induce or regulate different pathways [12,13]. GPXs are key players in many biological processes, such as fertility, anti-inflammatory- and anti-carcinogenesis associated routes [11,14,15]. In contrast to the animal GPXs, the plant isoenzymes are selenium-independent monomeric proteins containing cysteine in their catalytic site, therefore show lower peroxidase activity than the animal enzymes and most of them prefer thioredoxin (TRX) as electron donor rather than GSH [16,17,18]. Attacha et al. [19] proposed using the GPX-like (GPXL) nomenclature for the Arabidopsis thaliana isoforms to avoid any confusion resulting from protein names. Now we apply the GPXL abbreviation in all referees to the earlier published information about plant glutathione peroxidases.

The transcription levels of plant GPXLs were increased by various environmental stress conditions and numerous scientific reports provided evidence for their importance in different stress responses ([18] and references therein, [20,21,22,23]). Overexpression of GPXL genes led to higher tolerance against abiotic stress factors in tomato [24,25], rice [26], tobacco [27], A. thaliana [28,29,30] and Salvia miltiorrhiza [31]. Moreover, GPXLs are involved in development and growth of plants. Passaia et al. [32] compared phenotypes of 4-week-old Arabidopsis Atgpxl mutants and observed that although their shoot was largely similar to wild type plants, minor differences were found in the number of rosette leaves and lateral roots in the Atgpxl2, Atgpxl3, Atgpxl7, and Atgpxl8 mutants. Connections among the AtGPXLs and phytohormones such as auxin, abscisic acid, strigolactone hormones were suggested, thereby demonstrating the importance of AtGPXLs in the hormone-mediated regulation especially of lateral root development. These authors also proposed that GPXLs may be required to mediate glutathione and reduced thioredoxin functions in roots that impact on lateral root production or growth [32]. GSH and TRX are the main components of cellular redox homeostasis and are important for the processes that determine the root architecture [33]. In the case of acute shortage of glutathione more TRX is used as electron donor compared to the GSH, suggesting a possible connection between the GSH and TRX pools [34].

The involvement of the glutathione (γ-Glu-Cys-Gly) and GSH-related antioxidant enzymes in the redox homeostasis and signalling came to the forefront in last decades [18,20,35,36,37,38,39]. Glutathione is classically considered as one of the main antioxidants with low molecular weight, playing crucial role in ROS detoxification in cooperation with ascorbate in the Foyer-Halliwell-Asada pathway and as an electron donor to diverse antioxidant enzymes [40]. Reactive oxygen- and nitrogen species, ROS-processing enzymes, and other molecules with antioxidant activity, along with their redox states, all contribute to the overall redox homeostasis of the cell [41]. Now it is widely accepted that the ratio of reduced and oxidized glutathione (GSH/GSSG) and the glutathione half-cell reduction potential (E_GSSG/2GSH), which depends on the absolute glutathione concentration and the [GSH]:[GSSG] ratio, are effective markers of the redox homeostasis [42,43,44]. Enhanced ROS production temporarily shifts the redox potential to more oxidised values [45], which was proven in vivo by genetically encoded redox probes [46,47].

Genetically encoded biosensors are fluorescent proteins that change their properties specifically in response to a physiological parameter [48]. Among them, the roGFP2 (reduction-oxidation sensitive green fluorescent protein 2) is applied the most frequently for examining the redox state in plants, e.g., [19,34,48,49,50]. In this redox sensor specific amino acids of wild-type fluorescent protein that influence chromophore structure have been replaced and two Cys residues were introduced. The redox status of the Cys pair is dependent on GSH/GSSG redox state, and the rapid equilibration between them is facilitated by a linked glutaredoxin (GRX) molecule [46,51,52]. The oxidation of roGFP induces a change in the protonated state of the chromophore: it increases fluorescence intensity at the excitation band associated with the protonated form (400 nm) and decreases fluorescence intensity at the excitation band corresponding to the anionic form (490 nm), enabling ratiometric fluorescence readout. From the fluorescence intensities detected at the two excitation wavelengths, we can calculate the glutathione-dependent redox potential [46,53].

Being thiol peroxidases that harbour redox-active Cys in their active site, GPXL proteins are very sensitive to oxidation. Using GSH and/or TRX as a reductant, the GPXLs may control the redox status of these main redox compounds, also it may be an important connection between them. They can modify the thiol/disulfide balance and activity of other proteins and transcription factors [20,23]. They not only were considered to function as redox sensors but also as a link between ROS and functional redox signaling [20,32,54,55]. As efficient ROS-processing antioxidant enzymes GPXLs are involved in control ROS level and shaping ROS gradients e.g., in the maintenance of stem cell niche or triggering differentiation in the shoot and root apical meristems (SAM and RAM, respectively), and in the proper zygote/embryo development [56].

In the present study, we have introduced the roGFP2 redox probe into the Atgpxl1-8 mutants to detect the redox status of seedlings. Measuring the glutathione redox potential was conducted on 6-day-old seedlings. Our aim was to analyse the involvement of AtGPXLs in the growth and development of roots of seedlings and their role in the maintenance of ROS levels in shoots and roots. Beside the evaluation of O₂^•- and the overall level of intracellular ROS content (total ROS) by fluorescent dyes, the vitality (metabolic activity) of in vitro grown roots of Atgpxl1-8 mutants were investigated under control conditions and after 14-day-long treatment with mild salt- (50 and 100 mM NaCl) and osmotic stresses (100 and 200 mM mannitol). It was concluded that AtGPXLs have a fine-tuning role in processing ROS in seedlings both under control and abiotic stress conditions. Although some other antioxidant mechanisms can be activated, especially in case of stronger stresses, these are not capable of suppressing the ROS levels and thus mutation in one of the AtGPXL genes may cause change in overall redox status.

2. Materials and Methods

2.1. Plant Material and Growth Conditions

Arabidopsis thaliana (L.) ecotype Columbia (Col-0) as a wild type and T-DNA insertional mutants of the eight glutathione peroxidase-like genes (AtGPXL1-8) were used. The T-DNA insertional mutant lines were obtained from the Nottingham Arabidopsis Stock Centre (NASC) (Atgpxl1 [AT2G25080]: SALK_128885C; Atgpxl2 [AT2G31570]: SALK_082445C; Atgpxl3 [AT2G43350]: SALK_071176C; Atgpxl4 [AT2G48150]: SAIL_623_F09; Atgpxl5 [AT3G63080]: SALK_076628C; Atgpxl6 [AT4G11600]: WiscDsLox321H10; Atgpxl7 [AT4G31870]: SALK_072007C; Atgpxl8 [AT1G63460]: SALK_127691C), and the seeds of homozygous genotypes were used. According to our previous results, Atgpxl1, -2, -5 and -6 are knockdown, while Atgpxl3, -4, -7 and -8 are knock out mutants [57].

Seeds were sown on ½ MS media [58] and were pre-grown in growth chamber (Fitoclima S 600 PLH, Aralab, Rio de Mouro, Portugal) at 21°C under 100 mmol m^-2 s^-1 photon flux density with a 10 h day and 14 h night photoperiod. The redox potential was measured in 6-day-old seedlings. To estimate the growth parameters and ROS levels under different abiotic stresses, the 4-day-old seedlings were transferred to ½ MS media supplemented with 50 mM, 100 mM NaCl, or 100 mM, 200 mM mannitol in square Petri plates and placed vertically in the same chamber for 2 weeks to check differences in long term responses.

2.2. Analysis of the Redox Potential by Ratiometric Measurements of the roGFP2 Fluorescent Probe

To determine the redox potential of seedlings, a fusion protein of human glutaredoxin1 (GRX1) and roGFP2 protein with cytosolic localization (c-roGFP2) was used [43,59]. This vector, kindly provided by Prof. Dr. A. Meyer, had been introduced into wild type (Col-0) Arabidopsis earlier as was published in (Horváth et al. 2019). The redox sensor was introduced into Atgpxl1-8 mutants by crossing and homozygous mutants were established. Fluorescence measurements were performed with a confocal laser scanning microscope (Olympus Fluoview FV1000, Olympus Life Science Europe GmbH, Hamburg, Germany) in the cotyledon and the apical meristematic region (proximal meristem) of 6-day-old seedlings. Excitation wavelengths were 405 nm and 488 nm, the fluorescence was detected between the 505-530 nm emission wavelength [60].

The glutathione redox potential (E_GSH) was calculated using the formula by [46]:

$$E_{GSH}=E_{roGFP2}^0-\left({\displaystyle \frac{\mathrm{2,303}\mathrm{*}R\mathrm{*}T}{z\mathrm{*}F}}\right)\mathrm{*}{log}_{10}{\displaystyle \frac{1-{OxD}_{roGFP2}}{{OxD}_{roGFP2}}}$$

where: E⁰_roGFP2, the midpoint potential of roGFP2 is -272 mV at 30◦C pH 7; R: the gas constant (8.315 J K−1 mol−1); T: the absolute temperature (298.15 K); z: the number of transferred electrons (2); F: the Faraday constant (96485 C mol⁻¹). The degree of oxidation (OxD) of roGFP2 was:

$${OxD}_{roGFP2}={\displaystyle \frac{R-R_{red}}{\left(\frac{I_{488red}}{I_{488ox}}\right)\mathrm{*}\left(R_{ox}-R\right)+\left(R-R_{red}\right)}}$$

where: R is the ratio of excitation at 405/488 nm; R_red: the ratio of fully reduced form, using 10 mM dithiothreitol; R_ox: the ratio of the fully oxidized form, using 20 mM H₂O₂; I_488ox: the intensity at 488 nm for the fully oxidized form; I_488red: the intensity at 488 nm for the fully reduced form.

2.3. Measurements of Root Parameters

Root length and number of lateral roots of the Col-0 and Atgpxl1-8 mutants were analysed using ImageJ (Fiji) software [61] after scanning the square plates as was described earlier [30]. Lateral roots were counted, and lateral root density (LRD) was calculated by dividing the number of visible lateral roots by the primary root length [32].

2.4. Detection of the Vitality, Superoxide Radical and Total ROS Levels in Roots and Leaves

A. thaliana seedlings were incubated for 15 min in 3 mL 10 μM fluorescein diacetate (FDA) staining solution (prepared in 10/50 mM MES/KCl buffer, pH 6.15) to determine cell vitality [60], then washed four times with MES/KCl before detection.

For visualisation of total ROS level, 2’,7’-dichlorodihydrofluorescein diacetate (H₂DCFDA) fluorescent dye was used. Seedlings were incubated in 3 mL 50 μM H₂DCFDA staining solution (prepared in 50 mM sodium phosphate buffer, pH 7.5) for 30 min and washed twice with the buffer before detection [62].

Dihydroethidium (DHE) in Tris–HCl (10 mM, pH 7.4) buffer was used to detect superoxide (O₂^•-) radical level in roots (with the available filter the detection in the shoots was not possible). Seedlings were incubated in darkness for 30 min with 3 mL 10 μM DHE and then the samples were washed twice with buffer before detection [60].

Zeiss Axiowert 200M microscope (Carl Zeiss, Jena, Germany) equipped with a high-resolution digital camera (Axiocam HR, HQ CCD, Carl Zeiss, Jena, Germany) and filter set 10 (excitation 535–585 nm, emission 600–655 nm) for FDA and H₂DCFDA, the filter set 9 (exc.: 450–490 nm, em.: 515–∞ nm) for DHE was used to detect fluorescence in plants.

The intensities of fluorescein-2’,7’-dichlorofluorescein- and oxyethidium fluorescence were quantified on digital images using Axiovision Rel. 4.8 software in the proximal meristem of the roots in a circle with 50 μm radius [60], or in the middle of the leaves in a circle with 150 μm radius. The measurements were performed in three independent experiments (n ≥ 15) with the same microscopic settings.

2.5. Correlation Analysis

To determine the relationship between the investigated variable parameters, we calculated Pearson’s correlation coefficients [57]. The values of the correlation coefficient varied between +1 and −1. When the value is around +1 or −1, it indicates a close positive or negative relationship between the variables, respectively. As the correlation coefficient value approximates 0, the relationship between the two variables will become weaker.

2.6. Statistical Analysis

Experiments were carried out at least two times. Redox potential values presented here represent mean with standard deviation (± SD, n ≥ 3), the other data are mean with standard error (±SE), n ≥ 45. Statistical analysis was carried out with SigmaPlot 12.0 software (SigmaPlot, Milano, Italy). After analysis of variance (ANOVA), Duncan’s multiple comparisons were performed. Means were significantly different if p ≤ 0.05.

3. Results

3.1. The Redox Status of 6-Day-Old Atgpxl Seedlings Are More Oxidized Compared to the Wild Type

The determination of the redox potential, using the redox sensitive fluorescent probe (roGFP2) in vivo, uncovered more oxidized redox state in most of the untreated mutants (Figure 1, Table 1). In leaves, the highest oxidized redox status was found in gpxl1 and -7 seedlings, and ca. 25 mV difference in the redox status was detected between these mutants and the wild type plants. In roots, more than 20 mV increase was detected in case of gpxl1 and -4 mutants. Interestingly, in the cotyledons the redox potential of gpxl4 and -5 mutants was comparable to the wild type plants, moreover in the roots the E_GSH of gpxl3, -5, and -7 did not differ significantly from the Col-0 (Table 1). This may indicate that other mechanisms might complement the default of some AtGPXL isoenzymes in mutants, or the coding enzymes might have no relevant function in the maintenance of the redox homeostasis and ROS-processing in seedlings.

3.2. Several Atgpxl Mutants Had Longer Primary Roots Than the Wild Type Controls, and 50 mM NaCl Promoted the Growth of All Genotypes Compared to the Untreated Col-0 Primary Roots

The involvement of AtGPXLs in plant growth under control and abiotic stress responses was investigated after exposing the 4-day-old Atgpxl mutants to salt- and osmotic stresses for two weeks. Measuring the length of the primary roots of untreated 18-day-old seedlings revealed that the maximal root length of the Atgpxl1 and -3 was decreased, while the Atgpxl2, -5 and -7 did not differ significantly from the Col-0, but the growth of Atgpxl4, -6 and -8 roots was increased. Counting the number of visible lateral roots uncovered that Atgpxl5, -6 and -8 mutants had more lateral roots than the wild type plants, but the lateral root density was higher only at the Atgpxl5 roots (Figure 2).

Different concentrations of NaCl had inverse effect on the lengthwise growth of the wild type roots: 50 mM salt concentration promoted, while 100 mM inhibited root growth. Among the mutants, the root growth in the Atgpxl6, -7 and -8 mutants were increased compared to the Col-0. Interestingly, with the exception of Atgpxl3, all mutants had longer roots in the presence of 100 mM NaCl than the wild type plants. However, the number of lateral roots differed significantly from the salt stressed wild type roots only in the 100 mM NaCl treated Atgpxl3 roots and the density of lateral roots was significantly lower in the 100 mM NaCl treated Atgpxl1 mutant (Figure 2 (a)panel).

The used isosmotic mannitol concentrations decreased both the total length of primary roots and the number of lateral roots in the wild type seedlings after the 14-day-long treatment, but the density of lateral root was not affected (Figure 2 (b) panel). Similarly to the effect of 100 mM NaCl treatment, the shortest root length among the mutants was measured at 200 mM treated Atgpxl3, but all the other mutants had longer roots in the presence of 200 mM mannitol than the wild type plants. The roots of Atgpxl2, -7 and -8 mutants grew better than Col-0 in both mannitol treatments. Although significant differences were found only at two mutants in the number of lateral roots, compared to the wild type plants, but the lateral root density of Atgpxl1, -2, -3 and -5 was lower after applying 200 mM mannitol, while in Atgpxl6, -7 and -8 this parameter decreased after both osmotic stress treatments (Figure 2 (b) panel).

3.3. Slightly Elevated ROS Levels Decreased the Vitality of the Mutant Roots

Using dihydroethidium allows us to investigate the superoxide (O₂^•-) levels in roots. Detection of the DHE fluorescence in the non-treated 18-day-old seedlings revealed elevated O₂^•- levels in the roots of Atgpxl2, -3, -4, -5 and -7 seedlings under control conditions (Figure 3 and Figure 4). Supplementation of the media with 50 or 100 mM NaCl increased the O₂^•- level in wild type roots by 37 and 59%, respectively, after 14 days. These treatments usually resulted in moderated increase in the fluorescence intensity of mutants. Exceptions are the effect of the 50 mM NaCl on Atgpxl1 and the 100 mM NaCl on Atgpxl5 and -6 roots, where ca. 185, 210 and 195% O₂^•- accumulation was detected, respectively, as shown in Figure 3 (a) panel as relative percentages to the untreated Col-0 control. On the contrary, the isosmotic 100 mM mannitol did not change the O₂^•- content of the wild type roots and, except the elevation in Atgpxl6 roots, this stress did not significantly influence the level of this ROS in the mutants. However, in the presence of 200 mM mannitol the O₂^•- content of Col-0 increased by ca. 50%. Due to the applied osmotic stress, the O₂^•- accumulation was similar in most the mutants as well, but it was even higher in the Atgpxl1 and -5 roots (Figure 3 (b) panel, Figure 4), which indicates specific roles of the encoded AtGPXLs in the processing of O₂^•.

Detecting the total ROS levels in the roots by H₂DCFDA staining revealed that untreated Atgpxl mutants did not accumulate significantly higher amount of ROS compared to the wild type plants and, except the decreased vitality of Atgpxl1 after the osmotic stress treatment, their vitality was similar to the Col-0 roots (Figure 3, Figure 5 and Figure 6).

50 and 100 mM NaCl elevated the total ROS levels of the Col-0 plants by 55 and 39%, respectively, while applying 100- or 200 mM mannitol increased them by ca. 35%, but these changes did not prove to be significant. The vitality of the wild type roots decreased by 26 and 71% due to 50 and 100 mM NaCl treatments, respectively, and the isosmotic 100- and 200 mM mannitol caused ca. 55 and 70% decrease, respectively. In the mutants, the investigated parameters changed rather similarly to the wild type, with a few exceptions, such as the higher ROS accumulation in Atgpxl7 roots after treatment with 100 mM NaCl, in Atgpxl6 after applying 100 mM mannitol, and in Atgpxl2, -3 and -4 roots after 200 mM mannitol treatments (Figure 3, Figure 5 and Figure 6).

3.4. The Leaves of Untreated Mutant Seedlings Accumulated More Total ROS and Most of Them Showed Less Metabolic Activity

Comparison of the ROS levels of the 18-day-old untreated leaves revealed that the mutation of a single AtGPXL gene caused elevated dichlorofluorescein fluorescence in all investigated T-DNA insertional mutants’ shoots compared to Col-0. Parallelly, the vitality detected by FDA decreased in most of the mutants (exceptions are Atgpxl1 and Atgpxl2).

In the presence of 50 and 100 mM NaCl, the total ROS level of the Col-0 increased by 56 and 127%, while the vitality decreased by 20 and 40%, respectively. There was no significant difference between the ROS level of Atgpxl1-8 mutants and Col-0 cotyledons after applying 50 mM NaCl treatment, but the vitality of the Atgpxl2-8 mutants was lower than in the wild type. Interestingly, the 100 mM NaCl treatment did not result in any difference in the ROS levels and in the vitality of seedling’s shoots, moreover the ROS content in Atgpxl8 shoots were even lower than in other genotypes (Figure 7 (a) panel), suggesting that other ROS processing mechanisms were induced and might ensure similar ROS levels.

The 14-day-long 100- and 200 mM mannitol treatments elevated the ROS level in the wild type by ca. 80% simultaneously, meanwhile the vitality decreased by 34 and 46%, respectively. Generally, the ROS and the vitality in the mutants’ leaves changed similarly after mannitol treatment as in the presence of isosmotic NaCl concentrations, but more ROS accumulated in the leaves of Atgpxl3, -4 and -8 seedlings in the case of 200 mM mannitol. The vitality decreased more in the Atgpxl2, -5 and -8 mutants compared to the wild type plants after applying 100 mM mannitol, but there was no difference between the metabolic activity of the seedlings after treatments with higher concentration of mannitol (Figure 7 (b) panel).

3.5. Correlations Between the Measured Parameters Strengthen the Involvement of AtGPXLs in ROS Homeostasis

Since the determination of E_GSH was performed on 6-day-old seedlings (when the lesser amount of cell layers allow more precise detection of ratiometric roGFP2 fluorescence), but the other parameters were measured after two weeks of stress treatments, the correlation analysis was performed on data including the growth parameters, ROS levels and metabolic activity of 18-day-old seedlings. Our results indicated a very strong positive correlation between the shoot and root vitality, which were in negative correlation with ROS levels in shoots and roots in all plants (Figure 8 A, B). However, the correlation between the ROS levels and the vitality in the two organs is not that strong, even absent in the case of Atgpxl1, -2 and -6 (Figure 8 C, D). This analysis revealed a stronger negative correlation between the ROS levels of wild type shoot and development of lateral roots, than between the ROS of roots proximal meristematic region and growth of lateral roots. These correlations altered in the Atgpxl mutants: in some cases they were absent (Atgpxl1, -5 and -7), while in some other mutants (Atgpxl4, -6 and -8) they became stronger (Figure 8 E-H).

4. Discussion

4.1. The ROS-Processing Roles of AtGPXLs in the Mutants Were Not Substituted Completely Either Under Control or Abiotic Stress Conditions

Several reports provided direct or indirect evidence that plant GPXLs may have multiple roles particularly during abiotic stresses. Chen et al. [63] reported that tomato GPXLs in tobacco plants functioned as cytoprotectors preventing Bax-, heat-, salt- and abiotic stress-induced cell death. Beside protection of biological membranes by the reduction of lipid peroxides GPXLs may preserve proteins and DNA in case of oxidative stress [16,17,28,63]. It was proven that A. thaliana plants overexpressing AtGPXL8 had less oxidized proteins and -nucleotides under oxidative stress [28]. It was suggested that their role is more important in the detoxification of peroxides other than H₂O₂ [16,17], or even, they might have a crucial role in redox signaling [54]. However, the spectrophotometrically measured H₂O₂ level in our previous experiments in Atgpxl5 mutant plants was higher both under control conditions and in the presence of 100 mM NaCl than that in the wild type, while in OX-AtGPXL5 plants it was on the same level as in untreated Col-0 wild type plants [30]. Our present results confirm that AtGPXLs are involved in the maintenance of the ROS homeostasis and mutation in one of their coding genes can negatively affect the vitality of seedlings. Interestingly, the highest differences among the O₂^•- and total ROS levels of wild type and Atgpxl1-8 mutants were found usually under control conditions or after applying the lower concentrations of salt and osmotic stress (50 mM NaCl and 100 mM mannitol, respectively). This indicates a rather special role of AtGPXLs, since probably other ROS processing mechanisms are responsible for maintaining the vitality of these mutants under more severe stresses.

It is worth noting that the produced H₂O₂ can be reduced by catalase and other peroxidases as well. In the maintenance of ROS homeostasis, by the elimination of H₂O₂ and organic hydroperoxides, the function of Arabidopsis GPXLs shows some overlaps with peroxiredoxins and glutathione transferases (GSTs) with peroxidase activity [18,64]. Measuring the non-enzymatic antioxidants in the 6-week-old hydroponically grown Atgpxl mutants revealed altered ascorbic acid contents in shoots and increased GSH levels, particularly in the roots [57]. Moreover, it is also possible that the deficiency of one GPXL isoenzyme could be compensated by other GPXLs, as was found in the case of Arabidopsis tau group GSTs [65]. Even so, after two weeks of 100 mM NaCl or 200 mM mannitol treatments, the elevated ROS level and decreased metabolic activity in several mutant seedlings indicates their importance as antioxidant isoenzymes which have specific roles in ROS processing or redox signaling under control conditions and mild stress.

4.2. AtGPXLs Have Specific Functions in Tissues and Organs

Phylogenetic analysis of Arabidopsis GPXL protein sequences revealed that the 8 isoenzymes can be grouped to four pairs that show stronger similarity to each other than to the rest of the family. These are the AtGPX1 and AtGPX7, AtGPXL2 and AtGPXL3, AtGPX4 and AtGPX5, AtGPX6 and AtGPX8 [21,66,67]. Detailed analysis of their subcellular localization showed that AtGPXL1 and -7 are chloroplastic enzymes, AtGPXL2 and -8 exhibited cytosolic-nuclear localization, AtGPXL3 is a luminal protein that can be bound to the ER and Golgi membranes, AtGPXL4 and -5 are anchored to the plasma membrane, while AtGPXL6 can be found mainly in mitochondria [19]. Interestingly, their involvement in different growth or developmental processes was also reported and their expression depends on tissues and developmental stages [e.g., 18,21,23, and references therein]. Analysis of microarray data has shown that AtGPXL2, -3 and -8 were highly upregulated during germination, while AtGPXL1, -4, -5, -6 and -7 were downregulated. In the seedling stage AtGPXL3 showed the highest expression in the shoot apical meristem, but a rather high expression was reported in the case of AtGPXL1, -2 and -7 or -8 in the shoot tip, cotyledon and generally in the leaf tissues [21]. In the roots, the AtGPXL8 was expressed on the highest level in the radicle, but the transcript amounts of AtGPXL2, -3, -5 and -6 were also rather high both in the seedling stage and later on. However, their expression patterns depended on the root zones [21]. Focusing on the growth of roots, in line with the findings in the literature, the length of primary roots and/or number of lateral roots has changed significantly in the untreated Atgpxl2, -5, -6 and -8 mutants in the present experiments (Figure 2). Our correlation analysis revealed higher negative correlation between the ROS levels, and the numbers and growth of lateral roots compared to the found correlation in the wild type plants (Figure 7). However, this correlation was significantly disturbed in most Atgpxl mutants. Because auxin is the key hormone of the lateral root initiation, this may further confirm the relationship between GPXL proteins and auxin.

It was reported that H₂O₂ has a role in root hair formation and inhibition of root growth, while O₂^•- is necessary for root elongation [8,68]. The main sources of O₂^•- production are the electron transport systems and different enzymatic mechanisms, like the activity of NADPH oxidase, xanthene dehydrogenase and aldehyde oxidase [69]. Although the O₂^•- might convert to H₂O₂ even spontaneously at low pH, the role of SOD isoenzymes has high importance in the process, and they can be regarded as the first line of the antioxidant defence [70] or can act as a member of the ROS signaling route. Interestingly, it was suggested that GPXL2 can be part of a protein complex involved in O₂^•- conversion that can be activated under stress [21,71], but as far as we know, this theory has not been proven yet. Even if probably there are no direct connections between the GPXL enzymes and the O₂^•-, the O₂^•- concentrations were significantly higher in Atgpxl2 and in more than half of the mutants in our experiments than in the wild type plants. This difference was noticeable not just under control conditions, but in the case of certain mutants (Atgpxl5, -6), during stress as well.

4.3. AtGPXls Are Related to Redox Signaling

Increased ROS production temporarily shifts the redox balance toward a more oxidized state, which can be interpreted either as an active acclimation (adaptation) response or as damage [72]. Several reports indicated the involvement of GPXLs in redox signaling mechanisms, even more, their mediating role in the crosstalk between GSH and TRX pathways was suggested [2,32,34,54,55,73]. There is increasing evidence that these redox components are involved in healthy plant growth, development, successful organogenesis, and regeneration of cultured cells [55,64,74,75,76,77]. The Arabidopsis root meristemless1 (rml1) mutant, which has a defect in GSH biosynthesis [78], is unable to maintain the root apical meristem, however the apical meristem in its shoot is not affected, possibly because of the TRX-dependent control [79]. Characterization of the rml1 mutant discovered altered expression of several hundred genes, among them numerous encode redox-related proteins, such as glutaredoxins (GRXs), h-type thioredoxins (TRXhs) and GPXLs [34].

Thiol redox biochemistry is considered to have a fundamental role in cellular processes, even in the growth of cells. GSH can alter the expression of genes through modulation of redox state of proteins and transcription factors and is thought to be the main redox regulator, but this mechanism is yet poorly understood [80]. Glutathione limits the lifetime of the oxidative signals because of its antioxidant function, and the maintenance of a high GSH level and GSH/GSSG ratio are important in proper function of the cells and organs [81]. In general, the GSH-dependent redox system is regarded to be a house-keeping redox system, and E_GSH as a relative constant parameter in the plant cytosol [48]. Interestingly, measuring the GSH and GSSG amounts in 6-week-old hydroponically grown Atgpxl mutants revealed elevated GSH content, and a more reduced redox state compared to the wild type [57]. Our other results indicated that the key mechanism of AtGPXL5 might be the modification of the redox status [30]. Overexpression of AtGPXL5 increased the amount of total GSH, thus the glutathione redox potential became more negative than in the wild type. These changes were statistically significant under control conditions only in the shoots, but 24-h treatment with 100 mM NaCl caused significant differences in the roots as well. However, the elevated peroxide levels of Atgpxl5 mutant compared to wild type indicated that the AtGPXL5 has a role in fine-tuning the ROS homeostasis [30]. Among the 6-day-old untreated Atgpxl1-8 seedlings we found the highest (ca. 25 mV) changes in the E_GSH in the leaves of gpxl1 and -7 seedlings, demonstrating their importance in maintenance of the redox homeostasis.

It is widely accepted, that the environment-responsive developmental plasticity is linked to ROS, enzymatic and non-enzymatic antioxidants, and are in strict relationship with hormonal control in the development of the roots [82,83,84,85]. It was reported that the altered ROS levels and redox state influenced the growth of proximal meristem of 7-day-old Arabidopsis roots [49]. Using the roGFP1 redox sensor Jiang et al. [47] demonstrated that the redox profile of 3-9-day-old Arabidopsis roots shifted toward the more oxidized state due to 50-150 mM NaCl treatments at the early stage of stress, but it was re-established after 6 days. Moreover, these authors described the connection between the changes in EGSH and auxin transport using in vivo redox probes [47].

Changes in the redox profile of the root tissues, caused by salt treatments, influenced the patterns of both PIN1 and PIN2 auxin efflux carriers and the AUX1 influx carrier, therefore changing the root meristem size [47]. Since GPXLs are involved in forming the redox status of plants, they may also affect the hormone transport, metabolism and signaling as well. The relevance of AtGPXL7 in hormone-mediated root development, especially in lateral root development, was also demonstrated by 1-naphtaleneacetic acid and synthetic strigolactone treatments [32]. The interaction of auxin and ROS in plant growth and development has a role under salt stress conditions to trigger dynamic responses [86,87]. The elevated ROS production negatively modulates root growth by regulating the PIN-mediated auxin polar transport in the roots, and it affects root meristem activity [87]. Auxin may alter root growth and development through modulating the redox balance e.g., by regulation of APX activity [86,88], and enhanced redox metabolism was reported in the auxin receptor mutant tir1afb2 that showed enhanced tolerance to salt stress [89].

In our present study elevated ROS levels were found after both the NaCl and mannitol treatments especially in shoots, and the correlation analysis revealed higher negative correlation between ROS levels and the number and growth of lateral roots, instead of the growth parameters of the primary roots in the wild type plants. However, this correlation was significantly altered in the mutants (Figure 7). The relationship between the activation of GPXLs enzymes and the triggered hormonal changes needs further investigation. Our previous experiments uncovered the phenotypical effects of the mutation and overexpression of AtGPXL5 under NaCl stress, also the interaction between AtGPXLs and the ethylene signaling was demonstrated [90]. Our present correlation analysis indicates the primary role of ROS-dependent effects of AtGPXLs in altering the development of seedlings and in response to abiotic stresses. The isoenzyme-specific pattern highlighted that AtGPXLs having distinct interactions with other pathways are integrated members of signaling network.

5. Conclusions

AtGPXLs play a fine-tuning role in regulating ROS levels in seedlings under both control conditions and after two weeks of abiotic stress. Their significance under normal and mild stress conditions is likely linked to redox signaling. While other antioxidant mechanisms may be activated under severe stress, they are insufficient to fully suppress ROS levels. Consequently, mutations in an AtGPXL gene can alter the overall cellular redox status. Hormones such as auxin may also exert their effects in a redox-dependent manner, with AtGPXLs potentially contributing to this process. Our correlation analysis reinforces the role of GPXL enzymes in supporting plant growth and development by modulating ROS levels in an organ-specific manner. The observed differences in root phenotype and metabolic activity of Atgpxl mutant plants reflect the complex interplay between AtGPXLs, redox homeostasis, and plant hormones.