# Responses to Hypoxia and Endoplasmic Reticulum Stress Discriminate the Development of Vitreous and Floury Endosperms of Conventional Maize (Zea mays) Inbred Lines

Mathieu Gayral, Khalil Elmorjani, Michèle Dalgalarrondo, Sandrine Balzergue, Stephanie Pateyron, Marie Helene Morel, Sylvie Brunet, Laurent Linossier, Caroline Delluc, Benedicte Bakan, et al.

► To cite this version:

Mathieu Gayral, Khalil Elmorjani, Michèle Dalgalarrondo, Sandrine Balzergue, Stephanie Pateyron, et al. Responses to Hypoxia and Endoplasmic Reticulum Stress Discriminate the Development of Vitreous and Floury Endosperms of Conventional Maize (Zea mays) Inbred Lines. Frontiers in Plant Science, 2017, 8, pp.557. 10.3389/fpls.2017.00557. hal-01602776

HAL Id: hal-01602776

https://hal.science/hal-01602776

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

# Responses to Hypoxia and Endoplasmic Reticulum Stress Discriminate the Development of Vitreous and Floury Endosperms of Conventional Maize (Zea mays) Inbred Lines

Mathieu Gayral¹, Khalil Elmorjani¹, Michèle Dalgalarrondo¹, Sandrine M. Balzergue²,³†,
Stéphanie Pateyron²,³, Marie-Hélène Morel⁴, Sylvie Brunet⁵, Laurent Linossier⁵,
Caroline Delluc⁶, Bénédicte Bakan¹ and Didier Marion¹*

¹ Biopolymers, Interactions, Assemblies, Institut National de la Recherche Agronomique, Nantes, France, ² POPS (transcriptOmic Platform of iPS2) Platform, Centre National de la Recherche Scientifique, Institut of Plant Sciences Paris Saclay, Institut National de la Recherche Agronomique, Université Paris-Sud, Université Evry, Université Paris-Saclay, Orsay, France, ³ Institute of Plant Sciences Paris-Saclay, Paris Diderot, Sorbonne Paris-Cité, Orsay, France, ⁴ Agropolymer Engineering and Emerging Technologies, Institut National de la Recherche Agronomique, Montpellier, France, ⁵ Limagrain Céréales Ingrédients, Riom, France, ⁶ Limagrain Europe, Chappes, France

OPEN ACCESS

**Edited by:**
Soren K. Rasmussen,
University of Copenhagen, Denmark

**Reviewed by:**
Bryan C. Gibbon,
Florida A&M University, USA
Hardy Rolletschek,
Institute of Plant Genetics and Crop
Plant Research (LG), Germany

*Correspondence:
Didier Marion
didier.marion@inra.fr

†Present Address:
Sandrine Balzergue,
L'institut de Recherche en Horticulture
et Semences, Institut National de la
Recherche Agronomique,
AGROCAMPUS-Ouest, Université
d'Angers, Beaucouzé, France

**Specialty section:**
This article was submitted to
Crop Science and Horticulture,
a section of the journal
Frontiers in Plant Science

Received: 10 January 2017

Accepted: 28 March 2017

Published: 13 April 2017

**Citation:**
Gayral M, Elmorjani K, Dalgalarrondo M, Balzergue SM, Pateyron S, Morel M-H, Brunet S, Linossier L, Delluc C, Bakan B and Marion D (2017) Responses to Hypoxia and Endoplasmic Reticulum Stress Discriminate the Development of Vitreous and Floury Endosperms of Conventional Maize (Zea mays) Inbred Lines. Front. Plant Sci. 8:557.
doi: 10.3389/fpls.2017.00557

Major nutritional and agronomical issues relating to maize (Zea mays) grains depend on the vitreousness/hardness of its endosperm. To identify the corresponding molecular and cellular mechanisms, most studies have been conducted on opaque/floury mutants, and recently on Quality Protein Maize, a reversion of an opaque2 mutation by modifier genes. These mutant lines are far from conventional maize crops. Therefore, a dent and a flint inbred line were chosen for analysis of the transcriptome, amino acid, and sugar metabolites of developing central and peripheral endosperm that is, the forthcoming floury and vitreous regions of mature seeds, respectively. The results suggested that the formation of endosperm vitreousness is clearly associated with significant differences in the responses of the endosperm to hypoxia and endoplasmic reticulum stress. This occurs through a coordinated regulation of energy metabolism and storage protein (i.e., zein) biosynthesis during the grain-filling period. Indeed, genes involved in the glycolysis and tricarboxylic acid cycle are up-regulated in the periphery, while genes involved in alanine, sorbitol, and fermentative metabolisms are up-regulated in the endosperm center. This spatial metabolic regulation allows the production of ATP needed for the significant zein synthesis that occurs at the endosperm periphery; this finding agrees with the zein-decreasing gradient previously observed from the sub-aleurone layer to the endosperm center. The massive synthesis of proteins transiting through endoplasmic reticulum elicits the unfolded protein responses, as indicated by the splicing of bZip60 transcription factor. This splicing is relatively higher at the center of the endosperm than at its periphery. The biological responses associated with this developmental stress, which control the starch/protein balance, leading ultimately to the formation of the vitreous and floury regions of mature endosperm, are discussed.

Keywords: maize, endosperm, vitreousness, transcriptome, hypoxia, unfolded protein response (UPR)
# INTRODUCTION

Endosperm texture is an important quality trait of maize (Zea mays L.), determined by the proportions of vitreous and floury endosperm. Thus, the identification of the physicochemical mechanisms leading to vitreousness is fundamental in terms of both nutritional and agronomical issues. This was mainly achieved by studying different opaque/floury mutants, especially opaque2 (o2), and Quality Protein Maize (QPM), where vitreous endosperm is restored by introgressing "o2 modifier" genes (mo2) within an o2 background (Gibbon and Larkins, 2005). These studies highlighted the close relationships between zeins, starch, and maize vitreousness. Zeins, the storage proteins, are synthesized in the endoplasmic reticulum (ER) and stored in protein bodies. Numerous floury/opaque mutants are affected in zein biosynthesis, that leads to defect in the structure and shape of protein bodies (Wu and Messing, 2010). Indeed, o2 encodes a defective transcription factor (TF) that, in its wild-type form, regulates the expression of several genes during endosperm development. Such mutations result in a decrease in zein content and smaller protein bodies (Schmidt et al., 1987, 1990). *Floury2*, *floury4*, *Defective endosperm-B30*, and *Mucronate* mutants display mutations in α-zeins or γ-zein, leading to unshaped protein bodies and low zein content (Coleman et al., 1997; Kim et al., 2004, 2006; Wang et al., 2014). *Opaque1* and *floury1* mutants are impacted in the ER machinery, leading to disrupted protein body shape without a significant decrease in zein content (Holding et al., 2007; Wang G. et al., 2012). Starch biosynthesis is also involved in endosperm vitreousness. Indeed, *opaque5* (o5), and *shrunken4* (sh4) mutants that are affected in starch biosynthesis also display a floury endosperm (Yu et al., 2001; Myers et al., 2011). Finally, the recovery of a vitreous endosperm in QPM is related to an increase in amylose content and shorter amylopectin branching (Gibbon et al., 2003; Salazar-Salas et al., 2014) as well as an increase in 27 kDa γ-zein content and well-shaped protein bodies (Wu et al., 2010).

However, opaque mutations and *mo2* genes have pleiotropic effects in regard to both the number and differential expression level of up- and down-regulated genes when compared with the wild type (Hunter et al., 2002). Recently, characterization of flint and dent maize inbred lines, used in breeding programs, revealed gradients of proteins and starch closely related to the vitreous/floury transition (Gayral et al., 2016). In mature endosperm, the decrease in protein content from the periphery to the center of the endosperm is mainly attributed to α-zeins. Amylose content follows the same distribution, leading to a starch crystallinity gradient. This starch crystallinity is also affected by higher amounts of starch-bound lipids in the vitreous area (Gayral et al., 2015). These gradients are in agreement with the centrifugal development of endosperm from mitotic cell division to programmed cell death (PCD). In maize, PCD starts at the center of the endosperm (the developing floury region) to finally affect the entire endosperm, including the periphery (the developing vitreous region; Young and Gallie, 2000a,b). The origin of endosperm cell death has not yet been established; nevertheless, it is known that developing maize endosperm endures different stresses. Specifically, endosperm faces hypoxic

conditions from 10 days after pollination (DAP; Rolletschek et al., 2005). In fact, oxygen deficiency leads to different responses, which result in plant survival or cell death (Bailey-Serres et al., 2012). Characterization of the molecular processes associated with endosperm development is essential to assist in breeding of conventional maize hybrids for improved vitreousness adapted to food uses. In particular, in the context of climate changes, close relationships may be expected between certain genes expressed during endosperm development and those expressed during abiotic (heat and drought) stress. Therefore, a transcriptomic approach combined with the analyses of key metabolites was performed on maize inbred lines. The results showed that the development of vitreous and floury endosperm is clearly associated with significant differences in the responses to hypoxia and ER stress through the regulation of storage protein biosynthesis and energy metabolism.

# EXPERIMENTAL PROCEDURES

## RNA Extraction

For RNA extraction we used 10 grains from the same ear of a flint and a dent inbred line grown in the Limagrain research station (Mons, France) and harvested at 15 and 20 DAP. The inbred lines L3 (flint) and L5 (dent) were chosen from a previous study where 13 flint and dent inbred lines were analyzed for their protein, starch, and lipid contents in relation to endosperm vitreousness (Gayral et al., 2015, 2016). Two biological repetitions were carried out with grains from different ears. The grains were hand-dissected. In a first time, the middle third of the grain was isolated with a scalpel blade, then an endosperm piece of tissue was extracted from the central endosperm (which will develop in the floury region) and the peripheral endosperm (which will develop in the vitreous region) using a small cookie cutter (Ø = 0.9 mm), and immediately poured into liquid nitrogen before grinding. Total RNA extraction for microarray and RT-PCR (reverse-transcription PCR) analysis was performed using the RNeasy Plant Mini Kit (Qiagen, Les Ulis, France).

## Statistical Analysis of Microarray Data

Experiments were designed with the Genomics Networks team of the Institute of Plant Sciences (IPS2). For each array, the raw data comprised the logarithm of median intensities of feature pixels at wavelengths of 635 nm (red) and 532 nm (green). For each array, a global intensity-dependent normalization using the loess procedure (Yang et al., 2002) was performed to correct the dye bias. The differential analysis was based on averaging the log-ratios over the duplicate probes and over the technical replicates. Hence, the number of available data for each gene equals the number of biological replicates and is used to calculate the moderated *t*-test (Smyth, 2004). Under the null hypothesis, no evidence that the specific variances vary between probes was highlighted by Limma analysis, and consequently the moderated *t*-statistic was assumed to follow a standard normal distribution. To control the false discovery rate (FDR), adjusted *p*-values found using the optimized FDR approach described by Storey and Tibshirani (2003) were calculated. We considered the probes with an adjusted *p*-value ≤ 0.05 to be differentially expressed. Analysis
was performed using R software. The function SqueezeVar of the library Limma was used to smooth the specific variances by computing empirical Bayesian posterior means. The library kerfdr (kernel-based local fdr) was used to calculate the adjusted *p*-values.

## Data Deposition

Microarray data from this article were deposited in the international repository GEO, Gene Expression Omnibus (Edgar et al., 2002; http://www.ncbi.nlm.nih.gov/geo/), accession no. GSE88993. All steps of the experiment, from growth conditions to bioinformatic and statistical analyses, were detailed in CATdb (Complete *Arabidopsis* Transcriptome database; Gagnot et al., 2008; http://tools.ips2.u-psud.fr/CATdb/; Project: 4plex_Maize_2015_02) according to the “Minimum Information About a Microarray Experiment” standards.

## Annotation and Gene Ontology

Genes were annotated by sequence homology with *A. thaliana* or *O. sativa* according to MaizeGDB (Maize Genetics and Genomics Database; http://www.maizegdb.org; Sen et al., 2009).

Identification of GO (gene ontology) categories significantly enriched ($p \le 0.05$) was based on “plant GO slim” and done using online AgriGO software (http://bioinfo.cau.edu.cn/agriGO; Du et al., 2010).

## qRT-PCR

Total RNA (4 µg) was reverse-transcribed using SuperScript III and PolyT primer (Invitrogen, ThermoFisher Scientific, Courtaboeuf, France). qPCR (qualitative PCR) quantification was performed with CFX Connect (Bio-Rad, Marnes-la-Coquette, France) and Power SYBR Green PCR Master Mix (Applied Biosystems, ThermoFisher Scientific, Courtaboeuf, France). The relative gene expression level was normalized to ubiquitin (GRMZM2G109977), and was defined as $\Delta Cq = 2^{-(CqGene - CqUbi)}$. The data represent the mean of two technical replicates for two biological repetitions ($\pm SD$). All primers used in the study are listed in Table S1.

## Metabolite Extraction and Quantification

For metabolite extraction, we used five grains (30 DAP) from the same ear. Three biological repetitions were carried out from different ears. Grains were hand-dissected, and samples were freeze-dried and ground in liquid nitrogen. Extraction was conducted on 20 mg of sample with 500 µL of CH3OH containing 0.2 mM norleucine and 0.4 mM adonitol as an internal standard. After 15 min of stirring, 250 µL of CHCl3 and 500 µL of water were added. After stirring and centrifugation, the aqueous phase was collected and used for metabolite quantification. Aliquots (100 µL) were dried in a SpeedVac concentrator (Savant, ThermoFisher Scientific, Courtaboeuf, France) and derivatized for 2 h at 37°C in 50 µL 20 mg.mL⁻¹ methoxyamine in pyridine, then for 30 min at 37°C after the addition of 50 µL of BSTFA (N, O-bis (trimethylsilyl) trifluoroacetamide). Metabolites were analyzed by GC-FID/MS (gas chromatography with flame ionization detection and mass spectrometry) using Clarus SQ8C (Perkin Elmer, Villebon sur

Yvette, France) equipped with a DB-5MS capillary column (30 m × 0.25 mm × 0.25 µm; Agilent Technologies, Les Ulis, France), and identified according to mass spectra and retention time compared with standard solution. Metabolites were normalized and quantified according to adonitol.

For quantification of free amino acids, aliquots (200 µL) were dried in a SpeedVac concentrator, and derivatized in 20 µL H2O/CH3CH2OH/N(CH2CH3)3/PITC (phenyl isothiocyanate) (1/7/1/1) for 10 min at room temperature. Samples were then dried and solubilized in 100 µL 2 mM-Na2HPO4 pH 7.4, with 5% of CH3CN. Free amino acids were analyzed by reversed phase HPLC (reversed-phase high-performance liquid chromatography) using a Pico-Tag column (3.9 × 300 mm; Waters, Guyancourt, France) on an Alliance HT apparatus equipped with a UV detector (254 nm; Waters, Guyancourt, France). Elution (1 mL.min⁻¹) was performed at 40°C using a linear gradient from 100% eluent A (70 mM-CH3ČCOONa pH 6.45 with 2.5% CH3CN) to 100% eluent B (CH3CN/CH3OH/H2O (45/15/40)). The quantity of amino acid was normalized to the norleucine internal standard and quantified in regard to the calibration curve obtained by injection of standard amino acids (20,088, Pierce, ThermoFisher Scientific, Courtaboeuf, France) derivatized under the same conditions.

## RESULTS

### Transcriptome Specifications of Central and Peripheral Endosperm during Grain Filling

Maize grain development starts with the double fertilization of the female gamete and the central cells, which leads to the development of the embryo and the endosperm respectively. Upon 10 DAP, protein and starch begin to accumulate in the endosperm. To ascertain whether transcriptome modifications are related to vitreousness, we analyzed the RNAs extracted from central and peripheral endosperm that will become at maturity the floury and the vitreous areas, respectively. From here on, we will refer to these regions as maize endosperm center (MEC) and maize endosperm periphery (MEP). MEC and MEP transcriptomes were analyzed on a flint and a dent inbred line harvested at 15 and 20 DAP. By considering only genes up-regulated in MEC, i.e., 5,518, or in MEP, i.e., 5,167, at least in a developing stage and in a genetic background (flint or dent), gene ontology (GO) analysis revealed enrichment in “metabolic process,” “cellular process,” “cellular metabolic process,” and “biosynthetic process” in MEP, whereas in MEC we notably observed enrichments in “response to stress,” “response to stimulus,” and “response to abiotic stimulus” (Figure 1A). These results suggest that MEC and MEP follow their own development program, especially with respect to the stronger stress responses in MEC.

From 15 to 20 DAP developing stages, at least 9369 genes were found to be differentially expressed when MEC and/or MEP of flint and/or dent inbred lines were compared. Considering genes significantly regulated in MEC and MEP, we did not observe strong GO enrichment pattern differences (Figure S1A).
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | MEP and MEC transcriptome pattern. (A) Biological process GO categories significantly enriched in genes only up-regulated in MEC and MEP in almost a tested condition. (B) Annotated genes related to stress differentially expressed between MEC and MEP at 15 DAP and 20 DAP in flint and dent maize.</figcaption></figure>

Enriched GO categories are mainly related to metabolic and biosynthetic processing (Figure S1A), in agreement with the massive synthesis of storage compounds. However, we only observed GO enrichment of “response to stress” in MEC, whereas “generation of precursor metabolites and energy” and “multicellular organismal process” were enriched in MEP. Therefore, MEC and MEP follow a close developmental process, and differences in stress responses observed between MEC and MEP persist throughout grain filling.

To further characterize the relationships between developmental stress responses and maize vitreousness, we focused on genes that showed significant differences in expression levels between MEC and MEP in flint and dent endosperm at 15 and 20 DAP: 634 genes were up-regulated in MEP at 15 and 20 DAP in flint and dent inbred lines, and 525 genes were annotated, while in MEC, 776 genes were up-regulated and 602 were annotated. Differentially expressed genes involved in stress management are shown in **Figure 1B**. In MEP, the expression of genes involved in oxido-reduction (including responses to oxidative burst), ER-associated degradation (ERAD) and abiotic/biotic stress responses was clearly activated. In MEC, we also observed the expression of genes involved in vesicular trafficking, PCD, autophagy, hypoxia, and regulation of transcription (**Figure 1B**); these data suggested a higher stress response in this area at 15 and 20 DAP.

To further investigate the temporal modifications linked to vitreousness, we focused on genes showing significant differences in MEP and/or MEC in both flint and dent lines between 15 and 20 DAP. 378 genes were significantly up-regulated from 15 to 20 DAP in MEP and/or MEC in the two inbred lines, and 258 of these genes were annotated. On the other hand, 348 genes were found to be significantly down-regulated (from 15 to 20

DAP) in MEP and/or MEC in the two inbred lines, and of these genes, 280 were annotated. If we consider the genes involved in stress responses, their expression was mainly up-regulated at 20 DAP (Figure S1B). In fact, almost the same genes were found to be up-regulated in MEP and MEC, with the exception of those implicated in ERAD, essentially up-regulated in MEC. Between 15 and 20 DAP, PCD genes were found to be up-regulated in both MEC and MEP (Figure S1B), confirming that PCD increases throughout endosperm development as previously shown (Young and Gallie, 2000a,b).

## Starch and Zein-Associated Transcriptome Signatures

In regard to the genes encoding enzymes involved in core starch biosynthesis, ADP-glucose pyrophosphorylase, starch synthases, branching, and debranching enzymes, we did not observe any coherent differential expression. Indeed, certain genes were found to be up-regulated in MEP, while their homologs were found to be up-regulated in MEC (Figure S2). These findings are probably due to the fact that starch deposition follows gradients from the upper to the basal part of the endosperm (Dochlert, 1990) as well as from the periphery to the center (Gayral et al., 2016). This fact could also be related to the regulation of starch biosynthesis at the post-translational level through phosphorylations and formation of multienzyme complexes (Hennen-Bierwagen et al., 2009). In fact, starch biosynthesis starts with ADP-glucose a product provided by the sucrose synthase (SUS) pathway. Indeed, ectopic overexpression of a potato SUS in maize induces an increase in ADP-glucose and starch contents (Li et al., 2013). Three *sus* genes are strongly expressed in maize endosperm: *sus1* (GRMZM2G152908), *sus2* (GRMZM2G318780), and *sus-sh1* (GRMZM2G089713;
Figure S2). Previous studies suggest that SUS1 and SUS-SH1 are involved in starch and cell wall polysaccharide synthesis, respectively (Chourey et al., 1998). In contrast with invertases, SUS is an energy-conserving pathway for sucrose breakdown since it uses one PPi molecule and produces one UTP molecule. UTP can be used to phosphorylate fructose- and subsequent entry of fructose-6-phosphate in the glycolysis pathway (Geigenberger, 2014). *Sus1* (GRMZM2G152908) and *sus2* (GRMZM2G318780) are up-regulated in MEP (Figure S2), and display the highest differential expression among all genes that discriminate developing MEP from developing MEC. This is consistent (i) with the localization of the intense corresponding activity at the periphery of developing maize endosperm (Wittich and Vreugdenhil, 1998) (ii) the formation of SUS1-SUS2 hetero-oligomers (Duncan et al., 2006) and (iii) the up-regulation of *sus1* in hypoxic conditions (Zeng et al., 1998). Any differential expression is observed for *sus-sh1* (GRMZM2G089713) that encodes SUS-SH1 (Figure S2) in agreement with the fact that the expression of this gene is not up-regulated by hypoxia but by anoxia, a huge depletion of oxygen (Zeng et al., 1998). However, up-regulation of *sus1* and *sus2* in MEP is contradictory to (i) slightly lower starch contents in the vitreous than in the floury endosperm (Gayral et al., 2015) and (ii), as shown below, to up-regulation of genes markers of hypoxia in MEC. Altogether these results suggest that the management of hypoxia differs between MEP and MEC.

In contrast with starch deposition, protein deposition is significantly higher in the periphery of the endosperm (Gayral et al., 2016) than in its longitudinal axis, i.e., from the upper to the lower part of the endosperm (Dochlert, 1990). Furthermore, expression of zein genes is essentially regulated at the transcriptional level (Zhang et al., 2015; Qiao et al., 2016). Indeed, transcriptomic analysis revealed stronger zein expression in MEP at 15 and 20 DAP, especially for α and δ-zeins (Figure 2). These results agree with those of mRNA *in-situ* hybridization experiments, which revealed stronger α-zein expression at the periphery of developing endosperm (Woo et al., 2001). Furthermore, we observed a global up-regulation of genes specifying zeins at 20 DAP (Figure 2), in accordance with epigenetic regulation of α-zein expression in endosperm at 15–23 DAP (Locatelli et al., 2009). These results are consistent with the protein decrease observed previously from the sub-aleurone layer to the endosperm center (Gayral et al., 2016).

## ER-Stress, UPR and Endosperm Development

Down-regulation of genes specifying zeins in MEC is associated with up-regulation of genes involved in stress responses (Figure 1). When checked, chaperones, protein degradation (ERAD), vesicular trafficking, autophagy, and the transcription process are up-regulated in MEC (Figure 1B). All these mechanisms are generally induced in response to ER stress, leading to the corresponding signaling cascade named the unfolded protein response (UPR) and suggest stronger UPR in MEC than in MEP. In plants, two UPR pathways have been described. The first UPR pathway involves bZIP17/bZIP28, i.e.,

two membrane-associated TFs that are exported from ER to Golgi apparatus en route to the nucleus; the second pathway involves IRE1, i.e. a dual-protein kinase, and the bZIP60 TF. When unfolded or misfolded proteins accumulate in the ER, IRE1 is activated and the *bzip60* mRNA is spliced (Howell, 2013). The spliced form of *bzip60* leads to mRNA encoding a TF with a canonical nuclear targeting signal (Li et al., 2012; Wang B. et al., 2012). In addition to promoting UPR genes, bZIP60 up-regulates its own expression. In fact, microarray analysis showed a higher expression of *bzip60* in MEC than in MEP for both flint and dent maize and at both 15 and 20 DAP (Table S2). To characterize the UPR level in MEC and MEP, we followed *bzip60* splicing during endosperm development by quantifying the unspliced and spliced forms of this TF. RT-PCR analyses revealed spliced mRNA in both MEC and MEP at 15 and 20 DAP for both flint and dent maize (Figure 3A), attesting to the fact that UPR took place in developing endosperm. In flint, and especially in dent, inbred lines, qRT-PCR analyses globally showed a higher expression level at 20 DAP (Figure 3). These results reflect the increase in stress during endosperm development, as shown in Figure S1. In terms of total *bzip60* and unspliced *bzip60* expression levels, a higher level was observed in MEC, confirmed by microarray analysis (Table S2). When compared, splicing of *bzip60* is more pronounced in MEC than in MEP, confirming a stronger UPR in this region. This difference

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Spatiotemporal evolution of UPR during endosperm development. (A) *bzip60* (GRMZM2G025812) splicing was detected by RT-PCR in MEP and MEC at 15 and 20 DAP. Ubi: RT-PCR of *ubiquitin* (GRMZM2G109977). Arrow: spliced mRNA. Relative expression of *bzip60* (B), unspliced *bzip60* (C) and spliced *bzip60* (D) measured by qRT-PCR in MEP and MEC at 15 and 20 DAP. Relative expression level was normalized to *ubiquitin* (GRMZM2G109977). Asterisks indicate significant difference (paired t-test; \*p < 0.05; \*\*p < 0.01).</figcaption></figure>

was found to be more acute at 20 DAP. It is interesting to note that ER stress was previously associated with higher content of triacylglycerols (TAG) in maize endosperm of the *floury2* mutant

(Shank et al., 2001), in accordance with higher amounts of these storage lipids in the floury regions of conventional maize (Gayral et al., 2015). The higher TAG deposition in MEC relative to that in MEP is also related to major changes in the proteins that cover and stabilize oil droplets. Indeed, while an oleosin-like protein is highly up-regulated in MEP (GRMZM2G480954), an ortholog of *Arabidopsis* lipid droplet-associated proteins (LDAP) (GRMZM2G150367) is one of the most up-regulated proteins in MEC. These LDAPs are associated with stress responses and are mostly expressed in vegetative tissues (Laibach et al., 2015). Therefore, ER stress and UPR could be responsible for TAG accumulation in the floury regions of maize endosperm. Taken together, these results reveal a permanent ER stress in the endosperm, with a stronger UPR in MEC than in MEP.

## From Endosperm Stress Responses to Autophagy and PCD

UPR evolution is associated with spatiotemporal up-regulation of autophagy and PCD (Figure 1B and Figure S1), suggesting a possible relationship between these processes, especially in MEC. Such a relationship controls the pro-life and pro-death evolution of cells and is commonly encountered in all eukaryotes (Suzuki et al., 2017). Indeed, it is worth noting that the gene specifying BAX inhibitor-1, an ER-resident transmembrane protein that attenuates PCD owing to ER-stress in *Arabidopsis* (Watanabe and Lam, 2008), is up-regulated in MEC (Table S2). In the same way, the gene specifying Bcl-2-associated athanogene 1 (BAG1) is up-regulated in MEC (Table S2). In *Arabidopsis*, BAG1 acts as a co-chaperone in various processes related to environmental stress and PCD (Nawkar et al., 2016). The DAD1 (defender against cell death) gene, also localized in the ER, is another regulator of cell death. *dad1* is well expressed in both MEC and MEP, with a slight tendency to be more up-regulated in MEC (Table S2). DAD1 is a glycosyltransferase that is known to mediate, in *Arabidopsis*, N-glycosylation of nascent proteins and to suppress cell death induced by abiotic stresses (Danon et al., 2004). From 15 DAP, expression of the genes specifying these proteins, in maize endosperm, as well as their up-regulation in MEC, confirm the early PCD observed in this area. Nine metacaspases have been identified in *Arabidopsis*, AtMC1 to AtMC9, and their homologs, i.e., ZmMC1 to ZmMC9, were found in maize (Minina et al., 2013). Though these proteins differ from their animal counterparts, plant metacaspases are markers of PCD (Lam and del Pozo, 2000; Tsiatsiani et al., 2011; Fagundes et al., 2015). Microarray analyses revealed that two type II metacaspases are differentially expressed in maize endosperms: *mc7* is up-regulated in MEC, while *mc9* is up-regulated in MEP (Table S2), suggesting different PCD pathways in these two areas. Indeed, in Norway spruce, a close relationship was established between a type II metacaspase, autophagy and the subsequent PCD pathway (Minina et al., 2013).

Autophagy is a ubiquitous complex process allowing the degradation of cell organelles and the recycling of these components. Many genes encoding autophagy-related proteins (ATG proteins) have been identified (Li et al., 2015), most of them being expressed during endosperm development, suggesting that
the autophagic recycling machinery is fully active (Table S2). Remarkably, although differential expression of these genes between MEC and MEP is not systemically evidenced, some genes are up-regulated in MEP as *atg4b*, while *atg6a* and *atg8c* are up-regulated in MEC. In fact, most *atg8* genes, i.e., *atg8b*, *atg8c*, *atg8d*, and *atg8e* tend to be up-regulated in MEC (Table S2). ATG8 proteins are essential in the autophagic process since their lipidation by phosphatidylethanolamine is the actual *primum movens* of the autophagic process. ATG12, which participates in ATG8 lipidation, is also essential for autophagy (Suzuki et al., 2017), and consequently *atg12a* is up-regulated in MEC (Table S2). Expression of all *atg* genes confirms that an autophagic process is implicated in maize endosperm development (Chung et al., 2009), and suggests that a stronger autophagy takes place in MEC. These results agree with spatiotemporal relationships between UPR, autophagy and PCD in this endosperm area and show that UPR corresponds with the spatiotemporal development of PCD in maize endosperm.

## Hypoxia Response in Maize Endosperm
### Expression of Hypoxia-Responsive Genes

Transcriptomic analyses revealed an up-regulation of genes implicated in hypoxic responses in MEC (Figure 1B). Numerous TF homologs to *Arabidopsis* hypoxia-induced TF are up-regulated in MEC. Notably, VII ethylene response factors (VII ERFs), that plays a major role in the response to low oxygen in *Arabidopsis* (Bailey-Serres et al., 2012). VII ERF homologs are found in maize, especially the homolog of RAP2.12/RAP2.2, which is up-regulated in MEC (Table S2). RAP2.12/RAP2.2 induce the expression of different hypoxia-responsive genes as *alcohol dehydrogenase* (*adh*) and *pyruvate decarboxylase* (*pdc*; Bailey-Serres et al., 2012). These hypoxia-induced enzymes were indeed found to be up-regulated in MEC (Table S2).

## Soluble Sugars and Free Amino Acids

Analysis of soluble sugars revealed high levels of sucrose in flint and dent endosperm (Figure 4A) consistent with their import from maternal tissue in order to produce starch. In fact, these sucrose levels decrease from the basal to the upper endosperm (Doehlert and Kuo, 1990), following exactly the opposite gradient observed for starch accumulation (Doehlert, 1990). Therefore, the higher amount of soluble sugars in MEC is obviously in total discordance with the higher amounts of starch observed in the floury endosperm (Gayral et al., 2015). Conversely, the lower amounts of almost all free amino acids found in MEP, when compared with those in MEC (Figure 4B), are quite consistent with the significant level of protein synthesis in MEP, leading to a higher protein accumulation in the vitreous region of mature endosperm (Gayral et al., 2016). In flint and dent maize, our results showed high accumulation of alanine, glutamine, glutamate, and serine in whole endosperm, in full accordance with the hypoxic conditions already described (António et al., 2016). In fact, the production of alanine contributes to the preservation of carbon derived from glycolysis in competition with ethanol fermentation under hypoxic conditions. Consequently, the higher alanine amount

found in MEC (Figure 4B) supports the hypothesis of a stronger hypoxia (António et al., 2016) affecting this endosperm region.

## Differential Regulation of Respiratory and Fermentative Metabolisms in MEC and MEP

Hypoxia is directly linked to changes in energy metabolism, particularly to a switch from an oxidative to a fermentative metabolism (Zabalza et al., 2009). To investigate whether this switch is involved in the formation of vitreous endosperm, we analyzed the expression of genes implicated in the glycolytic and fermentative pathways in MEP and MEC (Figure 5 and Table S3). In almost all tested conditions (in flint and dent endosperm at 15 and 20 DAP), the results revealed an up-regulation of glycolytic enzymes from fructose-6-P to pyruvate in MEC (Figure 5). Moreover, *pdc* and *adh* are clearly up-regulated in this area at 15 and 20 DAP. Ethanol fermentation thus ensures the NAD (nicotinamide adenine dinucleotide) supply required for glycolysis, allowing energy production under hypoxic conditions. The sorbitol dehydrogenase (SDH) pathway can also generate NAD during conversion of fructose to sorbitol. The gene specifying SDH is slightly up-regulated in MEC, which correlates with accumulation of sorbitol in this region of the endosperm (Figure 4A and Figure S2) and its up-regulation by oxygen depletion (de Sousa et al., 2008). Furthermore, in MEC, we also observed an up-regulation of *alanine aminotransferase* (*ala-t*), which fully concurs with the higher alanine content of the floury endosperm (Figure 4).

It has been shown that pyrophosphate-fructose-6-phosphate1-phosphotransferase (PFP) and pyruvate orthophosphate dikinase (PPDK) activities are related to vitreousness recovery in QPM (Guo et al., 2012), the latter may potentially play a role in starch/protein balance (Mechin et al., 2007; Prioul et al., 2008a). During glycolysis, fructose-6-P is phosphorylated by phosphofructokinase (PFK) or PFP. Interestingly, *pfp* is up-regulated in MEC, whereas *pfk* is slightly up-regulated in MEP (Figure 5). Actually, PFK uses ATP whereas PFP uses PPI, allowing thus to save ATP under hypoxic conditions. These results fully concur with the stronger hypoxic conditions in MEC. Downstream of glycolysis, the expression of genes specifying pyruvate kinase (PK) and PPDK both producers of pyruvate from phosphoenolpyruvate (PEP), did not have the same transcriptomic pattern (Figure 5). Whereas *pk* was rather up-regulated in MEC, *ppdk* was slightly up-regulated in MEP, which correlates with QPM vitreousness adaptation. In the next step, pyruvate is catalyzed by pyruvate dehydrogenase (PDH). Conversely to the previous step of glycolysis, *pdh* was found to be weakly up-regulated in the MEP of dent maize and remarkably stable in flint maize. This result confirmed the close relationship already observed between acetyl-CoA and oxygen availability (Rolletschek et al., 2005).

All these results show that both respiratory and fermentative pathways co-exist in MEC and MEP, but in MEC, glycolysis is rather oriented toward fermentation and alanine production, whereas in MEP, glycolysis is rather directed to the TCA cycle to ensure ATP production. In the same way, the gene specifying phosphoenolpyruvate carboxylase (PEPC), which catalyzes the decarboxylation of pyruvate to oxaloacetate
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Metabolite profiling of MEP and MEC. Soluble sugars (A) and free amino acids (B) of MEP and MEC in dent and flint maize at 30 DAP. Asterisks indicate significant differences (paired t-test; *p* < 0.05).</figcaption></figure>

(OAA), is up-regulated in MEP, whereas the gene specifying phosphoenolpyruvate carboxykinase (PEPCK), which can catalyze the reverse reaction, is up-regulated in MEC (Figure 5). In view of the fact that citrate synthesis from PEP via PEPC is active in maize kernel (Jeanneau et al., 2002), the expression pattern of *pepc* and *pepck* also suggests a higher fermentative metabolism in MEC. In MEP, a part of the pyruvate is probably phosphorylated, leading to PEP in order to produce PPI, an energy source substitute of ATP under low oxygen conditions. The produced PEP is then probably transferred to TCA cycle through PEPC. In MEC however, OAA could be carboxylated to generate PEP, which could be used for the fermentative

metabolism requirements. Interestingly, both *pdc* and *pepck* are strongly up-regulated between 15 and 20 DAP (Figure 5), suggesting a stronger fermentative metabolism at 20 DAP, which fully concurs with the stronger hypoxic conditions. All these results suggest that a fermentative metabolism starts earlier in MEC than in MEP.

## Differential Regulation of TCA Cycle and GABA Shunt in MEC and MEP

The first step in the TCA cycle is the production of citrate catalyzed by citrate synthase (CTS). The gene specifying CTS is up-regulated in MEC in dent maize. Subsequently,
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Metabolic switch between MEP and MEC. Transcriptome pattern of genes involved in glycolysis, fermentation and TCA cycle in MEP (V) and MEC (F) at 15 and 20 DAP in flint and dent inbred lines. Right box lines represent MEP/MEC $log_2$ expression ratio. Left box lines represent 20/15 DAP $log_2$ expression ratio. For multigenic families, the results represent the mean of all genes. Genes ID and $log_2$ ratio of all genes are listed in Table S3. Blue lines, amino acid synthesis pathway; PFK, phosphofructokinase; PFP, pyrophosphate-fructose-6-phosphate-1-phosphotransferase; ALD, fructose-bisphosphate aldolase; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; PGK, phosphoglycerate kinase; PGM, phosphoglycerate mutase; ENO, enolase; PPDK, pyruvate orthophosphate</figcaption></figure>

(Continued)
<caption>**FIGURE 5 | Continued**</caption>

dikinase; PK, pyruvate kinase; PDH, pyruvate; PEPC, phosphoenolpyruvate carboxylase; PEPCK, phosphoenolpyruvate carboxykinase; PDC, pyruvate decarboxylase; ADH, alcohol dehydrogenase; ALA-T, alanine aminotransferase; CTS, citrate synthase; ACO, aconitase; IDH, isocitrate dehydrogenase; OGDC, oxoglutarate dehydrogenase complex; SCS, succinyl-CoA synthetase; SDH, succinate dehydrogenase; FUM, fumarate hydratase; MDH, malate dehydrogenase; ME, malic enzyme; GOGAT, glutamate synthase; GAD, glutamate decarboxylase; GABA-T, γ-aminobutyric acid trans-aminase; SSADH, succinate-semialdehyde dehydrogenase.

aconitase (ACO) and isocitrate dehydrogenase (IDH) produce α-ketoglutarate using citrate. In flint and dent inbred lines, the results revealed the up-regulation of *aco* and *idh* expression in MEP (Figure 5). During the TCA cycle, α-ketoglutarate is decarboxylated by the 2-oxoglutarate dehydrogenase complex (OGDC), leading to succinyl CoA and NADH, the succinyl CoA being subsequently hydrolyzed to succinate by succinyl-CoA synthase (SCS). The *ogdc* gene does not display significant differential expression levels in MEP and MEC, while the *scs* gene is up-regulated in the MEC of dent endosperm (Figure 5). Succinate dehydrogenase complex (SDH, also known as complex II of the mitochondrial electron transport chain) catalyzes the oxidation of succinate to fumarate, which is then hydrated, via fumarase (FUM), leading to the production of malate. *sdh* and *fum* genes are both clearly up-regulated in MEP (Figure 5). All these results support the hypothesis that maintaining a functional TCA cycle is probably a key step in the formation of vitreous endosperm. Moreover, the up-regulation of *sdh* increases both the net output of NADH and FADH as substrates for the electron transport chain and the direct input of electrons into this chain via SDH itself. These results support the idea of a metabolic switch when MEC and MEP are compared, where the latter maintains a significant respiratory metabolism.

Under hypoxic conditions, α-ketoglutarate is also used by glutamate synthase (GOGAT) and glutamate decarboxylase (GAD), which respectively produce glutamate and GABA via the GABA shunt pathway, leading to accumulation of glutamate and GABA (António et al., 2016). The genes specifying GOGAT and GAD are up-regulated in MEC (Figure 5), which concurs with the higher amount of glutamate in this tissue area (Figure 4). Hypoxia leads to a pH decrease (Felle, 2005) since protons are consumed during GABA production. The up-regulation of *gad* in MEC also supports stronger hypoxia in this area. The GABA shunt leads ultimately to succinate production via γ-aminobutyric acid trans-aminase (GABA-T) and succinic semialdehyde dehydrogenase (SSADH). The genes specifying these enzymes are clearly up-regulated in MEP (Figure 5). These results suggest that in MEP, the GABA shunt leads to the production of succinate and NADH production necessary to maintain the TCA cycle and respiration, respectively. These differences in the regulation of GABA shunt enzymes testify in favor of differences in hypoxia tolerance between MEP and MEC.

DISCUSSION

After their differentiation, maize endosperm cells start a
cellular program allowing the storage of nitrogen (zeins) and
carbon (starch). Deposition of these storage compounds follows
gradients from the periphery to the center of the endosperm
(Gayral et al., 2016) closely related to cell size and ploidy

gradient generated by mitosis and endoreplication (Dante et al., 2014). Remarkably, in conventional maize the transition from a floury to a vitreous endosperm also follows a centrifugal gradient (Gayral et al., 2016). In this work, the metabolic specificities of developing central and peripheral endosperms were characterized and useful data were collected in order to delineate the molecular mechanisms that lead to the formation of floury and vitreous endosperm. This was achieved by using differential transcriptomics of MEP and MEC, completed by analysis of some key metabolites. Finally, by performing this work on a flint and a dent maize inbred line, it was possible to target only specific metabolic features of cells that will form, later in the course of development, the vitreous and floury regions of the endosperm independently from the genetic background.

Transcriptomic analyses revealed significant differences between MEP and MEC that persist throughout endosperm development. Therefore, MEP and MEC display different development kinetics in agreement with the corresponding gradient of cell division and growth during endosperm development. Most of these differences were observed in both flint and dent maize inbred lines despite the fact that their genetic backgrounds are significantly different (Unterseer et al., 2016).

These differences are clearly associated with how the
endosperm endures various stresses. Hypoxia is probably the
primum movens of endosperm stress responses. Indeed, hypoxia,
previously characterized during maize endosperm development
(Rolletschek et al., 2005), was confirmed during the present
work, through the expression of different hypoxia-induced
genes. In particular, hypoxia induces an oxidative stress, i.e.,
the production of reactive oxygen species (ROS). ROS act
as signaling molecules and lead to strong modifications of
energy metabolism to ensure energy supply (Fukao and Bailey-
Serres, 2004). This situation generally leads to a metabolic
switch from respiration to fermentation in the endosperm of
maize (Prioul et al., 2008b), rice (Xu et al., 2008), and wheat
(Zhen et al., 2016) during grain filling process. In the present
work, a spatial metabolic adaptation between MEP and MEC
is suggested when oxygen deprivation is marked throughout
the endosperm (Rolletschek et al., 2005). Our results suggest
that this adaptation is different between the periphery and the
center of maize endosperm. Indeed, in MEP sus and GABA
shunt up-regulation allows sufficient glycolysis and TCA cycle
activities to sustain the ATP and NADH synthesis while, in
MEC, fermentative and sorbitol dehydrogenase pathways are
clearly up-regulated, ensuring thus a sufficient supply of NAD for
glycolysis to allow energy production under hypoxic conditions.
Consequently, differential transcriptomic results suggest that
both respiration and fermentation occur in both MEC and
MEP, but with different relative intensities. Indeed, we can make
the assumption that respiration could prevail in MEP allowing
sufficient production of ATP for the massive synthesis of storage proteins while fermentation could prevail in MEC limiting thus ATP production and protein synthesis. In that scenario, the lower ATP concentration observed in the periphery than in the central endosperm (Rolletschek et al., 2005) should be due to a rapid ATP consumption.

In regard to this metabolic adaptation, a special focus can be addressed here concerning the regulation of the expression of the three SUS isoforms, i.e., SUS1, SUS2, and SUS-SH1. These three genes are highly expressed but *sus1* and *sus2* expressions are up-regulated in MEP while *sus-sh1* expression remains stable. This constant expression of *sus-sh1* throughout the starchy endosperm, confirms previous experiments showing that *sh1* is up-regulated by anoxia (Zeng et al., 1998) and that developing endosperm is not anoxic (Rolletschek et al., 2005). In regard to up-regulation of hypoxia-related genes in MEC, the up-regulation of *sus1* and *sus2* in MEP is rather paradoxical since *sus1* is up-regulated by hypoxia (Zeng et al., 1998). However, it was previously shown that the oxygen concentration is already low in MEP and similar to that found in MEC, i.e., 1.4% of atmospheric saturation (Rolletschek et al., 2005). This means that (i) up-regulation of hypoxia-related genes in MEP reflects differences in the management of hypoxia rather than in the level of oxygen depletion and (ii) another signaling mechanism induces *sus* up-regulation. Actually, it was shown that *sus1* can be also up-regulated by sucrose (Koch et al., 1992) and, therefore, in the peripheral region of the maize endosperm, this could be due to the sucrose flux from phloem. Up-regulation of these two isoenzymes from 15 to 20 DAP could account for the lower sucrose accumulation in MEP than in MEC observed in late developing endosperm stage (30 DAP), in agreement with a higher SUS activity in MEP than in MEC (Wittich and Vreugdenhil, 1998). The higher sucrose degradation in MEP could supply the glycolysis-TCA pathway through the production of fructose in order to sustain the ATP production required for protein synthesis, particularly elevated in this area. This is in accordance with the fact that sucrose up-regulates both the expression of SUS (Koch et al., 1996) and storage proteins (Hattori et al., 1990). An increase of SUS-dependent sucrose breakdown leads to production of UDP-glucose, a substrate that could enter cell wall polysaccharide and starch biosynthesis pathways. Previous studies suggested that SUS1 is rather involved in starch biosynthesis pathway while SUS-SH1 is rather involved in cell wall biosynthesis (Chourey et al., 1998). Furthermore, the ectopic overexpression of a potato SUS in maize induces an increase of ADP-glucose and starch contents as well as of amylose/amylopectin mass ratio (Li et al., 2013). These findings are contradictory with the slightly lower starch contents in the vitreous endosperm when compared with the floury one but are in full agreement with the higher amylose content observed in vitreous maize endosperm (Gayral et al., 2015). Previous proteomic and transcriptomic analyses of glycolysis enzymes revealed the central role of pyruvate metabolism in the formation of a vitreous endosperm (Mechin et al., 2007; Prioul et al., 2008a). Indeed, in MEP, a part of the synthesized pyruvate is probably phosphorylated by PPDK to produce PEP and PPI, leading potentially to an indirect inhibition of starch synthesis.

This inhibition could account for the slightly lower amount of starch in the peripheral vitreous than in the central floury mature endosperm. However, up to now the roles of SUS2, of the corresponding SUS1-SUS2 hetero-oligomers (Duncan et al., 2006) and the complexes between SUS-SH1 and enzymes of the amyloplast starch synthesis pathway (Hennen-Bierwagen et al., 2009) are unknown and need to be finely explored in relation with starch and protein syntheses. Finally, our results emphasize the significance of the regulation of the sugar metabolism and especially of the energy metabolism linked to hypoxia in defining the protein/starch balance and *in fine* endosperm vitreousness.

The second type of stress response disclosed during this study is related to ER stress. In eukaryote cells, ER stress and the resulting response (UPR) regulate ER protein synthesis and degradation to limit the accumulation of unfolded/misfolded protein in the lumen (Howell, 2013; Unterseer et al., 2016). The expression of *bzip60* gene, as well as its splicing, a clear hallmark of such response, confirms that a higher UPR is activated in MEC. This process may impact transcription and protein translation, contributing hence to lower zein levels in MEC, which correlates with lower zein accumulation in floury endosperm (Gayral et al., 2016). UPR could also impact lipid metabolism, and especially, as observed in *floury2* (Shank et al., 2001) could be responsible for the higher amounts of storage lipids in the floury regions of conventional maize (Gayral et al., 2015). This was established in previous works showing that UPR is activated in wild-type maize endosperm, being stronger in many opaque mutants (Hunter et al., 2002; Morton et al., 2015). Moreover, ER stress is often associated with ROS (Blokhina and Fagerstedt, 2010) thus promoting a strong up-regulation of genes associated with oxido-reduction metabolism in MEC. Under these conditions, a higher level of hypoxia responses and ROS in MEC aggravates significantly the UPR process. In plants, ER stress also promotes autophagy (Liu et al., 2012; Pérez-Martín et al., 2014; Yang et al., 2016), up-regulated during maize endosperm development (Chung et al., 2009). MEC is characterized by stronger hypoxia, oxidative burst, and UPR. Moreover, in MEC, UPR is associated with precocious PCD. Although, there is no evidence linking UPR to PCD in maize endosperm, *o2* mutants which have stronger UPR display nevertheless a premature PCD (Young and Gallie, 2000a).

To summarize, the formation of endosperm vitreousness is linked to the capacity of the endosperm to respond to a strong developmental stress that is closely linked to oxygen depletion, regulating thus deposition of protein and starch within the maize endosperm. Among the responses to this developmental stress that controls starch/protein balance, in depth investigations should be done, in particular on the spatio-temporal regulation of SUS expression, oxidative and fermentative metabolisms and on the ER stress response, i.e., UPR, in relation with the interaction between both low oxygen and sucrose signaling.

## AUTHOR CONTRIBUTIONS

DM, MHM, SB, LL, and CD contributed to the design of the whole project; SB, LL, and CD produced and characterized the
plant material; MG and MD performed metabolite extraction and analyses; MG and KE performed PCR and qPCR analyses; SMB and SP designed and conducted microarray measurements; MG, SMB, and SP performed statistical analyses and gene annotation of microarray data; all co-authors contributed to the interpretation of data. MG, DM, KE, and BB wrote the manuscript

## FUNDING

This work was supported by grants from the Fond unique interministériel (FUI GranoFlakes, F1204004C). The platform POPS is supported by the Laboratory of

Excellence Saclay Plant Sciences (LabExSPS) (ANR-10-LABX-0040-SPS).

## ACKNOWLEDGMENTS

We thank Jean-Pierre Segain and Christian Bonnet for their assistance and advice regarding qRT-PCR analyses.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00557/full#supplementary-material

## REFERENCES

António, C., Päpke, C., Rocha, M., Diab, H., Limami, A. M., Obata, T., et al. (2016). Regulation of primary metabolism in response to low oxygen availability as revealed by carbon and nitrogen isotope redistribution. *Plant Physiol.* **170**, 43–56. doi: 10.1104/pp.15.00266

Bailey-Serres, J., Fukao, T., Gibbs, D. J., Holdsworth, M. J., Lee, S. C., Licausi, F., et al. (2012). Making sense of low oxygen sensing. *Trends Plant Sci.* **17**, 129–138. doi: 10.1016/j.tplants.2011.12.004

Blokhina, O., and Fagerstedt, K. V. (2010). Oxidative metabolism, ROS and NO under oxygen deprivation. *Plant Physiol. Biochem.* **48**, 359–373. doi: 10.1016/j.plaphy.2010.01.007

Chourey, P. S., Taliercio, E. W., Carlson, S. J., and Ruan, Y. L. (1998). Genetic evidence that the two isozymes of sucrose synthase present in developing maize endosperm are critical, one for cell wall integrity and the other for starch biosynthesis. *Mol. Gen. Genet.* **259**, 88–96.

Chung, T., Suttangkakul, A., and Vierstra, R. D. (2009). The ATG autophagic conjugation system in maize: ATG transcripts and abundance of the ATG8-lipid adduct are regulated by development and nutrient availability. *Plant Physiol.* **149**, 220–234. doi: 10.1104/pp.108.126714

Coleman, C. E., Clore, A. M., Ranch, J. P., Higgins, R., Lopes, M. A., and Larkins, B. A. (1997). Expression of a mutant α-zein creates the *floury2* phenotype in transgenic maize. *Proc. Natl. Acad. Sci. U.S.A.* **94**, 7094–7097. doi: 10.1073/pnas.94.13.7094

Danon, A., Rotari, V. I., Gordon, A., Mailhac, N., and Gallois, P. (2004). Ultraviolet-C overexposure induces programmed cell death in *Arabidopsis*, which is mediated by caspase-like activities and which can be suppressed by caspase inhibitors, p35 and Defender against Apoptotic Death. *J. Biol. Chem.* **279**, 779–787. doi: 10.1074/jbc.M304468200

Dante, R. A., Sabelli, P. A., Nguyen, H. N., Leiva-Neto, J. T., Tao, Y., Lowe, K. S., et al. (2014). Cyclin-dependent kinase complexes in developing maize endosperm: evidence for differential expression and functional specialization. *Planta* **239**, 493–509. doi: 10.1007/s00425-013-1990-1

de Sousa, S., del Giúdice Paniago, M., Arruda, P., and Yunes, J. A. (2008). Sugar levels modulate sorbitol dehydrogenase expression in maize. *Plant Mol. Biol.* **68**, 203–213. doi: 10.1007/s11103-008-9362-0

Dochlert, D. C. (1990). Distribution of enzyme activities within the developing maize (*Zea mays* L.) kernel in relation to starch, oil and protein accumulation. *Physiol. Plant* **78**, 560–567. doi: 10.1111/j.1399-3054.1990.tb05242.x

Doehlert, D. C., and Kuo, T. M. (1990). Sugar metabolism in developing kernels of starch-deficient endosperm mutants of maize. *Plant Physiol.* **92**, 990–994. doi: 10.1104/pp.92.4.990

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

Duncan, K. A., Hardin, S. C., and Huber, S. C. (2006). The three maize sucrose synthase isoforms differ in distribution, localization, and phosphorylation. *Plant Cell Physiol.* **47**, 959–971. doi: 10.1093/pcp/pcj068

Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. *Nucleic Acids Res.* **30**, 207–210. doi: 10.1093/nar/30.1.207

Fagundes, D., Bohn, B., Cabreira, C., Leipelt, F., Dias, N., Bodanese-Zanettini, M. H., et al. (2015). Caspases in plants: metacaspase gene family in plant stress responses. *Funct. Integr. Genomics.* **15**, 639–649. doi: 10.1007/s10142-015-0459-7

Felle, H. H. (2005). pH regulation in anoxic plants. *Ann. Bot.* **96**, 519–532. doi: 10.1093/aob/mci207

Fukao, T., and Bailey-Serres, J. (2004). Plant responses to hypoxia-is survival a balancing act? *Trends Plant Sci.* **9**, 449–456. doi: 10.1016/j.tplants.2004.07.005

Gagnot, S., Tamby, J. P., Martin-Magniette, M. L., Bitton, F., Taconnat, L., Balzergue, S., et al. (2008). CATdb: a public access to *Arabidopsis* transcriptome data from the URGV-CATMA platform. *Nucleic Acids Res.* **36**, 986–990. doi: 10.1093/nar/gkm757

Gayral, M., Bakan, B., Dalgalarrondo, M., Elmorjani, K., Delluc, C., Brunet, S., et al. (2015). Lipid partitioning in maize (*Zea mays* L.) endosperm highlights relationships among starch lipids, amylose, and vitreousness. *J. Agric. Food Chem.* **63**, 3551–3558. doi: 10.1021/acs.jafc.5b00293

Gayral, M., Gaillard, C., Bakan, B., Dalgalarrondo, M., Elmorjani, K., Delluc, C., et al. (2016). Transition from vitreous to floury endosperm in maize (*Zea mays* L.) kernels is related to protein and starch gradients. *J. Cereal Sci.* **68**, 148–154. doi: 10.1016/j.jcs.2016.01.013

Geigenberger, P. (2014). "Adaptation of storage metabolism to oxygen deprivation," in *Low-Oxygen Stress in Plants: Oxygen Sensing and Adaptive Responses to Hypoxia*, eds T. J. van Dongen and F. Licausi. (Vienna: Springer Vienna), 223–244.

Gibbon, B. C., and Larkins, B. A. (2005). Molecular genetic approaches to developing Quality Protein Maize. *Trends Genet.* **21**, 227–233. doi: 10.1016/j.tig.2005.02.009

Gibbon, B. C., Wang, X., and Larkins, B. A. (2003). Altered starch structure is associated with endosperm modification in Quality Protein Maize. *Proc. Natl. Acad. Sci. U.S.A.* **100**, 15329–15334. doi: 10.1073/pnas.2136854100

Guo, X., Ronhovde, K., Yuan, L., Yao, B., Soundararajan, M. P., Elthon, T., et al. (2012). Pyrophosphate-dependent fructose-6-phosphate 1-phosphotransferase induction and attenuation of Hsp gene expression during endosperm modification in Quality Protein Maize. *Plant Physiol.* **158**, 917–929. doi: 10.1104/pp.111.191163

Hattori, T., Nakagawa, S., and Nakamura, K. (1990). High-level expression of tuberous root storage protein genes of sweet potato in stems of plantlets grown *in vitro* on sucrose medium. *Plant Mol. Biol.* **14**, 595–604. doi: 10.1007/BF00027505

Hennen-Bierwagen, T. A., Lin, Q., Grimaud, F., Planchot, V., Keeling, P. L., James, M. G., et al. (2009). Proteins from multiple metabolic pathways associate with starch biosynthetic enzymes in high molecular weight complexes: a model for regulation of carbon allocation in maize amyloplasts. *Plant Physiol.* **149**, 1541–1559. doi: 10.1104/pp.109.135293
Holding, D. R., Otegui, M. S., Li, B., Meeley, R. B., Dam, T., Hunter, B. G., et al. (2007). The maize flouryl1 gene encodes a novel endoplasmic reticulum protein involved in zein protein body formation. *Plant Cell* 19, 2569–2582. doi: 10.1105/tpc.107.053538

Howell, S. H. (2013). Endoplasmic reticulum stress responses in plants. *Annu. Rev. Plant Biol.* 64, 477–499. doi: 10.1146/annurev-arplant-050312-120053

Hunter, B. G., Beatty, M. K., Singletary, G. W., Hamaker, B. R., Dilkes, B. P., Larkins, B. A., et al. (2002). Maize opaque endosperm mutations create extensive changes in patterns of gene expression. *Plant Cell* 14, 2591–2612. doi: 10.1105/tpc.003905

Jeanneau, M., Vidal, J., Gousset-Dupont, A., Lebouteiller, B., Hodges, M., Gerentes, D., et al. (2002). Manipulating PEPC levels in plants. *J. Exp. Bot.* 53, 1837–1845. doi: 10.1093/jxb/erf061

Kim, C. S., Gibbon, B. C., Gillikin, J. W., Larkins, B. A., Boston, R. S., and Jung, R. (2006). The maize Mucronate mutation is a deletion in the 16-kDa gamma-zein gene that induces the unfolded protein response. *Plant J.* 48, 440–451. doi: 10.1111/j.1365-313X.2006.02884.x

Kim, C. S., Hunter, B. G., Kraft, J., Boston, R. S., Yans, S., Jung, R., et al. (2004). A defective signal peptide in a 19-kD alpha-zein protein causes the unfolded protein response and an opaque endosperm phenotype in the maize De-B30 mutant. *Plant Physiol.* 134, 380–387. doi: 10.1104/pp.103.031310

Koch, K. E., Nolte, K. D., Duke, E. R., McCarty, D. R., and Avigne, W. T. (1992). Sugar levels modulate differential expression of maize sucrose synthase genes. *Plant Cell* 4, 59–69. doi: 10.1105/tpc.4.1.59

Koch, K. E., Wu, Y., and Xu, J. (1996). Sugar and metabolic regulation of genes for sucrose metabolism: potential influence of maize sucrose synthase and soluble invertase responses on carbon partitioning and sugar sensing. *J. Exp. Bot.* 47, 1179–1185. doi: 10.1093/jxb/47

Laibach, N., Post, J., Twyman, R. M., Gronover, C. S., and Prüfer, D. (2015). The characteristics and potential applications of structural lipid droplet proteins in plants. *J. Biotechnol.* 201, 15–27. doi: 10.1016/j.jbiotec.2014.08.020

Lam, E., and del Pozo, O. (2000). Caspase-like protease involvement in the control of plant cell death. *Plant Mol. Biol.* 44, 173–184. doi: 10.1007/978-94-010-0934-8_13

Li, F., Chung, T., Pennington, J. G., Federico, M. L., Kaeppler, H. F., Kaeppler, S. M., et al. (2015). Autophagic recycling plays a central role in maize nitrogen remobilization. *Plant Cell* 27, 1389–1408. doi: 10.1105/tpc.15.00158

Li, J., Baroja-Fernández, E., Bahaji, A., Muñoz, F. J., Ovecka, M., Montero, M., et al. (2013). Enhancing sucrose synthase activity results in increased levels of starch and ADP-glucose in maize (*Zea mays* L.) seed endosperms. *Plant Cell Physiol.* 54, 282–294. doi: 10.1093/pcp/pcs180

Li, Y., Humbert, S., and Howell, S. H. (2012). ZmbZIP60 mRNA is spliced in maize in response to ER stress. *BMC Res. Notes* 5:144. doi: 10.1186/1756-0500-5-144

Liu, Y., Burgos, J. S., Deng, Y., Srivastava, R., Howell, S. H., and Bassham, D. C. (2012). Degradation of the endoplasmic reticulum by autophagy during endoplasmic reticulum stress in *Arabidopsis*. *Plant Cell* 24, 4635–4651. doi: 10.1105/tpc.112.101535

Locatelli, S., Piatti, P., Motto, M., and Rossi, V. (2009). Chromatin and DNA modifications in the Opaque2-mediated regulation of gene transcription during maize endosperm development. *Plant Cell* 21, 1410–1427. doi: 10.1105/tpc.109.067256

Mechin, V., Thevenot, C., Le Guilloux, M., Prioul, J. L., and Damerval, C. (2007). Developmental analysis of maize endosperm proteome suggests a pivotal role for pyruvate orthophosphate dikinase. *Plant Physiol.* 143, 1203–1219. doi: 10.1104/pp.106.092148

Minina, E. A., Filonova, L. H., Fukada, K., Savenkov, E. I., Gogvadze, V., Clapham, D., et al. (2013). Autophagy and metacaspase determine the mode of cell death in plants. *J. Cell Biol.* 203, 917–927. doi: 10.1083/jcb.201307082

Morton, K. J., Jia, S., Zhang, C., and Holding, D. R. (2015). Proteomic profiling of maize opaque endosperm mutants reveals selective accumulation of lysine-enriched proteins. *J. Exp. Bot.* 67, 1381–1396. doi: 10.1093/jxb/erv532

Myers, A. M., James, M. G., Lin, Q., Yi, G., Stinard, P. S., Hennen-Bierwagen, T. A., et al. (2011). Maize opaque5 encodes monogalactosyldiacylglycerol synthase and specifically affects galactolipids necessary for amyloplast and chloroplast function. *Plant Cell* 23, 2331–2347. doi: 10.1105/tpc.111.087205

Nawkar, G. M., Maibam, P., Park, J. H., Woo, S. G., Kim, C. Y., Lee, S. Y., et al. (2016). *In silico* study on *Arabidopsis* BAG gene expression

in response to environmental stresses. *Protoplasma* 254, 409–421. doi: 10.1007/s00709-016-0961-3

Pérez-Martín, M., Pérez-Pérez, M. E., Lemaire, S. D., and Crespo, J. L. (2014). Oxidative stress contributes to autophagy induction in response to endoplasmic reticulum stress in *Chlamydomonas reinhardtii*. *Plant Physiol.* 166, 997–1008. doi: 10.1104/pp.114.243659

Prioul, J. L., Mechín, V., and Damerval, C. (2008a). Molecular and biochemical mechanisms in maize endosperm development: the role of pyruvate-Pi-dikinase and opaque2 in the control of C/N ratio. *C. R. Biol.* 331, 772–779. doi: 10.1016/j.crvi.2008.07.019

Prioul, J. L., Méchin, V., Lessard, P., Thévenot, C., Grimmer, M., Chateau-Joubert, S., et al. (2008b). A joint transcriptomic, proteomic and metabolic analysis of maize endosperm development and starch filling. *Plant Biotechnol. J.* 6, 855–869. doi: 10.1111/j.1467-7652.2008.00368.x

Qiao, Z., Qi, W., Wang, Q., Feng, Y., Yang, Q., Zhang, N., et al. (2016). ZmMADS47 regulates zein gene transcription through interaction with opaque2. *PLoS Genet.* 12:e1005991. doi: 10.1371/journal.pgen.1005991

Rolletschek, H., Koch, K., Wobus, U., and Borisjuk, L. (2005). Positional cues for the starch/lipid balance in maize kernels and resource partitioning to the embryo. *Plant J.* 42, 69–83. doi: 10.1111/j.1365-313X.2005.02352.x

Salazar-Salas, N. Y., Pineda-Hidalgo, K. V., Chavez-Ontiveros, J., Gutierrez-Dorado, R., Reyes-Moreno, C., Bello-Pérez, L. A., et al. (2014). Biochemical characterization of QTLs associated with endosperm modification in Quality Protein Maize. *J. Cereal Sci.* 60, 255–263. doi: 10.1016/j.jcs.2014.04.004

Schmidt, R. J., Burr, F. A., Aukerman, M. J., and Burr, B. (1990). Maize regulatory gene opaque-2 encodes a protein with a leucine-zipper motif that binds to zein DNA. *Proc. Natl. Acad. Sci. U.S.A.* 87, 46–50. doi: 10.1073/pnas.87.1.46

Schmidt, R. J., Burr, F. A., and Burr, B. (1987). Transposon tagging and molecular analysis of the maize regulatory locus opaque2. *Science* 238, 960–963. doi: 10.1126/science.2823388

Sen, T. Z., Andorf, C. M., Schaeffer, M. L., Harper, L. C., Sparks, M. E., Duvick, J., et al. (2009). MaizeGDB becomes “sequence-centric.” *Database* 2009:bap020. doi: 10.1093/database/bap020

Shank, K. J., Su, P., Brglez, I., Boss, W. F., Dewey, R. E., and Boston, R. S. (2001). Induction of lipid metabolic enzymes during the endoplasmic reticulum stress response in plants. *Plant Physiol.* 126, 267–277. doi: 10.1104/pp.126.1.267

Smyth, G. K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.* 3, Article3. doi: 10.2202/1544-6115.1027

Storey, J. D., and Tibshirani, R. (2003). Statistical significance for genomewide studies. *Proc. Natl. Acad. Sci. U.S.A.* 100, 9440–9445. doi: 10.1073/pnas.1530509100

Suzuki, H., Osawa, T., Fujioka, Y., and Noda, N. N. (2017). Structural biology of the core autophagy machinery. *Curr. Opin. Struct. Biol.* 43, 10–17. doi: 10.1016/j.sbi.2016.09.010

Tsatsiani, L., Van Breusegem, F., Gallois, P., Zavialov, A., Lam, E., and Bozhkov, P. V. (2011). Metacaspases. *Cell Death Differ.* 18, 1279–1288. doi: 10.1038/cdd.2011.66

Unterseer, S., Pophaly, S. D., Peis, R., Westermeier, P., Mayer, M., Seidel, M. A., et al. (2016). A comprehensive study of the genomic differentiation between temperate Dent and Flint maize. *Genome Biol.* 17:137. doi: 10.1186/s13059-016-1009-x

Wang, B., Zheng, J., Liu, Y., Wang, J., and Wang, G. (2012). Cloning and characterization of the stress-induced bZIP gene ZmbZIP60 from maize. *Mol. Biol. Rep.* 39, 6319–6327. doi: 10.1007/s11033-012-1453-y.

Wang, G., Qi, W., Wu, Q., Yao, D., Zhang, J., Zhu, J., et al. (2014). Identification and characterization of maize *floury4* as a novel semidominant opaque mutant that disrupts protein body assembly. *Plant Physiol.* 165, 582–594. doi: 10.1104/pp.114.238030

Wang, G., Wang, F., Wang, G., Wang, F., Zhang, X., Zhong, M., et al. (2012). Opaque1 encodes a myosin XI motor protein that is required for endoplasmic reticulum motility and protein body formation in maize endosperm. *Plant Cell* 24, 3447–3462. doi: 10.1105/tpc.112.101360

Watanabe, N., and Lam, E. (2008). BAX inhibitor-1 modulates endoplasmic reticulum stress-mediated programmed cell death in *Arabidopsis*. *J. Biol. Chem.* 283, 3200–3210. doi: 10.1074/jbc.M706659200
Wittich, P. E., and Vreugdenhil, D. (1998). Localization of sucrose synthase activity in developing maize kernels by *in situ* enzyme histochemistry. *J. Exp. Bot.* 49, 1163-1171. doi: 10.1093/jxb/49.324.1163

Woo, Y. M., Wang-Nan Hu, D., Larkins, B. A., and Jung, R. (2001). Genomics analysis of genes expressed in maize endosperm identifies novel seed proteins and clarifies patterns of zein gene expression. *Plant Cell* 13, 2297-2317. doi: 10.1105/tpc.13.10.2297

Wu, Y., Holding, D. R., and Messing, J. (2010). γ-zeins are essential for endosperm modification in Quality Protein Maize. *Proc. Natl. Acad. Sci. U.S.A.* 107, 12810-12815. doi: 10.1073/pnas.1004721107

Wu, Y., and Messing, J. (2010). RNA interference-mediated change in protein body morphology and seed opacity through loss of different zein proteins. *Plant Physiol.* 153, 337-347. doi: 10.1104/pp.110.154690

Xu, S. B., Li, T., Deng, Z. Y., Chong, K., Xue, Y., and Wang, T. (2008). Dynamic proteomic analysis reveals a switch between central carbon metabolism and alcoholic fermentation in rice filling grains. *Plant Physiol.* 148, 908-925. doi: 10.1104/pp.108.125633

Yang, X., Srivastava, R., Howell, S. H., and Bassham, D. C. (2016). Activation of autophagy by unfolded proteins during endoplasmic reticulum stress. *Plant J.* 85, 83-95. doi: 10.1111/tpj.13091

Yang, Y. H., Dudoit, S., Luu, P., Lin, D. M., Peng, V., Ngai, J., et al. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. *Nucleic Acids Res.* 30:e15. doi: 10.1093/nar/30.4.e15

Young, T. E., and Gallie, D. R. (2000a). Programmed cell death during endosperm development. *Plant Mol. Biol.* 44, 283-301. doi: 10.1007/978-94-010-0934-8_4

Young, T. E., and Gallie, D. R. (2000b). Regulation of programmed cell death in maize endosperm by abscisic acid. *Plant Mol. Biol.* 42, 397-414. doi: 10.1023/A:1006333103342

Yu, Y., Mu, H. H., Wasserman, B. P., and Carman, G. M. (2001). Identification of the maize amyloplast stromal 112-kD protein as a plastidic starch phosphorylase. *Plant Physiol.* 125, 351-359. doi: 10.1104/pp.125.1.351

Zabalza, A., Van Dongen, J. T., Froehlich, A., Oliver, S. N., Faix, B., Gupta, K. J., et al. (2009). Regulation of respiration and fermentation to control the plant internal oxygen concentration. *Plant Physiol.* 149, 1087-1098. doi: 10.1104/pp.108.129288

Zeng, Y., Wu, Y., Avigne, W. T., and Koch, K. E. (1998). Differential regulation of sugar-sensitive sucrose synthases by hypoxia and anoxia indicate complementary transcriptional and posttranscriptional responses (1998). *Plant Physiol.* 116, 1573-1583. doi: 10.1104/pp.116.4.1573

Zhang, Z., Yang, J., and Wu, Y. (2015). Transcriptional regulation of zein gene expression in maize through the additive and synergistic action of *opaque2*, prolamine-box binding factor, and *o2* heterodimerizing proteins. *Plant Cell* 27, 1162-1172. doi: 10.1105/tpc.15.00035

Zhen, S., Dong, K., Deng, X., Zhou, J., Xu, X., Han, C., et al. (2016). Dynamic metabolome profiling reveals significant metabolic changes during grain development of bread wheat (Triticum aestivum L.). *J. Sci. Food Agric.* 96, 3731-3740. doi: 10.1002/jsfa.7561

**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 © 2017 Gayral, Elmorjani, Dalgarrondo, Balzergue, Pateyron, Morel, Brunet, Linossier, Delluc, Bakan and Marion. 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) or licensor 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.