Over-Expression of Dehydroascorbate Reductase Improves Salt Tolerance, Environmental Adaptability and Productivity in Oryza sativa

Abiotic stress induces reactive oxygen species (ROS) generation in plants, and high ROS levels can cause partial or severe oxidative damage to cellular components that regulate the redox status. Here, we developed salt-tolerant transgenic rice plants that overexpressed the dehydroascorbate reductase gene (OsDHAR1) under the control of a stress-inducible sweet potato promoter (SWPA2). OsDHAR1-expressing transgenic plants exhibited improved environmental adaptability compared to wild-type plants, owing to enhanced ascorbate levels, redox homeostasis, photosynthetic ability, and membrane stability through cross-activation of ascorbate-glutathione cycle enzymes under paddy-field conditions, which enhanced various agronomic traits, including root development, panicle number, spikelet number per panicle, and total grain yield. dhar2-knockdown plants were susceptible to salt stress, and owing to poor seed maturation, exhibited reduced biomass (root growth) and grain yield under paddy field conditions. Microarray revealed that transgenic plants highly expressed genes associated with cell growth, plant growth, leaf senescence, root development, ROS and heavy metal detoxification systems, lipid metabolism, isoflavone and ascorbate recycling, and photosynthesis. We identified the genetic source of functional genomics-based molecular breeding in crop plants and provided new insights into the physiological processes underlying environmental adaptability, which will enable improvement of stress tolerance and crop species productivity in response to climate change.


Introduction
Rice is a major staple crop worldwide [1], and although rice yield is compromised by land erosion, degradation, and climate change, the demand for rice production is steadily increasing. Of late, the climate change caused by global warming has intensified. The Here, we used gene expression profiling and genetic network to analyze the expression patterns in OsDHAR1-overexpressing plants under conditions of environmental stress, and to identify the target genes of stress-related transcription factors and potential cis-acting DNA elements by combining the expression data. The profiles of plant responses to environmental stresses are expected to enable the identification of regulators that will be useful in biotechnological approaches aimed at improving stress tolerance, in addition to providing new insights into the physiological processes involved in the mechanisms underlying environmental adaptability.

Vector Construction and Agrobacterium-Mediated Rice Transformation
Full-length cDNA encoding Oryza sativa L. japonica cytosolic dehydroascorbate reductase (OsDHAR1) (accession no. AY074786) was cloned from rice leaves using reversetranscription polymerase chain reaction (RT-PCR). The primers, OsDHAR-FC and OsDHAR-RC were used for PCR (Table S1). The cDNA-designated as OsDHAR1-was subsequently cloned into the NcoI/KpnI site of the pBT1 vector, downstream of the SWPA2 promoter. The SWPA2::OsDHAR1::Tnos terminator was cloned into the HindIII site of pCAMBIA3300 ( Figure 1A) and subsequently designated as pOsDHAR1. The promoter and OsDHAR1 were sequenced to confirm that the open reading frame was intact, i.e., there were no frame shifts or nucleotide conversions. Finally, the pOsDHAR1 binary vector was introduced into Agrobacterium strain LBA 4404, which was then used to transform rice calli from the scutellum of the mature seeds of O. sativa L. japonica 'Ilmi variety,' as described previously [25,26]. underlying stress tolerance and responses of higher plants, but also for improving the stress tolerance of crops by gene manipulation.
Here, we used gene expression profiling and genetic network to analyze the expression patterns in OsDHAR1-overexpressing plants under conditions of environmental stress, and to identify the target genes of stress-related transcription factors and potential cis-acting DNA elements by combining the expression data. The profiles of plant responses to environmental stresses are expected to enable the identification of regulators that will be useful in biotechnological approaches aimed at improving stress tolerance, in addition to providing new insights into the physiological processes involved in the mechanisms underlying environmental adaptability.

Vector Construction and Agrobacterium-Mediated Rice Transformation
Full-length cDNA encoding Oryza sativa L. japonica cytosolic dehydroascorbate reductase (OsDHAR1) (accession no. AY074786) was cloned from rice leaves using reversetranscription polymerase chain reaction (RT-PCR). The primers, OsDHAR-FC and OsDHAR-RC were used for PCR (Table S1). The cDNA-designated as OsDHAR1-was subsequently cloned into the NcoI/KpnI site of the pBT1 vector, downstream of the SWPA2 promoter. The SWPA2::OsDHAR1::Tnos terminator was cloned into the HindIII site of pCAMBIA3300 ( Figure 1A) and subsequently designated as pOsDHAR1. The promoter and OsDHAR1 were sequenced to confirm that the open reading frame was intact, i.e., there were no frame shifts or nucleotide conversions. Finally, the pOsDHAR1 binary vector was introduced into Agrobacterium strain LBA 4404, which was then used to transform rice calli from the scutellum of the mature seeds of O. sativa L. japonica 'Ilmi variety,' as described previously [25,26]. stress-inducible promoter (SWPA2p), the OsDHAR1 coding regions, and the Tnos terminator, which was derived from the nopaline synthase gene. The bar gene, which encodes the herbicide bialaphos (Basta) and whose expression is driven by the 35S promoter, was used as a selection marker. (B) Selection of salt-tolerant homozygous transgenic lines. Three independent transgenic lines (TR1, TR2, and TR3) were identified upon screening 4-week-old seedlings using 100 mM NaCl. (C) Visualization of flanking T-DNA sequencing amplicon. Electrophoresis was performed to determine the copy number and insertion position of the integrated OsDHAR1 transgene. Flanking T-DNA sequencing was performed as described in SI Materials and Methods. The PCR amplicon digested with the indicated restriction enzymes is represented and amplicons marked with arrows were sequenced (date not shown). (D) Phenotypes of transgenic and WT seedlings under conditions of salt stress. Upper panel, 4-week-old WT and transgenic plants in normal state; lower panel, 4-week-old WT and transgenic plants recovered after exposure to 100 NaCl for 25 days. Phenotype analysis was performed at least three times and represented as a representative image. (E) Confirmation of OsDHAR1 expression using semi-quantitative RT-PCR and (F) western blotting. Signal intensity of WT plants was due to endogenous DHAR1 expression. tubulin (Tub) and corresponding protein were used as loading controls. (G) Total fresh weight of recovered plants after exposure to 100 NaCl for 25 days. WT, wild-type rice; TR, transgenic rice.

Plant Materials and Growing Conditions
The Ilmi variety of O. sativa L. japonica was used as the host for the homologous overexpression of OsDHAR1. Approximately 20 independent T 0 transgenic rice plants were screened using Basta. Among the transformants, 13 independent T 1 transgenic plants grown from T 1 seeds were selected through Basta and PCR genotyping. Four independents homozygous T 2 plants were identified using the same process, and one independent transgenic line overlapped with another independent transgenic line. Three independent transgenic lines (TR1, TR2, and TR3) were selected for subsequent experiments. For genotyping, both transgenic and wild-type (WT) plants were cultivated in natural paddy fields located at the Gunwi campus (36 • 24 N, 128 • 53 E, 95 m a.s.l.) of Kyungpook National University, Korea, from June to October. Seeds were germinated at 32 • C in a controlled growth chamber, transplanted into soil pots (12.5 cm diameter), and maintained at 28-32 • C under a 16-h light/8-h dark cycle. For the salt stress assay, 4-week-old seedlings were treated with 100 mM NaCl for 25 days. Leaves were harvested when seedlings were exposed to 100 mM NaCl for 4 days or grown for ∼4 weeks after transplantation in paddy fields. The harvested leaves were immediately soaked in liquid nitrogen and stored at −70 • C until subsequent use. In addition, a T-DNA insertion mutant line (3A-03259), a putative dehydroascorbate reductase gene (OsDHAR2, Os06g0232600) knockdown mutant, was provided by the Salk Institute Genomic Analysis Laboratory (http://signal.salk.edu/ cgi-bin/RiceGE, 18 November 2019). Homozygous lines were identified as described for the DHAR1-overexpressing transgenic plants.

Copy Number and Location of the OsDHAR1 Transgene
Determination of the copy number and location of the OsDHAR1 transgene consisted of two steps. First, genomic DNA was digested with restriction enzymes and ligated to an adaptor, and second, the fragments were subjected to PCR amplification using two primer sets that are specific to the OsDHAR1 transgene and adaptor regions. In the present study, genomic DNA (500 ng) was digested with 2 U of MspI and HincII and ligated using 5 U of T4 DNA ligase (Takara Bio) in 20-micoliter reactions; the reaction mixtures contained T4 DNA ligase buffer and 50 pmol of adaptors and were incubated for 1 h at 37 • C. The first PCR was conducted in 20-microliter reactions that contained 0.5 pmol of each primer (ADA1 with either LB1 or RB1, and ADA2 with either LB2 or RB2) (Table S1) and 1 µL digestion/ligation product, with the following conditions: an initial denaturation step at 95 • C for 5 min, 20 cycles of 30 s at 94 • C and 1 min at 72 • C, and a final elongation step at 72 • C for 10 min. A second PCR was conducted with 5 µL of the first PCR product under the following conditions: an initial denaturation step at 94 • C for 5 min; 40 cycles of 30 s at Antioxidants 2022, 11, 1077 5 of 20 94 • C, 30 s at 60 • C, and 1 min at 72 • C; and a final elongation step at 72 • C for 10 min. The amplified products were separated on a 1% agarose gel, purified using a HiYield Gel/PCR DNA Extraction Kit (Real Biotech, Taipei, Taiwan), and sequenced using LB2 and RB2 primers (Table S1).

Measurement of the Ascorbate Content
The total AsA (tAsA), AsA, and DHA contents were spectrophotometrically determined as described previously [27]. Briefly, ∼0.2 g of frozen leaf sample was ground with a mortar and pestle and vortexed vigorously for 10 min after adding 2 mL of 5% (v/v) m-phosphoric acid. The homogenate was centrifuged at 12,000× g for 20 min and the clear supernatant was collected. The tAsA content was determined in a reaction mixture containing 100 µL crude extract, 500 µL 150 mM KH 2 PO 4 buffer (pH 7.4; 5 mM EDTA), and 100 µL 10 mM dithiothreitol, which reduced DHA to AsA. The reaction mixture was incubated at 25 • C for 10 min and supplemented with 100 µL of 0.5% (w/v) N-ethylmaleimide (NEM) to remove excess dithiothreitol. The AsA content was assayed by replacing dithiothreitol and NEM with 200 µL of deionized H 2 O as described previously. Color was developed in both reaction mixtures by adding 400 µL of 10% (w/v) trichloroacetic acid, 400 µL of 44% (v/v) o-phosphoric acid, 400 µL of α,α-dipyridyl prepared in 70% (v/v) ethanol, and 200 µL of 30% FeCl 3 . The reaction mixture was incubated at 40 • C for 1 h and the absorbance was quantified at 525 nm. The DHA concentration was estimated by subtracting the AsA content from the tAsA content.

Enzyme Activity Assay
To prepare the crude protein extracts for the MDHAR and DHAR assays, leaf samples (0.2 g fresh weight) were homogenized at 4 • C in 1 mL of extraction buffer (50 mM Tris-HCl (pH 7.5), 3 mM MgCl 2 , 1 mM EDTA, 1 mM PMSF, and protease inhibitor cocktail (Sigma-Aldrich, St. Louis, MI, USA)) with a mortar and pestle. The homogenate was centrifuged at 15,000× g for 20 min at 4 • C and the clear supernatant was used for enzymatic assays. To prepare crude protein extracts for APX and GSH reductase (GR) assays, leaf samples (0.2 g) were extracted in 1 mL of extraction buffer (mM potassium phosphate (pH 7.0), 7 mM β-ME, and EDTA-free protease inhibitor cocktail (Roche, Basel, Switzerland)). In addition, crude protein extract for the APX assay was prepared in the presence of 2 mM AsA. All enzyme assays were performed by monitoring the corresponding reaction mixtures at 25 • C using a UV-1800 UV-VIS spectrophotometer (Shimadzu, Kyoto, Japan). MDHAR activity was assayed in 1 mL reaction mixture containing 50 mM potassium phosphate (pH 7.2), 0.2 mM NADH, 2 mM AsA, 1 U AsA oxidase, and crude extract. The absorbance was measured at 340 nm, the activity was calculated using an absorbance coefficient of 6.2 mM −1 cm −1 ; and one unit of MDHAR activity was defined as the amount of enzyme that oxidized 1 nmol of NADH per min at 25 • C [28]. DHAR activity was assayed in 1 mL reaction mixture containing 50 mM potassium phosphate (pH 7.0), 0.1 mM EDTA, 0.5 mM DHA, 2.5 mM GSH, and crude extract. The absorbance was measured at 265 nm, the activity was calculated using an absorbance coefficient of 14.6 mM −1 cm −1 [29], and one unit of DHAR activity was defined as the amount of enzyme that produced 1 nmol of AsA per min at 25 • C. GR activity was assayed in 1 mL reaction mixture containing 50 mM potassium phosphate (pH 7.6), 1.7 mM EDTA, 1.0 mM GSSG, 0.1 mM NADPH, and crude extract. Absorbance was measured at 340 nm [30]. APX activity was assayed in 1 mL reaction mixtures containing 50 mM sodium phosphate (pH 7.0), 0.2 mM EDTA, 0.5 mM AsA, 0.2 mM EDTA, 2.5 mM H 2 O 2 , and crude extract (50 µg equivalent). Enzyme activity was recorded as the decrease in absorbance at 290 nm for 1 min, and the amount of oxidized AsA was calculated using an extinction coefficient of 2.8 mM −1 cm −1 [31]. All enzymatic activities in the transgenic plants were calculated relative to those in the WT plants grown under normal conditions; enzyme activities observed in the WT plants were set as 100%.

Redox State Analysis
Leaf samples (0.2 g fresh weight) were homogenized at 4 • C in 1 mL extraction buffer (20 mM sodium phosphate (pH 7.2), 1 mM PMSF, and EDTA-free protease inhibitor cocktail (Roche)), using a mortar and pestle. The homogenate was centrifuged at 15,000× g for 20 min at 4 • C and the clear supernatant was collected. Intracellular hydroperoxide levels were determined based on ferrous ion oxidation in the presence of xylenol orange, a ferric ion indicator. Briefly, crude extract (50 µL) was added to 950 µL of FOX reagent (100 µM xylenol orange, 250 µM ammonium ferrous sulfate, 100 mM sorbitol, and 25 mM sulfuric acid). The mixture was incubated at 25 • C for 30 min and then centrifuged to remove the flocculated material before measuring the absorbance at 560 nm [32]. Lipid peroxidation was spectrophotometrically determined by measuring malondialdehyde or thiobarbituric acid-reactive substances, as reported previously [33]. Briefly, 1.0 mL crude extract was mixed with 2.0 mL of TCA-thiobarbituric acid-HCl reagent (15% thiobarbituric acid, 0.375% TCA, and 0.25 N HCl), mixed thoroughly, and heated in a boiling water bath (95 • C) for 15 min. After cooling, the supernatant was centrifuged at 12,000× g for 10 min, and the absorbance of the sample was measured at 535 nm and compared with that of a control sample. Malondialdehyde concentration in the sample was calculated using an extinction coefficient of 1.56 × 10 5 M −1 cm −1 [34]. For 3,3 -diaminobenzidine staining, the leaves of whole rice plants were sprayed with 50 µM methyl viologen solutions containing 0.1% Tween-20 and incubated for 12 h, after which the leaves were cut, placed in 3,3 -diaminobenzidine (1 mg mL −1 )-HCl (pH 3.8), incubated for 8 h at 28 • C, and then cleared in boiling ethanol (96%) for 10 min before photographing [35].

Ion Leakage and Chlorophyll Fluorescence Analysis under Conditions of MV Treatment
Five leaves were collected from the TR and WT rice plants grown for 4 weeks in the natural paddy field, after which ten leaf discs were cut (approximately 1 cm). The leaf discs were then incubated in water supplemented with 10 µM methyl viologen (MV) in the dark at 25 • C for 12 h to allow the MV to diffuse into the leaf discs. After pre-incubation, leaf discs were placed in a growth chamber (16 h light: 8 h dark) at 25 • C for 120 h. The conductivity of the suspension was measured using a conductance meter (455C; Isteck Co., Ltd., Seoul, Korea) before and after autoclaving at 121 • C for 15 min to release all electrolytes. Relative ion leakage was expressed as a percentage of total conductivity [36]. MV-treated leaves of TR and WT plants were dark-adapted for 15 min. Chlorophyll fluorescence was measured based on the photosynthetic activity of photochemical yield (Fv/Fm), which represents the maximum yield of the photochemical reaction of photosystem II at room temperature, using a Handy Fluor Cam FC 1000-H (Photon Systems Instruments Ltd., Brno, Czech Republic) [37]. Fv/Fm values were measured as described previously [38].

Microarray-Based Gene Expression Profiling
Leaf samples (pooled from five plants per family) were harvested from WT and TR1 plants that had been grown for ∼30 days under natural paddy field conditions. TR1 plants were used as representative transgenic lines. Total RNA was extracted from rice leaves using the TRIzol reagent (Sigma-Aldrich) according to the manufacturer's instructions. Following homogenization, 1 mL of the solution was transferred to a 1.5 mL Eppendorf tube and centrifuged at 12,000× g for 10 min at 4 • C to remove any insoluble material. The supernatant containing the RNA was collected, mixed with 0.2 mL of chloroform, and centrifuged at 12,000× g for 15 min at 4 • C. The upper phase containing the RNA, was transferred to a new tube, and RNA was precipitated by mixing the solution with isopropyl alcohol (0.5 mL) and recovered by centrifuging the sample at 12,000× g for 10 min at 4 • C. The RNA pellet was washed briefly with 1 mL of 75% ethanol by centrifuging at 7500× g for 5 min at 4 • C. Finally, the RNA pellet was dissolved in DNase-free water and RNA quality and quantity were assessed using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Gene expression was analyzed using GeneChip Rice Genome Arrays (4 × 44 K; Affymetrix), which (each) comprised over 54,000 probe sets, and Antioxidants 2022, 11, 1077 7 of 20 11 pairs of oligonucleotide probes were synthesized in situ on the arrays. Biotinylated cRNA was prepared from 100 ng of total RNA, according to the standard Affymetrix protocol. Following fragmentation, 15 µg of RNA was hybridized for 16 h at 45 • C using a GeneChip Rice Genome Array. GeneChips were washed and stained in the Affymetrix Fluidics Station 450 and scanned using the Affymetrix GeneChip Scanner 3000 7G. The data were analyzed using robust multi-array analysis-with default Affymetrix settings-together with the global scaling normalization method. The trimmed mean target intensity of each array was arbitrarily set at 100, and both the normalized and log-transformed intensity values were analyzed using GeneSpring GX 13.0 (Agilent Technologies). The fold-change filters used defined the up-and downregulated genes as those that are expressed at levels that were at least 200% less than 50% of the expression of the corresponding genes in the controls. Hierarchical clustering was performed using GeneSpring GX 13.0 (Agilent Technologies), with the clustering algorithm being the Euclidean distance and average linkage. The genetic network of the differentially expressed genes was framed using String (http://apps.cytoscape.org/apps/stringapp, 18 December 2020) and ClueGo (http: //apps.cytoscape.org/apps/cluego, 18 December 2020).

Estimation of Agronomic Traits in Paddy Fields
To examine the yield parameters of the transgenic rice plants, three independent T 3 or T 4 homozygous lines and their respective controls (WT plants) were transplanted into natural paddy fields located at the Kyungpook National University campus in 2020 (T 3 ) and 2021 (T 4 ). A completely randomized block design was employed with two replicates, and each plot comprised 12 plants, with two rows for each line, 0.2 m between plants, and 0.3 m between plots. When the plants reached maturity and the grains ripened, yield components were scored by digging up 8 rice plants per plot to a radius of 15 cm. After washing the plant roots, total plant weight (g), stem weight (g), and root weight were measured and number of panicles per hill and number of spikelets per panicle were counted. The filling rate (%) of the rice plants was calculated relative to 100% of WT rice plants, which were defined had expressed the degree of the filling rate as panicles and spikelets under field conditions. To measure the total grain weight (g), and 1000-grain weight (g), whole seeds per plant were collected and weighed, and then 1000 seeds were counted and weighed additionally. Unfilled and filled grains were separated, randomly counted, and weighed. The following agronomic traits were scored: total plant weight (g), culm weight (g), root weight (g), number of panicles per hill, number of spikelets per panicle, filling rate (%), total grain weight (g), and 1000-grain weight (g).

Data Analysis
All experiments were performed with at least three independent repetitions, and the results are expressed as mean ± standard deviation (SD). All biochemical results of the transgenic plants were calculated relative to those of WT plants grown under normal conditions, which were defined as 100%.

Development of SWPA2::OsDHAR1 Transgenic Rice Plants
To assess the stress tolerance of OsDHAR1-overexpressing plants, we produced 20 lines (T 0 generation) of transgenic plants carrying the SWPA2::OsDHAR1 transgene construct. Fifteen seeds (T 1 seeds) per primary transgenic line were cultivated in a paddy field and their progenies (T 2 seeds) were harvested.
In the present study, three independent salt-tolerant homozygous transgenic families (TR1, TR2, and TR3) were identified ( Figure 1B). In each of these families, only a single SWPA2::OsDHAR1 transgene copy was inserted into the rice genome ( Figure 1C; Table S2). In addition, all salt-tolerant plants harbored the SWPA2::OsDHAR1 transgene, which was stably inherited by the following generation. To investigate the mechanism by which OsDHAR1 overexpression increases salt tolerance, the phenotypes and gene expression of the transgenic plants were examined under salt stress. There was no difference in the phenotypes of the transgenic and WT plants grown under normal conditions (0 day; Figure 1D). Under conditions of salt stress, most of the WT plants withered after 25 days, but the transgenic plants were relatively unaffected (25 days; Figure 1D). The survival rate of the transgenic plants was higher than that of the WT plants. In addition, the transgenic and WT plants exhibited clear phenotypic differences with respect to their ability to recover from salt damage. More specifically, the transgenic plants recuperated quickly after exposure to salt stress, whereas the WT plants did not ( Figure 1D), and the fresh weight of the transgenic plants was approximately 2-fold greater than that of the WT plants ( Figure 1G). We constructed dhar2-knockdown transgenic lines (dhar2-1 and dhar2-2) and further analyzed the response of the mutant to salt stress. As a result, the dhar2-knockdown mutants showed hypersensitivity to salt stress ( Figure S3D). Homologous overexpression of OsDHAR1 was confirmed by measuring the levels of OsDHAR1 mRNA and protein in 4-week-old seedlings exposed to 100 mM NaCl for 4 days. The analysis revealed that OsDHAR1 was expressed at greater levels in the transgenic plants than in the WT plants ( Figure 1E), and western blotting confirmed that the salt-stressed transgenic plants expressed higher levels of the OsDHAR1 protein ( Figure 1F). On the other hand, the dhar2-knockdown mutants showed weak DHAR expression compared to in WT plants ( Figure S3C). Thus, our results confirmed that the expression of OsDHAR1 is related to stress, and the OsDHAR1 transgene is effectively translated into a functional DHAR protein, as expected.

Performance of OsDHAR1-Expressing Transgenic Plants under Salt Stress
We investigated, at least partially, the OsDHAR1-mediated mechanism underlying the enhanced salt tolerance of transgenic plants by challenging healthy plants with salt stress and then analyzing the biochemistry of the leaf samples. In the following experiments, TR1 and TR2 plants, which were better than TR3 plants in Figure 1, were used as representatives of three independent transgenic plants. The activity of MDHAR, APX, and GR, which are involved in the AsA-GSH cycle, was greatly increased by salt stress in the transgenic plants (TR1 and TR2) (compared to WT plants), even though the WT plants themselves exhibited a moderate increase in enzyme activity. Further, APX and GR exhibited the highest and lowest activities, respectively, in transgenic plants ( Figure 2A). Under normal conditions, the enzyme activity of the transgenic plants was slightly elevated, possibly because of transgene insertion.
To evaluate the metabolic effect of OsDHAR1 overexpression, the AsA pool of leaves from transgenic (TR1 and TR2) and WT plants was compared. Although salt stress increased the AsA levels in both transgenic and WT plants, the AsA levels of the transgenic plants were higher than those of the WT plants. Meanwhile, the DHA level of the transgenic plants was distinctly decreased compared to that of the WT plants, which subsequently resulted in an increased AsA/DHA ratio. Indeed, under conditions of salt stress, the AsA/DHA ratio of the transgenic plants was ∼1.67-fold higher than that of the WT plants, whereas, only a slight difference was observed under normal conditions ( Figure 2B). These results indicated that the DHAR activity of the OsDHAR1-expressing transgenic plants is increased, and the dhar2 mutants, in contrast, showed reduced DHAR activity under salt stress conditions ( Figure S3E). Based on these results, hydroperoxide and malondialdehyde levels were measured to determine whether the increased AsA pool of the transgenic plants affected the redox state. Under salt stress, the hydroperoxide and malondialdehyde contents of the transgenic plants were approximately 1.5-fold lower than those of the WT plants ( Figure 2C), and there was no difference under normal conditions. Ion leakage of the transgenic plants was also lower than that of the WT plants under conditions of salt stress ( Figure 2D). Therefore, our results indicate that OsDHAR1 overexpression improves redox homeostasis, the AsA pool, and membrane permeability in transgenic rice plants by cross-activating antioxidant enzymes, including MDHAR, APX, GR, and DHAR. To evaluate the metabolic effect of OsDHAR1 overexpression, the AsA pool of lea from transgenic (TR1 and TR2) and WT plants was compared. Although salt stress creased the AsA levels in both transgenic and WT plants, the AsA levels of the transge plants were higher than those of the WT plants. Meanwhile, the DHA level of the tra genic plants was distinctly decreased compared to that of the WT plants, which sub quently resulted in an increased AsA/DHA ratio. Indeed, under conditions of salt str the AsA/DHA ratio of the transgenic plants was ~1.67-fold higher than that of the W plants, whereas, only a slight difference was observed under normal conditions ( Fig  2B). These results indicated that the DHAR activity of the OsDHAR1-expressing tra genic plants is increased, and the dhar2 mutants, in contrast, showed reduced DHAR tivity under salt stress conditions ( Figure S3E). Based on these results, hydroperoxide a malondialdehyde levels were measured to determine whether the increased AsA poo the transgenic plants affected the redox state. Under salt stress, the hydroperoxide a malondialdehyde contents of the transgenic plants were approximately 1.5-fold low than those of the WT plants ( Figure 2C), and there was no difference under normal c ditions. Ion leakage of the transgenic plants was also lower than that of the WT pla under conditions of salt stress ( Figure 2D). Therefore, our results indicate that OsDHA overexpression improves redox homeostasis, the AsA pool, and membrane permeabi in transgenic rice plants by cross-activating antioxidant enzymes, including MDHA APX, GR, and DHAR.

Agronomic Traits of OsDHAR1-Expressing Transgenic Plants
We postulated that OsDHAR1 overexpression could protect transgenic plants from other types of environmental stresses because the OsDHAR1-transgenic plants exhibited greater salt tolerance owing to their improved AsA pool and redox state. To test this hypothesis, we transplanted three independent T 3 or T 4 homozygous transgenic (TR1-TR3) plants and WT plants into paddy fields using a completely randomized design. The plants were grown to the seed maturation stage, and agronomic parameters were scored from the vegetative to the reproductive stage. Transgenic plants performed better than the WT plants during both the vegetative and post-flowering harvest stages (Figurse 3B and S1C). These observations prompted us to examine other yield components. In 2020 and 2021, the total plant weight, culm weight, root weight, number of panicles per hill, number of spikelets per panicle, filling rate, total grain weight, and 1000-grain weight were 20.7 and 19.6%, 17.4 and 17.9%, 12.7 and 12.3%, 17.2 and 20.3%, 14.9 and 16.8%, 11.0 and 9.6%, 18.2 and 20.4%, and 3.7 and 5.0% greater, respectively, in transgenic plants than in WT plants ( Figure 3C). The rice grain yield factors TGW of the transgenic plants was higher than that of the WT plants and appeared to result from an increased number of spikelets per panicle and 1000-grain weight of the transgenic plants. Therefore, the number of panicles per hill in transgenic plants increased by ∼18.7% compared to that in WT plants, and as a result, the total grain weight of transgenic plants increased by ∼19.3%. There was also no difference in the agronomic traits of the transgenic plants grown over the two years, despite obvious differences in the mean temperature, precipitation, humidity, and total sunshine ( Figure 3A). On the contrary, the field-grown dhar2 mutant plants exhibited reduced agronomic traits ( Figure S4). More specifically, the total plant weight, culm weight, root weight, number of panicles per hill, number of spikelets per panicle, filling rate, total grain weight, and thousand seeds weight of the dhar2 mutant plants was ∼48.0, 30.8, 71.6, 19.3, 17.5, 35.0, 51.3, and 44.6% lower, respectively, than those of the WT plants ( Figure S4C). Among these traits, the reduced total biomass (total plant weight and root weight) and poor seed maturation of dhar2 mutant plants distinctly diminished rice grain yield (total grain weight and hundred seeds weight). These results indicated that OsDHAR1 overexpression significantly enhanced the rice grain yield and total biomass in response to various environmental stresses under paddy-field conditions. of the WT plants and appeared to result from an increased number of spikelets per panicl and 1000-grain weight of the transgenic plants. Therefore, the number of panicles per hi in transgenic plants increased by ~18.7% compared to that in WT plants, and as a resul the total grain weight of transgenic plants increased by ~19.3%. There was also no diffe ence in the agronomic traits of the transgenic plants grown over the two years, despit obvious differences in the mean temperature, precipitation, humidity, and total sunshin ( Figure 3A). On the contrary, the field-grown dhar2 mutant plants exhibited reduced ag ronomic traits ( Figure S4). More specifically, the total plant weight, culm weight, roo weight, number of panicles per hill, number of spikelets per panicle, filling rate, total grai weight, and thousand seeds weight of the dhar2 mutant plants was ~48.0, 30.8, 71.6, 19.3 17.5, 35.0, 51.3, and 44.6% lower, respectively, than those of the WT plants ( Figure S4C Among these traits, the reduced total biomass (total plant weight and root weight) an poor seed maturation of dhar2 mutant plants distinctly diminished rice grain yield (tota grain weight and hundred seeds weight). These results indicated that OsDHAR1 overex pression significantly enhanced the rice grain yield and total biomass in response to var ous environmental stresses under paddy-field conditions.

Effect of OsDHAR1 Overexpression on Overall Gene Expression
As OsDHAR1-expressing transgenic plants exhibited increased grain yield and biomass under paddy-field conditions, we performed gene expression profiling to identify the mechanisms underlying OsDHAR1-mediated environmental tolerance. TR1 plants were used because in this line, the SWPA2::OsDHAR1 transgene was inserted into an intergenic region. We found that ∼404 genes were upregulated, and ∼525 genes were downregulated (Figure 4; Tables S3 and S5). The upregulated genes were related to amino acid (glycine, serine, and threonine) metabolism, glycosphingolipid biosynthesis, porphyrin and chlorophyll metabolism, cell growth, transport, and translation; genes related to lipid metabolism were highly expressed (Table S4). In contrast (Table S6), the downregulated genes were related to nitrogen metabolism, amino acid (phenylalanine, tyrosine, and tryptophan) biosynthesis, organic metabolism, organic biosynthesis, cellular metabolism, and cellular biosynthesis. regulated (Figure 4; Tables S3 and S5). The upregulated genes were related to amino acid (glycine, serine, and threonine) metabolism, glycosphingolipid biosynthesis, porphyrin and chlorophyll metabolism, cell growth, transport, and translation; genes related to lipid metabolism were highly expressed (Table S4). In contrast (Table S6), the downregulated genes were related to nitrogen metabolism, amino acid (phenylalanine, tyrosine, and tryptophan) biosynthesis, organic metabolism, organic biosynthesis, cellular metabolism, and cellular biosynthesis. There was also a distinct difference in the expression of genes involved in biological regulation and biosynthetic processes in WT and TR1 plants. Briefly, the genes, Os03g0115800, Os11g0206000, and Os03g0115800 encoding conserved hypothetical proteins were most highly expressed, the expression of OsDHAR1 was increased by 2.4-fold, and the expression of the transcription factors HSfA2a, HSfA2d, TAGL 12, AT103, WRKY10, and Os_F0727 was upregulated in the TR1 plants. The following were the upregulated gene categories: pathogenesis-related protein, planthopper resistance, isoflavone reductase homolog IRL, heat shock proteins and molecular chaperonins, root-specific protein (RCc2), plant secretory proteins, rapid alkalization factor 2, tubulin (OsTUB6), actin depolymerizing factor (OsADF3), annexin p35, photosystems, phytochrome-interacting protein, molybdopterin biosynthesis CNX2 protein, sodium/hydrogen exchanger protein, abiotic stress responsive protein (OsASR2), heavy metal transport/detoxification protein, heat-and acid-stable phosphoprotein, drought-induced protein, ABA/WDS induced protein, multidrug-resistance associated protein, jasmonate ZIM-domain protein (OsJAZ2), methyl chloride transferase which is harmless to ozone layer 1 (OsHOL1), aquaporin, glutaredoxins (OsGRX4 and OsGRX29), peroxidase (Prx59), MDHAR (OsMDHAR1), AsA peroxidases (OsAPX1 and OsAPX8), GR1, and AsA biosynthesis (L-galactono-1,4-lactone dehydrogenase) ( Figure S5; Table S3). The highly expressed genes were networked into categories that corresponded to the regulation of signal transduction, glycosphingolipid biosynthesis, hydrolase activity, protein complex subunit organization, zinc ion transmembrane transporter activity, chlorophyll envelope, response to ROS, water-soluble vitamin biosynthesis, redox homeostasis, and isoprenoid biosynthesis (Figures 5 and 6; Table S4).   Supplementary Table S4.

Discussion
High salinity, which inhibits plant growth and photosynthesis in crop plants, is one of the key factors that are most detrimental to rice productivity, especially as soil salinization increases every year owing to anthropgenic activities and climate change. Approximately 6.5% of the world's total area and ∼20% of the irrigated lands have been severely salinized [39,40]. Therefore, there is a need to improve salt tolerance in crop plants.
In the present study, OsDHAR1-overexpressing transgenic rice plants exhibited improved salt tolerance by enhancing their AsA pool, redox homeostasis, membrane stability, and ion homeostasis following ion leakage. Transgenic plants achieve this through co-activation of enzymes (MDHAR, APX, and GR) involved in the AsA-GSH cycle (Figures 1 and 2). The physiological and molecular mechanisms of salt tolerance in various plant species have been widely reported [41], and 'omics'-based plant molecular breeding has been widely employed to increase crop tolerance to salt stress [39,40]. Rice plants also exhibit various adaptive mechanisms to overcome salt damage, including regulation of seedling vigor, and efficient transport of salt ions through the root system [39]. For example, the specific overexpression of ion transporters and pump-responsive genes that encode Na + /H + antiporter (SOS1, SOD2, nhaA, and NHX), Na + transporter (HKT1), Na + /K + antiporter (HKT2), Na + ATPase (ENA), and vacuolar pyrophosphatase (vacuolar H + -PPase) has been shown to improve salt tolerance in a variety of transgenic crop plants by increasing biomass (shoot and root) production, germination rate, and proline content by enhancing ion homeostasis [41]. Taken together, our findings suggest that OsDHAR1 plays an important role in salt tolerance by improving the AsA pool and redox and ion homeostasis, indicating that DHAR expression in WT plants is insufficient to overcome salt stress.
The transgenic plants exhibited a better phenotype than the WT plants under salt stress (Figure 1), which was due to the improvement in redox homeostasis and the AsA pool through the cross-activation of MDHAR, APX, and GR ( Figure 2). We hypothesized that high salinity induces ROS-induced oxidative stress. Salt-induced oxidative stress has been recognized as an important source of cellular damage in plants [42], and ROS scavenging depends on both enzymatic and non-enzymatic systems. Among these systems, the most important non-enzymatic antioxidant is AsA, which functions as a reductant for peroxidases [43,44]. In the present study, the greater AsA levels in the transgenic plants resulted from increased MDHAR and DHAR activities (Figure 2A,B). Using AsA as an electron donor, APX catalyzes the reduction of H 2 O 2 , which enhances redox homeostasis and proteostasis by minimizing oxidative damage through the buffering of stress-induced ROS [44]. During exposure to salt stress, GR activity also differed between transgenic and WT plants. The enzyme activity was ∼1.3-fold higher in the transgenic plants than in the WT plants (Figure 2A), indicating that MDHAR and GR were more efficient in converting oxidized AsA (MDHA or DHA) to reduce AsA.
In contrast, the dhar2-knockdown transgenic lines, dhar2-1 and dhar2-2, which exhibited lower DHAR activity, were sensitive to salt stress ( Figure S3). According to previous reports, homologous overexpression of DHAR (CrDHAR1; Cre10.g456750) enhances the growth and survival of Chlamydomonas reinhardtii through the improvement of redox state, photosynthesis efficiency, GSH and AsA pools, and GR and APX activities, which subsequently improve resistance to ROS-induced oxidative stress [45]. Previous studies have shown that increasing AsA concentrations through the overexpression of AsA biosynthesis genes enhances stress tolerance in transgenic tobacco, Arabidopsis, and lettuce [46], and the overexpression of ROS detoxification-responsive genes, namely APX, GST, SOD, MDHAR, and CAT, increases salt tolerance in transgenic tobacco and rice by augmenting photosynthetic efficiency and plant growth [41]. In contrast, several Arabidopsis DHAR-knockout mutants (dhar1, 2, and 3) that exhibited decreased tAsA and increased DHA were more sensitive to bright light and oxidative stress than the WT plants. The mutants possessed lower levels of reduced AsA, MDHAR, and DHAR activities and higher levels of hydroperoxide, which indicates that AsA can limit ROS scavenging under abiotic stress [44,47]. CrDHAR1-knockdown microRNA lines with reduced DHAR expression and AsA recycling ability were sensitive to light and oxidative stress [45].
The salt-tolerant barnyard grass Echinochloa crus-galli L. depends on higher activities of guaiacol peroxidase, APX, and GR than WT plants when exposed to 200 mM NaCl, whereas salt-sensitive plants with reduced enzyme activities exhibited slower growth rates, higher lipid peroxidation rates, and higher levels of ROS in their leaves [42]. Moreover, in saltstressed Vicia faba seedlings, exogenous supplementation of AsA was reported to alleviate oxidative stress [48,49], which indicates that steady-state AsA content was insufficient for effective antioxidative defense under salt stress and that its biosynthesis and recycling were necessary [42]. Taken together, the results of the present study suggest that enhanced AsA recycling via elevated DHAR activity and the combination of enzymes involved in the AsA-GSH pathway increases salt tolerance and that OsDHAR1 can be used to facilitate the development of improved rice cultivars with enhanced AsA content.
During the vegetative to reproductive stages of growth and development, rice plants grow in paddy fields under a variety of environmental conditions, including temperature, drought, flooding, ozone, UV radiation, nutrient depletion, air pollution, and photoperiod, which can decrease productivity [50]. To overcome the factors that limit rice productivity, researchers have concentrated on improving the yield of superior phenotypes using upand downstream regulatory genes [51,52]. Among these genes, the expression levels of OsDHAR1, which are a downstream regulatory gene, are higher in the three independent homologous TR rice lines than in WT rice plants in paddy fields at transcriptional and translational levels, as a result, OsDHAR1 overexpressing TR plants have been shown to affect environmental adaptability and productivity in rice grown in natural paddy fields, which are frequently exposed to abiotic and biotic stressors ( Figure S1). In the present study, transgenic rice plants overexpressing OsDHAR1 exhibited better phenotypes than WT plants in both the vegetative ( Figure S1A,C) and reproductive ( Figure 3B) stages. The major physiological differences resulted from the improvement of redox state, photosynthetic ability, and membrane permeability, which was achieved by buffering ROS-induced oxidative damage through the cross-activation of MDHAR, APX, GR, and DHAR in transgenic plants ( Figures S1 and S2). The increased tAsA in the transgenic plants resulted from the induction of the AsA biosynthesis gene GalDH. Augmented MDHAR and DHAR activities increased the AsA/DHA ratio by elevating and reducing AsA and DHA levels, respectively, in transgenic plants, which then enhanced the AsA pool via AsA recycling ( Figure S1E,F). Plants that are susceptible to abiotic and biotic stresses are generally depleted of AsA. In contrast, stress-resistant plants exhibit increased AsA levels and a steady-state or decreased oxidation status, depending on whether DHA is increased along with AsA or remains constant [21,53], as in the present study. Therefore, changes in the AsA pool following the AsA/DHA ratio can be used as an indicator of stress response ability in the field.
In the present study, field-grown transgenic rice plants that overexpressed OsDHAR1 possessed significantly improved agronomic traits, including rice grain yield and total biomass (root and shoot), compared to WT plants ( Figure 3). The roots and shoots of the transgenic plants were longer and healthier than those of the WT plants, which contributed to enhanced total biomass. The total plant weight, culm weight, root weight, number of panicles per hill, number of spikelets per panicle, filling rate, total grain weight, and 1000-grain weight of the transgenic plants were ∼20. 1, 17.6, 12.5, 18.7, 15.3, 10.3, 19.3, and 4.3% greater, respectively, than those of the WT plants during 2020 and 2021, even though the environmental conditions varied significantly between the two years ( Figure 3C). In particular, the increase in grain yield was caused by increased panicle architecture (number of panicles per hill), spikelet density (number of spikelets per panicle), and 1000-grain weight-based total grain weight ( Figure 3C). In contrast, field-grown salt-sensitive dhar2 plants exhibited poor phenotypes during both the vegetative and reproductive stages, and as a result, the dhar2 mutants exhibited reduced agronomic traits ( Figure S4). Specifically, the total plant weight, culm weight, root weight, number of panicles per hill, number of spikelets per panicle, filling rate, total grain weight, and 1000-grain weight of the dhar2 plants were ∼44.0, 28.4, 72.7, 20.1, 19.0, 36.7, 57.5, and 48.4% lower, respectively, than those of the WT plants ( Figure S4C). Among these traits, the reduced total biomass (total plant weight and root weight) and poor seed maturation of dhar2 mutants distinctly diminished rice grain yield (total grain weight and 1000-grain weight). Indeed, AsA content is correlated with seed development and seedling growth [46]. AsA is present in its reduced state during early embryonic development but undergoes progressive oxidation so that DHA is more prevalent by the time the seeds reach maturity [46,54,55]. Thus, DHA produced during seed germination must be rapidly recycled into AsA.
In previous studies, field-grown transgenic rice plants that overexpress cell rescuerelated genes have exhibited improved growth, development, and yield. More specifically, plants that overexpressed OsLSK1 (large spike S-domain receptor-like kinase 1), OsLRK1 (leucine-rich repeat receptor-like kinase 1), OsEBS1 (enhancing biomass and spikelet number), AtSUC2 (A. thaliana phloem-specific sucrose transporter), and CaPLA1 (Capsicum annuum phospholipiase A1) have been reported to possess more tillers, larger spikelet hulls, longer panicles, more branches per panicle, heavier grains, and greater cellular proliferation than WT plants, thereby improving total grain yield [56][57][58][59]. Meanwhile, the overexpression of OsAGSW1 (ABC1-like kinase related to grain size and weight) has been shown to significantly increase grain size, grain weight, grain filling rate, and 1000-grain weight by regulating the number of external parenchyma cells and the development of vascular bundles, which results in wider and longer spikelet hulls [60]. In addition, elevated expression of OsNRT2.3b (nitrate transporter), OsSPL14 (squamosal promoter binding protein-like 14), OsAMT (ammonia transporter), and OsPTR9 (peptide transporter/nitrate transporter 1 homologue) in transgenic rice enhances pH homeostasis by increasing the uptake of ions (N, Fe, and P), improving grain yield, and improving nitrogen use efficiency [61][62][63][64]. In contrast, the repression of OsCKX2, which encodes cytokinin oxidase 2, increases tiller number and rice yield, whereas constitutive overexpression reduces tiller number and growth in a gene dose-dependent manner [65]. Taken together, it is apparent that genes involved in signaling pathways and transport drive changes in yield-associated traits, whereas antioxidant genes such as OsDHAR1 do not. In addition, our findings also suggest that OsDHAR1 affects functional tolerance to environmental stresses and consequently leads to improvement in growth, seed maturation (shape and size), and productivity of rice plants grown in paddy fields.
Plant growth depends on the synergistic interactions of genes or proteins, and because the yield potential of crop plants is determined by the complex interactions of these factors, we investigated the expression profiles and network maps of a large number of genes in paddy fields [22]. Microarray-based gene expression profiling revealed that OsDHAR1 mediates the environmental adaptability of transgenic rice plants grown under paddy field conditions, despite the underlying mechanism being complex. In the OsDHAR1expressing TR1 plants, genes involved in a wide range of functional categories, including cell cycle, transcription factors, amino acid (serine and threonine) metabolism, glycosphingolipid biosynthesis, photosynthesis process, secretory proteins, cell rescue systems (redox and ion homeostasis, and proteostasis), antioxidants (isoflavone and AsA biosynthesis), and metal detoxification, were upregulated. The elevated genes were networked to the gene ontology corresponding to the regulation of signal transduction, glycosphingolipid biosynthesis, hydrolase activity, protein complex subunit organization, membrane transporter activity, chlorophyll envelope, response to ROS, water-soluble vitamin biosynthesis, redox homeostasis, and isoprenoid biosynthesis (Figures 5 and 6; Table S4). Genes interacting with OsDHAR1, OsGSTU42, GalDH, and OsAPX1 were highly expressed in TR1 plants ( Figure S8).
In addition, a wide range of genes encoding heat shock proteins, molecular chaperone proteins and their cofactors, and genes related to AsA biosynthesis and regeneration were also upregulated in TR1 plants (Table S3). AsA, which is the most abundant antioxidant in a variety of plant species, plays a major role in regulating redox status along with other antioxidants [43]. In addition, AsA also stimulates cell division during embryo development by increasing the proportion of cells progressing through the G1-to-S transition, whereas DHA blocks cell cycle progression and regulates the cell cycle and flowering time [66][67][68]. Accordingly, the suppression of DHAR expression in transgenic tobacco reduces the rates of leaf expansion, delays flowering time, and reduces foliar dry weight [46,54]. High expression of chaperonin-related genes can increase proteostasis by minimizing ROS-induced oxidative damage, such as protein carbonylation, which inhibits and inactivates enzymes, and by enhancing the refolding of aggregated proteins or folding of synthesized proteins under environmental conditions. Moreover, AsA, which functions as an essential cofactor for collagen synthesis enzymes, is linked to AsA-dependent protein folding. High concentrations of DHA enhance the formation of disulfide bonds (oxidative protein folding) in secretory proteins. DHAR harboring protein disulfide isomerase or chaperone-like activity can enhance protein refolding and DHA reduction under oxidative condition [69].
On the other hand, the expression of genes involved in nitrogen metabolism, amino acid (phenylalanine, tyrosine, and tryptophan) biosynthesis, apoptotic process, cell adhesion, flavonoid biosynthesis, and auxin-metabolic pathway were downregulated in TR1 plants. The downregulated genes were networked to the gene ontology corresponding to the MAPK signaling pathway, (secondary active) membrane transporter activity, methyltransferase activity, aromatic amino acid family biosynthetic process, and phenylalanine ammonia-lyase activity ( Figure S7; Tables S5 and S6). Interestingly, the OsDHAR1overexpressing TR1 plants regulated auxin-related genes, which could lead to rice grain yield via panicle and spikelet development, and improved leaf architecture. Accordingly, BG1 (big grain 1)-expressing transgenic rice plants exhibited significantly improved plant biomass, seed weight, and grain yield via regulation of auxin transport. However, the knockdown of BG1 decreased auxin sensitivity and grain size [70].
The total genetic network indicated that the differentially expressed genes formed a complex network and provided important information about the mechanisms underlying environmental adaptability (Figures S9 and S10). Although the mechanism is theoretically explained, it is not known. Recently, gene expression profiling was conducted for specific organs, tissues, and growth stages of field-grown WT rice (O. sativa L. japonica). In Table S7, this field-based profiling provides baseline information for the functional characterization of genes and reveals critical developmental and physiological transitions involved in growth potential [22]. Taken together, the mutual cooperative network of the differentially expressed genes in OsDHAR1-expressing TR1 plants promotes root growth and development, nutrient uptake, ion homeostasis, senescence regulation, pathogen resistance, photosynthesis optimization, redox homeostasis, proteostasis, and lipid metabolism, which together accelerate agronomic performance and productivity by regulating auxin metabolism.

Conclusions
The present study characterized the functions of OsDHAR1 with respect to the alleviation of abiotic stress, especially under paddy field conditions. We found that OsDHAR1-expressing transgenic plants were more tolerant to environmental stress, including high salinity, and we attributed the OsDHAR1-mediated environmental adaptability to increased ROS-neutralizing ability. The observed ROS-neutralizing ability was achieved by enhancing the AsA pool, redox and ion homeostasis, proteostasis, lipid-based membrane stability, photosynthesis capacity, and auxin transport through mutual organic networks between the genes, which improved various agronomic traits, including plant growth (root development), biomass, panicle and tiller number, spikelet per panicle, seed maturation, and grain yield. Furthermore, our findings regarding the molecular, genetic, and phenotypic bases for the responses of rice to abiotic and biotic stresses will facilitate the application of molecular breeding, genetic engineering, and other integrated approaches for developing stress-tolerant crop cultivars.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/antiox11061077/s1, Figure S1: Environmental adaptability of OsDHAR1-expressing transgenic plants in a natural paddy field, Figure S2: Redox state and photosynthetic capacity of OsDHAR1-expressing transgenic plants, Figure S3: Characterization of dhar2knockdown mutants, Figure S4: The poor the agronomic parameters of the dhar2 plants grown in the paddy field, Figure S5: Gene ontology of the genes upregulated in transgenic rice (TR1) plants, Figure S6: Gene ontology of the genes downregulated in transgenic rice (TR1) plants, Figure S7: Genetic network of the genes downregulated in transgenic rice (TR1) plants, Figure S8: Genetic interaction network of OsDHAR1 and other genes in transgenic rice (TR1) plants, Figure S9: Genetic network of up-and downregulated genes in transgenic rice (TR1) plants, Figure S10: Gene to gene network of differentially expressed genes, Table S1: Oligonucleotides used in the present study, Table S2: Copy number and location of the OsDHAR1 transgene, Table S3: Genes upregulated in paddy field-grown transgenic rice (TR1) plants, Table S4: ClueGO classification of the genes upregulated in transgenic (TR1) plants, Table S5: Genes downregulated in paddy field-grown transgenic rice (TR1) plants, Table S6: ClueGO classification of genes downregulated in transgenic (TR1) plants, Table S7: Functional classification of genes differentially expressed in transgenic rice (TR1) plants.