Enhanced Drought Stress Tolerance by the Arbuscular Mycorrhizal Symbiosis in a Drought-Sensitive Maize Cultivar Is Related to a Broader and Differential Regulation of Host Plant Aquaporins than in a Drought-Tolerant Cultivar

The arbuscular mycorrhizal (AM) symbiosis has been shown to improve maize tolerance to different drought stress scenarios by regulating a wide range of host plants aquaporins. The objective of this study was to highlight the differences in aquaporin regulation by comparing the effects of the AM symbiosis on root aquaporin gene expression and plant physiology in two maize cultivars with contrasting drought sensitivity. This information would help to identify key aquaporin genes involved in the enhanced drought tolerance by the AM symbiosis. Results showed that when plants were subjected to drought stress the AM symbiosis induced a higher improvement of physiological parameters in drought-sensitive plants than in drought-tolerant plants. These include efficiency of photosystem II, membrane stability, accumulation of soluble sugars and plant biomass production. Thus, drought-sensitive plants obtained higher physiological benefit from the AM symbiosis. In addition, the genes ZmPIP1;1, ZmPIP1;3, ZmPIP1;4, ZmPIP1;6, ZmPIP2;2, ZmPIP2;4, ZmTIP1;1, and ZmTIP2;3 were down-regulated by the AM symbiosis in the drought-sensitive cultivar and only ZmTIP4;1 was up-regulated. In contrast, in the drought-tolerant cultivar only three of the studied aquaporin genes (ZmPIP1;6, ZmPIP2;2, and ZmTIP4;1) were regulated by the AM symbiosis, resulting induced. Results in the drought-sensitive cultivar are in line with the hypothesis that down-regulation of aquaporins under water deprivation could be a way to minimize water loss, and the AM symbiosis could be helping the plant in this regulation. Indeed, during drought stress episodes, water conservation is critical for plant survival and productivity, and is achieved by an efficient uptake and stringently regulated water loss, in which aquaporins participate. Moreover, the broader and contrasting regulation of these aquaporins by the AM symbiosis in the drought-sensitive than the drought-tolerant cultivar suggests a role of these aquaporins in water homeostasis or in the transport of other solutes of physiological importance in both cultivars under drought stress conditions, which may be important for the AM-induced tolerance to drought stress.

The arbuscular mycorrhizal (AM) symbiosis has been shown to improve maize tolerance to different drought stress scenarios by regulating a wide range of host plants aquaporins. The objective of this study was to highlight the differences in aquaporin regulation by comparing the effects of the AM symbiosis on root aquaporin gene expression and plant physiology in two maize cultivars with contrasting drought sensitivity. This information would help to identify key aquaporin genes involved in the enhanced drought tolerance by the AM symbiosis. Results showed that when plants were subjected to drought stress the AM symbiosis induced a higher improvement of physiological parameters in drought-sensitive plants than in drought-tolerant plants. These include efficiency of photosystem II, membrane stability, accumulation of soluble sugars and plant biomass production. Thus, drought-sensitive plants obtained higher physiological benefit from the AM symbiosis. In addition, the genes ZmPIP1;1, ZmPIP1;3, ZmPIP1;4, ZmPIP1;6, ZmPIP2;2, ZmPIP2;4, ZmTIP1;1, and ZmTIP2;3 were down-regulated by the AM symbiosis in the drought-sensitive cultivar and only ZmTIP4;1 was up-regulated. In contrast, in the drought-tolerant cultivar only three of the studied aquaporin genes (ZmPIP1;6, ZmPIP2;2, and ZmTIP4;1) were regulated by the AM symbiosis, resulting induced. Results in the drought-sensitive cultivar are in line with the hypothesis that down-regulation of aquaporins under water deprivation could be a way to minimize water loss, and the AM symbiosis could be helping the plant in this regulation. Indeed, during drought stress episodes, water conservation is critical for plant survival and productivity, and is achieved by an efficient uptake and stringently regulated water loss, in which aquaporins participate. Moreover, the broader and contrasting regulation of these aquaporins by the AM symbiosis in the INTRODUCTION Crop adaptation to new environments is of crucial importance, especially in a climate change scenario. In order to secure food production in the future, efforts need to be directed to understand the mechanisms of plant adaptation and tolerance to abiotic stresses like water shortage, as these events are expected to intensify in coming years (Elliott et al., 2014). Plants cope with drought stress by recruiting drought avoidance and/or drought tolerance mechanisms, which include osmotic adjustment, regulation of stomatal conductance and photosynthesis, production of antioxidant and scavenger compounds or regulation of water uptake and flow in their tissues (Ruiz-Lozano et al., 2012b;Candar-Cakir et al., 2016). Maize is a primary food crop, even more important than other cereals such as rice or wheat since 2012 (Min et al., 2016). The impact of drought on productivity of rice, wheat, and maize will become of capital importance, as these crops represent the 50% of total consumed calories in most populated regions (Lobell et al., 2008).
Maize is fairly susceptible to drought stress, especially in the reproductive phase, experiencing important decreases in yields under drought stress in different world regions (Daryanto et al., 2016). Indeed, maize requires more water at the later vegetative and reproductive stages that at seedlings stage, but at the early crop establishment phase, water stress also influences seedlings adaptation and their grain yield potential, because of premature flowering and a longer anthesis-silk interval (Cao and Wj, 2004;Min et al., 2016) Despite the amount of information about crop responses to water deficit, our knowledge about the mechanisms originating drought tolerance in maize seedlings is still restricted (Min et al., 2016). Previous studies of drought tolerance in maize have shown that tolerant cultivars enhanced antioxidant activity, presented lower lipid peroxidation, improved accumulation of osmolytes and turgor adjustment, maintained photosynthetic activity and regulated aquaporin genes (Anjum et al., 2016;Min et al., 2016).
In this context, the symbiosis of arbuscular mycorrhizal (AM) fungi with plant roots has been shown to be helpful to tolerate and overcome water stress episodes in different plant species (Gholamhoseini et al., 2013;Chitarra et al., 2016), including maize (Boomsma and Vyn, 2008;Bárzana et al., 2014Bárzana et al., , 2015. Authors have previously reported that AM-plant association leads to better plant antioxidant activity, osmotic regulation and root hydraulic properties (Ruiz-Lozano et al., 2012a,b). Also, AM inoculated plants generally present a higher level of photosynthetic pigments, enhanced chlorophyll fluorescence parameters and net photosynthetic rate (Yooyongwech et al., 2016), as well as, a different hormone regulation compared to control plants (Aroca et al., 2008a,b).
In maize, the improvement of physiological plant status of AM inoculated plants when subjected to drought stress has been related to a better uptake of soil nutrients and water, reduced oxidative damage, enhanced root water transport capacity, or facilitated switching between apoplastic and cell-to-cell water transport pathways (Boomsma and Vyn, 2008;Bárzana et al., 2012Bárzana et al., , 2015. Furthermore, the establishment of the AM symbiosis originates extensive morphological alterations in plant root cells, in order to accommodate the presence of an endophytic symbiont, with most of these changes concerning cytoplasmic or vacuolar membranes (Krajinski et al., 2000). Thus, it is not surprising that AM plants may present different pattern of membrane proteins such as aquaporins, candidate proteins to be involved in the exchange of nutrients and water between both organisms (Uehlein et al., 2007;Maurel and Plassard, 2011;Bárzana et al., 2014). Aquaporins are small membrane intrinsic proteins located in different cell membranes and constitute a highly diverse protein family in plants, with at least 30 isoforms in most higher plants. They transport water but some of them can also facilitate the membrane diffusion of other relevant molecules for the plant such as CO 2 , silicon, boron, urea, or ammonia (Li et al., 2014). Recently, oxygen has also been shown to be transported by several Nicotiana tabacum aquaporins, with NtPIP1;3 as the most promising one, which points to the significance of pore-mediated O 2 transport for respiration and opens new perspectives for aquaporins roles in plant physiology (Zwiazek et al., 2017). Each aquaporin isoform often contributes, in concert with other isoforms, to several physiological functions. Thus, their numerous functions in plant growth and development seem to be essential but not well understood yet (Chaumont and Tyerman, 2014;Li et al., 2014;Afzal et al., 2016). However, their role in the maintenance of water homeostasis in the whole plant and in the stress responses has been well established (Chaumont and Tyerman, 2014;Afzal et al., 2016), affecting the radial water flow through the cell-to-cell pathway, which is predominant under conditions of low transpiration such as under drought stress (Steudle and Peterson, 1998). To this regard, it is also remarkable that several aquaporin genes have been found to be AM-responsive in numerous plant species (Krajinski et al., 2000;Aroca et al., 2007;Guether et al., 2009;Bárzana et al., 2014;Chitarra et al., 2016;He et al., 2016;Liu et al., 2016).
There are 36 different aquaporin isoforms in maize (Chaumont et al., 2001). In a recent study, 16 out of these 36 maize aquaporins, belonging to the four maize aquaporin subfamilies (PIPs, TIPs, NIPs, and SIPs), were found to be regulated by the AM fungus R. irregularis (Bárzana et al., 2014). The expression of these proteins varies according to the severity of the stress and depends on the duration of the water shortage period (Bárzana et al., 2014). Essentially, these results highlight the complex regulation of these proteins in the presence of AM symbiosis and their putative role in drought alleviation (Bárzana et al., 2014). Previous studies have provided evidences that the beneficial effects of the AM symbiosis on plant stress tolerance are generally larger in plants sensitive to the imposed stress than in tolerant ones, or under more limiting growing conditions (Subramanian et al., 1995;Subramanian and Charest, 1997;Gianinazzi et al., 2010;Bonneau et al., 2013;Yooyongwech et al., 2016). This has been emphasized also for maize plants (Boomsma and Vyn, 2008). Thus, the above approach can be combined with the use of drought-sensitive and drought-tolerant cultivars for comparative analyses (Rigano et al., 2014;Zhang et al., 2016) and for identification of key aquaporins whose expression is altered by the AM symbiosis in the sensitive cultivar to render it more tolerant (Subramanian and Charest, 1997;Yooyongwech et al., 2016). The present study deals with the hypothesis that aquaporin regulation by the AM symbiosis plays a significant role in the improvement of host plant tolerance to drought stress. Under such situation, aquaporin modulation mediated by the AM symbiosis could lead to improvements of the use of soil water and mineral resources, resulting in higher drought tolerance. The objective was to highlight the differences in aquaporin regulation by comparing the effects of AM symbiosis on root aquaporin gene expression and plant physiology in two maize cultivars with contrasting drought sensitivity. This information would help to identify key aquaporin genes involved in the enhanced drought tolerance by the AM symbiosis. A similar approach has been followed to study aquaporins involved in stomatal gating in rice plants (Vinnakota et al., 2016). Moreover, the present work deeps on the role of aquaporins in drought tolerance and their regulation by AM fungi.

Experimental Design and Statistical Analysis
The experiment consisted of a factorial design with two factors: (1) inoculation treatment, with non-inoculated control plants (C) and plants inoculated with the AM fungus Rhizophagus irregularis, strain EEZ 58 (Ri); (2) water regime, so that one half of the plants were cultivated under well-watered conditions (WW) throughout the entire experiment and the other half of the plants were subjected to drought stress for 12 days before harvest (DS). In addition, two maize cultivars with contrasting tolerance to drought stress were used. One cultivar was sensitive to drought (PR34B39) and the second was tolerant to drought (PR34G13). The different combinations of these factors gave a total of four treatments for the sensitive cultivar and four treatments for the tolerant cultivar. Ten replicates were used for each treatment, giving a total of 80 plants.
Within each maize cultivar, data were subjected to analysis of variance (ANOVA) with inoculation treatment, water regime and inoculation treatment-water regime interaction as sources of variation. Post hoc comparisons with the Duncan's test were used to find out differences between groups. Within each water regime, drought-sensitive and drought-tolerant cultivars were also compared by means of Duncan's test. The expression of the AM fungal aquaporins was analyzed by means of Student's T-test.
Maize (Zea mays L.) seeds from a drought-sensitive (PR34B39) and a drought-tolerant (PR34G13) cultivar were provided by Pioneer Hi-Bred, Spain (DuPont Pioneer Corporation). Seeds were pre-germinated on moist sand for 5 days and then transferred to pots filled with 1250 g of the soil/sand mixture described above.
Mycorrhizal inoculum was bulked in an open-pot culture of Z. mays L. and consisted of soil, spores, mycelia and infected root fragments. The AM fungus was Rhizophagus irregularis (Schenck and Smith), strain EEZ 58. Ten grams of inoculum with about 60 infective propagules per gram (according to the most probable number test), were added to appropriate pots at sowing time. Non-inoculated control plants received the same amount of autoclaved mycorrhizal inoculum together with a 3 ml aliquot of a filtrate (<20 µm) of the AM inoculum in order to provide a general microbial population free of AM propagules.

Growth Conditions
The experiments were carried out under greenhouse conditions with temperatures ranging from 19 to 25 • C, 16/8 light/dark period, a relative humidity of 50-60% and an average photosynthetic photon flux density of 800 µmol m −2 s −1 , as measured with a light meter (LICOR, Lincoln, NE, United States, model LI-188B). Plants were cultivated for a total of 9 weeks.
Soil moisture was measured with the ML2 ThetaProbe (AT Delta-T Devices Ltd., Cambridge, United Kingdom). Water was supplied daily to maintain soil at 100% of field capacity during the first 6 weeks after sowing. The 100% soil water holding capacity corresponds to 22% volumetric soil moisture measured with the ThetaProbe, as determined experimentally in a previous experiment using a pressure plate apparatus. Then, half of the plants were allowed to dry until soil water content reached 60% of field capacity (1 day needed), while the other half were maintained at field capacity. At this stage AM and non-AM plants of both genotypes had comparable size. The 60% of soil water holding capacity corresponds to 7% volumetric soil moisture measured with the ThetaProbe (also determined experimentally with a pressure plate apparatus in a previous assay). The soil water content was daily measured with the ThetaProbe ML2 before rewatering (at the end of the afternoon), reaching a minimum soil water content around 55% of field capacity in the drought-stressed treatments. This water deficit treatment resulted in severe drought stress for maize plants, as evidenced by the decrease in stomatal conductance and efficiency of the photosystem II. The amount of water lost was added to each pot in order to keep the soil water content at the desired levels of 7% of volumetric soil moisture (Porcel and Ruiz-Lozano, 2004). Plants were maintained under such conditions for 12 additional days before harvesting.

Biomass Production and Symbiotic Development
At harvest (8 weeks after sowing) the shoot and root system of five replicates per treatment were separated and the dry weight (DW) measured after drying in a forced hot-air oven at 70 • C for 2 days.
The percentage of mycorrhizal fungal colonization in maize plants was estimated by visual observation according to Phillips and Hayman (1970). The extent of mycorrhizal colonization was calculated according to the gridline intersect method (Giovannetti and Mosse, 1980) in five replicates per treatment.

Stomatal Conductance
Stomatal conductance was measured 2 h after the onset of photoperiod with a porometer system (Porometer AP4, Delta-T Devices Ltd., Cambridge, United Kingdom) following the user manual instructions. Stomatal conductance measurements were taken in the second youngest leaf from eight different plants of each treatment.

Photosynthetic Efficiency
The efficiency of photosystem II was measured with FluorPen FP100 (Photon Systems Instruments, Brno, Czechia), which allows a non-invasive assessment of plant photosynthetic performance by measuring chlorophyll a fluorescence. FluorPen quantifies the quantum yield of photosystem II as the ratio between the actual fluorescence yield in the light-adapted state (FV ) and the maximum fluorescence yield in the lightadapted state (FM ), according to Oxborough and Baker (1997). Measurements were taken in the second youngest leaf of eight different plants of each treatment.

Membrane Electrolyte Leakage
Leaf electrolyte leakage (EL) was determined in six plants per treatment. Leaf samples were washed with deionized water to remove surface-adhered electrolytes. The samples were placed in closed vials containing 10 mL of deionized water and incubated at 25 • C on a rotary shaker (at 100 rpm) during 3 h, and the electrical conductivity of the solution (L 0 ) was determined using a conductivity meter (Metler Toledo AG 8603, Switzerland). Samples were then placed at −80 • C for 2 h. Subsequently, tubes were incubated again at room temperature under smoothly agitation and the final electrical conductivity (L f ) was obtained after 3 h under these conditions. The EL was defined as follows: [(L 0 − L water )/(L f − L water )] × 100, where L water is the conductivity of the deionized water used to incubate the samples.

Oxidative Damage to Lipids
Lipid peroxides were extracted by grinding 500 mg of fresh leaf tissues with and ice-cold mortar and 6 ml of 100 mM potassium phosphate buffer (pH 7). Homogenates were filtered through one Miracloth layer and centrifuged at 15,000 × g for 20 min. The chromogen was formed by mixing 200 ml of supernatants with 1 ml of a reaction mixture containing 15% (w/v) trichloroacetic acid (TCA), 0.375% (w/v) 2-thiobarbituric acid (TBA), 0.1% (w/v) butyl hydroxytoluene, 0.25 N HCl and by incubating the mixture at 100 • C for 30 min (Minotti and Aust, 1987). After cooling at room temperature, tubes were centrifuged at 800 × g for 5 min and the supernatant was used for spectrophotometric reading at 532 nm. Lipid peroxidation was estimated as the content of 2-thiobarbituric acid-reactive substances (TBARS) and expressed as equivalents of malondialdehyde (MDA) according to Halliwell and Gutteridge (1985). The calibration curve was made using MDA in the range of 0.1-10 nmol. A blank for all samples was prepared by replacing the sample with extraction medium, and controls for each sample were prepared by replacing TBA with 0.25 N HCl. In all cases, 0.1% (w/v) butyl hydroxytoluene was included in the reaction mixtures to prevent artifactual formation of 2-TBARS during the acid-heating step of the assay.

Total Soluble Sugars Accumulation
At harvest, total soluble sugars (TSS) were extracted from 1 g fresh leaf tissues in 100 mM potassium phosphate buffer for TSS. Soluble sugars were analyzed by 0.025 mL of plant extract reacting with 3 ml freshly prepared anthrone [200 mg anthrone + 100 ml 72% (w:w) H 2 SO 4 ] and placed in a boiling water bath for 10 min according to Irigoyen et al. (1992). After cooling, the absorbance at 620 nm was determined in a spectrophotometer Hitachi U-1900 (Hitachi Corporation, Japan). The calibration curve was made using glucose in the range of 0.2 to 0.4 mg/ml.

Hydrogen Peroxide Content
Hydrogen peroxide content was determined by Patterson's method (Patterson et al., 1984;Aroca et al., 2003), with slight modifications as described previously by Aroca et al. (2003). Five hundred milligrams of fresh leaf tissues were homogenized in a cold mortar with 5 ml 5% (w/v) TCA containing 0.1 g of activated charcoal and 1% (w/v) PVPP. The homogenate was centrifuged at 18,000 × g for 10 min. The supernatant was filtered through a Millipore filter (0.22 mm). A volume of 1.2 ml of 100 mM potassium phosphate buffer (pH = 8.4) and 0.6 ml of the colorimetric reagent were added to 130 ml of the supernatant. The colorimetric reagent was freshly made by mixing 1:1 (v/v) 0.6 mM potassium titanium oxalate and 0.6 mM 4-2 (2-pyridylazo) resorcinol (disodium salt). The samples were incubated at 45 • C for 1 h and the absorbance at 508 nm was recorded. The blanks were made by replacing leaf extract by 5% TCA.

Root Hydraulic Conductivity (Lo)
Eight weeks after sowing the sap flow rate (Jv) and Lo were measured on detached roots exuding under atmospheric pressure for 2 h (Aroca et al., 2007). Osmotic root hydraulic conductivity (Lo) was calculated as Lo = Jv/ , where Jv is the exuded sap flow rate and the osmotic potential difference between the exuded sap and the nutrient solution where the pots were immersed. These measurements were carried out 3 h after the onset of light.

Quantitative Real-Time RT-PCR
Total RNA was isolated from maize roots harvested at noon 8 weeks after sowing and kept at −80 • C, by a phenol/chloroform extraction method followed by precipitation with LiCl (Kay et al., 1987). The RNA was subjected to DNase treatment and reversetranscription using the QuantiTect Reverse Transcription Kit (Qiagen), following the instructions provided by manufacturer. To rule out the possibility of a genomic DNA contamination, all the cDNA sets were checked by running control PCR reactions with aliquots of the same RNA that have been subjected to the DNase treatment but not to the reverse transcription step.
The expression of the group of maize aquaporins previously selected as regulated by the AM symbiosis (Bárzana et al., 2014) was studied by real-time PCR by using iCycler system (Bio-Rad, Hercules, CA, United States) adjusting protocols to optimize the PCR reaction to each gene. The primer sets used to amplify each aquaporin gene were designed in the 3 and 5 untranslated regions of each gene in order to avoid unspecific amplification of the different aquaporin genes (Hachez et al., 2006;Bárzana et al., 2014). The specificity of amplicons was checked with a heat dissociation protocol (from 70 to 100 • C), after the final PCR cycle. The efficiency of the primer sets was evaluated with the software Bio-Rad iQ5 (version 2.1.97.1001) by analyzing the ratio Ct/fluorescence at four-six independent points of PCR curves (Ramakers et al., 2003), giving values between 90 and 98%. The sequences of primers used for the aquaporin and constitutive genes are those described in Bárzana et al. (2014). Standardization was carried out based on the expression of the best-performing reference gene under our growing conditions. Thus, aquaporin expression levels were normalized according to the elongation factor 1 (gi:2282583).
The fungal aquaporin genes GintAQP1, GintAQPF1, and GintAQPF2 were also analyzed using the primers and conditions described previously (Aroca et al., 2009;. Standardization was carried out based on the expression of the fungal elongation factor 1a gene in each sample. The relative abundance of transcripts was calculated by using the 2 − ct method (Livak and Schmittgen, 2001). Real-time PCR measurements were carried out in at least three independent RNA samples per treatment, with the threshold cycle (Ct) determined in duplicate. Negative controls without cDNA were used in all PCR reactions.

AM Root Colonization and Plant Biomass
Arbuscular mycorrhizal inoculated plants from drought-sensitive cultivar, PR34B39, had an average of 54% of mycorrhizal root length, with no significant differences due to the water treatment ( Figure 1A). In the case of the drought-tolerant cultivar, PR34G13, mycorrhizal root length was 50%, also with no significant differences due to the water treatment. Uninoculated maize plants did not exhibit AM root colonization ( Figure 1A).
Shoot dry weight (SDW) in the sensitive line was similar for both AM and non-AM plants when cultivated under well-watered conditions. Drought stress decreased FIGURE 1 | (A) Percentage of mycorrhizal root length, (B) shoot dry weight (SDW), and (C) root dry weight (RDW) in two maize genotypes differing in drought tolerance and inoculated or not with an arbuscular mycorrhizal (AM) fungus. Data represents the means of five values ± SE for Mycorrhization and RDW and 5 values ± SE for SDW. Different letter indicates significant differences between treatments (p < 0.05) based on Duncan's test for sensitive (uppercase) and tolerant (lowercase) genotypes. Asterisks indicate significant differences between drought-sensitive and drought-tolerant genotypes within each watering regime, according to Duncan's test.
SDW by 37% in non-AM plants but only by 17% in AM plants ( Figure 1B). When subjected to drought stress AM plants produced 35% more SDW than non-AM plants ( Figure 1B). In the drought-tolerant cultivar, no effect of the AM symbiosis on SDW was observed either under well-watered conditions or under drought stress. Drought stress decreased SDW by 17 and 22% in non-AM and AM plants, respectively ( Figure 1B). In any case, under drought stress conditions, significant differences in SDW between drought-sensitive and drought-tolerant non-AM plants were observed, with the latter growing 41% more than the former ( Figure 1B).
Drought and AM inoculation had a similar effect on root dry weight (RDW) as in SDW for both PR34B39 and PR34G13 lines ( Figure 1C). In non-AM plants drought decreased significantly RDW in the sensitive genotype, and AM plants produced higher root biomass under drought stress conditions only in the sensitive cultivar ( Figure 1C). In contrast, under well-watered conditions, AM plants enhanced RDW only in the tolerant cultivar.  (gs) and (B) photosystem II efficiency in the light-adapted state ( Fv/Fm ) in two maize genotypes differing in drought tolerance and inoculated or not with an AM fungus. Data represents the means of eight values ± SE. Different letter indicates significant differences between treatments (p < 0.05) based on Duncan's test for sensitive (uppercase) and tolerant (lowercase) genotypes. Asterisks indicate significant differences between drought-sensitive and drought-tolerant genotypes within each watering regime, according to Duncan's test.

Stomatal Conductance (gs) and Efficiency of Photosystem II
The stomatal conductance (gs) of drought-sensitive cultivar was enhanced by the AM symbiosis under well-watered conditions (36% of increase) but not under water deficit. In the droughtsensitive cultivar drought did not significantly affect this parameter (Figure 2A). The drought-tolerant cultivar showed enhanced gs by the AM symbiosis both under well-watered conditions (27%) and under drought stress conditions (143%) (Figure 2A). However, drought decreased this parameter as compared to well-watered conditions. This decrease was 69% in non-AM plants and 41% in AM plants (Figure 2A). Under well-watered conditions both AM and non-AM plants exhibited higher gs values in the drought-tolerant cultivar than in the drought-sensitive one (Figure 2A). In contrast, under drought stress conditions, non-AM plants from the drought-tolerant cultivar had lower gs values than the corresponding droughtsensitive ones.
The light-adapted maximum quantum yield of PSII primary photochemistry ( F/Fm ) in plants from drought-sensitive cultivar was affected by drought stress in the non-AM plants only, which reduced this parameter by 42% (Figure 2B). In contrast, in the AM plants, no significant effect was observed. In the case of the drought-tolerant cultivar, the F/Fm was enhanced by the AM symbiosis both under well-watered conditions (13% of increase) and under drought stress conditions (36% of increase) ( Figure 2B). In this cultivar, drought stress also reduced this parameter (by 19%) in non-AM plants only (Figure 2B). Under drought stress conditions, significant differences in F/Fm between drought-sensitive and drought-tolerant non-AM plants were observed, with the latter having values 30% higher than the former (Figure 2B).

Membrane Electrolyte Leakage
The membrane EL was reduced by the AM symbiosis in drought-sensitive plants, both under well-watered conditions (50% of decrease) and under drought stress conditions (67% of decrease) (Figure 3A). The imposed drought stress increased this parameter by 58% but only in non-AM plants. In the droughttolerant cultivar the membrane EL increased by drought stress only in non-AM plants (by 279%), while AM plants did not increase this parameter as a consequence of drought ( Figure 3A).
The EL values were higher in non-AM drought-sensitive plants than in non-AM drought-tolerant ones, both under well-watered and under drought stress conditions ( Figure 3A).

Oxidative Damage to Lipids (MDA)
The AM symbiosis reduced the oxidative damage to lipids measured as MDA equivalents in the drought-sensitive cultivar regardless of the water regime ( Figure 3B). In contrast, drought stress did not significantly affect this parameter either in the AM or in the non-AM plants ( Figure 3B). In the drought-tolerant cultivar, the oxidative damage to lipids was not significantly affected by the AM symbiosis or by the drought stress imposed ( Figure 3B).

Total Soluble Sugars
The leaf TSS concentration was significantly increased by the AM symbiosis in both maize cultivars, but only under drought stress conditions ( Figure 3C). Plants cultivated under well-watered conditions did not alter their TSS content as consequence of the AM symbiosis ( Figure 3C).

Accumulation of Hydrogen Peroxide
The accumulation of hydrogen peroxide was significantly affected by the AM symbiosis only in the drought-sensitive cultivar, increasing the values in AM plants under drought stress conditions (Figure 3D). Under drought stress conditions, hydrogen peroxide accumulation was higher in non-AM drought-tolerant plants than in non-AM drought-sensitive ones ( Figure 3D).

Root Hydraulic Conductivity (Lo)
In the drought-sensitive cultivar root hydraulic conductivity (Lo) was strongly reduced by drought, but this reduction reached 95% in non-AM plants and 73% in AM plants (Figure 4). Thus, under drought stress conditions AM plants exhibited enhanced Lo values by five-fold when compared to non-AM plants (Figure 4).
In the drought-tolerant cultivar, AM symbiosis reduced Lo by 33% under well-watered conditions but increased it by 82% under drought stress conditions. In this cultivar, drought stress also reduced significantly this parameter (Figure 4). Thus non-AM plants decreased Lo by 81% as consequence of drought. The decrease was by 49% in AM plants (Figure 4). Under drought ) concentration in two maize genotypes differing in drought tolerance and inoculated or not with an AM fungus. Data represents the means of six values ± SE for EL and three values ± SE for MDA, TSS, and H 2 O 2 . Different letter indicates significant differences between treatments (p < 0.05) based on Duncan's test for sensitive (uppercase) and tolerant (lowercase) genotypes. Asterisks indicate significant differences between drought-sensitive and drought-tolerant genotypes within each watering regime, according to Duncan's test.
FIGURE 4 | Osmotic root hydraulic conductivity (Lo) in two maize genotypes differing in drought tolerance and inoculated or not with an AM fungus. Data represents the means of four values ± SE. Different letter indicates significant differences between treatments (p < 0.05) based on Duncan's test for sensitive (uppercase) and tolerant (lowercase) genotypes. Asterisks indicate significant differences between drought-sensitive and drought-tolerant genotypes within each watering regime, according to Duncan's test. stress conditions, Lo values were significantly higher (by 360%) in non-AM drought-tolerant plants than in non-AM droughtsensitive ones (Figure 4).

Expression of Maize and Fungal Aquaporins
We analyzed the expression of 16 maize aquaporins previously shown to be regulated by the AM symbiosis under drought stress conditions (Bárzana et al., 2014). One of these genes, ZmTIP4;2, was not detected in any of the two genotypes, likely because of its low expression level. Besides, ZmNIP1;1 was only detected in the sensitive genotype, but its expression was also very low, and it was not possible to detect any modification due to mycorrhization (Supplementary Figure S1). When analysing the expression patterns in both maize cultivars some of these aquaporin genes were not affected by the AM symbiosis or by the drought stress in the drought-tolerant cultivar (ZmPIP1;2, ZmPIP1;4, ZmTIP1;2, ZmNIP2;2, and ZmSIP2;1) (Supplementary Figure  S1). The Figure 5 shows the expression data of the aquaporin genes that are affected by the AM symbiosis and/or drought stress in both maize cultivars or at least in the drought-sensitive cultivar.
The expression of ZmPIP1;1 in the drought-sensitive cultivar was unaltered by the AM symbiosis under well-watered conditions. However, its expression was enhanced by drought stress by 108% in non-AM plants, while in AM plants its expression did not change as consequence of drought. Thus, under drought stress conditions, the expression of ZmPIP1;1 gene was 77% lower in AM than in non-AM plants ( Figure 5A). On the contrary, in the drought-tolerant cultivar drought induced ZmPIP1;1 expression by 60% in AM plants only (Figure 5A). Under drought stress conditions, the expression of ZmPIP1;1 was significantly higher in non-AM drought-sensitive plants than in non-AM drought-tolerant ones ( Figure 5A).
In the drought-sensitive cultivar the expression of ZmPIP1;3 gene was reduced under well-watered conditions by 70% due to mycorrhization ( Figure 5B). In the same way, the exposition to FIGURE 5 | Expression of ZmPIP1;1 (A), ZmPIP1;3 (B), ZmPIP1;6 (C), ZmPIP2;2 (D), ZmPIP2;4 (E), ZmTIP1;1 (F), ZmTIP2;3 (G), ZmTIP4;1 (H), and ZmNIP2;1 (I), in two maize genotypes differing in drought tolerance and inoculated or not with an AM fungus. Values in the Y -axis represent the expression levels in relative units. Data represents the means of three values ± SE. Different letter indicates significant differences between treatments (p < 0.05) based on Duncan's test for sensitive (uppercase) and tolerant (lowercase) genotypes. Asterisks indicate significant differences between drought-sensitive and drought-tolerant genotypes within each watering regime, according to Duncan's test.
drought stress reduced the expression of this gene by 65% in non-AM plants, reaching expression values similar to those in AM plants. AM plants showed unaltered expression levels under well-watered and drought stress conditions ( Figure 5B). In the drought-tolerant cultivar AM and non-AM plants showed no significant differences in ZmPIP1;3 expression levels under wellwatered and under drought stress conditions. Drought stress only reduced the expression of this gene in AM plants as compared to non-AM plants under well-watered conditions ( Figure 5B). Under drought stress conditions, significant differences in the expression of ZmPIP1;3 gene between non-AM drought-sensitive and drought-tolerant plants were observed, being higher in the latter than in the former ( Figure 5B).
ZmPIP1;6 was down-regulated by the AM symbiosis under well-watered conditions in the drought-sensitive cultivar, showing 60% of inhibition as compared to non-AM plants ( Figure 5C). Drought stress inhibited the expression of this gene in non-AM plants, while in AM plants no further inhibition was observed, as compared to well-watered conditions. In the drought-tolerant cultivar, the AM symbiosis up-regulated by 150% the expression of ZmPIP1;6 under well-watered conditions ( Figure 5C). However, when plants were subjected to drought stress such up-regulation was avoided, reaching similar values than non-AM plants. No changes in gene expression due to water regime were observed in non-AM plants for this gene. Non-AM plants exhibited higher expression levels of ZmPIP1;6 gene in the drought-sensitive cultivar than in the drought-tolerant one, regardless of water regime. In contrast, under well-watered conditions, AM plants had significantly higher expression levels in the drought-tolerant cultivar than in the drought-sensitive one ( Figure 5C).
The expression of ZmPIP2;2 in the drought-sensitive cultivar was significantly reduced by mycorrhization only under drought stress conditions (reduction by 62%) (Figure 5D). Under well-watered conditions this decrease was not significant. The expression of this gene in AM plants subjected to drought was also significantly lower (by 64%) than in non-AM plants under well-watered conditions. In the case of the drought-tolerant cultivar the behavior was different since AM plants up-regulated this gene by 71% under well-watered conditions. In contrast, drought stress inhibited the expression of this gene in AM plants by 58% as compared to well-watered counterparts ( Figure 5D).
In AM plants cultivated under well-watered conditions the expression of ZmPIP2;2 was higher in the drought-tolerant cultivar than in the drought-sensitive one. The opposite was observed in non-AM plants when cultivated under drought stress conditions ( Figure 5D).
The mRNA level of ZmPIP2;4 was reduced by 72% by mycorrhization in the drought sensitive cultivar when cultivated under well-watered conditions (Figure 5E). The expression of this gene was not further inhibited by drought stress in AM plants, while in non-AM plants it was reduced by 59%. In the drought-tolerant cultivar the expression of ZmPIP2;4 did not show significant differences due to mycorrhization or water regime ( Figure 5E).
In the drought-sensitive cultivar the expression of ZmTIP1;1 gene was significantly affected by the AM symbiosis only under drought stress conditions, reducing its expression levels by 63% in AM plants as compared to non-AM ones (Figure 5F). Drought stress itself did not significantly affect the expression of this gene in both AM and non-AM plants. In the drought-tolerant cultivar the expression of ZmTIP1;1 was unaltered by mycorrhization or water regime (Figure 5F). Under drought stress conditions, the expression of ZmTIP1;1 was significantly higher in AM droughttolerant plants than in AM drought-sensitive plants ( Figure 5F).
The expression of ZmTIP2;3 in the drought sensitive cultivar was inhibited by mycorrhization when cultivated under wellwatered conditions, with a reduction of 77% ( Figure 5G). The expression of this gene was not further inhibited by drought stress. In the drought-tolerant cultivar the expression of ZmTIP2;3 was unaltered by mycorrhization or water regime ( Figure 5G).
The ZmTIP4;1 expression was up-regulated under wellwatered conditions in the drought-sensitive cultivar as consequence of AM root colonization, with an increase in expression levels by 122% ( Figure 5H). However, the drought stress reduced the expression of this gene by 58%, reaching similar expression levels than non-AM plants. In the case of the drought-tolerant cultivar, the expression of ZmTIP4;1 gene in non-AM plants was low and it was induced by the AM symbiosis both under well-watered (by 210%) and under drought stress conditions (by 310%) (Figure 5H). Non-AM plants had higher ZmTIP4;1 expression levels in the drought-sensitive cultivar than in the drought-tolerant one, regardless of water regime. Moreover, under well-watered conditions, AM plants also exhibited significantly higher ZmTIP4;1 expression in the drought-sensitive cultivar than in the drought-tolerant one ( Figure 5H).
In the drought-sensitive cultivar the expression of ZmNIP2;1 was only affected by drought stress, which reduced its expression in both non-AM (by 68%) and AM plants (by 87%) (Figure 5I). In the drought-tolerant cultivar similar data were observed, with a reduction of gene expression by drought in non-AM plants (by 90%) and in AM plants (by 79%) (Figure 5I).
The expression of GintAQP1 was slightly induced by drought stress in the drought-sensitive cultivar ( Figure 6A). The expression of this gene was significantly higher in the droughttolerant cultivar under well-watered conditions, but it resulted considerably inhibited (by 80%) by drought stress in this cultivar. The gene GintAQPF1 resulted similarly inhibited by drought stress in both cultivars ( Figure 6B). However, the expression of this gene was lower than that of the other two fungal genes. Finally, the expression of GintAQPF2 resulted unaltered by drought stress in both maize cultivars ( Figure 6C).

DISCUSSION
This study highlights the divergent responses to AM symbiosis of two maize genotypes differing in drought tolerance: PR34G13, a drought-tolerant cultivar, and PR34B39, a drought-sensitive cultivar (DuPont Pioneer Corporation). It particularly focused on the differential regulation of root aquaporins by the AM symbiosis under well-watered and drought stress conditions and its impact on plant performance. We also featured the influence of such factors on plant growth as well as on traits showing the effects of drought and AM symbiosis on plant physiology.

AM Effects on Plant Physiological Status
The AM fungal root colonization in both genotypes exceeded 50%, being not significantly affected by the drought stress treatment, probably due to its limited duration of only 12 days. The AM symbiosis has been previously reported to enhance drought tolerance of host plants (Augé, 2001). In the present study, the beneficial effect of the AM fungus was firstly observed in plant biomass production. Indeed, plant biomass production is an integrative index of plant performance under stressful conditions and the efficiency of the AM symbiosis has often been measured in terms of host plant biomass improvement (Ruiz-Lozano et al., 2012a). Droughted AM plants from the sensitive genotype presented higher SDW and RDW compared to noninoculated plants. Contrariwise, no enhancement of SDW and RDW was observed in the case of the drought-tolerant genotype, highlighting genotype-depending responses to AM inoculation (Subramanian et al., 1995;Subramanian and Charest, 1997;Gianinazzi et al., 2010;Yooyongwech et al., 2016). Anyway, water deficit negatively affected growth in both maize cultivars, but to a lesser extent in the drought-tolerant genotype.
Many of the physiological responses of plants to drought stress are directed toward the control of transpiration, of root hydraulic conductivity and of osmotic adjustment . Stomatal closure is a conserved mechanism in both maize genotypes studied, regardless of AM inoculation. A recent metaanalysis of 460 studies revealed that even if AM-inoculated C3 plants usually show higher gs values, C4 plants featured increases in gs of around 12% (Augé et al., 2015). In agreement with this, AM symbiosis increased gs in both genotypes, especially in the case of the tolerant genotype under drought conditions. However, no differences were found in gs values in the drought-sensitive genotype subjected to the water stress. This could be probably related to the larger SDW of these plants with the consequent increased transpiring area, or to the fact that drought-sensitive plants had generally lower gs values than the drought-tolerant ones. Nonetheless, it is noteworthy that the maize gs response to fungal inoculation could be described as inconsistent, ranging from unaltered to increased by AM (Boomsma and Vyn, 2008).
The alteration of plant physiology by AM is also confirmed with the better efficiency of Photosystem II, a highly sensitiveto-drought component from the plant photosynthetic apparatus (Ma et al., 1995). The highest effect of the AM symbiosis was shown under drought in the drought-sensitive genotype, with enhanced performance of PSII by 72% as compared to 36% enhancement in the drought-tolerant cultivar. This indicates that photochemical apparatus of droughted AM plants did not lose functionality in light conversion, that is the proportion of the light absorbed by chlorophyll associated with PSII, to reaction centers (Maxwell et al., 2000), as it was reported in other species under several stresses (Hajiboland et al., 2010;Porcel et al., 2015;Yooyongwech et al., 2016).
The percentage of membrane EL, an estimation of cell membrane stability, has been postulated as a good indicator of the tolerance to water stress (Ortiz et al., 2015). Accordingly, non-AM drought-tolerant plants had lower EL values than the corresponding drought-sensitive ones. In addition, whereas in the case of the tolerant genotype droughted AM inoculated plants showed steady state levels, in sensitive plants AM symbiosis helped to stabilize the membranes both under well-watered and under drought stress conditions. In this sense, the higher membrane stability is often related to lower MDA levels (Abid et al., 2016) accumulated as a result of lipid peroxidation. These results are in agreement with previous studies where MDA production was reduced by AM fungi . Furthermore, as expected, it is remarkable the similarity of results between EL and MDA concentration.
Plants need to maintain root osmotic potential below soil osmotic potential to take-up water. Previous studies have demonstrated that the AM fungi improve the plant osmotic adjustment by accumulation of different compounds (sugars, proline, free amino acids, etc.) (Bheemareddy and Lakshman, 2011;Sheng et al., 2011). This regulation by the AM symbiosis has been proposed as a mechanism allowing plants to grow under water stress (Ruiz-Lozano, 2003). In leaves of droughted plants, AM plants increased TSS in both genotypes, although to a lesser extent in the tolerant cultivar, suggesting an increased osmotic adjustment in AM plants during drought. The key effect of AM on sugar accumulation has been often reported under drought conditions (Wu et al., 2007;Zhang et al., 2014;Yooyongwech et al., 2016) as it is also shown here in maize plants from both sensitive and tolerant genotypes.
In this study, when plants were subjected to drought stress the AM symbiosis induced a higher improvement of physiological parameters in drought-sensitive plants than in drought-tolerant plants. These include efficiency of photosystem II, membrane stability, accumulation of soluble sugars and shoot and root DWs. Thus, drought-sensitive plants obtained greater physiological benefit from the AM symbiosis.

AM Regulation of Root Hydraulic Properties
Osmotic root hydraulic conductivity (Lo) can be considered as an estimation of water flow via the cell-to-cell pathway, and is highly related to the activity or density of water channels in the plasma membrane (Tyerman et al., 1999). A reduction in Lo is usually reported in plants exposed to water deprivation Aroca et al., 2012) probably as a mechanism for preventing water loss. This fact is consistent with our results, as a sharp drop in root hydraulic conductivity was observed in both genotypes when submitted to water stress. However, under drought stress the drought-tolerant genotype maintained a higher Lo values by 360% as compared to drought-sensitive genotype. Interestingly, AM increased Lo under drought compared to control plants in both genotypes, and this enhancement is in accordance with previous studies on AM plants under drought (Porcel et al., 2005;Bárzana et al., 2014;Sánchez-Romera et al., 2016). The increase of Lo in AM plants could be related to an increased expression of plant or fungal aquaporins (Sánchez-Romera et al., 2016). However, fungal aquaporins seem not to be involved in such increase since one gene was unaltered by drought, another gene was inhibited considerably in both maize cultivars, and the third one was only slightly induced in the drought-sensitive cultivar, but inhibited in the drought-tolerant one. The lack of a clear correlation between Lo and aquaporin gene expression ( Supplementary  Table S1) suggest that, the increase of Lo in AM plants may be due to additional mechanisms such as increased abundance and/or activity of the plants aquaporins due to post-translational modifications of these proteins (Chaumont and Tyerman, 2014) or to changes in density or size of plasmodesmata in AM roots (Blee and Anderson, 1998). Indeed, symplastic movement of water via plasmodesmata may also contribute significantly to Lo values (Galmés et al., 2007).
Aquaporin abundance in root cortex cells may alter Lo, especially during water shortage (Maurel et al., 2015), where aquaporins are thought to be regulated for the maintenance of the adequate water balance (Jang et al., 2007a,b). Among them, PIPs were proved to contribute to the adaptation of plants to drought episodes, also contributing to rehydration of the whole plant after water shortage . In addition to that, transcriptome analysis of drought tolerant and sensitive RILs in maize suggested that down-regulation of aquaporins is a mechanism contributing to the drought tolerance by upholding tissue turgor over longer time than drought-sensitive line (Min et al., 2016).
In the present study, 16 maize aquaporins previously shown to be regulated by the AM symbiosis under different drought scenarios (Bárzana et al., 2014) were analyzed to check a possible differential regulation by the AM symbiosis in two maize cultivars with contrasting drought sensitivity. We first observed that there were differences in the expression of several of the studied aquaporins between the drought-sensitive and the drought-tolerant genotypes. But these differences depended on the water regime and also on the presence or absence of the AM fungus. In the sensitive genotype, a general down-regulation of aquaporins by the AM symbiosis, under drought and/or wellwatered conditions (ZmPIP1;1, ZmPIP1;3, ZmPIP1;4, ZmPIP1;6, ZmPIP2;2, ZmPIP2;4, ZmTIP1;1, and ZmTIP2;3) was featured ( Figure 5 and Supplementary Figure S1). Similar result was also found in maize by Bárzana et al. (2014) and in other plant species Chitarra et al., 2016). However, AM regulation of aquaporins in the drought-tolerant genotype was weaker, and only three aquaporins (ZmPIP1;6, ZmPIP2;2, and ZmTIP4;1) were found to be altered. It is noteworthy that these three aquaporins were even up-regulated under well-watered conditions, which is an opposite behavior than in the sensitive genotype, similar to results reported by Liu et al. (2013) or Vinnakota et al. (2016) in two rice varieties and two Malus species with contrasting drought sensitivity. Also, upland rice and lowland rice with different responses to drought were compared to study the role of aquaporins in drought resistance and authors found important differences in PIP aquaporin transcriptional regulation in both types of rice (Lian et al., 2006).
During drought stress episodes, water conservation is critical for plant survival and productivity, and is achieved by an efficient uptake and stringently regulated water loss, in which aquaporins participate (Vinnakota et al., 2016). Our results in the drought-sensitive cultivar are in line with the hypothesis that down-regulation of aquaporins under water deprivation could be a way to minimize water loss, and the AM symbiosis could be helping the plant in this regulation. Through downregulation of aquaporin expression, roots from the droughtsensitive plants may be preventing drought damages by reducing water flow through cell membranes and upholding tissue turgor as a response to the soil water deficit Min et al., 2016). Indeed, dehydration avoidance during drought stress is a consequence of a tight balance between stomatal movements, root water uptake capacity and water distribution along plant tissues Ionenko et al., 2012). Nevertheless, the drought-tolerant genotype may not need this adjustment as other naturally occurring mechanisms such as deeper root development, improved turgor adjustment and photosynthetic efficiency or altered hormonal levels (Min et al., 2016) protected this genotype from the damage produced by drought.
It is also remarkable that under drought stress conditions ZmPIP1;1, ZmTIP1;1, and ZmPIP2;2 were downregulated by AM only in the drought-sensitive genotype. Among these three aquaporin genes, ZmTIP1;1 is the most expressed TIP in maize (Chaumont et al., 2001) and, besides water, it has the capacity to transport different compounds (urea, ammonia, boron, H 2 O 2 ) (Bárzana et al., 2014). ZmPIP2;2 showed a high water permeability (Pf) when expressed in Xenopus laevis oocytes (Bárzana et al., 2014). Thus, such tight regulation makes sense with a fine control of water balance in roots. Moreover, the specific regulation of these aquaporins by the AM symbiosis in the drought-sensitive cultivar point out a putative role of these three aquaporins in the AM-induced tolerance to drought stress, being possible targets for future studies.
In this sense, it must be taken into account that plant aquaporins can transport water, but also many other physiological substrates such as urea, glycerol, boric acid, silicic acid, hydrogen peroxide or gaseous molecules such as carbon dioxide, ammonia, or oxygen (Heinen et al., 2014;Li et al., 2014;Zwiazek et al., 2017). Among the different plant aquaporin subfamilies, NIPs is a versatile group with high diversity of substrates and a broad range of subcellular localization patterns (Maurel et al., 2008). Regulation of NIP genes by the AM symbiosis has been shown in different plant species such as Medicago truncatula (Uehlein et al., 2007), Lotus japonicus (Giovannetti et al., 2012), Zea mays (Bárzana et al., 2014) or Solanum lycopersicum (Chitarra et al., 2016). MtNIP1 had putative plasma membrane localization and was induced by mycorrhization. LjNIP1 was expressed in the inner membrane system of arbuscule-containing cells and could transport water. ZmNIP1;1 was shown to transport glycerol as well as silicon, while ZmNIP2;2 could transport silicon. LeNIP3;1 was overexpressed in AM tomato plants subjected to drought stress. Altogether, their transport capacities and localizations suggest that the regulation of NIP genes by the AM symbiosis could be involved in cell turgor regulation and in the exchange of water and solutes between both symbionts (Uehlein et al., 2007;Giovannetti et al., 2012;Bárzana et al., 2014;Chitarra et al., 2016), which may be of physiological importance to cope with drought stress.
Given the diversity of substrates that can be transported by plant aquaporins, those isoforms regulated by the AM symbiosis may have a role in regulation of leaf and root hydraulics, stomatal movement, nutrient uptake and translocation along plant tissues, carbon fixation or signaling processes. In this context, regulation of aquaporins having urea or ammonium transport capacity suggests that these aquaporins could be involved in the fungusbased nitrogen nutrition of the host plants or in plant nitrogen mobilization and metabolism (Bárzana et al., 2014), as was also proposed for ectomycorrhizal fungi (Dietz et al., 2011). Indeed, in the AM symbiosis, ammonium is suggested to be the major nitrogen compound transferred to the host plant, with urea playing a role as an intermediate solute (Tian et al., 2010). Studies by Gustavsson et al. (2005) suggested that export of plantderived glycerol may be important for symbiotic fungi. Thus, the regulation of plant aquaporins which can transport glycerol (i.e., ZmNIP1;1 and ZmTIP4;1) may be important for the AM symbiosis or for the plant-fungus interaction under drought stress conditions. Similarly, the regulation by the AM symbiosis of aquaporins with boron and/or silicon transport capacity could have structural functions in maize plants. Hydrogen peroxide is one of the most abundant reactive oxygen species continuously produced in the metabolism of aerobic organisms. As oxidant molecule, it reacts with various cellular targets causing cell damage, while at low concentration it acts as a signal molecule, controlling different essential processes in plants (Bienert et al., 2006). Thus, aquaporins with H 2 O 2 transport capacity such as ZmTIP1;1 could play a key role in the detoxification of excess H 2 O 2 generated under stress conditions, or in signaling events mediated by H 2 O 2 (Bárzana et al., 2014). That means that elucidating the in vivo transport capacities of the aquaporins regulated by the AM symbiosis is required to understand the role of these proteins in the AM-induced drought tolerance.

CONCLUSION
In summary, under water limiting conditions AM plants enhanced maize growth, especially in the case of the drought sensitive cultivar as reflected by the larger biomass (shoots and roots) accumulation. This beneficial effect of the AM symbiosis was linked to a better efficiency of PSII, to the higher membrane stability and to lower lipid peroxidation.
It is noteworthy that ZmPIP1;1, ZmPIP1;3, ZmPIP1;4, ZmPIP1;6, ZmPIP2;2, ZmPIP2;4, ZmTIP1;1, ZmTIP2;3, and ZmTIP4;1 gene expression was regulated by the AM symbiosis in the drought-sensitive cultivar, while in the drought-tolerant cultivar only ZmPIP1;6, ZmPIP2;2, and ZmTIP4;1 genes were regulated by the AM symbiosis. In the drought-sensitive cultivar, the genes ZmPIP1;1, ZmPIP2;2, and ZmTIP1;1 were down-regulated by the AM symbiosis when the plants were subjected to drought stress. Moreover, in this cultivar the genes ZmPIP1;3, ZmPIP1;4, ZmPIP1;6, ZmPIP2;4, ZmTIP2;3 were also down-regulated when the plants grew under well-watered conditions and only ZmTIP4;1 was up-regulated. In the drought tolerant cultivar the three genes regulated by the AM symbiosis were indeed up-regulated under well-watered conditions and ZmTIP4;1 was in addition up-regulated under drought stress. Thus, the broader and contrasting regulation of these aquaporins by the AM symbiosis in the drought-sensitive than the droughttolerant cultivar suggests a role of these aquaporins in water homeostasis or in the transport of solutes of physiological importance in both cultivars under drought stress conditions, which may be important for the AM-induced tolerance to drought stress. Grondin et al. (2016) found recently that a differential regulation of PIP aquaporins in six rice varieties was related to the drought stress tolerance of these varieties. Further research on the in vivo transport capacities by these aquaporins is needed to understand the specific role of these proteins in the AM-induced drought tolerance.

AUTHOR CONTRIBUTIONS
GQ and GE initiated the experiment and performed all the experimental measurements. JR-L designed the experiment and analyzed all the results. GQ, GE, and JR-L wrote the manuscript. RA and FC discussed the results and commented on the manuscript, and all the authors approved the final version of the article.

FUNDING
This work was financed by Ministerio de Economía y Competitividad (Project AGL2014-53126-R). GQ is a FPI Ph.D. Grant Holder.

ACKNOWLEDGMENTS
We thank DuPont Pioneer Corporation for providing us the maize seeds. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01056/ full#supplementary-material TABLE S1 | Pearson correlation coefficients between Lo and expression of the different maize aquaporin genes in a drought-sensitive and a drought-tolerant genotype.