Hypoxia-Induced Aquaporins and Regulation of Redox Homeostasis by a Trans-Plasma Membrane Electron Transport System in Maize Roots

In plants, flooding-induced oxygen deficiency causes severe stress, leading to growth reduction and yield loss. It is therefore important to understand the molecular mechanisms for adaptation to hypoxia. Aquaporins at the plasma membrane play a crucial role in water uptake. However, their role during hypoxia and membrane redox changes is still not fully understood. The influence of 24 h hypoxia induction on hydroponically grown maize (Zea mays L.) was investigated using an oil-based setup. Analyses of physiological parameters revealed typical flooding symptoms such as increased ethylene and H2O2 levels, an increased alcohol dehydrogenase activity, and an increased redox activity at the plasma membrane along with decreased oxygen of the medium. Transcriptomic analysis and shotgun proteomics of plasma membranes and soluble fractions were performed to determine alterations in maize roots. RNA-sequencing data confirmed the upregulation of genes involved in anaerobic metabolism, biosynthesis of the phytohormone ethylene, and its receptors. Transcripts of several antioxidative systems and other oxidoreductases were regulated. Mass spectrometry analysis of the plasma membrane proteome revealed alterations in redox systems and an increased abundance of aquaporins. Here, we discuss the importance of plasma membrane aquaporins and redox systems in hypoxia stress response, including the regulation of plant growth and redox homeostasis.


Introduction
Flooding is a major global problem and occurs more often due to climate change. Flooding induces low-oxygen environments (hypoxia, anoxia), leading to reduced oxygen uptake and cell respiration in the roots. Contrarily, increasing levels of CO 2 and ethylene were observed [1][2][3] whereby an increase in ethylene is also caused by an upregulation of its biosynthesis [4]. Associated with hypoxia is the overproduction of reactive oxygen species (ROS) that causes lipid peroxidation, protein and nucleic acid oxidation, and, furthermore, cell death [5]. Stress adaptations of plants that are tolerant or resistant to submergence generated two survival strategies: (i) The "low-oxygen escape strategy" (LOES) facilitates oxygen supply by an ethylene-induced shoot elongation, which is believed to be more efficient under shallow long-term flooding. (ii) The "low-oxygen quiescence strategy" (LOQS) maintains energy consumption without shoot elongation which is used by plants coping with deep short-term flooding [6]. Ethylene plays an important role in abiotic stress response. On one side, as part of the LOQS, ethylene triggers the development of adventitious roots [7]. It also stimulates cellulases and pectinases to break down cell walls and therefore form aerenchyma [8,9]. Both adventitious roots and aerenchyma enable the oxygen supply in deficient root zones. On the other side, ethylene triggers an increase in shoot elongation [10,11]. This enables plants to avoid the submergence as part of the LOES. To additionally ensure ATP generation under oxygen deficiency, plants switch to the less efficient anaerobic pathway [11][12][13].
The oxidative stress in low-oxygen environments is a result of an imbalance of ROS production and scavenging [14]. In both, peroxidases (Prx) play a crucial role [15][16][17]. The majority of soluble class III peroxidases were downregulated by hypoxia stress, whereas plasma membrane (PM)-bound peroxidases were induced for membrane protection and cell wall remodeling [18,19]. Other antioxidants such as superoxide dismutase (SOD) and catalase (CAT), as well as thioredoxin (TRX), enzymes of the Foyer-Halliwell-Asada pathway, and L-ascorbate, are involved in the detoxification of ROS [5,[20][21][22]. In PM, not only ROSproducing and -scavenging proteins but also transport proteins are located. Aquaporins (AQPs) belong to the major intrinsic proteins (MIPs) [23]. They consist of six transmembrane spanning helices with the N-and C-terminus at the cytosol, five loops, and two highly conserved Asn-Pro-Ala (NPA) motifs [24,25]. Based on sequence similarities, plasma membrane intrinsic proteins (PIPs), as a subgroup of MIPs, are further divided into the subclasses PIP1 and PIP2 [26]. It was shown that PIP1 proteins remain in the endoplasmic reticulum after expression, whereas PIP2 are localized at the PM [27]. The central function of AQPs in plant growth was shown via a four-fold decreased maize root elongation rate after treatment with AQP-inhibiting mercurial compounds [28]. Aquaporins facilitate water transport in either direction [29]. Proteins of the PIP2 subgroup seem to be more efficient as water channels than PIP1 [30,31]. Some PIP1 only act as water channels after forming heterotetramers with PIP2 monomers (e.g., ZmPIP1;2 with ZmPIP2;1) [27,32]. Studies showed that PIPs can also transport small organic molecules, for example, CO 2 [31,33,34] or oxygen [23]. Their role in abiotic stress is increasingly studied. The water channel activity is either downregulated by dephosphorylation under drought stress or protonation under flooding stress [35]. The closure of AQPs by protonation of a fully conserved His residue under flooding stress was observed for SoPIP2;1 [35]. Therefore, regulation of AQP gating by cytosolic pH (acidosis) in anoxic stress has been demonstrated [36]. Evidence for redox regulation was given by oxidative gating of water channels in maize via a decrease in hydraulic conductivity after H 2 O 2 treatment [37,38]. At least for AQP8, redox regulation by persulfidation of Cys-53 was demonstrated [39]. Aquaporin-8 transports H 2 O 2 across the PM and is reversibly gated during cell stress, modulating signal strength and duration.
Regarding other PM proteins, physiological functions of the flavocytochrome b family, e.g., respiratory burst oxidase homologs (Rboh) in flooding stress and ROS signaling, were explored [40,41]. Furthermore, a constitutive transmembrane redox system, the so-called standard system, was demonstrated in animal and plant PM [42][43][44][45][46]. Evidence for an electron transfer from cytosolic NAD(P)H:quinone oxidoreductases via vitamin K to a cytochrome b561 was presented [44,47]. The physiological functions of the cytochrome b561 protein family and other PM redox systems still need further elucidation. In addition, the regulation of AQPs and PM redox systems upon hypoxia is not fully understood.
The present study introduced an oil-based hydroponic setup to investigate molecular changes by hypoxia in maize roots using an integrative transcriptomic and proteomic approach. We discuss the importance of AQPs and PM redox systems in hypoxia stress response.

Plant Material and Growth Conditions
Maize caryopses (Zea mays L. cv. Gelber Badischer Landmais, Saatenunion, Hannover, Germany) were cultivated as described before [19]. First, caryopses were sterilized with 3% hydrogen peroxide for 10 min, then washed and soaked in deionized water for at least 4 h with regularly changing the water, then sterilized again. Trays (410 × 300 × 85 mm) were sterilized with 70% ethanol and filled with wetted germination tissue. About 90 caryopses were placed separately onto and covered with wetted tissue. Finally, the trays were covered with aluminum foil for dark incubation at 26 • C for four days. The germinated seedlings were transferred to 9 L plastic boxes filled with hydroponic culture medium (5. 25  The culture medium was replaced after ten days. For 14 days, the medium was oxygenated by KOH-washed air (compressor type LK60, OSAGA, Glandorf, Germany). Then, hypoxia stress was induced by stopping the oxygenation and preventing further oxygen supply by adding 500 mL of commercially available rapeseed oil to the culture medium. This resulted in decreased concentrations of oxygen in the hydroculture medium after 24 h of stress induction (21% to 3.5 ± 0.5% oxygen). Control plants were continuously supplied with air (21% oxygen). For both treatments, the pH thereby stayed stable at 6.8 ± 0.5 to 6.3 ± 0.7 within 24 h. After 24 h, the roots were harvested between 9 and 10 a.m. (CET). Adhered oil was removed from the roots by washing with 0.1% Triton X-100 for 15-30 s and finally with deionized water.

Ethylene and CO 2
About five maize roots (3-4 g fresh weight) of stressed and non-stressed plants (n = 3 biological replicates) were pooled and placed in sealed glass tubes (40 mL volume) for 2 h. With a Hamilton syringe (Chromatographie Service GmbH, Langerwehe, Germany), 500 µL of each gas phase was injected twice (n = 2 technical replicates) into a gas chromatograph (Shimadzu GC-14A with C-R5A chromatograph, Kyoto, Japan; 2 m Porapek column with mesh 80-100, column temp.: 50 • C, injector temp.: 80 • C, detector temp.: 80 • C, carrier gas: N 2 with 30 mL·min −1 flow rate, combustion gas: compressed air with 300 mL·min −1 flow rate and H 2 with 30 mL·min −1 flow rate). The peak areas for ethylene and CO 2 were normalized to a reference measurement with an internal standard (10 ppm reference solution in air) and to the fresh weight of the maize roots. Standard deviation and Student's t-test were used to determine significant changes between control and stressed samples.
2.4. RNA Analysis 2.4.1. RNA Sequencing RNA isolation was done as described before [19]. The whole root systems, including primary and lateral roots, of five maize plants were harvested after 24 h of hypoxia induction and pooled for each of the three biological replicates (control and stressed), and then ground with liquid nitrogen to a very fine powder (about 0.3 mg fresh weight). This was used for total RNA isolation with the NucleoSpin ® RNA Plant and Fungi Kit (Macherey-Nagel, Düren, Germany). For RNA sequence (RNA-Seq) analysis, a final step of ethanol precipitation with a 1/10th volume of 3 M sodium acetate pH 5.2 and 3 volumes of absolute ethanol (>99.5%) was added before delivering the samples to Macrogen Inc. (Seoul, Republic of Korea) for further analyses. Quality control of the samples was performed by Macrogen using agarose gel electrophoresis and an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Six high-quality RNA samples (RNA Integrity Number > 7) were used for cDNA library construction and sequencing using Illumina Sequencing. The following was used as the reference gene: "B73 RefGen_v4". Available online: https://www.ncbi.nlm.nih.gov/assembly/GCF_000005005.2/ (accessed on 1 April 2019). Expression profiles were represented as the read count and normalized based on the transcript length and depth of coverage. The fragments per kilobase of transcript per million fragments mapped (FPKM) were used for log2 transformation and quantile normalization. Genes with more than 2-fold expression and p ≤ 0.05 were considered as differentially expressed genes (DEG).

Quantitative Reverse Transcription Polymerase Chain Reaction (RT-qPCR)
To determine an earlier transcriptional response, nine representative redox and PIP2 proteins with increased abundance in PM after 24 h hypoxia were selected and analyzed after 6 (7.8 ± 1% oxygen) and 12 h (6.4 ± 1.3% oxygen) of treatment using RT-qPCR as described before [19]. The concentration and purity of the isolated total RNA were determined with a Nanodrop spectrophotometer (Fisher Scientific GmbH, Schwerte, Germany). The concentration of total RNA ranged from 200-360 ng·µL −1 with an absorbance ratio of 2.13 to 2.19 at 260/280 nm. The cDNA was prepared from 100 ng of total RNA with the First Strand cDNA Synthesis Kit (Fisher Scientific GmbH, Schwerte, Germany) according to the manufacturer s protocol. A set of primers were designed (Eurofins Genomics Germany GmbH, Ebersberg, Germany, Table 1). Efficiencies of the primers were checked with a 1:2 to 1:200 dilution series by RT-qPCR (5 min 95 • C, 40 cycles of 10 s 95 • C, and 30 s 60 • C terminating in a melting curve from 65 to 95 • C with 0.5 • C s −1 steps), using Blue S Green qPCR Mix Separate ROX (Biozym Scientific GmbH, Hessisch Oldendorf; Germany) and the CFX 96 Cycler (CFX96 Touch system, Bio-Rad, Munich, Germany). The efficiencies of the primers were above 90%. As a housekeeping gene, Zea mays translational elongation factor EF-Tu (zmtufM, AF264877.1, Q9FUZ6) was used. For RT-qPCR, 1 µg of 1:2 diluted cDNA was used. For statistical analysis, RT-qPCR was performed twice for three biological replicates of each treatment and compared to the housekeeping gene (normalized expression ∆∆Cq). Using the CFX manager software version 3.1 (Bio-Rad, Hercules, CA, USA), the expression (relative to zero) was calculated for each stressed sample. Standard deviation and Student s t-test were used to determine significant changes in the expressions.

Preparation of Subcellular Fractions
Subcellular fractions were prepared as described before [19]. Maize roots were washed (3 mM KCl, 0.5 mM CaCl 2 , 0.125 mM MgSO 4 ) and homogenized (0.25 M sucrose, 50 mM HEPES, 5 mM Na 2 -EDTA, pH 7.5 supplied with 1 mM dithiothreitol and 1% polyvinylpolypyrrolidone) using a Waring blender 7011HS (Waring, Stamford, CT, USA). The homogenate was filtered through a nylon net (125 µM mesh, Hydro-Bios, Kiel, Germany) and 1 mM phenylmethylsulfonyl fluoride was added. This homogenate was centrifuged (10,000× g for 10 min at 4 • C and 48,000× g for 30 min at 4 • C, Avanti J-E centrifuge, rotor types JA-14 and JA-25.50, Beckman Coulter, Krefeld, Germany) which resulted in a supernatant with mainly soluble components and a microsomal pellet. Proceeding from the pellet, PMs were isolated with the aqueous polymer two-phase partitioning using 36 g phase systems [49]. Therefore, the pellet was resolved in phase buffer (0.25 M sucrose, 5 mM KCl, 5 mM phosphate buffer, pH 7.8) and applied to the 27 g phase mixture to give a 36 g phase system (0.25 M sucrose, 5 mM phosphate buffer, pH 7.8, 5 mM KCl, 6.5% Dextran T500, 6.5% polyethylene glycol 3350). The phase system was mixed by inverting the tube 20-25 times and then centrifuged at 1,500× g for 5 min at 4 • C (Heraeus Christ, Hanau, Germany). The upper phase was transferred to a fresh lower phase and the steps were repeated five times. The final PEG-rich upper phase was centrifuged at 105,000× g at 4 • C for 45 min (Optima XPN-80 ultracentrifuge, Beckman Coulter, Krefeld, Germany) and washed twice with buffer (0.25 M sucrose, 50 mM HEPES, pH 7.0) to remove the PEG. The PM pellet was resuspended in this washing buffer and stored at −76 • C. The proteins of the 48,000× g soluble fraction were precipitated with 90% saturated (662 g/L) ammonium sulfate at 4 • C overnight, then pelleted at 15,000× g at 4 • C for 20 min (Avanti J-E centrifuge, Beckman Coulter, Krefeld, Germany) and resolved (0.25 M sucrose, 50 mM HEPES, pH 7.0) for storage at −76 • C. Proteins were quantified with a modified Bradford assay with 0.01% Triton X-100 using Bovine Serum Albumin Standard (PierceTM, Thermofisher Scientific GmbH, Dreieich, Germany) for calibration [50].

Tetramethylbenzidine Oxidase Activity
The activity of haem and copper-containing proteins was determined in soluble and PM fractions of maize roots of at least four biological replicates per treatment. For this, a solution with 4.7 mM 3,3 ,5,5 -tetramethylbenzidine (TMB, PanReac AppliChem, Darmstadt, Germany), 30% methanol in 50 mM sodium acetate buffer pH 5.0) was prepared prior to use [adapted from 49]. The assay contained 790 µL 25 mM sodium acetate buffer pH 5.0, 100 µL 0.3% H 2 O 2 (AppliChem, Darmstadt, Germany), 100 µL TMB solution, and 10 µL protein sample with an average of 0.6 µg soluble protein and 4.5 µg PM protein, respectively. The turnover of TMB was measured photometrical (UV-1800 spectrophotometer, Shimadzu Cooperation, Kyoto, Japan) for 2 min at 652 nm (ε 652 nm = 39 mM −1 ·cm −1 ). The buffer and the two substrates served as a reference. Standard deviation and Student s t-test were used to determine significant changes between control and stressed samples.

Shotgun Proteomics of the Plasma Membrane and Soluble Fraction
Sample preparation (at least three biological and two technical replicates of control and stressed samples) and further mass spectrometry (MS) analysis were done as described before [19]. First, PMs (100 µg) were washed in 2 mL (250 mM sucrose, 50 mM HEPES, 150 mM KCl, 0.01% Triton X-100) for 30 min, then pelleted for 1 h at 16,060× g at 4 • C (Biofuge Fresco, Heraeus, Hanau, Germany). The pellets were incubated in 200 µL solubilization buffer (125 mM Tris pH 6.5, 2% sodium dodecyl sulfate, 5% β-mercaptoethanol, 6 M urea) for 1 h at room temperature and centrifuged again at 16,060× g at 4 • C for 60 min. Proteins of the resulting supernatant (PMs) and proteins of the soluble 48,000 g fraction (100 µg) were precipitated with 1.8 mL methanol:chloroform (4:1) at −20 • C overnight, centrifuged, and washed three times in 0.5 mL pure methanol with centrifugation at 16,060× g at 4 • C for 20 min. The pellets were air-dried and digested in 50 µL buffer (200 mM NH 4 CO 3 pH 8.5, 8 M urea, 10% acetonitrile) with 0.2 µg lysine C (New England Biolabs GmbH, Frankfurt am Main, Germany) at 37 • C for 16-18 h. Then, the samples were diluted 1:3 in 10% acetonitrile, and 10 µL trypsin beads (Poroszyme bulk immobilized trypsin, Applied Biosystem, Waltham, USA) were added and incubated at 37 • C for 16-24 h. The digestion was stopped by adding three times the volume of 0.3% heptafluorobutyric acid and the trypsin beads were removed by multiple centrifugations. Peptides were washed and dried using ZipTips (Agilent Technologies, Santa Clara, CA, USA) in accordance with the manufacturer s protocol.
The peptides were dissolved in 2% acetonitrile and 0.1% formic acid. For each sample, 1 µg was applied randomly on a C18 column (15 cm, 50 mm, column, PepMapR RSLC, 2 µM particle size, Thermo Scientific, Bremen, Germany) for separation during a 90 min gradient at a flow rate of 300 nL min −1 . Measurements were done on a QExactive (Thermo Fisher Scientific, Bremen, Germany) with the following settings: full scan range 350-1800 m/z, max 20 MS2 scans (activation type CID), repeat count 1, repeat duration 30 s, exclusion list size 500, exclusion duration 60 s, charge state screening enabled with a rejection of unassigned and +1 charge states, minimum signal threshold 500. Proteins were identified and quantified using a UniprotKB FASTA download for Zea mays (UP000007305, downloaded July 2017 for PM and downloaded November 2021 for soluble fraction) and the software MaxQuant v1.6.5.0 with the following parameters: first search peptide tolerance 20 ppm, main search tolerance 4.5 ppm, ITMS MS/MS match tolerance 0.6 Da. A maximum of three of the following two variable modifications were allowed per peptide: oxidation of methionine and acetylation of the N-term. A maximum of two dismissed cleavages were tolerated. The best retention time alignment function was determined in a 20 min window. Identifications were matched between runs in a 0.7 min window. An FDR cutout at 0.01 (at Peptide Spectrum Match and protein level) was set with a reversed decoy database. A minimum of seven amino acids was required for the identification of peptides and at least two peptides were required for protein identification. The resulting data matrix was filtered so that there are label-free quantification (LFQs) values in more than four replicates (biological and/or technical replicates) of at least one of the treatments. Missing values were corrected with COVAIN [51]. The LFQ intensities of the replicates were averaged for statistical analyses. Standard deviation and Student s t-test were used to determine significant changes between control and stressed samples.

Analyses of Shotgun Proteomics of the PMs and RNA Sequencing Data
Using the UniProt accessions of the shotgun proteomics data, the gene identification numbers (Entrez_Gene_IDs) were identified. Gene-set enrichment (GO) analyses of the RNA-Seq and shotgun data of the PMs were performed using DAVID v6.8 (for Database for Annotation, Visualization and Integrated Discovery; Available online: https://david.ncifcrf.gov/ (accessed on 2 July 2021) [52,53] with a Fisher s test and Bonferroni correction related to the categories "biological process", "molecular function", and "cellular component". DAVID tool also provides the KEGG (for Kyoto Encyclopedia of Genes and Genomes) pathway, INTERPRO, and SMART domains (for Simple Modular Architecture Research Tool). Prediction of the transmembrane helices was performed with TMHMM Server v2.0 (Available online: http://www.cbs.dtu.dk/services/TMHMM/ (accessed on 3 June 2021)) [54]. Protein parameters (molecular weight and pI) were predicted with ProtParam (Available online: https://web.expasy.org/protparam/ (accessed on 3 June 2021)) [55]. Promoter analyses were done with NewPLACE (Available online: https://www.dna.affrc.go.jp/PLACE/?action=newplace (accessed on 3 June 2021)) [56]. Expression profiles that were obtained from the NimbleGen microarray, provided by [57], were extracted from the Maize eFP Browser.
The proteomics data are available via ProteomeXchange with identifiers PXD028987 (PM) and PXD033386 (soluble fraction).
Genes of the anaerobic pathway (ldh, pdc, adh) were found in RNA-Seq analysis (Supplementary Table S1). Two of the four ldh genes showed a more than 2-fold increase (gene IDs 100193279 and 100282503). The three pdc genes were all more than 5-fold upregulated, with one of them even differentially expressed (gene ID 542651). Fourteen adh genes were found, of which two genes were 6-and 9-fold upregulated (Zmadh1, gene ID 542363 and Zmadh2, gene ID 542364). Two putative adh superfamily genes were differentially downregulated by 4.5-and 2.8-fold (gene IDs 100383565 and 100272357), respectively.
The cytochrome bc1 complex and the cytochrome oxidase 2 are key enzymes of the aerobic pathway. The four found genes were all downregulated, one significantly by 1.8-fold (gene ID 100383274). An alternative oxidase was not regulated on the transcript level (Supplementary Table S1). In the shotgun data of the soluble fraction, one ADH protein (ZmADH, A0A1D6GEX5) was found with a 5-fold decreased abundance (Supplementary Table S2).

RNA Sequence Analysis
The quality control report of the RNA-Seq analysis revealed more than 90% of the bases with quality over a Phred score Q30 for each sample, where the base call accuracy is 99.9% (Supplementary Figure S1). In total, 49,607 genes were determined. After quantile normalization, 29,703 genes were analyzed regarding up-or downregulation. From these normalzsed genes, 1461 genes were differentially expressed (DEG; 975 down-and 486 upregulated) based on 2-fold change and a p-value of p < 0.05 of the comparison pair (stressed versus control) ( Figure 3A). The hierarchical clustering of the DEGs is shown in the heat map ( Figure 3B).

Gene Ontology Analysis of Differentially Regulated Genes
Gene ontology analysis with DAVID revealed that the DEGs were only partially mapped (Supplementary Table S3A). The top 10 terms of each category for all DEGs are shown in the graph (Figure 4, Supplementary Table S4A-C). The main biological processes were "response to oxidative stress" (GO:0006979), "hydrogen peroxide catabolic process" (GO:0042744), and "plant-type cell wall organisation" (GO:0006270), which indicated stress adaptations. The main molecular functions were "heme binding" (GO:0020037), "transporter activity" (GO:0005215), "peroxidase activity" (GO:0004601), and "water channel activity" (GO:0015250), which were also related to these processes. The nucleus (GO:0005634) was the main cellular component, followed by the extracellular region, PM, and cytosol. In addition, the main domains found with INTERPRO were haem peroxidase (IPR010255, IPR002016) and AQP-related (IPR000425, IPR023271). Other domains regarding histones were determined with INTERPRO and SMART. The main domain in SMART was Leucinerich repeats (LRR).
Regarding the L-ascorbate biosynthesis pathway, seven genes were found including the last two enzymes L-galactono-1,4-lactone dehydrogenase (gene ID 100381436) and L-galactose dehydrogenase (gene ID 732776) which were downregulated. In addition, GDP-D-mannose 3 ,5 -epimerase (gene IDs 100192102 and 100280325) and one of the three found phosphomannomutases (gene ID 100272867) were downregulated.
Genes of seven nucleobase-ascorbate transporters were either not or downregulated. Genes of 13 Casparian strip-like proteins (CASPL) were found, with 3 more than 2-fold downregulated (gene IDs 100281882, 100282577, and 100285037) and 1 upregulated (gene ID 100285613). Genes of six aluminum activated malate transporters (ALMTs) were found in RNA-Seq analysis, with one (gene ID 103627532) significantly 1.55-fold upregulated. Eleven genes of the receptor-like kinases (RLK) were found, with two that were upregulated.
The expression of 35 aquaporin genes, including pips, tonoplast intrinsic proteins, nodulin26-like intrinsic proteins, and small basic intrinsic proteins, were determined. Almost all maize PIPs decreased significantly by 2-fold, excluding zmpip1;5, zmpip2;1, and zmpip2;2 where the gene expression showed no significant decrease between controls and hypoxia-stressed samples ( Table 2). The gene expression of zmpip2;7 was not detected.

Gene Expression and Promotor Analysis of PM-Bound Redox Systems and PIPs
The gene expression of PM-bound proteins revealed a relative strong expression in primary roots for zmlox1, zmlox3, zmlox4, zmfqr1, and zmdohcytb561 ( Figure 5A). The minimal expression was determined for zmair12, the maximal expression for zmdohcytb561. In addition, pips showed a relative strong expression in primary roots. The minimal expression was determined for zmpip2;1, the maximal expression for zmpip1;5.
The analyses of cis-regulatory elements ( Figure 5B) concluded a regulation by abiotic factors (anaerobic, hypoxia, and oxygen) including all pips. In addition, motifs for binding sites of ethylene responsive factors (ERF) were found (GCC box, ERE, and DRE), in spite of the fact that the motif for the hormonal regulator ethylene was not determined for zmdohcytb561, zmair12, and zmlox1 genes, and only in four of the 12 pips (zmpip1;1, zmpip1;5, zmpip2;1, and zmpip2;7). A motif found in Arabidopsis thaliana ADH gene promoter (Adh; S000133; CCACGTGG) was detected in zmpdil1-1 and zmpdil1-2, zmlox4, and zmsku5 and for zmpip1;5 only. Motifs regarding oxidative stress were determined in zmlox1, zmlox2, zmlox3, and zmsku5, but not for the pips.

RT-qPCR
The expression analysis with RT-qPCR ( Figure 6   . Regulation of PM-bound redox systems and PIPs by hypoxia stress. Expression analysis of selected genes using SYBRGreen RT-qPCR with total RNA of three non-stressed (white columns) and three hypoxia-stressed (dark gray columns) maize roots either after 6 or 12 h after treatment. Shown are (A) six redox system-related genes and (B) three representative pip genes. Expression was normalized to zmtufM as a housekeeping gene using Biometra CFX96 software. Significances were calculated with Student s t-test and marked with p < 0.05 * and p <0.01 ** for three biological and two technical replicates per sample.

Shotgun Analysis
In the shotgun analysis of the PMs, a total of 485 proteins were analyzed, of which 165 were filtered out (see methods). Data are available via ProteomeXchange. Protein immunoblots of the PMs showed an enrichment of the H + -ATPase and a lower amount of V-PPase and Cox2 signals compared to the corresponding microsomal fractions (Figure 7). In addition, marker analyses of the shotgun data revealed highly enriched PM preparations (Supplementary Table S5A). Plasma membrane-specific P-type ATPases (K7TX67, A0A1D6MV33, and A0A1D6DVJ7) were found in all PM samples, but with higher amounts in stressed samples. Vacuolar V-ATPases (C0PHC0, B6UHI4, and A0A1D6JW70) were detected but with less amount than H + -ATPase and a reduction in stressed samples. Vacuolar V-PPases and mitochondrial cytochrome c oxidase were not detected. Proof of luminal-binding protein 2 BiP2 (P24067) and calreticulin-2 (A0A1D6EN43) suggested typical contamination by the endoplasmic reticulum. The Golgi marker GTP-binding protein SAR1 (B7ZZP2), but not the Sec21p protein, was detected. Typical cell wall proteins (expansins) were not found. Cytosolic markers such as ADH, glyceraldehyde-3-phosphate dehydrogenase, or fructose-1,6-bisphosphatase were not observed in PM samples.

TMB Oxidation
Focusing on the subcellular fractions of maize roots, the TMB oxidase activity of putative haem and copper-containing proteins in the soluble fraction decreased from 27.00 ± 9.22 in the control to 23.56 ± 8.64 µmol·min −1 ·mg −1 in stressed samples ( Figure 9A). This was a decrease to 87% activity in stressed samples relative to the controls. In the PM fraction, the activity increased significantly by 2.9-fold from 2.84 ± 1.02 in control samples to 8.48 ± 2.64 µmol·min −1 ·mg −1 in stressed samples. PMs, the activities of haem and copper-containing proteins were determined spectrophotometrically using the substrates H 2 O 2 and 3,3 ,5,5 -tetramethylbenzidine (TMB). (B) Plasma membranes of control and 24 h hypoxia-stressed maize roots were used for mass spectrometry analyses. Labelfree quantification (LFQ) data were determined with MaxQuant software and used for further statistical analyses. Control (white columns) and stressed samples (gray columns) were shown. Error bars indicate standard deviation. Significances are indicated with asterisk (p < 0.05 *, p < 0.01 **, p < 0.001 ***) determined with Student s t-test for TMB oxidation (n > 4 biological replicates and n = 3 technical replicates) and for shotgun (n = 3 biological and n = 2 technical replicates per treatment).
Three of the PM-bound RLKs showed a more than 2-fold increased abundance (A0A1D6JEH7, A0A1D6DZ88, and A0A1D6IE45). None of the nucleobase-ascorbate transporters, Casparian strip-like proteins, or aluminum-activated malate transporters were found.

Shotgun of Soluble Fraction
In the shotgun analysis of the soluble fractions (Supplementary Table S2), ZmCat1, 5 ZmAPx and 33 ZmPrx were identified. ZmCat1, ZmAPx01, ZmAPx02, ZmAPx05, and 18 ZmPrx showed a more than 2-fold decreased abundance. Only ZmPrx99 showed a 16-fold increase in abundance. In addition, three ZmGPx and four cysteine peroxiredoxins with only one showing a 7-fold decreased abundance (Zm2CysPrx02, B4FM07) were determined. Three SODs were found, with one (ZmSODM4, B4F9H6) showing a 5-fold decreased abundance. Four disulfide isomerases, one thioredoxin (A0A804Q653), and a thioredoxin superfamily protein (C4IZH7), as well as one glutathione reductase (B4FWU6), showed a more than 2-fold decreased abundance. Two DHAR and three MDHAR were found, where the MDHAR showed a more than two-fold decreased abundance. Four of seven found MDHs showed a more than 2-fold increased abundance. In addition, one ZmWrbA (B6TFN1) was determined with a 5-fold decreased abundance.
One of the proteins of the L-ascorbate biosynthesis pathway (B4FBC2) and one RLK (A0A804NHW6) were found. The latter showed a strong decrease in abundance.

Discussion
In the present study, we demonstrated that the oil-based setup induced typical symptoms of flooding-like hypoxia stress in maize roots grown in hydroponics. We showed the induction of AQPs of the PIP2 subgroup and several PM-bound redox systems. Here, we discuss the functions of PM redox systems in root growth and redox homeostasis during hypoxia stress. Finally, the function of hypoxia-induced water channels and their regulation by PM redox systems were reconsidered.

The Oil-Based Setup Induced Typical Hypoxia Stress Symptoms
The efficiency of the oil-based setup for hypoxia stress induction was demonstrated by typical flooding symptoms. For 24 h of hypoxia stress, this includes (i) decreased concentrations of oxygen in the hydroculture medium, (ii) an elevated level of ethylene in maize roots ( Figure 2A) corresponding with an upregulated expression of ethylene biosynthesis genes (Supplementary Table S1), (iii) the formation of aerenchyma [19], and (iv) an increased ADH activity ( Figure 2B) corresponding with an upregulation of anaerobic genes (Supplementary Table S1). The function of ethylene-induced aerenchyma in the facilitation of gas transport and the expression of anaerobic genes as a switch from aerobic to anaerobic metabolism for ongoing ATP energy production was discussed elsewhere [6]. An elevated level of H 2 O 2 in maize roots and higher abundances of class III peroxidases with increased activity were determined in PM preparations under hypoxia stress [19].
Gene ontology analysis of RNA-Seq data revealed that metabolic and biosynthetic processes are downregulated on the gene level ( Figure 4) probably due to potential energy deficit under low-oxygen conditions. The main biological processes of the DEGs and the PM-bound proteins ("response to oxidative stress", "hydrogen peroxide catabolic process", and "plant-type cell wall organization") confirmed previous results [19] and indicated alterations in additional redox proteins. "Peroxidase activity" or "oxidoreductase activity" as the main molecular function and "haem-peroxidase domains" and "LOX do-mains" correlated well with these findings. Here, an increase in redox activity ( Figure 9A) and alterations in abundances of redox proteins at the PM were determined (Table 2). Additionally, "transporter activity" and "water channel activity" as molecular functions and "aquaporin-related domains" on gene and protein levels indicated an alteration in water channels.

Hypoxia-Induced Redox Systems Effected Root Architecture
A function of anaerobic genes and air12 in lateral root formation was shown [59,60]. The lower transcript level and protein abundance of ZmAir12 suggested that its role in lateral root formation may be negligible after 24 h of hypoxia stress ( Table 2). The result demonstrated that this protein is transcriptionally regulated. Alterations of ZmSKU5 and ZmLOX fit with changes in root architecture under hypoxia. SKU5 is involved in directional root growth [61]. Its functions in cell polar expansion and cell wall synthesis have been demonstrated [62]. Cell wall remodeling was supported by upregulation and higher abundances of PM-bound peroxidases under hypoxia. Additionally, upregulation of key enzymes of the phenylpropanoid pathway, higher abundances of dirigent proteins, and a fasciclin-like arabinogalactan were observed [19].
Lipoxygenases inhibit lateral root development [63] but increase the formation of primordia [64]. Lateral root formation seems to be linked to the cooperation of 9-LOX and 13-LOX [64]. Here, ZmLOX1, ZmLOX2, and ZmLOX3 showed increased protein abundances in PM that might be linked to primordia formation under hypoxia ( Table 2). The upregulated expressions of zmlox4, zmlox5, zmlox8, and zmlox10 (Supplementary Table S1) suggested an involvement in ongoing lateral root formation. ZmLOX3 has a function in normal plant growth, whereas an overexpression of zmlox5, zmlox8, and zmlox10 was shown to compensate for a missing zmlox3 expression [65]. Here, the upregulation of zmlox5, zmlox8, and zmlox10, accompanied by the zmlox3 expression and ZmLOX3 abundance, correlated with the observed hypoxia-induced adaptations in root growth [66].
Aquaporins accomplish hydraulic conductivity, which is essential for plant growth. Both the inhibition and mutation of AQPs reduced hydraulic conductivity [28,37,67] accompanied with an impaired root growth [68]. Function of zmpip1;3, zmpip2;1, zmpip2;2, zmpip2;4, and zmpip2;6 in root growth has been suggested [69]. Hypoxic and anoxic treatment also reduced hydraulic conductivity [70][71][72]. Under anoxic conditions, the cytosolic pH was decreased which regulates the gating of AQPs [35,36,73]. Contrarily, the overexpression of PIPs (e.g., scpip1 and ospip1;2) increased root growth [74,75]. The lower abundances or the absence of PIP1s in hydroponically grown maize roots ( Figure 9B) suggests a function in the development of primary roots. Contrary to this, higher abundances of PIP2 under hypoxia might regulate root hair development in hydroponics to compensate nutrient uptake. The mutation in atpip2;4 resulted in longer root hairs in the control and phosphatedeficient environment [76]. Additionally, suberization of the endodermis inhibits water flow through the apoplastic pathway. Suberized endo-and exodermal cells were determined in maize root cross-sections under hypoxia [19]. This effect requires AQPs with a high hydraulic conductivity, which was shown for PIP2 isoforms [30,31]. The abundances of ZmPIP2;1, ZmPIP2;4, ZmPIP2;5, and ZmPIP2;6 increased under hypoxia ( Figure 9B). It was demonstrated that ZmPIP2;5 was more abundant in cortical cells of the suberized root hair zone but not in the exodermis, suggesting a role in radial water transport, restricted by the Casparian strip [77,78]. The function of each PIP isoform in specific root cells and under hypoxia stress needs further evaluation.

Hypoxia-Induced PM Redox Systems Are Involved in ROS Production
The upregulation of zmrboh10 and the observed high levels of H 2 O 2 confirmed ROS production at the PM under hypoxia stress in maize roots [19]. Thus, the higher abundance of ZmPIP2;5 under hypoxia (Table 2) correlated with its function in H 2 O 2 transport [79]. High levels of ROS caused an imbalance of the cellular redox homeostasis. The RNA-Seq analysis revealed a differential regulation of several redox systems in maize roots under hy-poxia (Supplementary Table S1). It appears that the majority of antioxidant genes showed a lack-phase in their response to hypoxia with an upregulation after 48 h [80], possibly as a mechanism to ensure high H 2 O 2 levels for aerenchyma formation [19]. Genes of major antioxidant systems including those of the Foyer-Halliwell-Asada cycle were not or downregulated after 24 h; the only exceptions were zmapx04 and the zmsodm4. Downregulation of zmapx, zmcat, and zmsod and their decreased abundance in shotgun data of the soluble fraction correlated with the observed lower TMB activity of these fractions compared to controls ( Figure 9A). The observed upregulation of mitochondrial SOD and the ubiquinone biosynthetic pathway suggested ROS production by the respiratory chain [81]. The upregulation of the anaerobic pathway makes this suggestion unlikely (Supplementary Table S1). Transcripts of alternative oxidase were not regulated and complex III was downregulated (Supplementary Table S1). The observed upregulation of zmrboh10 and zmdiox01 correlated with their function in signal transduction during hypoxia [82,83].
Shotgun proteomics revealed alterations in abundances for several PM-bound redox systems by hypoxia ( Table 2). The increased TMB oxidation by PM fractions in the presence of H 2 O 2 ( Figure 9A) suggested higher abundances of haem and copper-containing proteins in stressed samples [84,85]. In accordance with this observation, an increase of members of the Cytb561 family (DoHCytb561) and the copper oxidase ZmSKU5 were detected in PM under hypoxia ( Table 2). The higher abundances of class III peroxidases confirmed the observed upregulation and increased activity in maize root PM under hypoxia [19]. In the shotgun experiment, the abundance of ZmRboh10 was below the limits of detection. Either the transcript was not translated or ZmRboh10 increased only in a few specialized cells for aerenchyma formation, as demonstrated for OsRbohH [19,86]. Thus, low-abundant members of the flavocytochrome b family (Rboh) and sulfhydryl-groups of membrane proteins may also account for the observed TMB oxidation [87].
Apart from the upregulation of zmrboh10, the increased abundance of ZmLOX1, ZmLOX2, and ZmLOX3 supported a production of ROS at the PM. The influence of the various LOX isoforms in oxidative stress-induced lipid peroxidation is still unclear. The higher activity and abundance of ZmLOX2 and ZmLOX9 together with increased malondialdehyde levels, a marker for lipid peroxidation, was shown for waterlogged maize and pepper leaves overexpressing calox1 [88,89]. Contrary to this, a lack of atlox1 and atlox5 increased the susceptibility to singlet oxygen and led to lipid peroxidation [90]. The association of ZmLOX1, ZmLOX2, and ZmLOX3 with the maize root PM and their enrichment under hypoxia supported the role in lipid peroxidation.
The turnover of glycolytically derived ATP produce protons as a source of intracellular acidosis [91]. Different proton-consuming metabolic pathways have been discussed for maize under hypoxia [92,93]. A higher abundance of the PM H + -ATPase was observed in hypoxia stressed samples ( Table 2). In addition, cytosolic MDH acts as a biochemical buffer system [94]. Two MDHs showed higher abundances in the PM of stressed samples ( Table 2). At least one of these enzymes showed a significant higher affinity for oxaloacetate compared to malate [95]. The consumption of protons by the reaction of PM-bound MDH may be part of the mechanism that counteract against intracellular acidosis [96]. A PM localized malate shuttle involved in pH control was postulated [97]. This hypothesis was supported by the observed differential regulation of five putative PM-bound aluminumactivated malate transporters (ALMT) by hypoxia, with ALMT10 significantly upregulated (Supplemental Table S1).

Hypoxia-Induced PM Redox Systems Are Involved in Membrane Protection and Redox Homeostasis
Oxidative stress caused increased activities of antioxidative systems [5]. The upregulation of zmapx04, zmsodm4 (Supplementary Table S1), and membrane-bound antioxidative systems in hypoxia-stressed samples was observed [19]. A function in membrane protection was already discussed for the PM-bound class III peroxidases ZmPrx03, ZmPrx24, ZmPrx81, and ZmPrx85 in maize roots under hypoxia [19]. Although physiological func-tions of PM redox systems are not completely understood, the data at hand supported a regulatory role in the apoplastic redox homeostasis [45,98]. This was demonstrated for members of the SKU5 similar protein family [99]. Interestingly, the abundance of a ZmDoHcytb561 increased measurably under 24 h of hypoxia (Table 2). It was suggested that AtDoHcytb561 regenerates apoplastic ascorbate by a transmembrane electron transfer and thereby regulates redox homeostasis during oxidative stress [45]. The importance of this reaction becomes obvious by the downregulation of cellular antioxidant systems, ascorbate biosynthesis, and ascorbate transporters during this early phase of hypoxia stress ( Table 2, Supplementary Table S1) which was also confirmed by the lower abundances of antioxidative systems in the soluble fractions (Supplementary Table S2).
The molecular mechanism for the maintenance of the apoplastic redox homeostasis under oxidative stress was reflected in the co-regulation of ROS producing PM redox systems (Rboh10 and LOX) and membrane protection by antioxidative systems such as PM-bound class III peroxidases, the DoHcytb561, and SKU5. Additionally, vitamin K 1 has been identified in maize and soybean PM [44,47] that acts as a mobile electron and proton carrier and as an effective antioxidant inside the membrane [100,101]. The upregulation of menG (Supplementary Table S1) and the higher abundance of the NAD(P)H:quinone reductase (Table 2) indicated a de novo synthesis of vitamin K 1 in root PM [102]. Further, the upregulation of the ubiquinone biosynthetic pathway (Supplementary Table S1) and downregulation of complex III further supported the function of vitamin K 1 in membrane protection. Both ubiquinone and vitamin K 1 were demonstrated to act as effective antioxidants [103]. During oxidative stress, NAD(P)H:quinone reductases recycle the vitamin K 1 pool inside the PM. Vitamin K 1 may also have the potential to reduce the DoHcytb561. At least a NAD(P)H-dependent cytochrome b reduction has been observed by naphthoquinones in maize root PM [44].

Redox Regulation of H 2 O 2 Transporting Aquaporins
The gene expression of almost all maize root zmpips decreased after 24 h of hypoxia (Table 2), which was also determined in hypoxic lateral and adventitious roots [104]. This downregulation might be related to energy deficits caused by the stress. Additionally, the increased ethylene level in maize roots after 24 h of hypoxia ( Figure 2A) might lead to reduced pip expression, as shown in rose leaves [105]. The missing gene expression of zmpip2;7 in maize roots is probably associated with a restriction to specific cell types beyond root tissue, as shown by [24] and [106].
Shotgun analysis revealed six ZmPIPs with ZmPIP2;1, ZmPIP2;2, ZmPIP2;4, ZmPIP2;5, and ZmPIP2;6 that increased in abundance after 24 h of hypoxia, two of them significantly (ZmPIP2;1 and ZmPIP2;4) ( Table 2). Other maize root PIPs with no or downregulated gene expression might have a protein abundance below the limits of detection. Interaction partners of AQPs for hetero-tetramers or complexes with other proteins were investigated [26,[107][108][109]. On the contrary, a re-localization of AtPIP2;1 from the PM to subcellular membranes was induced by H 2 O 2 [110] that revealed higher levels in hypoxia stressed maize roots [19]. The function of PIPs in water transport and, in this regard, the more efficient PIP2 subgroup was sufficiently studied [30,111]. Ethylene was found to be a direct regulator of the AtPIP's phosphorylation, leading to enhanced water channel activity [112]. The increased ethylene level (Figure 2A) might serve as an activator of the PIP2 subgroup that showed higher abundances in PM of hypoxia-stressed maize roots. In accordance with these results, upregulation of PtPIP2;4 and PtPIP2;5 has been observed in aspen (Populus tremuloides) roots under hypoxia [72]. Anyway, a function in the transport of water and dissolved gases (e.g., oxygen) under hypoxia stress in hydroponics might be significantly reduced. This hypothesis was supported partially by the downregulation of pip1;1 and pip1;2 and the absence of the proteins corresponding to these genes in the PM ( Table 2). The oxygen deficiency is probably counteracted with aerenchyma; the energy deficit with the switch to anaerobic metabolism.
Additionally, the interaction of AtDoHCytb561 with AtPIP1;2 and AtPIP2;1 has been shown under control and stress conditions [118]. The lower abundance of ZmPIP1;2 and the higher abundance of ZmPIP2;1 fit well with the higher abundance of ZmDoHcytb561 and its postulated function in the regulation of the apoplastic redox homeostasis ( Table 2). Besides ZmPIPs, class III peroxidases, ZmDoHcytb561, Rboh, and other PM redox systems have been detected in detergent-resistant membranes that present functional microdomains [119,120]. For AtPIP2;1, the interaction with CASPL and members of the RLK family have been observed [118,121]. At least one CASPL was upregulated in maize roots under 24 h of hypoxia (Supplementary Table S1). While CASPL was suggested to direct PIPs to microdomains with specialized functions, RLKs (RKL1 and Feronia) modulate PIP activity [121]. Three members of the RLK family showed higher abundances in maize root PM and two transcripts were upregulated under hypoxia (Supplementary Table S5B). At least RK20-1 has a function in fine-tuning cell growth, controlling apoplastic pH, and ROS production [122][123][124]. Moreover, the lipid composition of microdomains may have an additional effect on the gating of PIPs [125].
An early idea on the function of PM redox was the regulation of SH-groups of transmembrane proteins, e.g., transporters [42]. Currently, the occurrence of vitamin K 1 in PM is generally accepted and quinone perception via Leucin-rich repeat RLK has been demonstrated [44,47,126]. Vitamin K 1 directly effects the redox status of SH-groups and serves as an electron carrier for protein disulfide bond formation [127][128][129][130]. The later reaction is important for the assembly of protein complexes and redox regulation [131]. Further, H 2 O 2 can also cause intermolecular disulfide bonds [132,133].
Redox regulation of maize AQPs was suggested by [38], who observed an oxidative gating of water channels in maize roots by H 2 O 2 . Amino acid sequences of maize PIPs were highly conserved with three to four Cys residues (Supplementary Figure S3). For ZmPIP2;5, a conserved disulfide bond was demonstrated for Cys-79 between two monomers in loop A [134]. It was shown that Cys-79 was not a target of redox regulation [109]. Redox regulation, as found in AQP8 [39], may be also excluded, because the Cys-53 residue was not conserved in ZmPIPs (Supplementary Figure S3).
The prediction of the putative structure of ZmPIP2;5 and ZmPIP2;6 revealed a conserved Cys residue at transmembrane helix III that was oriented to the lipid bilayer ( Supplementary Figures S2 and S3). A comparable Cys residue has been shown in SoPIP2;1 by crystallography [35]. The higher level of H 2 O 2 [19] and the observed lipid peroxidation in maize roots by hypoxia [135] may cause oxidation of those SH-groups in ZmPIP2;1, ZmPIP2;2 ZmPIP2;4; ZmPIP2;5, or ZmPIP2;6 to sulfinic, sulfenic, or sulfonic acid groups. Vitamin K 1 , as an antioxidant acting inside the PM, can keep those SH-groups reduced.

Transcriptome to Proteome Regulation
After 24 h of hypoxia, transcript and protein abundance at the PM were mainly negatively correlated with a decreased gene expression and an increased protein abundance. A transcriptional positive regulation was indicated for 9.5% of the proteins listed in Table 2 including Air12, WrbA, Prx85, and the PIP1;5. A decrease in pip expression and PIP accumulation was observed under different stress conditions [136,137]. Nonetheless, gene expression usually proceeds early after stress induction. Therefore, genes with a negative correlation were selected and analyzed by RT-qPCR at earlier time points (Figure 6). The results after 6 and 12 h confirmed the downregulation of those genes by hypoxia. To explain the negative correlation, a deeper look into mRNA-to-protein translation processes is necessary [138][139][140]. Some data presented an increased protein abundance with no alteration in its gene expression even under stress [141,142]. In addition, a correlation between mRNA expression and protein product depends on significance [143]. A difference between RNA-Seq and RT-qPCR was observed for smaller genes with fewer exons and lower expression [144], which was demonstrated for the regulation of the PM-bound class III peroxidases under hypoxia stress [19]. In addition, major post-transcriptional regulations during hypoxia were suggested. About 40% of the maize genome were intron-containing genes that undergo alternative splicing [145]. This post-transcriptional modification has been demonstrated at least for PIP2;7 under salt stress [146]. Besides, alternative splicing regulates the targeting of MDH [147]. Altogether, post-transcriptional protein regulation under stress is still a field that needs broader attention. For the MDHs (Table 2), the predicted S-palmitoylation may guide the proteins to the PM. Other post-translational modifications, e.g., ubiquitination, alter the abundances of proteins via proteasomal degradation, which was shown for Na + /K + -ATPase under hypoxia in eukaryotic cell lines [148] and Arabidopsis root microsomes under osmotic stress [149]. An effect of hypoxia on decreased protein degradation also seems possible. Additionally, an alternative protein translation pathway involving hypoxia-induced transcription factors was discussed [138,150,151].

Conclusions
In the present study, an oil-based setup was established to analyze hypoxia-induced changes in maize roots. A combination of transcriptomics and proteomics revealed general molecular mechanisms in hypoxia-response and proteins at the PM that effect root architecture and maintain redox homeostasis under oxidative stress. Proteins involved in ROS production were either upregulated (zmrboh10) or accumulated (ZmLOX1, ZmLOX2, and ZmLOX3) in stressed samples. The simultaneous induction of the PIP2 subgroup supported H 2 O 2 transport across the PM, whereas induction of PM redox systems (Prx, DoHcytb561, NAD(P)H:quinone reductase) increased ROS scavenging, recycled apoplastic ascorbate and the vitamin K 1 pool that act as an antioxidant inside the lipid bilayer. The observed upregulation of vitamin K 1 biosynthesis supported this suggestion. Redox regulation of PIPs will need future elucidation.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/antiox11050836/s1. Table S1: RNA-Seq analyses; Table S2: Overview of selected shotgun and RNA-Seq data of 24 h hypoxia-stressed maize root soluble proteins; Table S3: Quantifying mapping results of DAVID; Table S4: Gene ontology enrichment analysis of differentially expressed genes using DAVID; Table S5: Overview of selected shotgun and RNA-Seq data of 24 h hypoxia-stressed maize root plasma membrane proteins; Table S6: Gene ontology enrichment analysis of PM proteins using DAVID; Figure S1: Quality control of RNA isolation and RNA-Seq data; Figure