# Dehydration Stress Memory: Gene Networks Linked to Physiological Responses During Repeated Stresses of Zea mays

Laetitia Virlouvet, Thomas J. Avenson, Qian Du, Chi Zhang, Ning Liu, Michael Fromm, Zoya Avramova, Sabrina E. Russo

► To cite this version:

Laetitia Virlouvet, Thomas J. Avenson, Qian Du, Chi Zhang, Ning Liu, et al.. Dehydration Stress Memory: Gene Networks Linked to Physiological Responses During Repeated Stresses of Zea mays. Frontiers in Plant Science, 2018, 9, 10.3389/fpls.2018.01058. [hal-02621904](https://hal.inrae.fr/hal-02621904)

HAL Id: hal-02621904

https://hal.inrae.fr/hal-02621904

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Dehydration Stress Memory: Gene Networks Linked to Physiological Responses During Repeated Stresses of Zea mays

Laetitia Virlouvet¹,², Thomas J. Avenson³, Qian Du⁴, Chi Zhang⁴, Ning Liu¹,⁵,
Michael Fromm¹, Zoya Avramova⁴ and Sabrina E. Russo⁴*

¹ Department of Agronomy and Horticulture, University of Nebraska, Lincoln, NE, United States, ² Institut Jean-Pierre Bourgin, INRA, AgroParisTech, CNRS, Université Paris-Saclay, Versailles, France, ³ LI-COR Biosciences, Lincoln, NE, United States, ⁴ School of Biological Sciences, University of Nebraska, Lincoln, NE, United States, ⁵ State Key Laboratory of Plant Genomics, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China

Stress memory refers to the observation that an initial, sub-lethal stress alters plants' responses to subsequent stresses. Previous transcriptome analyses of maize seedlings exposed to a repeated dehydration stress has revealed the existence of transcriptional stress memory in *Zea mays*. Whether drought-related physiological responses also display memory and how transcriptional memory translates into physiological memory are fundamental questions that are still unanswered. Using a systems-biology approach we investigate whether/how transcription memory responses established in the genome-wide analysis of *Z. mays* correlate with 14 physiological parameters measured during a repeated exposure of maize seedlings to dehydration stress. Co-expression network analysis revealed ten gene modules correlating strongly with particular physiological processes, and one module displaying strong, yet divergent, correlations with several processes suggesting involvement of these genes in coordinated responses across networks. Two processes key to the drought response, stomatal conductance and non-photochemical quenching, displayed contrasting memory patterns that may reflect trade-offs related to metabolic costs versus benefits of cellular protection. The main contribution of this study is the demonstration of coordinated changes in transcription memory responses at the genome level and integrated physiological responses at the cellular level upon repetitive stress exposures. The results obtained by the network-based systems analysis challenge the commonly held view that short-term physiological responses to stress are primarily mediated biochemically.

**Keywords:** gene co-expression networks, stress memory, drought, photosynthesis, chlorophyll fluorescence, non-photochemical quenching, stomatal conductance

## OPEN ACCESS

**Edited by:**
Bend Mueller-Roeber,
Universität Potsdam, Germany

**Reviewed by:**
Chuang Ma,
Northwest A&F University, China
Gaurav Sablok,
University of Helsinki, Finland

*Correspondence:
Sabrina E. Russo
srusso2@unl.edu

**Specialty section:**
This article was submitted to Plant Abiotic Stress, a section of the journal Frontiers in Plant Science

**Received:** 08 February 2018

**Accepted:** 28 June 2018

**Published:** 24 July 2018

**Citation:**
Virlouvet L, Avenson TJ, Du Q, Zhang C, Liu N, Fromm M, Avramova Z and Russo SE (2018) Dehydration Stress Memory: Gene Networks Linked to Physiological Responses During Repeated Stresses of Zea mays. Front. Plant Sci. 9:1058. doi: 10.3389/fpls.2018.01058

## INTRODUCTION

Understanding responses of plants to water limitation is an imperative for maintaining productivity in both agricultural and natural settings in the face of climate change (Porter and Semenov, 2005; Zhu et al., 2010; Lobell and Gourdji, 2012; Trenberth et al., 2014). These responses involve a range of physiological and biochemical mechanisms that have cost-benefit trade-offs
in a fluctuating environment (Kozlowski and Pallardy, 2002; Chaves et al., 2003). Although the closure of stomatal pores helps avoid water loss (Schulze, 1986), it decreases CO₂ and increases O₂ concentrations in the leaf airspace (Flexas et al., 2004, 2007). In addition to direct limitation of CO₂ availability to chloroplasts, drought stress also causes perturbations in photosynthetic metabolism and biochemistry that reduce carboxylation (Lawlor and Cornic, 2002; Chaves et al., 2003) and cause oxidative damage, which can ultimately lead to cell death (Chaves and Oliveira, 2004; Avenson et al., 2004; Baker et al., 2007; Lawlor and Tezara, 2009).

To prevent damage under hydraulic conditions unfavorable for carbon fixation, electron transfer can be shifted to alternative pathways (Ruuska et al., 2000) and protection involves a composite of processes collectively known as non-photochemical quenching (NPQ) (Müller et al., 2001). The predominant component of NPQ, energy dependent quenching, or *qE* (Crofts and Yerkes, 1994) harmlessly dissipates excessively absorbed energy as heat and is characterized by rapid reversibility driven by pH changes in the thylakoid lumen (Kramer et al., 1999; Müller et al., 2001; Mittler, 2002; Baier and Dietz, 2005). Regulation of *qE* has been shown to involve concentration changes in inorganic phosphate (*P*₁) (Takizawa et al., 2008) and to depend on ATP synthase in the chloroplast (Kohzuma et al., 2009). The integrated regulation of these processes is poorly understood, but is thought to be coordinated by hormones, such as abscisic acid and other metabolic and signaling pathways (Parent et al., 2009; Tardieu et al., 2010).

Plants often experience repeated sublethal water stresses, with intervening water-recovery periods over the course of a growing season (Slatyer, 1967). Their phenotypically plastic responses to such temporal variation presumably increase long-term survival and productivity despite repeated stress (Nicotra and Davidson, 2010; Cavanagh and Kubien, 2014). Several lines of empirical evidence support this conjecture. Pre-exposure to stressors alters responses to subsequent stressors (Bruce et al., 2007; Byun et al., 2014; Ding et al., 2014; To and Kim, 2014; Fleta-Soriano and Munné-Bosch, 2016).

Horticultural practices, such as drought hardening, enhance dehydration tolerance and photochemical efficiency under future unfavorable conditions by pre-exposing plants to sublethal stress (Turner, 2003; Byun et al., 2014). These observations suggest that plants have some form of “memory” that alters their responses to a subsequent stress (Bruce et al., 2007; Waters et al., 2008). Evidence of such memory has been demonstrated in *Arabidopsis thaliana* and *Zea mays*. Plants subjected to repeated cycles of dehydration stress, alternating with periods of full-watered recovery, exhibited both transcriptional and physiological memory responses, including reduced water loss rate, during a subsequent dehydration stress, compared to plants experiencing dehydration stress for the first time (Ding et al., 2012, 2014). However, it is not known whether other aspects of photosynthetic metabolism and biochemistry that change during dehydration stress also display memory, nor is it known how these responses are regulated by changes in gene expression to adjust the plant phenotype optimally to repeated cycles of stress and recovery.

While stress memory may be protective, there are likely costs associated with altered metabolism and lost opportunities for resource assimilation (Crisp et al., 2016). While such physiological adjustments are ultimately regulated by the expression of genes involved in dehydration stress response (Hayano-Kanashiro et al., 2009; Lu et al., 2011; Krishna-Reddy et al., 2014), we reason that, in a fluctuating environment, stress memory should act to avoid potentially catastrophic losses due to inadequate stress responses, but return to maximal functioning when favorable conditions resume. In repeated dehydration stress/rehydration cycles, two types of gene response groups were revealed in *A. thaliana* and *Zea mays* (Ding et al., 2012, 2013, 2014): genes that produced similar levels of transcripts in response to each stress, after returning to the initial pre-stressed levels during rehydration (non-memory genes) and genes that in a subsequent stress exhibited a transcriptional response significantly different from the response in the first stress (memory genes). It is not known how transcriptional memory is translated into an integrated physiological response to repeated dehydration stress. However, discovering these mechanisms is critical for understanding stress memory and identifying key genes and pathways to increase drought-resistance of agricultural crops.

Our framework (Figure 1) posits that the dehydration stress memory response may be thought of as a system of coordinated changes stimulated by an initial dehydration stress that precipitates altered gene expression that underlies, directly or indirectly, physiological responses to achieve a new state of homeostasis due to memory. Many biological processes are controlled by complex gene networks (D'haeseleer et al., 2000; Xiong and Zhu, 2001), and so our framework further posits that the activity of memory genes can be grouped into subsystems, each of which may respond to endogenously and exogenously driven changes in the plant to engender particular physiological responses. We used a systems approach to connect whole-genome transcriptional patterns to leaf physiological responses to repeated exposures to dehydration stress in maize. Specifically, our goals were: (1) to identify sets of genes that exhibit coordinated changes in gene expression under these conditions; (2) to establish whether physiological processes involved in photosynthesis and photoprotection display behavior consistent with memory when exposed to repeated dehydration stress with an intervening rehydration period; (3) to correlate transcriptional and physiological responses in order to identify potential candidates for key genes mediating dehydration stress memory; and (4) to identify putative functions of key genes that may point to particular physiological processes that are most strongly involved in dehydration stress memory. Co-expression network analysis was used to identify modules of co-expressed genes and to determine how they correlated with drought-related physiological parameters. Gene ontology analysis was used to examine whether presumed functions of genes in modules matched expectations based on the network correlations between gene expression and physiological parameters.
<figure><img src="image_3.png" /></figure>

# MATERIALS AND METHODS

Our materials and methods are described in detail in **Appendix S1** and summarized here.

## Gene Transcriptional Analysis

Data from the whole-genome transcriptome analyses of two biological replicates of maize plants exposed to three consecutive dehydration stresses (amounting to six whole-genome RNA seq analyses) were used for the network correlation analyses. These data were previously published in Ding et al. (2014). Genes that in the first stress (S1) and under repeated stresses (S2, S3) were similarly expressed were categorized as non-memory; genes that in S2 or S3 produced transcripts at significantly different levels

from the levels produced in S1 were categorized as memory, as described in detail below.

Bioinformatics analysis of transcriptome sequencing data of watered (W), S1, and S2 samples was performed for two biological replicates per sample. The distribution of raw and mapped reads for the two samples is available in Ding et al. (2014, Supplementary Files), and the quality of the replicates (for each of the three samples) was determined as described in detail therein. Genes were classified as being significantly differentially expressed when all three of the following conditions were met: $q \le 0.05$; $|log₂(fold change)| \ge 1$; and the FPKM-normalized expression value of at least one sample out of the two must be larger than the 25th percentile of all FPKM values within that sample. This additional check for significance (besides *q*-value
and fold change) filters out genes that have very low expression in both samples, despite having statistically significant differential expression. For all memory genes, the significant differences in expression were displayed between the S1 and S2 treatments, while the differences between S2 and S3 were not significant. These transcription patterns indicated the memory responses occurred after S1 and were fully displayed in S2. Consequently, the physiological responses measured here were in W, S1, R1, and S2, and transcriptional responses are reported for W, S1, R1, and S3.

## Classification of Stress Memory Response Types

According to our operational definition, a memory response implies that changes in gene expression or in a physiological parameter induced by repeated dehydration stresses (S2 or S3) are significantly different from responses displayed in the first stress (S1). Whole transcriptome data were from W, S1, and S3, whereas physiological data were collected at W, S1, R1, and S2. Transcript analysis of selected genes showed similar expression responses in S2 and S3, and so we considered these stages to be comparable for the transcript and physiological data. Among the 3740 memory genes in *Z. mays*, a subgroup of 2924 genes did not exhibit altered expression in S1 compared to the initial, non-stressed expression in W, but significantly changed (increased or decreased) expression in subsequent stresses. These genes were referred to as delayed memory response genes (Ding et al., 2014) and annotated as [>=-] or [>=+], where the [=] sign indicates the values in S1 were not statistically different from values in W; the second sign indicates the values in S2 were significantly different (lower [-] or higher [+]) than in S1. The other memory response-types are annotated as [+/+], [-/-], [+/-], and [-/+], where the first sign indicates significant difference in S1 versus W, and the second sign indicates significant difference in S2 versus S1. Non-memory responses are annotated as [+/=] or [-/=] indicating significant positive or negative response in S1 versus W but similar responses in S1 and S2.

## Co-expression Network and Gene Ontology Analyses

Co-expression network analysis was used to identify genes with coordinated transcriptional responses (modules). Among the total of 4986 *Z. mays* genes exhibiting memory and non-memory responses to repeated dehydration stress (Ding et al., 2014), 4710 genes displaying sufficiently high variation (range of FPKM-normalized read counts across a total of six samples for W, S1, and S2 with two biological replicates in each of the three treatments for a given gene >0.05) were identified, and their values were used to construct a correlation matrix using the R package, WGCNA (Langfelder and Horvath, 2008). Dynamic tree cutting was adopted to identify modules with minModuleSize of 30 (Langfelder and Horvath, 2007). Eigengenes were used to cluster all identified modules using Average Hierarchical Clustering analyses (Langfelder and Horvath, 2007). Pairwise distances between modules were calculated from the correlation

between eigengenes for each module as an estimate of similarity. A distance of 0.2 was used as a cut-off to group modules into meta-modules.

The correlation coefficients used to construct the co-expression network could be noisy due to our sample size of six for the transcriptomic data, although it has been demonstrated that overall the quality of the transcriptomic data set was good enough to generate reliable results, as shown in Ding et al. (2014). To increase the signal-to-noise ratio and reduce the false positive rate, during network construction, we used a high threshold to fill out genes with low expression levels or high fluctuations, and only genes with high correlation were considered to belong to the same module. This high threshold for data analysis reduced the false positives significantly and facilitates obtaining robust results.

Gene ontology (GO) (Ashburner et al., 2000; Bard and Rhee, 2004) and term enrichment analyses using AgriGO (Du et al., 2010) were performed to identify putative functions of dehydration stress memory genes in the network analysis. Genes sharing the same GO term are considered to be involved in similar biological processes, allowing us to identify the putative functions represented in each module. Statistical analysis was performed using Fisher's exact test with subsequent adjustment for multiple testing by calculating the false discovery rate (FDR) as implemented in AgriGO. The minimum number of mapping entries was set to 5, and the list of genes from the maize genome AGPv3 release 5a database was used as a reference<fnref n="1" />.

## Plant Growth and Treatments

All measurements were performed on 2-week old greenhouse-grown *Z. mays* L. (cultivar B73) seedlings. Handling of plants after removal from soil and dehydration stress-rehydration treatments were performed as previously described (Ding et al., 2014). Briefly, whole seedlings were air-dried for 90 min (first dehydration stress, S1) followed by 22 h of full rehydration with roots in water in covered trays (R1). R1 plants were again exposed to air-drying for 90 min (S2), followed by another rehydration treatment (R2). In order to avoid confounding wounding responses from cut tissue, different sets of plants were used for gas exchange and chlorophyll fluorescence, versus for the transcriptome, RWC, abscisic acid, and chlorophyll measurements, since the latter require leaf harvesting.

## Physiological Parameters

In each dehydration treatment stage, we measured several physiological parameters related to water stress (Table 1): photosynthetic gas exchange, pulse-amplitude-modulation (PAM) chlorophyll *a* fluorescence (Schreiber et al., 1995), foliar relative water content (RWC), and foliar abscisic acid (ABA) and chlorophyll concentrations. The rates of net photosynthesis (*P*<sub>N</sub>, µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup>), stomatal conductance (*g*<sub>s</sub>, mmol H<sub>2</sub>O m<sup>-2</sup> s<sup>-1</sup>), transpiration (*E*, mmol H<sub>2</sub>O m<sup>-2</sup> s<sup>-1</sup>), and PAM chlorophyll *a* fluorescence were measured simultaneously on the fourth leaf of seedlings at W, S1, R1, and S2 stages using a LI-6400XT infrared gas analyzer with a 6400-40 fluorescence

<footnote n="1">¹https://maizegdb.org/</footnote>
<table><caption>**TABLE 1** | Physiological parameters and their units estimated for Zea mays leaves during dehydration stress treatments.</caption><thead><tr><th>Parameters</th><th>Name</th><th>Units</th></tr></thead><tbody><tr><td>RWC</td><td>Relative water content</td><td>%</td></tr><tr><td>g<sub>s</sub></td><td>Stomatal conductance</td><td>mmol H<sub>2</sub>O m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>E</td><td>Transpiration</td><td>mmol H<sub>2</sub>O m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>C<sub>i</sub></td><td>Leaf internal CO<sub>2</sub> concentration</td><td>µmol CO<sub>2</sub> mol<sup>-1</sup></td></tr><tr><td>P<sub>N</sub></td><td>Photosynthetic rate</td><td>µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>A<sub>G</sub></td><td>Gross Photosynthetic rate</td><td>µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>ETR</td><td>Electron transport rate in PSII</td><td></td></tr><tr><td>φPSII</td><td>Quantum yield of photochemical energy conversion</td><td></td></tr><tr><td>F<sub>v</sub>′/F<sub>m</sub>'</td><td>Efficiency of the PSII open reaction centers</td><td></td></tr><tr><td>qP</td><td>Coefficient of photochemical quenching</td><td></td></tr><tr><td>φCO<sub>2</sub></td><td>Quantum yield of CO<sub>2</sub> assimilation</td><td></td></tr><tr><td>qE</td><td>Coefficient of energy-dependent quenching</td><td></td></tr><tr><td>φqE</td><td>Quantum yield of energy-dependent non-photochemical quenching</td><td></td></tr><tr><td>ABA</td><td>Foliar abscisic acid concentration</td><td>pmol g<sup>-1</sup> DW</td></tr><tr><td>Chlorophyll</td><td>Foliar chlorophyll concentration</td><td>mg g<sup>-1</sup> DW</td></tr></tbody></table>

chamber, with a uniform integrated LED light source and PAM fluorometer (LI-COR Inc., Lincoln, NE, United States). These parameters were simultaneously recorded, and a mean for each parameter for each of ten plants at each stage was estimated from 13 measurements at steady state, taken every 25 s over 5 min. In estimating chlorophyll fluorescence parameters (Baker et al., 2007), we used approaches described in the literature (Edwards and Baker, 1993; Lal and Edwards, 1996; Maxwell and Johnson, 2000; Kramer et al., 2004a). We assumed for φPSII (photosystem (PS) II operating efficiency) that the proportion of QA (reaction centers in the oxidized state) was between 0 and 1 (Genty et al., 1990) and, for Fv/Fm' (PSII maximum efficiency), that the proportion of QA was 1. We estimated two expressions of energy-dependent quenching based on different assumptions of PSII oxidation state, qE and φqE (Ahn et al., 2009), as well as the ratio, qE/ETR, a proxy expression describing a phenomenon termed qE sensitivity (Kramer et al., 2004b). We also calculated the composite parameters, φPSII/φCO₂ and ETR/A₉, which have been associated with drought responses (Edwards and Baker, 1993).

The foliar abscisic acid (ABA) assay was performed as in Liu et al. (2014). Four technical readings in each treatment were used to obtain a mean, which were averaged across four plants within the W, S1, R1, and S2 stages. Foliar chlorophyll content was measured at the W, S1, R1, and S2 stages. Leaf disks were immediately ground in liquid nitrogen, homogenized in 80% (v/v) acetone, and absorbance (A) was read at 663 and 645 nm. The total chlorophyll content (mg g⁻¹ DW) was calculated as [(8.02 × A663 nm) + (20.2 × A645 nm)]/DW, where DW is the dry weight of ground leaf material (Arnon, 1949). Two different sets of plants were analyzed on two dates, each of which had three replicate pools per treatment, with each pool containing the same tissue amount from five plants from a treatment.

## Analysis of Physiological Parameters and Correlation Between Transcriptional and Physiological Dehydration Stress Responses

To quantify changes in physiological parameters across dehydration treatments, we used Student's t-tests. Chlorophyll concentrations were analyzed using a linear mixed-effect model, with experiment date as a random intercept, and differences due to treatment analyzed using type III tests of fixed effects with degrees of freedom determined by the Satterthwaite method as implemented in R statistical software (R Core Development Team, 2011) using the packages *lme4* (Bates et al., 2014) and *lmerTest* (Kuznetsova et al., 2014).

The novel component of our study is that the previously published transcriptomic data (Ding et al., 2014) were integrated with the detailed plant physiological data using bioinformatic and gene expression network analyses in order to investigate how changes in gene expression correlated with changes in specific physiological processes during response to repeated drought stress. To do so, the eigengene of each module from the co-expression network analysis was correlated with the response levels across dehydration treatments for each physiological parameter as a measure of the level of coordination between genes in that module and physiological responses.

## RESULTS

### Gene Co-expression Networks

Among the 4986 maize genes responding to repeated dehydration stress (Ding et al., 2014), co-expression network analysis identified 11 gene modules; 10 of them were grouped into three clusters (meta-modules), each containing four, two, and four modules, plus one module (M11) that did not group with any other module and remained as a 'solo module' (Figure 2 and Supplementary Table S2). The co-expressed stress-responding genes in these differed broadly in predominant GO slim classification functions revealed in GO term analysis (Figure 3 and Supplementary Table S3).

### Meta-Module A Is Enriched in Genes Involved in Photosynthesis and Fluorescence Mechanisms

Meta-module A, represented by 2757 co-expressed genes, comprised four highly correlated modules, M1 (1588 genes), M4 (925 genes), M3 (205 genes), and M2 (39 genes) (Supplementary Table S2). Of these, 95% (2613 genes) exhibited a delayed memory response, including genes up- (1466 [=/+]) and genes down- (1147 [=/-]) regulated in S3. The remaining 5% (143 genes) showed other memory responses; none were non-memory genes. Meta-module A was enriched in genes encoding proteins implicated in electron transport and energy pathways, in pigment, carbohydrate and glycoside metabolic processes, and in response to biotic and abiotic stimuli including light (Figure 3 and Supplementary Table S3). While no enrichment in particular functions was identified for the
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Hierarchical clustering analysis of co-expressed genes in leaves of *Zea mays* seedlings in response to repeated dehydration stress and rehydration treatments. Three meta-modules (modules M1–M10) and one single module (module M11) were identified based on the dendrogram of 11 discovered modules identified by a hierarchical clustering analysis based on average linkage method for eigenvectors of these modules (eigengene). The height of branches is the distance between two corresponding eigenvectors, which is one minus their correlation coefficients. An eigenvector is the first principal component of the corresponding module, and the eigenvector characterizes the gene expression pattern for all genes in this given module. Two modules with lower distance in the dendrogram have higher similarity of gene expression patterns.

combined 244 genes of modules M3 and M2 (Supplementary Table S3), genes implicated in the CBB cycle, in light harvesting, electron transport, non-photochemical quenching, and in overall photosynthesis, were over-represented in modules M1 and M4 (Supplementary Table S3). The expression patterns of 94% of the M1 and M4 genes corresponded to delayed memory responses (Supplementary Table S2). Effectively, in module M1, 47 [>=-] delayed memory response genes were implicated in light harvesting, including proteins for PSI and PSII, chlorophyll A/B binding proteins, electron transport and cyclic electron flow proteins. In addition, 13 [>=-] delayed memory response genes encoding enzymes for the CBB cycle, including ribose-5-phosphate isomerase, glyceraldehyde-3-phosphate dehydrogenases A and B, phosphoglycerate kinases, phosphoribulokinase, and ribulose bisphosphate carboxylases were also clustered in M1 (Supplementary Table S2).

## Meta-Module B Is Enriched With House-Keeping Genes

Meta-module B containing modules M5 and M6 was defined by 357 (79%) co-expressed memory and 95 (21%) non-memory genes (Figure 2 and Supplementary Table S2). The meta-module was largely enriched in house-keeping genes and genes involved in signaling, cell wall organization, protein modification, and response to stimuli (Figure 3 and Supplementary Table S3).

## Meta-Module C Is Linked to Regulation of Diverse Processes

Meta-module C was defined by 1125 co-expressed genes clustered in modules M7, M8, M9, and M10 (Figure 2 and Supplementary Table S2). Module M7, containing 602 (93%) and 44 (7%) co-expressed non-memory and memory genes.

In general, meta-module C was enriched by genes implicated in transport, signaling, transcription, regulation of numerous biological processes and response to stimuli (Figure 3 and Supplementary Table S3).

## The Solo Module Is Enriched in Small Molecule Biosynthetic Process

The solo module (M11) grouped 241 co-expressed genes: 128 (52%) genes were categorized as memory, 36 (15%) genes as delayed memory, and 77 (32%) genes as non-memory response genes (Supplementary Table S2). The GO analysis indicated enrichment in small molecule biosynthetic process representing proteins involved in lipid metabolism, stress responses, and in producing intermediates of the C4 metabolic pathway, including allosteric inhibitors (Gonzalez et al., 1986) (Figure 3 and Supplementary Table S3).

## Changes in Physiological Parameters in Response to Repeated Dehydration Stress

Foliar RWC decreased significantly in S1 and S2, but returned to pre-stressed watered (W) levels after each rehydration (R1 and R2, Figure 4A), consistent with a non-memory ([-/=]) response. Foliar ABA also exhibited a non-memory response, but its [+/=] pattern was opposite to that of RWC (Figure 4B), as expected, since increases in ABA trigger closure of stomatal pores during dehydration stress (Lim et al., 2015).

The rates of stomatal conductance (gs) and transpiration (E) declined from W to S1, as expected, but did not return to levels observed in W after rehydration in R1, despite full restoration of RWC and a significant decline in ABA during R1 (Figures 4A,B). Both declined significantly to levels lower in S2 compared to S1
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Heatmap of gene ontology (GO) enrichment analysis of the meta-modules genes expressed in Z. mays under repeated dehydration stress and rehydration treatments. The GO terms enriched significantly (FDR < 0.01) for at least one meta-module have been grouped accordingly to the GO slim list established for *Arabidopsis thaliana*. See **Supplementary Table S3** for the entire list of the GO terms. The legend shows the color scaling with FDR values.</figcaption></figure>

(Figures 4D,E), exhibiting a stress memory response [—/—]. The net rate of photosynthetic CO₂ fixation ($P_N$) did not initially respond strongly to dehydration stress, and it was only by S2 that $P_N$ had declined to rates significantly lower than in W (Figure 4C), consistent with a delayed [=/—] memory response. The CO₂ concentration in the leaf internal airspace ($C_i$) declined precipitously from W to S1, as expected due to closure of stomatal pores (Figure 4F). During R1 and S2, $C_i$ remained at similar levels, displaying a non-memory [—/—] response, likely a result of the combined effect of reduced gs being buffered by reduced $P_N$ in S2.

All chlorophyll fluorescence parameters displayed delayed memory responses, as did the ratio of $qE/ETR$, a proxy expression for $qE$ sensitivity (Figures 5A–G). These parameters exhibited

their most dramatic changes from R1 to S2. The responses of $\varphi$PSII, ETR, and $qP$ were similar to that for $P_N$, in that they declined slightly from W to S1 and only became significantly different from levels at W by S2 ([=/-]) (Figures 5A,B,D), whereas $F_v'/F_m'$ and $\varphi$CO₂ displayed a more noticeable, but still gradual, decline from W to S1 (Figures 5C,F). The decline in $F_v'/F_m'$ and in $qP$ between R1 and S2 indicates enhanced NPQ that can account for the observed sharp decline in $\varphi$PSII. Consistent with this, parameters reflecting increases in the amount of energy dissipated via non-photochemical sinks, $qE$, $\varphi qE$, and $qE/ETR$ (Figures 5E,G) remained at similar levels from W to S1, but increased dramatically from S1 to S2 ([=/+]), whereas $qP$ declined. Foliar chlorophyll declined continuously from W to S2 (Figure 5H), displaying what may best be classified as a memory response ([-/—]). The composite parameters, $\varphi$PSII/$\varphi$CO₂ and ETR/$A_G$ also showed a delayed memory response ([=/+]) in that only levels at W and S2 were significantly different from each other, although they exhibited a slight increase from W to S1 that was not statistically significant (Figure 6).

## Changes in Gene Expression in Response to Repeated Dehydration Stress

Among the genes responding to repeated dehydration stress, a total of 164 genes encoding known dehydration stress-associated proteins were classified into four broad functional categories related to photosynthesis, pigments, stomatal regulation, and production of ABA. Most of them, 119 genes, belonged to the delayed memory category, 19 genes exhibited a memory response, and 26 were non-memory genes (Table 2).

Ninety four genes (79%) of the delayed memory response group encoded proteins associated with the light-dependent and carbon fixation reactions of photosynthesis (Table 2), including proteins implicated in electron transport in PSI and PSII (e.g., cytochrome $b_6f$ and quinone/ferredoxin), chlorophyll a/b binding, signaling and redox regulation, as well as proteins involved in the Calvin–Benson–Bassham (CBB) cycle (Table 2 and Supplementary Table S1). Delayed memory responses were displayed also by seven of the 13 genes associated with the metabolism of pigments and NPQ functions. In particular, two chlorophyll catabolic genes activated in S3 (delayed memory [=/+] genes), and two chlorophyll biosynthesis genes downregulated in S3 ([=/-] delayed memory genes) displayed patterns consistent with the decline in chlorophyll levels observed during the dehydration stresses (Figure 5G). Two genes implicated in the xanthophyll cycle, violaxanthin de-epoxidase (VDE, npq1) and zeaxanthin epoxidase (ZE, npq2) also displayed delayed memory, complementing each other's memory responses ([=/-] and [=/+], respectively; Supplementary Table S1). These changes in transcription suggest violaxanthin levels would increase, whereas zeaxanthin would decrease in subsequent stresses, creating an apparent paradox, as conversion of violaxanthin to zeaxanthin is involved in energy-dependent NPQ, which increased during dehydration stress. However, the gene coding for the $\gamma$ subunit of $CF_0CF_1$
<figure><img src="image_6.png" /></figure>

**FIGURE 4** | Changes in foliar relative water content (RWC), abscisic acid (ABA) concentration, and gas exchange parameters during successive dehydration stress and rehydration treatments in *Z. mays*. (A) RWC, (B) ABA concentration, (C) net photosynthesis rate (*P*<sub>N</sub>), (D) stomatal conductance (*g*<sub>s</sub>), transpiration rate (E), and (F) internal CO<sub>2</sub> concentration (*C*<sub>i</sub>) were measured under the well-watered condition (W), after the first (S1) and second (S2) dehydration stresses, and after rehydration (R1; R2 only for RWC). Values are means ± SE (*n* = 10, 10 and 4 plants per treatment for RWC, gas exchange parameters, ABA concentration, respectively). Different letters correspond to statistically significant differences between treatments (*P* < 0.05).

ATP synthase, displayed a [=-] memory response. Its reduced expression in S3, implying putative reductions in amounts of ATP synthase, would ultimately increase *qE* (Kohzuma et al., 2009) (*Supplementary Table S2*).

Seventeen genes (49%) of the 35 dehydration stress-responding genes implicated in the biogenesis, distribution, and movement of guard cells, including ion channel proteins that drive changes in stomatal apertures, belonged to the delayed memory category (Table 2). Of the remaining 18 genes, 13 exhibited non-memory and five exhibited other memory-type responses in S3 (Table 2). Of note is that members of the same gene families involved in ABA-dependent changes in stomatal aperture, displayed a broad spectrum of transcriptional responses, despite encoding proteins with similar structure and predicted biochemical functions (Supplementary Table S1). The increased transcription in S3 of the genes implicated in calcium-mediated signaling of ABA-dependent pathways and in ABA-signaled stomata closure (Mori et al., 2006; Vahisalu et al., 2008; Dreyer et al., 2012; Imes et al., 2013; Yu et al., 2014), suggested they would stimulate stomatal closing. In contrast, the expression of genes driving stomatal opening [i.e., inward potassium channel protein, KZM1 (Guo et al., 2002) and two nitrate transporters (Philippar et al., 2003) were significantly downregulated in S3].

Ten ABA metabolic genes were induced in S1; upon repeated stress, five displayed non-memory, and five displayed memory responses (Table 2). Among the latter, two genes encoding

ABA-anabolic and two genes encoding ABA-catabolic enzymes declined in expression (displaying a [+/-] response pattern), while one (for a catabolic abscisic acid 8'-hydrolase) was super-induced ([+/+]) (Table 2 and Supplementary Table S1). Since both ABA-synthesizing and ABA-degrading biochemical pathways were activated when plants experienced dehydration stress, the results suggest that the foliar ABA levels were determined by the combined effect of both ABA anabolism and catabolism. One ABA degrading gene displayed a [=/-] delayed memory response.

## Correlation Between Transcriptional and Physiological Responses

The strength of correlation between physiological parameters and the expression of genes in modules exhibited wide variation (Figure 7). Strong correlation of specific modules with particular physiological responses suggests that physiological changes involved in dehydration stress memory are associated with distinct networks of genes operating in a coordinated way.

Expression of genes in meta-module A correlated strongly with photosynthesis and chlorophyll fluorescence parameters (Figure 7). The constituent modules, however, differed in the strength of correlations with individual physiological parameters: module M1 correlated strongly with φPSII, ETR, P<sub>N</sub>, φqE, and F<sub>v</sub>'/F<sub>m</sub>', and less strongly with φCO<sub>2</sub>, qE, and qP; all correlations, except with qE and φqE, were positive. The M4, M3, and M2 modules also correlated with φPSII, ETR, and P<sub>N</sub>, as well as
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Changes in photosynthetic biochemistry and foliar chlorophyll concentration during successive dehydration stress and rehydration treatments in *Z. mays*. (A) φPSII, (B) electron transport rate (ETR), (C) Fv/Fm', (D) qP, (E) qE and φqE, (F) φCO₂, (G) qE/ETR, and (H) foliar chlorophyll concentration were measured under the well-watered condition (W), after the first (S1) and second (S2) dehydration stresses, and after rehydration (R1). Fluorescence parameters (A–G) were derived from the chlorophyll fluorescence measurements performed simultaneously with gas exchange measurements. See **Table 1** for abbreviations, definitions, and units of fluorescence parameters. Values are means ± SE (n = 10 and 6–7 plants per treatment for the chlorophyll fluorescence and concentration, respectively). Different letters correspond to statistically significant differences between treatments (*P* < 0.05).</figcaption></figure>
<figure><img src="image_8.png" /></figure>

**FIGURE 6** | Changes in chlorophyll fluorescence parameters during dehydration stress and rehydration treatments in *Z. mays*. (A) φPSII/φCO₂, and (B) electron transport rate (ETR)/gross photosynthetic rate (A₉). were measured under the well-watered condition (W), after the first (S1) and second (S2) dehydration stresses, and after rehydration (R1). Fluorescence parameters were derived from chlorophyll fluorescence measurements performed simultaneously with gas exchange measurements. See **Table 1** for abbreviations, definitions, and units of fluorescence parameters. Values are means ± SE (n = 10 plants per treatment). Different letters correspond to statistically significant differences between treatments (*P* < 0.05).

$qP$, but, in contrast to M1, the correlations were negative and weaker. However, no strong correlations of M5 and M6 modules of meta-module B with the physiological parameters measured here were found (Figure 7). Interestingly, only module M7 from meta-module C correlated negatively with $C_i$ and RWC and positively with ABA (Figure 7). Weaker correlations with $C_i$ and ABA were displayed by the M8 module, while no significant correlations with the responses studied here were found for modules M9 and M10. The solo module, M11, displayed negative correlations with stomatal conductance ($g_s$), transpiration (E), $\phi CO_2$, and chlorophyll (Chl) content, a positive correlation with $qE$, and a weaker positive correlation with $\phi qE$ (Figure 7). Thirty-one [-/-] and 97 [+/+] memory genes in the solo module correlated positively or negatively, respectively, with the [-/-] memory responses displayed by the $g_s$ and E parameters (Figures 4D,E).

DISCUSSION

Being sessile, plants cannot escape environmental variability, and so have evolved integrated, plastic responses to maintain fitness despite fluctuating biotic and abiotic stresses. Although these mechanisms are becoming better understood, much remains to be discovered regarding the role of stress memory. Genome-wide screening of dehydration stress-responding genes in combination with detailed physiological measurements revealed that widespread transcriptional changes are an important component of physiological dehydration stress memory. Our discovery challenges the commonly held view that short-term physiological responses to stress are primarily mediated biochemically. Moreover, we discovered that processes of known importance in

plant drought response displayed contrasting memory patterns: stomatal conductance declined with initial dehydration and was further reduced in the subsequent stress, whereas photosynthesis and NPQ showed delayed memory with significant changes only in the subsequent stress. Adjustments in translation may also occur, but, to our knowledge, ours is the first study to demonstrate that the sophisticated physiological memory responses to repeated dehydration stress are, apparently, driven by widespread transcriptional changes in genes implicated in diverse photosynthetic and photoprotective functions. While they clearly interact, the relative importance of expression-independent versus dependent mechanisms in dehydration stress memory is unknown, but likely depends on the timescale of fluctuations in water stress.

Dehydration-induced transcriptional changes were clearly organized in distinct modules of co-expressed genes and were correlated to varying degrees with photosynthetic and photoprotective responses to repeated stresses. It is important to note that the expression of genes that did not correlate with any physiological parameter measured here had mainly generalized, house-keeping-type functions. Most genes, however, showed strong, consistent correlations with particular physiological parameters, implying that some physiological responses involved in dehydration stress memory were associated with expression changes in distinct networks of interdependently operating genes. In contrast, the expression patterns of genes that displayed strong, but divergent, correlations with several physiological parameters suggested that they might be less integral to particular responses or might integrate signals across groups of co-responding genes. Although many of these genes have not been functionally annotated, their high transcriptional connectivity and correlations with particular physiological
<table><caption>**TABLE 2** | Dehydration stress response genes classified into four broad functional categories related to photosynthesis, pigments, stomatal regulation, and production of ABA.</caption><thead><tr><th rowspan="2">Category</th><th rowspan="2">No. genes</th><th colspan="8">Transcriptional Memory pattern</th></tr><tr><th>+[/+]</th><th>+[/-]</th><th>+[/=]</th><th>=[/+]</th><th>[-/-]</th><th>[-/+]</th><th>[-/=]</th><th>[=/-]</th></tr></thead><tbody><tr><td><strong>Photosynthesis</strong></td><td><strong>105</strong></td><td><strong>1</strong></td><td><strong>3</strong></td><td><strong>2</strong></td><td><strong>11</strong></td><td><strong>3</strong></td><td><strong>0</strong></td><td><strong>2</strong></td><td><strong>83</strong></td></tr><tr><td>Calvin-Benson-Basham Cycle</td><td>16</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>14</td></tr><tr><td>NADP-ME type</td><td>2</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>NAD-ME type</td><td>2</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>PEPC</td><td>2</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>PEPCK enzyme type</td><td>2</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>PEPC kinase</td><td>4</td><td>-</td><td>1</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Other photosynthesis proteins</td><td>4</td><td>-</td><td>-</td><td>1</td><td>2</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Photosystem I</td><td>11</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>11</td></tr><tr><td>Photosystem II</td><td>13</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>1</td><td>10</td></tr><tr><td>Cytochrome b6f</td><td>4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td></tr><tr><td>Quinone/ferredoxin</td><td>5</td><td>-</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td><td>2</td></tr><tr><td>Chlorophyll a/b binding protein</td><td>16</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>16</td></tr><tr><td>Other light photosynthesis proteins</td><td>5</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>4</td></tr><tr><td>Electron transport</td><td>8</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>8</td></tr><tr><td>Nonphotochemical regulation</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Photosynthesis regulation</td><td>10</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>9</td></tr><tr><td><strong>Pigment</strong></td><td><strong>13</strong></td><td><strong>0</strong></td><td><strong>2</strong></td><td><strong>3</strong></td><td><strong>3</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>1</strong></td><td><strong>4</strong></td></tr><tr><td>Chlorophyll biosynthesis</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>Chlorophyll degradation</td><td>4</td><td>-</td><td>1</td><td>-</td><td>2</td><td>-</td><td>-</td><td>1</td><td>-</td></tr><tr><td>Carotene</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Zeaxanthin biosynthesis</td><td>4</td><td>-</td><td>1</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Zeaxanthin/violaxanthin interconversion</td><td>2</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td><strong>Stomatal regulation</strong></td><td><strong>35</strong></td><td><strong>2</strong></td><td><strong>2</strong></td><td><strong>12</strong></td><td><strong>12</strong></td><td><strong>0</strong></td><td><strong>1</strong></td><td><strong>1</strong></td><td><strong>5</strong></td></tr><tr><td>Channel</td><td>6</td><td>-</td><td>1</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>Stomatal regulation</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Stomatal movement</td><td>28</td><td>2</td><td>1</td><td>9</td><td>12</td><td>-</td><td>1</td><td>1</td><td>2</td></tr><tr><td><strong>ABA</strong></td><td><strong>11</strong></td><td><strong>1</strong></td><td><strong>4</strong></td><td><strong>5</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>0</strong></td><td><strong>1</strong></td></tr><tr><td>ABA biosynthesis</td><td>6</td><td>-</td><td>3</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>ABA degradation</td><td>5</td><td>1</td><td>1</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td><strong>Total genes</strong></td><td><strong>164</strong></td><td><strong>4</strong></td><td><strong>11</strong></td><td><strong>22</strong></td><td><strong>26</strong></td><td><strong>3</strong></td><td><strong>1</strong></td><td><strong>4</strong></td><td><strong>93</strong></td></tr></tbody></table>

parameters allow us to map how they interact to mediate the complex phenotype of dehydration stress memory.

## Memory Responses Limiting Gas Exchange and Protecting the Photosynthetic Apparatus

Previous research examining dehydration stress memory revealed that the fresh weight of maize seedlings declined more slowly in subsequent, relative to the first, dehydration stresses (Ding et al., 2014). Water loss through transpiration is governed by ABA-mediated changes in stomatal apertures (Schulze, 1986; Lim et al., 2015). In *Aptenia cordifolia*, a CAM plant, ABA showed differences in double-stressed plants compared with plants challenged with drought for the first time (Fleta-Soriano et al., 2015). Despite the presence of multiple copies of ABA-metabolic genes in the maize genome, only half of them participated in dehydration stress responses,

and the observed transcription patterns suggested a complex interaction among ABA-metabolizing activities. The lack of super-induced ABA levels in successive stress cycles rules out a model wherein retention or increase of ABA was responsible for the super-induced transcript levels of other memory genes or for the memory responses of $g_s$ and E. The memory response of conductance is likely due to a mechanism maintaining reduced stomatal apertures even after recovery, as also found for *Arabidopsis* in an analogous memory response (Virlouvet and Fromm, 2015). Moreover, the recovery of conductance may decrease or require longer times as the severity of the dehydration stress increases (Miyashita et al., 2005; Efeoğlu et al., 2009).

Drought stress is typically described as resulting in diminished supplies of CO$_2$ to chloroplasts due to reduced stomatal conductance, thereby causing diminished carboxylative demand and energy production, as evidenced by concurrent decreases in CO$_2$ assimilation and ETR (Lal and Edwards, 1996;
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | Correlations among gene co-expression modules and physiological parameters in Z. mays exposed to repeated dehydration stress and rehydration treatments. Each row corresponds to a module, and each column to a physiology parameter. One cell indicates the correlation between the eigenvector of this given module and the corresponding physiology parameter. In one cell, the corresponding *p*-value is shown below the correlation coefficient. The colors of cells, changing gradually from blue to orange to red, indicate the significance of correlations among module eigenvectors and physiological parameters, with the legend showing the correspondence of colors to −log₁₀(*p*-values).</figcaption></figure>

Medrano et al., 2002). Such changes could result from purely biochemical mechanisms; however, the transcript and GO term analyses showed dramatic evidence for delayed memory responses of many genes involved in photosynthesis. The reduced stomatal conductance between R1 and S2, apparently, did not cause a decrease in either carboxylative or energy production capacities *per se*, that is, by reductions in chloroplastic levels of CO₂, since Cᵢ levels remained constant between R1 and S2. Instead, transcriptional changes, including diminished expression of CBB and electron transfer-associated proteins, were elicited that were associated with a coordinated decrease in carboxylative and energy producing capacities, as shown by as shown by the reduction in both Pₙ and *ETR*. The non-memory response of Cᵢ is thus likely due to compensatory changes in stomatal conductance and the coupled carboxylation and *ETR* reactions that coordinate CO₂ assimilation. The diminished carboxylative and energy-producing capacities were accompanied by transcriptional-level decreases in energy-capture capacity, as suggested by down-regulation of light harvesting (i.e., chlorophyll-binding) proteins and chlorophyll biosynthesis enzymology, ultimately reducing foliar chlorophyll.

A novel aspect of this diminished energy capture-capacity is the first evidence supporting the hypothesis that the multifaceted photoprotective responses involved in dehydration stress memory are explicitly controlled at the transcriptional level. To maximize CO₂ assimilation while minimizing photodamage in a fluctuating environment, plants must flexibly adjust *qE*, the predominant component of NPQ (Avenson et al., 2005), and this ability has recently been shown to increase crop productivity (Kromdijk et al., 2016). In studies of drought (Avenson et al., 2005; Kohzuma et al., 2009; Kanazawa et al., 2017), biochemical perturbation of the carboxylative reactions was accompanied by diminished *ETR*, and yet *qE* remained as, if not more, robust in comparison to control conditions. In effect, *qE* became more sensitive to *ETR* under biochemical perturbation of the carboxylative reactions. In a biochemical model for the mechanisms underlying drought-induced increases in *qE* sensitivity, chloroplastic ATP synthase emerged as a central player (Kramer et al., 2004a). This model, however, did not invoke any transcription-mediated phenomena. It proposed that the relative resistance of chloroplastic ATP synthase to proton efflux increases via kinetic constraints imposed by
inorganic phosphate (P₁) limitation (Takizawa et al., 2008) and/or diminished levels of ATP synthase (Kohzuma et al., 2009). The resulting reduced proton efflux lowers the pH of the thylakoid lumen causing: (1) the rate of electron transfer through the cytochrome b6f complex, and consequently *ETR*, to slow down so as to match the diminished carboxylative capacity and (2) the pH-dependent mechanisms controlling the magnitude of *qE* (i.e., protonation of Psbs and activation of VDE) to be enhanced. The net result is that lower fluxes of *ETR* are capable of maintaining significant levels of *qE*. However, the reduced levels of ATP synthase were not explicitly attributed to translational and/or transcriptional phenomena (Kohzuma et al., 2009).

Here, we demonstrate that this well-known response of dynamic *qE* sensitivity is a key component of dehydration stress memory. Moreover, we provide the first evidence that drought-induced *qE* sensitivity modulation involves gene transcriptional changes. Specifically, we observed reduced expression and a delayed memory response of the gene coding for the γ subunit of CF₀CF₁ ATP synthase, the inactivation of which has been shown to result in high non-photochemical quenching in *Arabidopsis* (Dal Bosco et al., 2004). Transcriptional changes in the maize genome were coordinated with attenuation of *ETR* and enhancement of *qE*. Therefore, it can be inferred that the relative resistance of ATP synthase increased between R1 and S2, likely due to reductions in the amount of the enzyme, thereby effecting the observed increase in *qE* sensitivity. It is of significant interest to determine whether the gene networks controlling levels of ATP synthase share regulatory effectors with other gene networks involved in the dehydration stress memory response.

Why photoprotective responses showed delayed memory and were not induced immediately in the first stress may be partly due to C4-mediated carbon-concentration (Rao and Dixon, 2016). Since Pₐ did not decline dramatically until S2, carboxylation was still occurring in S1, presumably with CO₂ previously concentrated in the bundle sheath cells. Further accumulation of CO₂ was presumably minimal, since stomatal conductance was reduced in S1 and R1, and was reduced even further in S2, when photoprotective responses increased. Comparisons with C3 species would be instructive, as their Pₐ may be more likely to display a memory response to repeated dehydration stress. Consistent with this prediction is the observation that the most dramatic difference in the nature of the dehydration stress memory genes in maize versus the C3 eudicot, *A. thaliana*, was that photosynthesis-related genes in *Arabidopsis* displayed [-/-] and [-/+], memory responses, whereas the homologous genes in maize belonged mainly to the delayed [-/=] memory category (Ding et al., 2014).

Changes in other chlorophyll fluorescence parameters give insights into the physiological integration of the dehydration stress memory response. The decrease in φPSII, and φCO₂, between R1 and S2 can be explained by the concomitant increase in *qE* and decrease in *qP*. The delayed memory pattern of these parameters suggests that closure of PSII reaction centers occurred only after a severe dehydration stress and/or after induction of dehydration stress memory. Consistent with our findings, a study of drought memory over the course of a growing season in a C3 grass species also found enhanced photoprotection of

doubly stressed plants (Walter et al., 2011). The coordinated decline in F_v/F_m', a measure of the efficiency of PSII if all reaction centers were oxidized (Baker et al., 2007), suggests that the maximal capacity of PSII electron transport was also affected by stress memory, possibly due to damage of PSII reaction centers (Genty et al., 1989; Tyystjärvi and Aro, 1996) or migration of light-harvesting antennae from PSII to PSI (Genty et al., 1990). The increases in *ETR/A*<sub>G</sub> and φPSII/φCO₂ suggest that by S2, more energy was being dissipated by non-carboxylative mechanisms (Edwards and Baker, 1993; Loriaux et al., 2013).

All of the photosynthesis and fluorescence parameters exhibited a delayed memory response pattern that correlated remarkably well with the transcriptional response patterns of implicated genes. Thus, meta-module A, representing 2757 co-expressed delayed memory response genes, encoded proteins involved in light harvesting, electron transport, non-photochemical quenching, and in overall photosynthesis including enzymes for the CBB cycle. Genes known to be involved in NPQ, ZE (npq2) and VDE (npq1) (Niyogi et al., 1998), were found in both the M1 and the solo modules, both of which correlated with *qE* and φqE. All else being equal, the increased expression of zeaxanthin epoxidase (npq2) and decreased expression of violaxanthin de-epoxidase (npq1) would presumably cause higher violaxanthin levels and lower zeaxanthin levels in subsequent stresses. Although zeaxanthin is necessary for NPQ (Demmig-Adams and Adams, 2000; Adams et al., 2006), it is not sufficient: both the pH gradient across the thylakoid membrane (Müller et al., 2001) and the PSII protein, Psbs (Li et al., 2000), are also necessary (Takizawa et al., 2008). We found no changes in the expression of Psbs. However, the apparent paradox of putatively lower npq1 but higher *qE* in S2, could be explained by a smaller pool of npq1 that is more highly activated by a reduced pH in the thylakoid lumen (Kramer et al., 1999). The relationship between *qE* and zeaxanthin is likely to be more complex under drought. In support, other studies have also found transcript levels of npq1 and npq2 not to correlate well with short-term responses in these pigments (North et al., 2005).

CONCLUSION

Current models posit that short-term physiological responses of plants to environmental stress are predominantly mediated by biochemical changes. Our findings, however, demonstrate that changes in gene expression are also an important component of the stress response mechanism and that short-term photosynthetic and photoprotective responses to repeated dehydration stresses are partially controlled at the transcriptional level. The dehydration stress memory response in maize involves multiple physiological mechanisms that act to protect the photosynthetic machinery while maximizing carbon gain under fluctuating water availability. Our study is novel in that we quantified these mechanisms in detail and in coordination with foliar transcriptomic data from multiply stressed/rehydrated seedlings, allowing for new insights into the organization of the transcriptional and physiological responses involved. Moreover, that a gene for chloroplastic ATP synthase displayed
downregulated delayed memory, to our knowledge, provides the first evidence that drought-induced modulation of *qE* sensitivity involves a transcriptional component. That both gene expression-independent and dependent mechanisms drive this and other responses is a novel revelation of how the photosynthetic and photoprotective processes occurring in dehydration stress memory responses are regulated.

Establishing whether genes grouped in a particular meta-module share specific transcription factors' (TFs) binding motifs might provide additional insights into their memory patterns. However, data on the TFs and their binding sights in the maize genome are currently sparse, which limits analyses of their potential enrichment within a module. Extrapolation from data regarding TF binding sites in other species is impractical, as the grass cistrome is highly divergent from that of *A. thaliana*, and species-specific differences in the response of photosynthesis genes may be significant, even between grass species (Kümpers et al., 2017). Although critical for the transcription of regulated genes, the transcriptional behavior of a TF, alone, cannot explain, or predict, their memory patterns. For example, the memory expression pattern of the TF MYC2 under repeated stress did not correlate with the non-memory expression pattern of its target marker gene RD22 (Liu et al., 2014), and the (non-memory) expression patterns of the TFs, ABFs, does not correlate with the super-induced memory patterns of their directly regulated RD29B or RAB18 genes (Virlouvet et al., 2014). The roles of the TFs in the memory transcription of their directly regulated targets may need to be established for each case.

Why some genes and responses in these networks showed memory versus delayed memory is unknown, but may reflect cost-benefit trade-offs. For example, fewer genes were involved in the memory response of stomatal regulation, compared to substantially more that were involved in delayed memory responses of photosynthetic and photoprotective processes. Avoiding water loss by changing regulation of relatively fewer genes involved in the memory response of stomatal regulation may have lower costs and immediate benefits, compared to altering expression of substantially more genes involved in photosynthesis, which may require more energy and also come at the cost of an altered photosynthetic apparatus and reduced carboxylation capacity.

Our study demonstrates that stress memory contributes to environmental acclimation, considered plastic and reversible, but, depending on the mechanisms regulating changes in gene expression, it could also contribute to heritable variation in stress responses (Hauser et al., 2011). However, the transcription memory of several tested memory genes lasted 5-7 days after removal of the dehydration stress (Ding et al., 2012), suggesting that the dehydration stress memory is short-term and reversible, unlikely to be preserved *trans*-generationally, and functioning as a mechanism for transcriptional responses during recurring bouts of drought. By temporarily altering expression of specific genes and related physiological parameters, it may facilitate an individual plant's acclimation and survival under repeated seasonal drought spells. Further studies will provide more understanding of

how gene regulatory and physiological networks interact to produce phenotypically integrated drought memory responses in plants, which is fundamental to developing agricultural crops tuned to the rapid environmental changes occurring across the globe.

## AVAILABILITY OF DATA

The physiological data used in this study are available from the corresponding author on reasonable request. The raw transcriptome sequence files for watered, S1, and S3 have been uploaded, together with gene expression result files, to NCBI's Gene Expression Omnibus under sequence number GSE48507.

## AUTHOR CONTRIBUTIONS

ZA, MF, and SR designed the study and supervised data collection. LV and NL collected the physiological data. SR, TA, and LV analyzed and interpreted the physiological data. ZA and LV conducted the GO-term analyses. CZ and QD conducted the network analyses. LV, ZA, SR, TA, and CZ wrote the manuscript.

## FUNDING

This study was supported by the University of Nebraska (UNL) Life Sciences Competitive Grants Program in Enhancing Interdisciplinary Teams and NSF award MCB-1121898. This work was also supported by the Nebraska Soybean Board (Award #1728) and the National Science Foundation (Awards OIA-1557417) to CZ. Part of this work was conducted while SR was on Faculty Development Leave from UNL.

## ACKNOWLEDGMENTS

We thank Jean-Jack Marco Riethoven for assistance with bioinformatic analyses and Sylvie Coursol for comments on the study.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01058/full#supplementary-material

**APPENDIX S1** | Supplementary materials and methods.

**TABLE S1** | List of genes in each module and their degrees in the co-expression network.

**TABLE S2** | Gene ontology term enrichment analysis for the genes from each meta-module and module of the co-expression network analysis.

**TABLE S3** | Classification of memory genes into four broad functional categories related to photosynthesis, pigments, stomatal regulation, and production of abscisic acid.
REFERENCES

Adams, W., Zarter, C., Mueh, K., Amiard, V., and Demmig-Adams, B. (2006). "Energy dissipation and photoinhibition: a continuum of photoprotection," in *Photoprotection, Photoinhibition, Gene Regulation, and Environment*, eds B. Demmig-Adams, W. Adams, and A. Mattoo (Dordrecht: Springer), 49–64. doi: 10.1007/1-4020-3579-9_5

Ahn, T. K., Avenson, T. J., Peers, G., Li, Z., Dall'Osto, L., Bassi, R., et al. (2009). Investigating energy partitioning during photosynthesis using an expanded quantum yield convention. *Chem. Phys.* **357**, 151–158. doi: 10.1016/j.chemphys. 2008.12.003

Arnon, D. I. (1949). Copper enzymes in isolated chloroplasts. polyphenoloxidase in *Beta vulgaris*. *Plant Physiol.* **24**, 1–15. doi: 10.1104/pp.24.1.1

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. *Nat. Genet.* **25**, 25–29. doi: 10.1038/75556

Avenson, T. J., Cruz, J. A., Kanazawa, A., and Kramer, D. M. (2005). Regulating the proton budget of higher plant photosynthesis. *Proc. Natl. Acad. Sci. U.S.A.* **102**, 9709–9713. doi: 10.1073/pnas.0503952102

Avenson, T. J., Cruz, J. A., and Kramer, D. M. (2004). Modulation of energy-
dependent quenching of excitons in antennae of higher plants. *Proc. Natl. Acad.
Sci. U.S.A.* **101**, 5530–5535. doi: 10.1073/pnas.0401269101

Baier, M., and Dietz, K.-J. (2005). Chloroplasts as source and target of cellular redox regulation: a discussion on chloroplast redox signals in the context of plant physiology. *J. Exp. Bot.* **56**, 1449–1462. doi: 10.1093/jxb/eri161

Baker, N. R., Harbinson, J., and Kramer, D. M. (2007). Determining the limitations and regulation of photosynthetic energy transduction in leaves. *Plant Cell Environ.* **30**, 1107–1125. doi: 10.1111/j.1365-3040.2007.01680.x

Bard, J. B. L., and Rhee, S. Y. (2004). Ontologies in biology: design, applications and future challenges. *Nat. Rev. Genet.* **5**, 213–222. doi: 10.1038/nrg1295

Bates, D., Maechler, M., Bolker, B., and Walker, S. (2014). *Ime4: Linear Mixed-Effects Models Using Eigen and S4. R Package Version 1.1-7*. Available at: http://cran.r-project.org/package=Ime4

Bruce, T. J., Matthes, M. C., Napier, J. A., and Pickett, J. A. (2007). Stressful “memories” of plants: evidence for possible mechanisms. *Plant Sci.* **173**, 603–608. doi: 10.1016/j.plantsci.2007.09.002

Byun, Y.-J., Koo, M.-Y., Joo, H.-J., Ha-Lee, Y.-M., and Lee, D.-H. (2014). Comparative analysis of gene expression under cold acclimation, deacclimation and reacclimation in Arabidopsis. *Physiol. Plant*. **152**, 256–274. doi: 10.1111/ ppl.12163

Cavanagh, A. P., and Kubien, D. S. (2014). Can phenotypic plasticity in Rubisco performance contribute to photosynthetic acclimation? *Photosynth. Res.* **119**, 203–214. doi: 10.1007/s11120-013-9816-3

Chaves, M. M., Maroco, J. P., and Pereira, J. S. (2003). Understanding plant responses to drought &#821221; from genes to the whole plant. * Funct. Plant Biol.* **30**, 239–264. doi: 10.1071/FP02076

Chaves, M. M., and Oliveira, M. M. (2004). Mechanisms underlying plant resilience to water deficits: prospects for water-saving agriculture. *J. Exp. Bot.* **55**, 2365–2384. doi: 10.1093/jxb/erh269

Crisp, P. A., Ganguly, D., Eichten, S. R., Borevitz, J. O., and Pogson, B. J. (2016). Reconsidering plant memory: Intersections between stress recovery, RNA turnover, and epigenetics. *Sci. Adv.* **2:e1501340**. doi: 10.1126/sciadv.150 1340

Crofts, A. R., and Yerkes, C. T. (1994). A molecular mechanism for qE-quenching. *FEBS Lett.* **352**, 265–270. doi: 10.1016/0014-5793(94)00976-7

Dal Bosco, C., Lezhneva, L., Biehl, A., Leister, D., Strotmann, H., Wanner, G., et al. (2004). Inactivation of the chloroplast ATP synthase γ subunit results in high non-photochemical fluorescence quenching and altered nuclear gene expression in *Arabidopsis thaliana*. *J. Biol. Chem.* **279**, 1060–1069. doi: 10.1074/ jbc.M308435200

Demmig-Adams, B., and Adams, W. W. (2000). Photosynthesis: harvesting sunlight safely. *Nature* **403**, 371–374. doi: 10.1038/35000315

D'haeseleer, P., Liang, S., and Somogyi, R. (2000). Genetic network inference: from co-expression clustering to reverse engineering. *Bioinformatics* **16**, 707–726. doi: 10.1093/bioinformatics/16.8.707

Ding, Y., Fromm, M., and Avramova, Z. (2012). Multiple exposures to drought 'train' transcriptional responses in Arabidopsis. *Nat. Commun.* **3**:740. doi: 10.1038/ncomms1732

Ding, Y., Liu, N., Virlouvet, L., Riethoven, J. J., Fromm, M., and Avramova, Z. (2013). Four distinct types of dehydration stress memory genes in *Arabidopsis thaliana*. *BMC Plant Biol.* **13**:229. doi: 10.1186/1471-2229-13-229

Ding, Y., Virlouvet, L., Liu, N., Riethoven, J.-J., Fromm, M., and Avramova, Z. (2014). Dehydration stress memory genes of *Zea mays*: comparison with *Arabidopsis thaliana*. *BMC Plant Biol.* **14**:141. doi: 10.1186/1471-2229-14-141

Dreyer, I., Gomez-Porras, J. L., Riaño-Pachón, D. M., Hedrich, R., and Geiger, D. (2012). Molecular evolution of slow and quick anion channels (SLACs and QUACs/ALMTs). *Front Plant Sci.* **3**:263. doi: 10.3389/fpls.2012.00263

Du, Z., Zhou, X., Ling, Y., Zhang, Z., and Su, Z. (2010). agriGO: a GO analysis toolkit for the agricultural community. *Nucleic Acids Res.* **38**, W64–W70. doi: 10.1093/nar/gkq310

Edwards, G. E., and Baker, N. R. (1993). Can CO2 assimilation in maize leaves be predicted accurately from chlorophyll fluorescence analysis? *Photosynth. Res.* **37**, 89–102. doi: 10.1007/BF02187468

Efeoğlu, B., Ekmekçi, Y., and Çiçek, N. (2009). Physiological responses of three maize cultivars to drought stress and recovery. *S. Afr. J. Bot.* **75**, 34–42. doi: 10.1016/j.sajb.2008.06.005

Fleta-Soriano, E., and Munné-Bosch, S. (2016). Stress memory and the inevitable effects of drought: a physiological perspective. *Front Plant Sci.* **7**:143. doi: 10.3389/fpls.2016.00143

Fleta-Soriano, E., Pintó-Marijuan, M., and Munné-Bosch, S. (2015). Evidence of drought stress memory in the facultative CAM, *Aptenia cordifolia*: possible role of phytohormones. *PLOS One* 10:e0135391. doi: 10.1371/journal.pone.0135391

Flexas, J., Bota, J., Loreto, F., Cornic, G., and Sharkey, T. D. (2004). Diffusive and metabolic limitations to photosynthesis under drought and salinity in C3 plants. *Plant Biol.* **6**, 269–279. doi: 10.1055/s-2004-820867

Flexas, J., Diaz-Espejo, A., Galmés, J., Kaldenhoff, R., Medrano, H., and Ribas-Carbo, M. (2007). Rapid variations of mesophyll conductance in response to changes in CO2 concentration around leaves. *Plant Cell Environ.* **30**, 1284–1298. doi: 10.1111/j.1365-3040.2007.01700.x

Genty, B., Briantais, J.-M., and Baker, N. R. (1989). The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. *Biochim. Biophys. Acta* **990**, 87–92. doi: 10.1016/ S0304-4165(89)80016-9

Genty, B., Wonders, J., and Baker, N. R. (1990). Non-photochemical quenching of Fo in leaves is emission wavelength dependent: consequences for quenching analysis and its interpretation. *Photosynth. Res.* **26**, 133–139. doi: 10.1007/ BF00047085

Gonzalez, D. H., Iglesias, A. A., and Andreo, C. S. (1986). Active-site-directed inhibition of phosphoenolpyruvate carboxylase from maize leaves by bromopyruvate. *Arch Biochem. Biophys.* **245**, 179–186. doi: 10.1016/0003- 9861(86)90203-1

Guo, F.-Q., Young, J., and Crawford, N. M. (2002). The nitrate transporter AtNRT1.1 (CHL1) functions in stomatal opening and contributes to drought susceptibility in Arabidopsis. *Plant Cell* **15**, 107–117. doi: 10.1105/tpc.006312

Hauser, M. T., Aufsatz, W., Jonak, C., and Luschnig, C. (2011). Transgenerational epigenetic inheritance in plants. *Biochim. Biophys. Acta* **1809**, 459–468. doi: 10.1016/j.bbagrm.2011.03.007

Hayano-Kanashiro, C., Calderon-Vazquez, C., Ibarra-Laclette, E., Herrera-Estrella, L., and Simpson, J. (2009). Analysis of gene expression and physiological responses in three Mexican maize landraces under drought stress and recovery irrigation. *PLoS One* **4:e7531**. doi: 10.1371/journal.pone.0007531

Imes, D., Mumm, P., Böhm, J., Al-Rasheid, K. A. S., Marten, I., Geiger, D., et al. (2013). Open stomata 1 (OST1) kinase controls R-type anion channel QUAC1 in Arabidopsis guard cells. *Plant J.* **74**, 372–382. doi: 10.1111/tpj.12133

Kanazawa, A., Ostendorf, E., Kohzuma, K., Hoh, D., Strand, D. D., Sato-Cruz, M., et al. (2017). Chloroplast ATP synthase modulation of the thylakoid proton motive force: Implications for photosystem I and photosystem II photoprotection. *Front. Plant Sci.* **8**:719. doi: 10.3389/fpls.2017.00719

Kohzuma, K., Cruz, J. A., Akashi, K., Hoshiyasu, S., Munekage, Y. N., Yokota, A., et al. (2009). The long-term responses of the photosynthetic proton circuit to drought. *Plant Cell Environ.* **32**, 209–219. doi: 10.1111/j.1365-3040.2008. 01912.x

Kozlowski, T. T., and Pallardy, S. G. (2002). Acclimation and adaptive responses of woody plants to environmental stresses. *Bot Rev.* **68**, 270–334. doi: 10.1663/ 0006-8101(2002)068[0270:AAAROW]2.0.CO;2
Kramer, D. M., Avenson, T. J., and Edwards, G. E. (2004a). Dynamic flexibility in the light reactions of photosynthesis governed by both electron and proton transfer reactions. *Trends Plant Sci.* 9, 349–357.

Kramer, D. M., Johnson, G., Kiirats, O., and Edwards, G. E. (2004b). New fluorescence parameters for the determination of QA redox state and excitation energy fluxes. *Photosynth. Res.* 79, 209. doi: 10.1023/B:PRES.0000015391.99477.0d

Kramer, D. M., Sacksteder, C. A., and Cruz, J. A. (1999). How acidic is the lumen? *Photosynth. Res.* 60, 151–163. doi: 10.1023/A:1006212014787

Krishna-Reddy, S., Liu, S., Rudd, J. C., Xue, Q., Payton, P., Finlayson, S. A., et al. (2014). Physiology and transcriptomics of water-deficit stress responses in wheat cultivars TAM 111 and TAM 112. *J. Plant Physiol.* 171, 1289–1298. doi: 10.1016/j.jplph.2014.05.005

Kromdijk, J., Glowacka, K., Leonelli, L., Gabilly, S. T., Iwai, M., Niyogi, K. K., et al. (2016). Improving photosynthesis and crop productivity by accelerating recovery from photoprotection. *Science* 354:857. doi: 10.1126/science.aai8878

Kümpers, B. M., Burgess, S. J., Reyna-Llorens, I., Smith-Unna, R., Boursnell, C., and Hibberd, J. M. (2017). Shared characteristics underpinning C4 leaf maturation derived from analysis of multiple C3 and C4 species of Flaveria. *J. Exp. Bot.* 68, 177–189. doi: 10.1093/jxb/erw488

Kuznetsova, A., Brockhoff, P., and Christensen, R. (2014). *lmerTest: Tests for Random and Fixed Effects for Linear Mixed Effect Models (lmer Objects of lme4 Package)*. R Package Version 2.0-6. Available at: http://cran.r-project.org/package=lmerTest

Lal, A., and Edwards, G. (1996). Analysis of inhibition of photosynthesis under water stress in the C4 species *Amaranthus cruentus* and *Zea mays*: electron transport, CO2 fixation and carboxylation capacity. * Funct. Plant Biol.* 23, 403–412.

Langfelder, P., and Horvath, S. (2007). Eigengene networks for studying the relationships between co-expression modules. *BMC Syst. Biol.* 1:54. doi: 10.1186/1752-0509-1-54

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. *BMC Bioinformatics* 9:559. doi: 10.1186/1471-2105-9-559

Lawlor, D. W., and Cornic, G. (2002). Photosynthetic carbon assimilation and associated metabolism in relation to water deficits in higher plants. *Plant Cell Environ.* 25, 275–294. doi: 10.1046/j.0016-8025.2001.00814.x

Lawlor, D. W., and Tezara, W. (2009). Causes of decreased photosynthetic rate and metabolic capacity in water-deficient leaf cells: a critical evaluation of mechanisms and integration of processes. *Ann Bot.* 103, 561–579. doi: 10.1093/aob/mcn244

Li, X.-P., Bjorkman, O., Shih, C., Grossman, A. R., Rosenquist, M., Jansson, S., et al. (2000). A pigment-binding protein essential for regulation of photosynthetic light harvesting. *Nature.* 403, 391–395. doi: 10.1038/35000131

Lim, C., Baek, W., Jung, J., Kim, J.-H., and Lee, S. (2015). Function of ABA in stomatal defense against biotic and drought stresses. *Int. J. Mol. Sci.* 16, 15251. doi: 10.3390/ijms160715251

Liu, N., Ding, Y., Fromm, M., and Avramova, Z. (2014). Endogenous ABA extraction and measurement from Arabidopsis leaves. *Bio Protocol* 4:e1257. doi: 10.21769/BioProtoc.1257

Lobell, D. B., and Gourdji, S. M. (2012). The influence of climate change on global crop productivity. *Plant Physiol.* 160, 1686. doi: 10.1104/pp.112.208298

Loriaux, S. D., Avenson, T. J., Welles, J. M., McDermitt, D. K., Eckles, R. D., Riensche, B. et al. (2013). Closing in on maximum yield of chlorophyll fluorescence using a single multiphase flash of sub-saturating intensity. *Plant Cell Environ.* 36, 1755–1770. doi: 10.1111/pce.12115

Lu, H.-F., Dong, H.-T., Sun, C.-B., Qing, D.-J., Li, N., Wu, Z.-K., et al. (2011). The panorama of physiological responses and gene expression of whole plant of maize inbred line YQ7-96 at the three-leaf stage under water deficit and re-watering. *Theor. Appl. Genet.* 123:943. doi: 10.1007/s00122-011-1638-0

Maxwell, K., and Johnson, G. N. (2000). Chlorophyll fluorescence—a practical guide. *J Exp Bot.* 51, 659–668. doi: 10.1093/jexbot/51.345.659

Medrano, H., Escalona, J. M., Bota, J., GulfAs, J., and Flexas, J. (2002). Regulation of photosynthesis of C3 plants in response to progressive drought: stomatal conductance as a reference parameter. *Ann. Bot.* 89, 895–905. doi: 10.1093/aob/mcf079

Mittler, R. (2002). Oxidative stress, antioxidants and stress tolerance. *Trends Plant Sci.* 7, 405–410. doi: 10.1016/S1360-1385(02)02312-9

Miyashita, K., Tanakamaru, S., Maitani, T., and Kimura, K. (2005). Recovery responses of photosynthesis, transpiration, and stomatal conductance in kidney bean following drought stress. *Environ. Exp. Bot.* 53, 205–214. doi: 10.1016/j.envexpbot.2004.03.015

Mori, I. C., Murata, Y., Yang, Y., Munemasa, S., Wang, Y.-F., Andreoli, S., et al. (2006). CDPKs CPK6 and CPK3 Function in ABA regulation of guard cell S-type anion- and Ca2+-permeable channels and stomatal closure. *PLoS Biol.* 4:e327. doi: 10.1371/journal.pbio.0040327

Müller, P., Li, X.-P., and Niyogi, K. K. (2001). Non-photochemical quenching. A response to excess light energy. *Plant Physiol.* 125, 1558. doi: 10.1104/pp.125.4.1558

Nicotra, A. B., and Davidson, A. (2010). Adaptive phenotypic plasticity and plant water use. *Funct. Plant Biol.* 37, 117–127. doi: 10.1071/FP09139

Niyogi, K. K., Grossman, A. R., and Björkman, O. (1998). Arabidopsis mutants define a central role for the xanthophyll cycle in the regulation of photosynthetic energy conversion. *Plant Cell.* 10:1121. doi: 10.1105/tpc.10.7.1121

North, H. M., Frey, A., Boutin, J.-P., Sotta, B., and Marion-Poll, A. (2005). Analysis of xanthophyll cycle gene expression during the adaptation of Arabidopsis to excess light and drought stress: changes in RNA steady-state levels do not contribute to short-term responses. *Plant Sci.* 169, 115–124. doi: 10.1016/j.plantsci.2005.03.002

Parent, B., Hachez, C., Redondo, E., Simonneau, T., Chaumont, F., and Tardieu, F. (2009). Drought and abscisic acid effects on aquaporin content translate into changes in hydraulic conductivity and leaf growth rate: A trans-scale approach. *Plant Physiol.* 149:2000. doi: 10.1104/pp.108.130682

Philippar, K., Büchsenschütz, K., Abshagen, M., Fuchs, I., Geiger, D., Lacombe, B., et al. (2003). The K+ Channel KZM1 mediates potassium uptake into the phloem and guard cells of the C4 grass *Zea mays*. *J. Biol. Chem.* 278, 16973–16981. doi: 10.1074/jbc.M212720200

Porter, J. R., and Semenov, M. A. (2005). Crop responses to climatic variation. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* 360, 2021–2035. doi: 10.1098/rstb.2005.1752

R Core Development Team (2011). *A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Rao, X., and Dixon, R. A. (2016). The differences between NAD-ME and NADP-ME subtypes of C4 photosynthesis: more than decarboxylating enzymes. *Front. Plant Sci.* 7:1525. doi: 10.3389/fpls.2016.01525

Ruuska, S. A., Badger, M. R., Andrews, T. J., and von Caemmerer, S. (2000). Photosynthetic electron sinks in transgenic tobacco with reduced amounts of Rubisco: little evidence for significant Mehler reaction. *J. Exp. Bot.* 51, 357–368. doi: 10.1093/jexbot/51.suppl_1.357

Schreiber, U., Hormann, H., Neubauer, C., and Klughammer, C. (1995). Assessment of photosystem II photochemical quantum yield by chlorophyll fluorescence quenching analysis. *Funct. Plant Biol.* 22, 209–220. doi: 10.1016/j.bbbio.2010.01.016

Schulze, E. (1986). Carbon dioxide and water vapor exchange in response to drought in the atmosphere and in the soil. *Annu. Rev. Plant Physiol.* 37, 247–274. doi: 10.1146/annurev.pp.37.060186.001335

Slatyer, R. (1967). *Plant-Water Relationships*. Cambridge, MA: Academic Press.

Takizawa, K., Kanazawa, A., and Kramer, D. M. (2008). Depletion of stromal Pi induces high 'energy-dependent' antenna exciton quenching (qE) by decreasing proton conductivity at CFO-CF1 ATP synthase. *Plant Cell Environ.* 31, 235–243. doi: 10.1111/j.1365-3040.2007.01753.x

Tardieu, F., Parent, B., and Simonneau, T. (2010). Control of leaf growth by abscisic acid: hydraulic or non-hydraulic processes? *Plant Cell Environ.* 33, 636–647. doi: 10.1111/j.1365-3040.2009.02091.x

To, T., and Kim, J.-M. (2014). Epigenetic regulation of gene responsiveness in Arabidopsis. *Front. Plant Sci.* 5:458. doi: 10.3389/fpls.2013.00548

Trenberth, K. E., Dai, A., van der Schrier, G., Jones, P. D., Barichivich, J., Briffa, K. R., et al. (2014). Global warming and changes in drought. *Nat. Clim. Chang* 4, 17–22. doi: 10.1038/nclimate2067

Turner, N. (2003). "Drought hardening and presowing seed hardening," in *Encyclopedia of Water Science*, ed. S. W. Trimble (Boca Raton, FL: CRC Press), 166–172.
Tyystjärvi, E., and Aro, E. M. (1996). The rate constant of photoinhibition, measured in lincomycin-treated leaves, is directly proportional to light intensity. *Proc. Natl. Acad. Sci. U.S.A.* 93, 2213–2218. doi: 10.1073/pnas.93.5.2213

Vahisalu, T., Kollist, H., Wang, Y.-F., Nishimura, N., Chan, W.-Y., Valerio, G., et al. (2008). SLAC1 is required for plant guard cell S-type anion channel function in stomatal signalling. *Nature* 452, 487–491. doi: 10.1038/nature06608

Virlouvet, L., Ding, Y., Fujii, H., Avramova, Z., and Fromm, M. (2014). ABA signaling is necessary but not sufficient for RD29B transcriptional memory during successive dehydration stresses in *Arabidopsis thaliana*. *Plant J.* 79, 150–161. doi: 10.1111/tpj.12548

Virlouvet, L., and Fromm, M. (2015). Physiological and transcriptional memory in guard cells during repetitive dehydration stress. *New Phytol.* 205, 596–607. doi: 10.1111/nph.13080

Walter, J., Nagy, L., Hein, R., Rascher, U., Beierkuhnlein, C., Willner, E., et al. (2011). Do plants remember drought? Hints towards a drought-memory in grasses. *Environ. Exp. Bot.* 71, 34–40. doi: 10.1016/j.envexpbot.2010.10.020

Waters, M., Moylan, E., and Langdale, J. (2008). GLK transcription factors regulate chloroplast development in a cell-autonomous manner. *Plant J.* 56, 432–444. doi: 10.1111/j.1365-313X.2008.03616.x

Xiong, L., and Zhu, J.-K. (2001). Abiotic stress signal transduction in plants: molecular and genetic perspectives. *Physiol. Plant.* 112, 152–166. doi: 10.1034/j.1399-3054.2001.1120202.x

Yu, Q., An, L., and Li, W. (2014). The CBL–CIPK network mediates different signaling pathways in plants. *Plant Cell Rep.* 33, 203–214. doi: 10.1007/s00299-013-1507-1

Zhu, X.-G., Long, S. P., and Ort, D. R. (2010). Improving photosynthetic efficiency for greater yield. *Annu. Rev. Plant Biol.* 61, 235–261. doi: 10.1146/annurev-plant-042809-112206

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Virlouvet, Avenson, Du, Zhang, Liu, Fromm, Avramova and Russo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.