# An Integrated “Multi-Omics” Comparison of Embryo and Endosperm Tissue-Specific Features and Their Impact on Rice Seed Quality

Marc M. Galland, Dongli He, Imen I. Lounifi, Erwann Arc, Gilles Clément, Sandrine Balzergue, Stéphanie Huguet, Gwendal Cueff, Béatrice Godin, Boris B. Collet, et al.

► To cite this version:

Marc M. Galland, Dongli He, Imen I. Lounifi, Erwann Arc, Gilles Clément, et al. An Integrated “Multi-Omics” Comparison of Embryo and Endosperm Tissue-Specific Features and Their Impact on Rice Seed Quality. Frontiers in Plant Science, 2017, 8, pp.1984. 10.3389/fpls.2017.01984. hal-01645144

HAL Id: hal-01645144

https://hal.science/hal-01645144

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>

# An Integrated “Multi-Omics” Comparison of Embryo and Endosperm Tissue-Specific Features and Their Impact on Rice Seed Quality

OPEN ACCESS

**Edited by:**
Glória Catarina Pinto,
University of Aveiro, Portugal

**Reviewed by:**
Koh Aoki,
Osaka Prefecture University, Japan
Monica Escandon,
University of Aveiro, Portugal

*Correspondence:
Loïc Rajjou
loic.rajjou@agroparistech.fr

†Present Address:
Marc Galland,
Department of Plant Physiology,
Swammerdam Institute for Life
Sciences, University of Amsterdam,
Amsterdam, Netherlands
Erwann Arc,
Institute of Botany, University of
Innsbruck, Innsbruck, Austria
Sandrine Balzergue,
Institut de Recherche en Horticulture
et Semences (Institut National de la
Recherche Agronomique,
AgroCampus Ouest, Université
d'Angers), Beaucouzé, France
Joseph Tran and Halima Morin,
Institute of Plant Sciences
Paris-Saclay, Institut National de la
Recherche Agronomique, Centre
National de la Recherche Scientifique,
Université Paris-Sud, Université
d'Evry, Université Paris-Diderot,
Sorbonne Paris-Cité, Université
Paris-Saclay, Orsay, France
Benoit Valot,
UMR 6249 Chrono-Environnement,
Haut de Chazal, Ambroise Paré,
Besançon, France

Marc Galland¹†, Dongli He², Imen Lounifi¹, Erwann Arc¹†, Gilles Clément¹, Sandrine Balzergue³†, Stéphanie Huguet³, Gwendal Cueff¹, Béatrice Godin¹, Boris Collet¹, Fabienne Granier¹, Halima Morin¹†, Joseph Tran¹†, Benoit Valot⁴† and Loïc Rajjou¹*

¹ JPB, Institut Jean-Pierre Bourgin (INRA, AgroParisTech, CNRS, Université Paris-Saclay), Saclay Plant Sciences (SPS), Versailles, France, ² Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, China, ³ IPS2, Institute of Plant Sciences Paris-Saclay (INRA, CNRS, Université Paris-Sud, Université d'Evry, Université Paris-Diderot, Sorbonne Paris-Cité, Université Paris-Saclay), POPS-Transcriptomic Platform, Saclay Plant Sciences (SPS), Orsay, France, ⁴ GQE-Le Moulon, Génétique Quantitative et Evolution (INRA Université Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay), PAPPSO-Plateforme d'Analyse Protéomique de Paris Sud-Ouest, Saclay Plant Sciences (SPS), Gif-sur-Yvette, France

Although rice is a key crop species, few studies have addressed both rice seed physiological and nutritional quality, especially at the tissue level. In this study, an exhaustive “multi-omics” dataset on the mature rice seed was obtained by combining transcriptomics, label-free shotgun proteomics and metabolomics from embryo and endosperm, independently. These high-throughput analyses provide a new insight on the tissue-specificity related to rice seed quality. Foremost, we pinpointed that extensive post-transcriptional regulations occur at the end of rice seed development such that the embryo proteome becomes much more diversified than the endosperm proteome. Secondly, we observed that survival in the dry state in each seed compartment depends on contrasted metabolic and enzymatic apparatus in the embryo and the endosperm, respectively. Thirdly, it was remarkable to identify two different sets of starch biosynthesis enzymes as well as seed storage proteins (glutelins) in both embryo and endosperm consistently with the supernumerary embryo hypothesis origin of the endosperm. The presence of a putative new glutelin with a possible embryonic favored abundance is described here for the first time. Finally, we quantified the rate of mRNA translation into proteins. Consistently, the embryonic panel of protein translation initiation factors is much more diverse than that of the endosperm. This work emphasizes the value of tissue-specificity-centered “multi-omics” study in the seed to highlight new features even from well-characterized pathways. It paves the way for future studies of critical genetic determinants of rice seed physiological and nutritional quality.

**Keywords:** rice, seed, endosperm, embryo, multi-omics, translation, starch, glutelins
# INTRODUCTION

Seeds of high nutritional and physiological quality are essential for the benefit of mankind. Seeds are also at the forefront for the preservation of biodiversity through the plant conservation strategies in seed banks (Wyse Jackson and Kennedy, 2009; Westengen et al., 2013; Hay and Whitehouse, 2017). The quality of seeds comprises physiological, ecological and nutritional traits for agriculture, agroecology, and agro-food system. In terms of botanical provenance, the mature seed of Angiosperms is a patchwork of maternal and filial tissues (Walbot and Evans, 2003; Olsen, 2004; Nowack et al., 2010; Lafon-Placette and Kohler, 2014). The double fertilization of the haploid egg cell and diploid central cell by the two haploid sperm cells gives rise to a diploid embryo and to a triploid endosperm respectively. In addition, the seed coat (testa) is composed of several cell layers coming from the mother plant ovule and ovary. In cereals (*Poaceae*), the embryo is composed of the embryonic axis surrounded by a single cotyledon (scutellum) and it will form the future seedling. The mature endosperm is composed of four differentiated regions: the endosperm transfer cell (ETC) region, the embryo-surrounding region (ESR), the aleurone layer (AL), and the starchy endosperm (SE) cells (Olsen, 2004). While the ETC and AL cells remain alive at the end of cereal seed development, most of the ESR and SE cells have undergone programmed cell death (PCD) with characteristic DNA laddering and organelle degradation.

From an evolutionary point of view, one hypothesis stipulates that the endosperm would have derived from a supernumerary embryo originating from double fertilization and that would have evolved into an embryo-supporting structure (Friedman, 1998). Subsequently, the endosperm has evolved multiple roles related to the embryo. First, the endosperm protects the embryo from atmospheric oxygen that eventually leads to the formation of hydroperoxides and cell death (De Giorgi et al., 2015). In addition, critical cross-talk between abscisic acid (ABA) and gibberellin (GA) regulating seed development, size, dormancy or storage breakdown during germination are also the results of endosperm—embryo interactions (Fincher, 1989; Penfield et al., 2006; Bethke et al., 2007; Folsom et al., 2014; Yan et al., 2014; Bassel, 2016). Still, so far, few reports have addressed seed quality issues in terms of functional tissue specialization.

During seed maturation, orthodox seeds acquire desiccation tolerance (DT) and storability (longevity) as defined by the ability to withstand extreme water loss (to values around 0.1 g water per gram of dry weight) and to survive in the dry state (Hoekstra et al., 2001; Alpert, 2005). Among the mechanisms

that promote DT, the formation of a glassy cytoplasm and the subsequent decrease in molecular mobility is positively correlated with seed longevity (Buitink and Leprince, 2008). During the late phase of seed maturation, the accumulation of late embryogenesis abundant (LEA) proteins, heat shock proteins (HSPs), antioxidants and non-reducing sugars all together contribute to glassy state establishment (Boudet et al., 2006; Farrant and Moore, 2011; Kaur et al., 2016; Sano et al., 2016). How different tissues cooperate to establish seed DT and storability is still unclear especially in cereals with a persistent endosperm. Among the few existing studies, a proteomic analysis on maize *viviparous5* (*vp5*) mutant showed that LEA and HSPs were affected in the ABA-deficient embryos of the *vp5* mutants contrarily to the endosperm (Wu et al., 2014). Furthermore, a consequence of lipid degradation during storage, detoxification of the lipid peroxidation by-product malondialdehyde (MDA) by the rice aldehyde dehydrogenase 7 (OsALDH7) proved to be essential for DT (Shin et al., 2009). The null rice mutants of OsALDH7 showed increased MDA resulting in reduced seed viability (Shin et al., 2009). Yet, the exact site of MDA generation during dry storage and seed aging remains unknown. The impact of active lipid degradation on rice seed quality was further reinforced by transgenic analysis of two rice lipoxygenases OsLOX2 and OsLOX3 since silencing and overexpression of these two enzymes acts in opposite directions on seed germination and longevity (Huang et al., 2014; Xu et al., 2015). Interestingly, the suppression of OsLOX3 expression in the rice endosperm improved resistance to seed aging (Xu et al., 2015). The endosperm outer cuticle layer permeability can also preserve the embryonic components from oxidation by atmospheric oxygen (De Giorgi et al., 2015). A recent paper showed that during rice seed aging, the endosperm and embryo were differentially affected by seed aging in particular regarding glycolytic enzymes that decreased in abundance in the endosperm while increasing in the embryo (Zhang et al., 2016). Last but not least, seed storage proteins can also buffer the oxidative stress caused by seed aging as shown for Arabidopsis cruciferins (Rajjou et al., 2008; Nguyen et al., 2015). Despite these relevant studies, the majority of the molecular determinants of both DT and seed longevity in different seed tissues are still to be established especially in the embryo of cereal crops.

The generation of energy, digestion of seed storage proteins (SSPs), carbohydrates and/or triacylglycerols (TAGs) proved to be of paramount importance to obtain highly vigorous seeds. Within a few hours of imbibition, the seed embryo is rapidly resuming respiration thanks to the presence of a functional electron transport chain in undifferentiated pro-mitochondria (Ehrenshaft and Brambl, 1990; Logan et al., 2001; Howell et al., 2006). Later on, the differentiation of these pro-mitochondria in fully functional mitochondria participates to full metabolic resumption through establishment of the tricarboxylic acid cycle (TCA) (Lawlor and Vince, 2014). In contrast, it is less clear how the AL cells produce the requested energy to synthesize the large amounts of α-amylases for starch degradation (Fincher, 1989). During barley seed development, the inner SE is mildly to severely hypoxic while the AL is not (Rolletschek et al., 2011). Still, the precise investigation of oxygen requirement and

<footnote>**Abbreviations:** AL, aleurone layer; ALDH, aldehyde dehydrogenase; AsA, ascorbate; DHA, dehydroascorbate; DT, desiccation tolerance; eIF, eukaryotic initiation factor; ESR, embryo surrounding region; ETC, endosperm transfer cells; HAI, hours after imbibition; HSP, heat shock protein; LEA, late embryogenesis abundant; LOX, lipoxygenase; MDA, malondialdehyde; MDHA, monodehydroascorbate; MSR, methionine sulfoxide oxidase; NBB, naphthol blue black; PA, periodic acid; PCD, programmed cell death; PIMT, protein-L-isoaspartate O-methyltransferase; PLD, phospholipase D; QTL, quantitative trait locus; RFO, raffinose family oligosaccharides; ROS, reactive oxygen species; SE, starchy endosperm; SSPs, seed storage proteins; TAGs, triacylglycerols; TCA, tricarboxylic acid cycle; TIC, total ion current; XIC, extracted ion current.</footnote>
consumption during germination in the different cereal seed tissues remains to be established.

The metabolic system for mRNA translation and protein synthesis is the most energy-requiring process in most organisms. Through the use of the translational inhibitor cycloheximide, it has been shown that translation of stored mRNAs is necessary and sufficient for both Arabidopsis and rice germination (Rajjou et al., 2004; Sano et al., 2012, 2013). On the other hand, DNA transcription by RNA polymerase II is necessary for seed vigor and seedling growth (Rajjou et al., 2004; Sano et al., 2012). In cereal seeds, the importance of mRNA translation in the AL during germination is crucial. Indeed, a strong synthesis and accumulation of α-amylases, in response to gibberellic acid, participates to starch mobilization that in turns fuel the germinating embryo with oligosaccharides (Fincher, 1989). Interestingly, the scutellum also provides α-amylases contributing to the amylolytic activity of the starch endosperm during germination (Subbarao et al., 1998) suggesting a close relationship between endosperm tissues and embryo. Nevertheless, only few studies have addressed the precise content of the cereal endosperm in terms of translational machinery especially compared to that of the cereal embryo that was historically used as a model cell-free system to translate mRNAs *in vitro* (Takai and Endo, 2010).

As a model species for cereals with a well-annotated genome (International Rice Genome Sequencing Project, 2005), rice seeds are widely studied by taking advantage of “omics” approaches (Koller et al., 2002; Tarpley et al., 2005; Xu et al., 2008; Jiao et al., 2009; Wang et al., 2010; Lee and Koh, 2011; Nguyen et al., 2012; Xue et al., 2012). Most studies on cereal seed focused either on the isolated embryo (Howell et al., 2009; Kim et al., 2009; Han et al., 2014) or on the whole seed (Yang et al., 2007). In contrast, a small number of works were performed on both the embryo and endosperm in the same experiment (Gallardo et al., 2007; Sreenivasulu et al., 2008).

Here, we used an integrated “multi-omics” approach combining transcriptomics, label-free quantitative shotgun proteomics and gas chromatography coupled to mass spectrometry (GC-MS)-based metabolomics on dry mature rice seeds from a reference rice cultivar (*Oryza sativa* ssp *japonica* cv Nipponbare). The present work was firstly aimed at comparing the compartmentalization of nutritionally relevant pathways between the endosperm and embryo (nutritional quality). These pathways could be further improved through metabolic engineering based on knowledge of the fine composition of both embryo and endosperm. Secondly, we highlighted genes potentially important for seed storability and germination (agricultural quality). Thirdly, a targeted study of the factors associated with the seed storage compounds (starch and proteins) were analyzed carefully. Altogether, these exhaustive datasets emphasize determinants of rice seed quality in a tissue-specific manner.

# METHODS

## Rice Biological Material

Dry mature rice seeds (*Oryza sativa* ssp *japonica* cv Nipponbare) were harvested in September 2011 at the “Centre Français du Riz”

(Mas du Sonnailler, Arles, France). At the lab, seeds were dehulled and dissected in one embryo (E0) and one endosperm (A0) fraction with a sharp scalpel. The white rice (SE) fraction was obtained thanks to a lab bench rice milling machine as previously described (Galland et al., 2014a). Dry weight was determined on 10 bulks comprising 10 rice seeds. The corresponding E0, A0, and SE were placed in a dry oven at 105°C for 48 h and weighted on a precision lab balance (XP204, Mettler-Toledo, France).

## Optical and Confocal Microscopy

Fixation, inclusion into historesin and cutting of 5 µm semi-thin sections were done exactly as previously described (Galland et al., 2014a). Proteins and complex carbohydrates were revealed using a Periodic Acid (PA)—Schiff /Naphthol Blue Black (NBB) staining. Samples semi-thin sections (5 µm) were first hydrolyzed 5 min in 1% periodic acid (w/v), rinsed with tap water and distilled water and then colored with Schiff’s reagent for 10 min in complete darkness. Subsequently, sample slices were washed with sulfurous water that contains 5% (w/v) sodium metabisulfite, 250 mM HCl and distilled water for 1 min before washing with tap water and distilled water. NBB staining was done using a pre-heated (65°C) working solution that contains 0.1% (w/v) NBB, 10% (v/v) acetic acid and distilled water in which samples were placed for 5 min. One final thorough washing was done in tap water and samples sections were finally placed in acetic acid (7%, v/v) for 1 min or more if PA-Schiff staining was too weak. Finally, samples sections were mounted between glass slides in glycerol (Histomount, National Diagnostics, UK) for visualization before imaging with Leica optical microscopy (Leica Zeiss Axioplan, Leica Microsystems, Wetzlar, Germany).

Neutral lipids were imaged by confocal scanning fluorescence microscopy (Leica TCS SP2, Leica Microsystèmes SAS, France) using the Nile red dye (Greenspan et al., 1985). For neutral lipid observation, 100 µm wide sections of rice dry mature seeds were cut using a vibrating blade microtome (Leica VT1000 S, Leica Microsystèmes SAS, France) in sterile distilled water. Then, sections were quickly put on a glass slide in 100 µl of a Nile Red solution that contained 0.1% Nile Red (w/v) in 50% glycerol. The cell walls were counterstained by adding 100 µl of a Calcofluor solution that contained 1% Fluorescent Brightener 28 (w/v) in a carbonate/bicarbonate buffer pH 9.2. For Nile Red imaging, 488 nm was used for excitation and emission was collected between 593 and 654 nm. For Calcofluor imaging, 405 nm was used for excitation and emission was collected between 412 and 483 nm.

## Metabolome Analysis by Gas Chromatography Coupled to Mass Spectrometry (GC-MS)

For each tissue (i.e., embryo, endosperm), metabolite samples were obtained starting from three replicates of 100 rice seeds manually dissected in embryo and endosperm fractions. 100 embryos and 100 endosperms were grinded with mortar and pestle under liquid nitrogen for the embryos and with a Cyclotec™ 1093 Sample Mill (FOSS, Hillerød, Danemark) for the endosperms. All samples were lyophilized and around 20 mg dry weight (DW) of each sample were placed in 2 ml Safelock Eppendorf tubes (Eppendorf AG, Hamburg, Germany).
All analysis steps including extraction, derivatization, analysis, and data processing were adapted from the original protocol described by Fiehn et al. (2008) and following the procedure described by Avila-Ospina et al. (2017). The extraction solvent was prepared by mixing water:acetonitrile:isopropanol at the volume ratio 2:3:3 allowing to extract metabolites with a broad range of polarities. For derivatization step, N-methyl-N-trimethylsilyl-trifluoroacetamide (MSTFA; Sigma-Aldrich) was used in silylation procedure of metabolites. Samples were analyzed on an Agilent 7890A gas chromatograph coupled to an Agilent 5975C mass spectrometer. Raw Agilent datafiles were converted in NetCDF format and analyzed with AMDIS (Automated Mass Deconvolution and Identification System; http://chemdata.nist.gov/mass-spc/amdis/). A home retention indices/ mass spectra library built from the NIST, Golm, and Fiehn databases and standard compounds was used for metabolites identification. Peak areas were then determined using the QuanLynx software (Waters, Milford, USA) after conversion of the NetCDF file in MassLynx format.

## RNA Isolation and Microarray Analyses

Total mRNAs were isolated from three replicates of 100 embryos and 50 endosperms and hybridizations on the Affymetrix GeneChip® Rice Genome Array (Affymetrix, Santa Clara, CA, USA) were performed as previously described (Galland et al., 2014a). To obtain presence/absence calls for each probe, we normalized the CEL files by the MAS5 algorithm (Affymetrix). The CEL files were then normalized with the GC-RMA algorithm using the "gcrma" library available from the R Bioconductor suite of open-source softwares (Huber et al., 2015). To determine differentially expressed genes in the embryo and endosperm transcriptomes, we performed a usual two group *t*-test that assumes equal variance between groups. The raw *P*-values were adjusted by the Bonferroni method. We considered a gene as differentially expressed if adjusted-value is < 0.01. To establish the Pearson correlation, we plotted the embryo against the endosperm normalized mean probe intensity. All raw CEL files are available from the Gene Expression Omnibus under the accession GSE43780 (for the embryo: GSM1071216, GSM1071217, GSM1071204; for the endosperm: GSM1071199, GSM1071201, GSM1071210).

## Proteome Analysis

### Protein Extraction and In-Gel Digestion

For embryo protein extraction, three replicates of 50 embryos were ground in liquid nitrogen using mortar and pestle. Then, total soluble proteins were extracted at room temperature in 400 µl thiourea/urea lysis buffer (7 M urea, 2 M thiourea, 6 mM Tris-HCl, 4.2 mM Trizma® base (Sigma-Aldrich, Lyon, France), 4% (w/v) CHAPS) supplemented with 50 µl of the protease inhibitor cocktail Complete Mini (Roche Diagnostics, France, Meylan, France). Then, 15 µl of dithiothreitol (DTT, 1 M, Sigma-Aldrich), 2 µl of DNase I (Roche Diagnostics) and 5 µl of RNase A (Sigma-Aldrich) were added to the sample. For endosperm protein extraction, three replicates of 5 endosperms were ground in liquid nitrogen using mortar and pestle. Then, total soluble proteins were extracted at room temperature in

1 ml thiourea/urea lysis buffer (same composition as above) supplemented with 35 µl of DTT, 2 µl DNaseI and 10 µl RNAse A. Finally, protein extracts were let to agitate for 2 h at 4°C. All samples were then centrifuged at 20,000 g at 4°C for 15 min. The resulting supernatant was submitted to a second clarifying centrifugation as above. The final supernatant was kept and protein concentrations in the various extracts were measured according to Bradford (1976) using Bovine Serum Albumin as a standard.

Twenty-five microgram of embryo and endosperm soluble protein extracts (n = 3 biological replicates) were subjected to SDS-PAGE analysis with 10% acrylamide (Figure S1). Each lane was systematically cut in 16 bands and submitted to in-gel digestion with the Progest system (Genomic Solution) according to a standard trypsin protocol. Gel pieces were washed twice by successive separate baths of 10% acetic acid, 40% ethanol, and acetonitrile. They were then washed twice with successive baths of 25 mM NH₄CO₃ and ACN. Digestion was subsequently performed for 6 h at 37°C with 125 ng of modified trypsin (Promega) dissolved in 20% methanol and 20 mM NH₄CO₃. The peptides were extracted successively with 2% trifluoroacetic acid (TFA) and 50% ACN and then with ACN. Peptide extracts were dried in a vacuum centrifuge and suspended in 20 µl of 0.05% TFA, 0.05% HCOOH, and 2% ACN.

## LC-MS/MS Analysis

Peptide separation by NanoLC was performed as described previously (Bonhomme et al., 2012). Eluted peptides were analyzed on-line with a Q-Exactive mass spectrometer (Thermo Electron) using a nano-electrospray interface. Peptide ions were analyzed using Xcalibur 2.1 with the following data-dependent acquisition parameters: a full MS scan covering 300–1,400 range of mass-to-charge ratio (m/z) with a resolution of 70,000 and a MS/MS step (normalized collision energy: 30%; resolution: 17,500). MS/MS Step was reiterated for the 8 major ions detected during full MS scan. Dynamic exclusion was set to 45 s. A database search was performed with X!Tandem (Craig and Beavis, 2004). Enzymatic cleavage was declared as a trypsin with two possible misscleavage. Cys carboxyamidomethylation was set to static modifications. Met oxydation was set as possible modifications. Precursor mass and fragment mass tolerance were 10 ppm and 0.02 Th, respectively. The 7th annotation of the Rice Genome Annotation Project database (Kawahara et al., 2013; 66,338 proteins) and a contaminant database (trypsin, keratins) were used. Only peptides with an *E*-value smaller than 0.1 were reported.

Peptide quantification was performed by extracted ion current (XIC) using MassChroQ software (Valot et al., 2011). A 5 ppm precision windows was set for XIC extraction. We eliminated the peptide ions not specific of a single protein and, since a peptide ion was detected several times in one biological sample, we summed the Total Ionic Current (TIC) area under peak corresponding to the same peptide ion. We also removed peptide ions that were not reliably detectable by keeping only peptide ions detected at least twice out of the three biological replicates. We obtained a final number of 34,179 and 11,824 peptide ions in the embryo and endosperm respectively corresponded to 2,099
and 786 non-redundant proteins. Since several peptide ions corresponded to the same protein, we summed the total peptide ions' TIC area to get the overall protein abundance and we then log₂-transformed this protein abundance.

## Gene Singular Enrichment Analysis

Gene Singular Enrichment Analysis were performed using the Gene Ontology analysis toolkit provided by the AgriGO web resource (Du et al., 2010) with the Affymetrix (transcriptome) or the corresponding tissue proteome (2099 embryo or 786 endosperm proteins) as backgrounds. The *p*-values generated by a classical hypergeometric overrepresentation test were adjusted by the Yekutieli False Discovery Rate.

## Phylogenetic Analysis

The amino acid sequence of the 12 known glutelins and the putative new glutelin, namely Glu-X, were retrieved from the Rice Genome Annotation Project (Ouyang et al., 2007). Protein sequences were aligned with Clustal Omega (Sievers et al., 2011) with allowed gaps and a distance matrix computed (BLOSUM62 matrix). The corresponding phylogenetic tree was built using the Neighbor-joining method.

## Measurement of Translational Activity

We imbibed three biological replicates of 20 isolated embryos and 5 embryoless endosperms in 4 ml of sterile distilled water with 50 µCi of [³⁵S]-Met (PerkinElmer) at 30°C during 24 h in the dark. Samples were placed on filter papers to remove excess water and

grinded with mortar and pestle using liquid nitrogen. Proteins were then extracted according to previously published protocols (Rajjou et al., 2006). To avoid measuring the non-specific incorporation of radioisotopes into contaminants, we purified the total soluble proteins. In addition, dead seeds (autoclaved seeds) were used as a negative control in order to measure the background level related to non-specific incorporation of [³⁵S]-Met. Finally, 10 µl of protein extracts were added to 5 ml of scintillation liquid cocktail [Ecolite(+), MP Biomedicals, France]. Radioactivity was finally measured (3 biological and 3 technical replicates) using a liquid scintillation analyzer (Tri-Carb 2810TR, PerkinElmer, MA, USA) set between 5 and 100 keV with 10 min integration per sample.

# RESULTS

## Tissue Anatomy of the Rice Mature Seed

The pericarp, testa and AL were removed using a lab-polishing machine to obtain the inner part of the rice seed i.e., the starchy endosperm (Galland et al., 2014a; **Figure 1A**). On a dry weight basis, the embryo represents only 2% of the whole dry seed with the endosperm accounting for the remaining 98% (**Figure 1B**). Within the endosperm, the starchy endosperm represents 87.8% of the endosperm dry weight (**Figure 1B**). To describe the anatomy and content of the rice seed tissues, 5 µm semi-thin sections of dry seeds were obtained and embedded in resin and stained with specific reagents (**Figures 1C-E**). The embryo appeared as a very cell-dense tissue rich in proteins with no

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

**FIGURE 1** | Description of the dry mature rice seed. (A) The dry mature rice seed (Oryza sativa ssp. japonica cv Nipponbare) is composed of several tissues including the embryo, endosperm (pericarp, testa, aleurone layer and starchy endosperm) and inner starchy endosperm. Scale bar, 1 mm. (B) Dry weight per seed of the isolated embryo (E0), starchy endosperm and the aleurone layer/pericarp tissue. The average percentage of each seed tissue is indicated (average % per seed) along with its standard-deviation (n = 10). The endosperm (A0) is the combination of the starchy endosperm (SE) and of the aleurone layer/pericarp tissue. (C-E) Proteins (blue) and complex carbohydrates (including starch, pink) were revealed using a Periodic Acid Schiff-Naphthol Blue Black staining on 5 µm historesin-embedded semi-thick sections. Stained sections of the embryo shoot apical meristem (C, longitudinal cut), embryo radicule (D, transversal cut) and endosperm dorsal side (E) were visualized by optic microscopy. Scale bars represent 100 µm in (C,D) and 25 µm in (E). Al, Aleurone layer; Cp, coleoptile; PI, plumule; P/T, pericarp/testa; Rad, radicle; Sc, scutellum; SE, starchy endosperm.
complex sugars detectable except in cell walls and around the radicle (Figures 1C,D). In contrast, the endosperm appeared as heterogeneous tissue displaying a marked differentiation between the inner starchy endosperm (rich in starch and storage proteins) and the living AL (visible nuclei and numerous protein bodies) (Figure 1E). Lipids were detected by Nile Red tissue staining in both rice endosperm and embryo (Figures 2A-C), with local enrichments in the aleurone/subaleurone layers (Figure 2B) and scutellum epidermis (Figure 2C). These cytological observations show the important degree of compartmentation within the dry mature rice seed.

## Metabolic Composition of the Rice Mature Seed Compartments

Metabolomic data were generated from both embryo and endosperm (Table S1). Thus, 124 unique metabolites were identified at least once either in the rice seed embryo or endosperm and most of them (i.e., 117) were detected in all seed compartments (Figure 3A). Indeed, we only identified six embryo-specific metabolites i.e., γ-tocopherol, feruloylquinic acid, maltotriose, adenosine-5-P and two galactinol isomer (m/z equal to 204 and 433, Table 1). In contrast, ascorbate was detected only in the endosperm (Table 1). A quantitative analysis was performed on the 117 common metabolites detected in both seed compartments. The abundances of the metabolites were normalized according to the dry weight of the embryo and

endosperm and a differential analysis performed that revealed 72 differentially accumulated metabolites ($p < 0.05$, Table 1). We found a strong correlation between the embryo and endosperm metabolite abundance per seed (Figure 3B). Thus, despite the differential ploidy and origin of the mature rice seed tissues, the composition in terms of primary metabolites is rather similar on a per seed basis. Nevertheless, most of these metabolites were more abundant in the endosperm such as unsaturated fatty acids: oleic acid (C18:1), linoleic acid (C18:2) and linolenic acid (C18:3) (Figure 2D). Only a few of them (e.g., raffinose, citrate, α-tocopherol, glucaric acid, digalactosylglycerol) were significantly more abundant in the embryo (Table 1).

## Analysis of Long-Lived Stored mRNAs in the Embryo and Endosperm

The transcriptome of both the embryo and endosperm was analyzed using the Affymetrix Rice Genome Array and a dedicated workflow (Table S2, Figure S2A, Jung et al., 2008). Ambiguous and “absent” probe sets were removed which gave a final number of 15,339 and 16,998 detectable probe sets in the endosperm and embryo, respectively (equivalent to 12,964 and 14,150 unique genes representing 33 and 36% of the total genes) (Table S2C). The existence of a large overlap with 14,227 probe sets commonly detected in both endosperm and embryo was highlighted by these data (Figure 3C). Furthermore, a high significant correlation ($R = 0.948, p < 0.001$) was found

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Precise localization of oil bodies and compartmentation of fatty acids between rice seed tissues. (A–C) Lipids were visualized by confocal microscopy on fresh 100 µm sections of mature rice seeds. Fresh sections were stained by both Nile Red (red channel) to monitor neutral lipids and Calcofluor (green channel) to reveal cell walls. Scale bars are 100 µm in (A), 50 µm in (B), and 10 µm in (C). Sc, scutellum; Endo, endosperm; Al, aleurone layer; Subal, subaleurone layer; ESR, embryo surrounding region. (D) Bar plots represent the mean seed equivalent lipid abundance ($n = 3$). Asterisks indicate significant ($p < 0.05$, False Discovery Rate) different means. A0, endosperm; E0, embryo.</figcaption></figure>
<figure><img src="image_5.png" /></figure>

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

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

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

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

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

**FIGURE 3** | Tissue-specificity and correlations between component tissue molecular apparatus. (A,C,E) Venn diagrams of the number of metabolites (A), probe sets (B), and proteins (C) detected in the dry embryo and endosperm. (B,D,F) Pearson correlations were calculated based on the metabolite level (B), probe signal (D) and protein abundance (F) in the embryo and endosperm. Protein and metabolites are plotted on a log₂ and log₁₀ scale respectively. Transcript (probe sets) intensities are log₂ transformed during normalization and are plotted on non-transformed axis. Pearson correlation coefficients are indicated on each graph along with their significance level (***p < 0.001).

between the normalized probe intensities in the embryo and endosperm transcriptomes (Figure 3D). Yet, some tissue-specific transcripts were detected with 2771 embryo-specific (E0) probe sets (corresponding to 2613 single genes) and 1112 endosperm-specific (A0) probe sets (corresponding to 903 single genes, Figure 3C). Altogether, these results suggest that the endosperm

transcriptome is comparable qualitatively and quantitatively to that of the embryo.

The biological roles for the genes that were strongly (superior to median) and differentially regulated (*p* < 0.01) between the endosperm and the embryo were then analyzed, which yielded to 787 and 1,921 probes (corresponding to 728 and 1746 single
<table><caption>**TABLE 1** | Tissue-specific and differentially accumulated metabolites between embryo and endosperm of the dry mature rice seed.</caption><thead><tr><th>Compound<sup>a</sup></th><th>Class</th><th>Endosperm/embryo log<sub>2</sub> ratio<sup>b</sup></th><th>FDR<sup>c</sup></th></tr></thead><tbody><tr><td>γ-tocopherol</td><td>Apolar</td><td>NA</td><td>Embryo</td></tr><tr><td>Maltotriose</td><td>Carbohydrate</td><td>NA</td><td>Embryo</td></tr><tr><td>Ascorbate</td><td>Organic acids</td><td>NA</td><td>Endosperm</td></tr><tr><td>Feruoylquinic acid</td><td>Organic acids</td><td>NA</td><td>Embryo</td></tr><tr><td>Adenosine-5-P</td><td>Purine_Pyrimidine</td><td>NA</td><td>Embryo</td></tr><tr><td>Galactinol_isomer_1</td><td>Sugar alcohol</td><td>NA</td><td>Embryo</td></tr><tr><td>Galactinol_isomer_2</td><td>Sugar alcohol</td><td>NA</td><td>Embryo</td></tr><tr><td>4-hydroxyproline</td><td>Amino acids</td><td>1.95</td><td>0.05</td></tr><tr><td>Alanine</td><td>Amino acids</td><td>0.88</td><td>0.04</td></tr><tr><td>beta-Alanine</td><td>Amino acids</td><td>1.76</td><td>0.03</td></tr><tr><td>3-cyano-alanine</td><td>Amino acids</td><td>3.31</td><td>0.03</td></tr><tr><td>Cystein</td><td>Amino acids</td><td>1.77</td><td>0.05</td></tr><tr><td>4-Aminobutyrate (GABA)</td><td>Amino acids</td><td>2.29</td><td>0.01</td></tr><tr><td>Glutamine</td><td>Amino acids</td><td>3.00</td><td>0.01</td></tr><tr><td>Isoleucine</td><td>Amino acids</td><td>2.30</td><td>0.04</td></tr><tr><td>Leucine</td><td>Amino acids</td><td>2.21</td><td>0.03</td></tr><tr><td>Lysine</td><td>Amino acids</td><td>1.32</td><td>0.00</td></tr><tr><td>Ornithine</td><td>Amino acids</td><td>0.84</td><td>0.03</td></tr><tr><td>Phenylalanine</td><td>Amino acids</td><td>2.57</td><td>0.03</td></tr><tr><td>Pyroglutamate</td><td>Amino acids</td><td>0.97</td><td>0.03</td></tr><tr><td>Serine</td><td>Amino acids</td><td>2.05</td><td>0.03</td></tr><tr><td>Threonine</td><td>Amino acids</td><td>1.59</td><td>0.04</td></tr><tr><td>Tryptophan</td><td>Amino acids</td><td>1.20</td><td>0.01</td></tr><tr><td>Valine</td><td>Amino acids</td><td>2.39</td><td>0.04</td></tr><tr><td>α-Tocopherol</td><td>Apolar</td><td>-2.55</td><td>0.03</td></tr><tr><td>β-Sitosterol</td><td>Apolar</td><td>2.57</td><td>0.00</td></tr><tr><td>Campesterol</td><td>Apolar</td><td>1.72</td><td>0.00</td></tr><tr><td>Monopalmitin</td><td>Apolar</td><td>3.83</td><td>0.03</td></tr><tr><td>Monostearin</td><td>Apolar</td><td>4.62</td><td>0.02</td></tr><tr><td>Stigmasterol</td><td>Apolar</td><td>3.18</td><td>0.00</td></tr><tr><td>N-acetylmannosamine</td><td>Other</td><td>2.49</td><td>0.02</td></tr><tr><td>Arabinose</td><td>Carbohydrate</td><td>4.38</td><td>0.02</td></tr><tr><td>Cellobiose</td><td>Carbohydrate</td><td>1.95</td><td>0.03</td></tr><tr><td>Galactose</td><td>Carbohydrate</td><td>5.64</td><td>0.02</td></tr><tr><td>Glucopyranose</td><td>Carbohydrate</td><td>4.30</td><td>0.02</td></tr><tr><td>Mannose</td><td>Carbohydrate</td><td>1.54</td><td>0.01</td></tr><tr><td>Melibiose</td><td>Carbohydrate</td><td>1.49</td><td>0.01</td></tr><tr><td>Ribose</td><td>Carbohydrate</td><td>5.07</td><td>0.01</td></tr><tr><td>Sedoheptulose</td><td>Carbohydrate</td><td>3.56</td><td>0.02</td></tr><tr><td>Trehalose</td><td>Carbohydrate</td><td>1.53</td><td>0.00</td></tr><tr><td>Xylose</td><td>Carbohydrate</td><td>4.42</td><td>0.02</td></tr><tr><td>Raffinose</td><td>Carbohydrate</td><td>-0.76</td><td>0.02</td></tr><tr><td>Sucrose</td><td>Carbohydrate</td><td>0.98</td><td>0.01</td></tr><tr><td>Linoleic acid</td><td>Fatty acids</td><td>4.76</td><td>0.01</td></tr><tr><td>Linolenic acid</td><td>Fatty acids</td><td>4.61</td><td>0.02</td></tr><tr><td>Caffeate</td><td>Organic acids</td><td>4.73</td><td>0.03</td></tr><tr><td>Citrate</td><td>Organic acids</td><td>-1.24</td><td>0.01</td></tr><tr><td>Erythronate</td><td>Organic acids</td><td>3.57</td><td>0.05</td></tr><tr><td>Fumarate</td><td>Organic acids</td><td>3.20</td><td>0.05</td></tr><tr><td>4-hydroxybutanoate</td><td>Organic acids</td><td>5.34</td><td>0.01</td></tr></tbody></table>

(Continued)

<table><caption>**TABLE 1** | Continued</caption><thead><tr><th>Compound<sup>a</sup></th><th>Class</th><th>Endosperm/embryo log<sub>2</sub> ratio<sup>b</sup></th><th>FDR<sup>c</sup></th></tr></thead><tbody><tr><td>Glucaric acid</td><td>Organic acids</td><td>-2.69</td><td>0.01</td></tr><tr><td>Gluconate</td><td>Organic acids</td><td>1.96</td><td>0.03</td></tr><tr><td>Glycerate</td><td>Organic acids</td><td>3.80</td><td>0.03</td></tr><tr><td>Glycolate</td><td>Organic acids</td><td>4.57</td><td>0.03</td></tr><tr><td>Malate</td><td>Organic acids</td><td>2.77</td><td>0.01</td></tr><tr><td>Maleate</td><td>Organic acids</td><td>3.90</td><td>0.01</td></tr><tr><td>Malonate</td><td>Organic acids</td><td>2.35</td><td>0.05</td></tr><tr><td>Nicotinate</td><td>Organic acids</td><td>2.11</td><td>0.03</td></tr><tr><td>Phosphate</td><td>Organic acids</td><td>1.15</td><td>0.01</td></tr><tr><td>Plpecolate</td><td>Organic acids</td><td>0.68</td><td>0.01</td></tr><tr><td>Salicylate</td><td>Organic acids</td><td>4.03</td><td>0.05</td></tr><tr><td>Sinapinate</td><td>Organic acids</td><td>3.90</td><td>0.03</td></tr><tr><td>Sinapinate-trans</td><td>Organic acids</td><td>4.01</td><td>0.01</td></tr><tr><td>Succinate</td><td>Organic acids</td><td>3.28</td><td>0.03</td></tr><tr><td>Threonate</td><td>Organic acids</td><td>1.50</td><td>0.00</td></tr><tr><td>Spermidine</td><td>Polyamine</td><td>2.76</td><td>0.01</td></tr><tr><td>Guanine</td><td>Purine_Pyrimidine</td><td>4.31</td><td>0.03</td></tr><tr><td>Uracil</td><td>Purine_Pyrimidine</td><td>3.62</td><td>0.04</td></tr><tr><td>Erythritol</td><td>Sugar alcohol</td><td>3.59</td><td>0.00</td></tr><tr><td>Glycerol</td><td>Sugar alcohol</td><td>4.19</td><td>0.03</td></tr><tr><td>Glycerol-2-P</td><td>Sugar alcohol</td><td>1.56</td><td>0.01</td></tr><tr><td>Glycerol-3-P</td><td>Sugar alcohol</td><td>0.56</td><td>0.01</td></tr><tr><td>Inositol_isomer_3</td><td>Sugar alcohol</td><td>2.15</td><td>0.03</td></tr><tr><td>Inositol_isomer_4</td><td>Sugar alcohol</td><td>1.59</td><td>0.01</td></tr><tr><td>Inositol_isomer_5</td><td>Sugar alcohol</td><td>2.66</td><td>0.03</td></tr><tr><td>Inositol-1-P</td><td>Sugar alcohol</td><td>3.15</td><td>0.03</td></tr><tr><td>Mannitol</td><td>Sugar alcohol</td><td>3.98</td><td>0.03</td></tr><tr><td>Digalactosylglycerol</td><td>Sugar alcohol</td><td>-0.61</td><td>0.00</td></tr><tr><td>Galactosylglycerol</td><td>Sugar alcohol</td><td>4.38</td><td>0.03</td></tr></tbody></table>

<sup>a</sup>See additional file 1 for compound identification details (m/z and retention index)

<sup>b</sup>The average seed metabolite abundance in each tissue (endosperm or embryo) was used to calculate the log<sub>2</sub> ratio. A positive log<sub>2</sub> ratio indicates an endosperm-favored metabolite accumulation (negative for embryo-favored log<sub>2</sub> ratio).

<sup>c</sup>For each metabolite, a Student t-test followed by a False Discovery Rate (FDR) correction was computed. This table contains only metabolites with a *p* < 0.05. Tissue-specific metabolites are indicated (e.g., "Embryo").

genes, respectively) that showed a preferential accumulation in the endosperm and embryo, respectively. Among the GO terms enriched in the 787 endosperm-favored probe sets, the "serine-type endopeptidase inhibitor activity" category, which contains known rice allergenic proteins (RAL2-5) associated with α-amylase or trypsin inhibitory functions, was found (AgriGO, Du et al., 2010, *p* < 0.05, **Figure S3**). Concerning the 1,921 probes with an embryo-favored expression, highly overrepresented GO terms related to "ribosome biogenesis," "translation," "rRNA binding," "ribosomal large and small subunit," and "structural constituent of ribosome" was detected (**Figure S4**). It seems therefore that a large proportion of the embryo long-lived mRNAs will serve as the basis for translation, a conserved and essential process for seed germination (Rajjou et al., 2012).
# Proteome Analysis of the Rice Dry Mature Seed

Regarding seed biology, post-transcriptional and translational regulations add a significant level of complexity as exemplified by studies on developing Arabidopsis and Medicago seeds where numerous examples of delays between mRNA accumulation and protein synthesis have been documented (Gallardo et al., 2007; Hajduch et al., 2010; Verdier et al., 2013). Thus, to complete the description of each rice seed compartment, a proteomic analysis was carried out. The total soluble proteins of embryo and endosperm tissues were extracted and these samples were subjected to a label-free quantitative shotgun proteomic analysis (Figure S1, Figure S2B, S5; Table S3). A total of 2212 single proteins were identified of which only 30.4% (673 proteins) were common to both compartments thereby revealing 1426 embryo-specific and 113 endosperm-specific proteins (Figure 3E). These results showed that the embryo proteome is much more diversified than the endosperm proteome. Furthermore, the abundance of the 673 common proteins is poorly correlated (Figure 3F) contrasting with what was observed at the transcriptome and metabolome levels (Figures 3B,D). For each of the 673 common proteins, the endosperm to embryo protein log₂ ratio was calculated and we found 76 and 267 proteins with an endosperm-favored or embryo-favored abundance respectively (log₂ ratios superior to the median i.e., 1.7 and -2.9 for the endosperm and embryo). By combining these proteins showing a tissue-favored profile with the tissue-specific proteins, we obtained a list of 189 endosperm and 1,693 embryo proteins that we subsequently analyzed using the AgriGO tool (Du et al., 2010). In the endosperm, several expected enriched GO terms were retrieved such as “carbohydrate metabolic process” or “plastid” related to starch biosynthesis (Figure S6). In addition, others interesting enriched GO terms such as “response to endogenous stimulus,” “cell wall,” and “vacuole” were highlighted. These functions are probably related to the developed vesicle trafficking occurring during endosperm development and PCD. In the embryo, a strong GO enrichment for biological processes related to “translation,” “embryonic development,” “post-embryonic development” but also “response to stress” was found (Figure S7). At the subcellular level, several key organelles related terms were also overrepresented such as those related to the “mitochondrion,” “plastid,” or “ribosome” (Figure S7). Since the distribution of the log₂ ratios was similar to a normal distribution (Shapiro-test *p* = 0.95, Figure S5G), a z-score analysis was also performed. This revealed 64 proteins as being differentially accumulated between the embryo and the endosperm (*p* < 0.05; Table S3D). Among these 64 differentially accumulated proteins, 47 are more abundant in the endosperm with several classical proteins classically found enriched in the endosperm such as SSPs (glutelins, prolamins, globulin) and starch biosynthesis enzymes (Table 2). Among the 10 glutelins commonly detected to both seed tissues, all of them were significantly (*p* < 0.05) more abundant in the endosperm (7–13 fold, Figure 4A). An additional protein (Os08g03410, called Glu-X), annotated as a putative glutelin, is preferentially accumulated in the embryo (11 fold, Figure 4A). Glu-X display high sequence homology to other proteins belonging to the glutelin GluA/B/C/D families

(Figures 4B,C). Furthermore, we were able to find one eukaryotic translation factor eIF4A-1 or the rice homolog of the MOTHER of FT and TFL1 as two endosperm-favored proteins not directly related to classical endosperm proteins (Table 2).

# Identification of Post-transcriptional and Translational Regulations in the Mature Seed

The strong seed tissue differentiation observed at the proteome level could be well connected to the delay between transcript and protein accumulation observed at the end of seed development (Gallardo et al., 2007; Hajduch et al., 2010; Arc et al., 2011; Verdier et al., 2013). We quantitatively compared the 673 proteins fold-change between seed tissues with that of their corresponding mRNA level. These 673 proteins matched with 672 non-redundant genes that were used to retrieve their corresponding probe sets among the 14,227 common probe sets (Figure 3C). The probe sets corresponding to the same gene were removed leading to 504 unique probe set/protein pairs (Table S4). Then, the endosperm to embryo log₂ ratios of these 504 probe set/protein pairs were compared (Figure 5). From this analysis, it was obvious that the majority of the 504 protein log₂ ratios were poorly correlated to their corresponding mRNA log₂ ratio. On the one hand, proteins with an endosperm-favored log₂ ratio (≥0) also have similar endosperm-favored mRNA log₂ ratios. On the other hand, the proteins favorably accumulated in the embryo (log₂ ratio < 0) show poorly concordant mRNA profiles. These results once again highlight the seed as a tissue with major strong post-transcriptional and translational regulations probably related to the presence of long-lived mRNAs and to the seed metabolism being brought to a halt (Galland and Rajjou, 2015).

# Measurement of Translational Activities from Isolated Embryos and Embryoless Endosperms

During cereal seed germination, both embryo and endosperm tissues actively transcribe RNA and translate both stored mRNAs and *de novo* synthetized mRNAs (Fincher, 1989). The presently used rice seeds exhibited high vigor (T₅₀ = 16 HAI; Gmax = 24 HAI). The isolated rice embryos were also capable of germination as evidenced by coleoptile emergence upon 24 h of imbibition (Figure 6A). Thus, the translational activities from isolated embryos and embryoless endosperms were assessed. As known for decades, the cereal embryo is an efficient system for *in vitro* translation of mRNAs (Takai and Endo, 2010) while the aleurone layer is actively synthetizing starch-degrading enzymes during germination (Fincher, 1989). Accordingly, considering a single seed, the endosperm had a slightly higher significant translational activity by comparison to the embryo (Figure 6B). By contrast, relatively to total protein content, translational activity is significantly higher in the embryo than in the endosperm (Figure 6C). This result indicates that protein synthesis is more active in the embryo. In support of this statement, the proteomic data revealed that the embryo is better equipped with proteins involved in the translation initiation machinery compared with the endosperm (Figure 6D).
<table><caption>**TABLE 2** | Proteins differentially accumulated between endosperm and embryo.</caption><thead><tr><th>Protein</th><th>Description</th><th>Endosperm to embryo log<sub>2</sub> ratio<sup>a</sup></th><th>p-value<sup>b</sup></th></tr></thead><tbody><tr><td>Os02g06410.1</td><td>Putative SNF-1 related protein kinase</td><td>-8.13</td><td>0.01</td></tr><tr><td>Os12g29400.1</td><td>GRAM domain containing protein, expressed</td><td>-8.02</td><td>0.02</td></tr><tr><td>Os05g39250.1</td><td>PEBP (phosphatidylethanolamine-binding) family protein</td><td>-7.67</td><td>0.02</td></tr><tr><td>Os05g49440.2</td><td>DUF1264 domain containing protein</td><td>-7.62</td><td>0.02</td></tr><tr><td>Os03g15960.1</td><td>17.9 kDa heat shock protein</td><td>-7.49</td><td>0.02</td></tr><tr><td>Os03g14180.1</td><td>26.7 kDa heat shock protein</td><td>-7.47</td><td>0.03</td></tr><tr><td>Os02g32860.1</td><td>Poly [ADP-ribose] polymerase 3</td><td>-7.46</td><td>0.03</td></tr><tr><td>Os03g19290.1</td><td>Chloroplastic outer envelope pore protein</td><td>-7.42</td><td>0.03</td></tr><tr><td>Os03g04410.1</td><td>Aconitase</td><td>-7.41</td><td>0.03</td></tr><tr><td>Os05g37330.1</td><td>60S acidic ribosomal protein</td><td>-7.38</td><td>0.03</td></tr><tr><td>Os02g51750.1</td><td>Annexin</td><td>-7.32</td><td>0.03</td></tr><tr><td>Os10g30150.1</td><td>Universal stress protein</td><td>-7.22</td><td>0.03</td></tr><tr><td>Os03g06360.1</td><td>Late embryogenesis abundant protein D-34</td><td>-7.18</td><td>0.03</td></tr><tr><td>Os01g50910.1</td><td>Late embryogenesis abundant protein LEA_4 domain</td><td>-6.87</td><td>0.04</td></tr><tr><td>Os05g44340.1</td><td>Chaperone protein ClpB1, heat shock protein 101</td><td>-6.86</td><td>0.04</td></tr><tr><td>Os08g41390.1</td><td>Putative 70 kDa peptidylprolyl isomerase</td><td>-6.74</td><td>0.04</td></tr><tr><td>Os10g17280.1</td><td>ATP synthase gamma chain</td><td>-6.71</td><td>0.04</td></tr><tr><td>Os07g42490.1</td><td>Sucrose synthase 3</td><td>8.80</td><td>0.00</td></tr><tr><td>Os07g11630.1</td><td>LTPL163 LTP family protein</td><td>7.94</td><td>0.00</td></tr><tr><td>Os12g14070.1</td><td>70 kDa heat shock-related protein</td><td>7.93</td><td>0.00</td></tr><tr><td>Os05g26350.1</td><td>13 kDa prolamin PROLM4</td><td>6.93</td><td>0.00</td></tr><tr><td>Os05g41970.1</td><td>19 kDa globulin</td><td>6.74</td><td>0.00</td></tr><tr><td>Os11g47520.1</td><td>Xylanase inhibitor protein 2</td><td>6.74</td><td>0.00</td></tr><tr><td>Os06g09450.1</td><td>Sucrose synthase 4</td><td>6.62</td><td>0.00</td></tr><tr><td>Os02g14600.1</td><td>Glutelin GluB-7</td><td>6.53</td><td>0.00</td></tr><tr><td>Os02g32660.1</td><td>Starch branching enzyme</td><td>6.53</td><td>0.00</td></tr><tr><td>Os07g11410.1</td><td>Seed allergenic protein RAG2</td><td>6.19</td><td>0.00</td></tr><tr><td>Os01g44220.1</td><td>ADP-glucose pyrophosphorylase</td><td>6.08</td><td>0.00</td></tr><tr><td>Os06g51084.1</td><td>Starch branching enzyme 1</td><td>6.07</td><td>0.00</td></tr><tr><td>Os07g10570.1</td><td>13 kDa prolamin PROLM25</td><td>5.98</td><td>0.00</td></tr><tr><td>Os02g15090.1</td><td>Glutelin GluD-1</td><td>5.84</td><td>0.00</td></tr><tr><td>Os06g04200.1</td><td>Granule-bound starch synthase 1 (Waxy)</td><td>5.81</td><td>0.00</td></tr><tr><td>Os05g33570.1</td><td>PPDK1 Pyruvate phosphate dikinase 1</td><td>5.79</td><td>0.00</td></tr><tr><td>Os07g11650.1</td><td>LTPL164 LTP family protein</td><td>5.46</td><td>0.01</td></tr><tr><td>Os07g11380.1</td><td>Seed allergenic protein RAG2</td><td>5.33</td><td>0.01</td></tr><tr><td>Os07g10580.1</td><td>13 kDa prolamin PROLM26</td><td>5.29</td><td>0.01</td></tr><tr><td>Os07g11360.1</td><td>Seed allergenic protein</td><td>5.27</td><td>0.01</td></tr><tr><td>Os06g22060.1</td><td>6-phosphofructokinase pyrophosphate dependent</td><td>5.11</td><td>0.01</td></tr><tr><td>Os02g25640.1</td><td>Glutelin GluC-1</td><td>5.03</td><td>0.01</td></tr></tbody></table>

(Continued)

<table><caption>**TABLE 2** | Continued</caption><thead><tr><th>Protein</th><th>Description</th><th>Endosperm to embryo log<sub>2</sub> ratio<sup>a</sup></th><th>p-value<sup>b</sup></th></tr></thead><tbody><tr><td>Os02g15169.1</td><td>Glutelin GluB-1b</td><td>4.88</td><td>0.01</td></tr><tr><td>Os07g11510.1</td><td>Seed allergenic protein</td><td>4.82</td><td>0.01</td></tr><tr><td>Os07g34520.1</td><td>Isocitrate lyase</td><td>4.79</td><td>0.01</td></tr><tr><td>Os02g16820.1</td><td>Glutelin GluB-5</td><td>4.77</td><td>0.01</td></tr><tr><td>Os07g11330.1</td><td>Seed allergenic protein</td><td>4.68</td><td>0.01</td></tr><tr><td>Os04g40660.1</td><td>MA3 domain containing protein</td><td>4.59</td><td>0.01</td></tr><tr><td>Os06g31070.1</td><td>13 kDa prolamin PROLM24</td><td>4.57</td><td>0.01</td></tr><tr><td>Os03g55090.1</td><td>Starch phosphorylase</td><td>4.54</td><td>0.01</td></tr><tr><td>Os03g16440.1</td><td>Translocon at outer-envelope membrane of chloroplast</td><td>4.51</td><td>0.01</td></tr><tr><td>Os03g31360.1</td><td>Glutelin GluA-3</td><td>4.34</td><td>0.02</td></tr><tr><td>Os04g33150.1</td><td>desiccation-related protein</td><td>4.14</td><td>0.02</td></tr><tr><td>Os10g35010.1</td><td>Translocon at inner-envelope membrane of chloroplast</td><td>4.00</td><td>0.02</td></tr><tr><td>Os02g15150.1</td><td>Glutelin GluB-2</td><td>3.79</td><td>0.03</td></tr><tr><td>Os10g26060.1</td><td>Glutelin GluA-2</td><td>3.69</td><td>0.03</td></tr><tr><td>Os06g48750.1</td><td>elf4A-1 Eukaryotic translation initiation factor 4A-1</td><td>3.56</td><td>0.03</td></tr><tr><td>Os02g50350.1</td><td>Dihydropyrimidine dehydrogenase</td><td>3.55</td><td>0.03</td></tr><tr><td>Os01g55690.1</td><td>Glutelin GluA-1</td><td>3.51</td><td>0.03</td></tr><tr><td>Os11g40530.1</td><td>LTPL162 LTP family protein</td><td>3.44</td><td>0.03</td></tr><tr><td>Os07g11900.1</td><td>13 kDa prolamin PROLM19</td><td>3.30</td><td>0.04</td></tr><tr><td>Os01g47410.1</td><td>Aspartic proteinase oryzasin-1</td><td>3.24</td><td>0.04</td></tr><tr><td>Os06g46284.1</td><td>Probable alpha-glucosidase</td><td>3.15</td><td>0.04</td></tr><tr><td>Os06g30370.1</td><td>OsMFT1 homologous to Mother of FT and TFL1</td><td>3.08</td><td>0.05</td></tr><tr><td>Os07g11310.1</td><td>LTPL166 LTP family protein</td><td>3.06</td><td>0.05</td></tr><tr><td>Os05g02060.1</td><td>Chloroplastic outer envelope pore protein</td><td>3.03</td><td>0.05</td></tr><tr><td>Os02g32030.1</td><td>Elongation factor</td><td>3.00</td><td>0.05</td></tr></tbody></table>

## DISCUSSION

### Origin and Roles of the Endosperm, a Key Tissue in Seed Biology

The supernumerary embryo origin of the endosperm is backed by several results in our study and that of others. As previously observed in the case of developing (Belmonte et al., 2013) or germinating Arabidopsis seeds (Penfield et al., 2006), the transcriptomes of the rice embryo and endosperm are highly similar both quantitatively and quantitatively (Figures 3C,D). In developing maize kernels (8 days after pollination), tissue-specific transcriptomics showed that the embryo transcriptome resembled that of aleurone cell layers but also that of other endosperm regions as Spearman correlations ranged from 0.73 (central starchy endosperm) to 0.80 (aleurone layer) (Zhan et al., 2015). In this latter study, the observed total number of maize endosperm-specific expressed genes (3140) is in the same order

<footnote><sup>a</sup>The average seed protein abundance in each tissue (endosperm or embryo) was used to calculate the log<sub>2</sub> ratio. A positive log<sub>2</sub> ratio indicates an endosperm-favored metabolite accumulation (negative for embryo-favored log<sub>2</sub> ratio).</footnote>

<footnote><sup>b</sup>P-value obtained from z-score analysis of endosperm to embryo log<sub>2</sub> ratios.</footnote>
<figure><img src="image_11.png" /></figure>

**FIGURE 4** | Identification of a putative new glutelin with a preferential embryo accumulation. (A) Boxplots of protein abundances (log₂ transformed) for the 10 glutelins that are commonly detected in the endosperm (A0, dark red) and the embryo (E0, blue). The Glu-X protein corresponds to a newly discovered glutelin which is the only glutelin more abundant in the embryo although not significantly. Significantly different means are indicated according to the z-score analysis with *p*-values inferior to 0.01 (**) or 0.05 (†). (B) Close-up view of the conserved cleavage site of glutelin precursor into acidic and basic subunits. The GluC-1 protein sequence does not exhibit the characteristic NG motif. (C) Phylogenetic tree built from the 12 glutelin protein sequences together with the newly putative Glutelin, namely Glu-X (neighbor joining, BLOSUM62 matrix)

of magnitude as ours (1112; **Figure 3C**). Similar conclusions can be drawn in for the embryo (2235 vs. 2771). Embryo and endosperm transcriptomes in developing and mature cereal caryopses could therefore be quite comparable to each other.

## Stored Reserves, Building Blocks, and Energy Source for Germination

### Starch

In rice, starch can account for nearly 70–85% of the total seed weight and is made in different proportions of linear amylose and ramified amylopectin (Bao et al., 2008). Historically, the degradation of starch during germination has been considered as a key mechanism that produced the oligosaccharides required for energy production during germination (Fincher, 1989). Here, the AMY3E α-amylase protein (Os08g36900) was among the most abundant embryo-specific proteins (Table S3C) while a probe signal for the AMY3E mRNA was detectable in both embryo and endosperm (Table S2B). It is worth noting that no amylase activity could be measured in dry rice seeds (Guglielminetti et al., 1995) suggesting that the preformed embryo AMY3E enzyme is not functional at this stage. In the endosperm, only one α-amylase protein was detected (Os01g51754) mildly ortholog of the *Vigna mungo* AMY1.1 enzyme (Table S3C). Finally, as reported before, no β-amylase protein was detected in rice seed compartments. Still, many β-amylase transcripts were found also suggesting that post-transcriptional and/or translational regulation is likely to occur for these enzymes.

Along with the presence of few starch-degrading enzymes, a number of starch biosynthesis enzymes were identified (**Figure 7A**). Starch biosynthesis originates from glucose-1-phosphate that is converted to ADP-glucose by ADP-glucose pyrophosphorylases (AGPases). ADP-glucose serves for α-glucan chain elongation by starch syntheses with participation of starch branching enzymes and starch debranching enzymes (Jeon et al., 2010). Most of the studies on rice starch biosynthesis have focused on the endosperm. Accordingly, a complete enzymatic machinery for starch biosynthesis and ADP-glucose transport (OsBT1-1) was detected in the endosperm (**Figure 7A**). While the presence of these enzymes is not surprising, the present work highlighted a complete starch synthesis enzymatic set already present in the quiescent embryo i.e., AGPS2a, SSIIIb, GBSSII and Pho2 indicating that plastid is the main site of embryonic ADP-glucose synthesis (**Figure 7B**). In contrast, as a characteristic feature of graminaceous, ADP-glucose synthesis can occurs in the endosperm thanks to a cytosolic ADP-glucose pyrophosphorylase (Beckles et al., 2001). Our results reveal a molecular specificity of each rice seed compartment for starch biosynthesis possibly related to starch synthesis in the embryo during germination (Han et al., 2014).

## Seed Storage Proteins (SSPs)

Glutelins are the major SSPs accumulated in rice corresponding to approximately 60–80% of the total proteins in the endosperm. During germination, these proteins are mobilized by proteases
<figure><img src="image_12.png" /></figure>

hereby releasing free amino acids that can readily be incorporated into new proteins requested for germination. Based on amino acid sequence similarity, previous studies established that 12 genes classified into GluA, GluB, GluC, and GluD families encode for rice glutelins (Kawakatsu et al., 2008). The present proteomic data evidenced 11 glutelin forms belonging to the GluA/B/C/D classification (Table S3C). Since the endosperm was separated from the embryo, it is remarkable to note that only two glutelin isoforms belonging to GluB class (i.e., Os02g15150.2 and Os02g15070.1) were specific of this storage tissue (Table S3C). In agreement with a storage function of the endosperm, among the 10 glutelins common to both seed tissues, all glutelins were significantly (p < 0.05) more abundant in the endosperm (7–13 fold, **Figure 4A**). Unexpectedly, one protein (Os08g03410), annotated as a putative glutelin, is strongly accumulated in the embryo (11 fold, **Figure 4A**). This protein, hereby named Glu-X, is not classified in the glutelin GluA/B/C/D families but nevertheless presents the characteristic asparagine/glycine (NG) cleavage site specifically recognized by the vacuole-processing enzyme OsVPE (Wang et al., 2009) to process glutelin precursors into the corresponding acidic and basic subunits (**Figure 4B**; Kumamaru et al., 2010). This Glu-X protein is not closely related to the other glutelin families and could have a distinctive function in the embryo (**Figure 4C**). For instance, in Arabidopsis, cruciferins seed storage proteins protect the embryo from oxidative stress during seed aging (Rajjou et al., 2008; Nguyen et al., 2015).

## Lipids

In most cereal seeds, lipids generally account for only 2–3% of the dry weight (Barthole et al., 2012). In seeds, lipids are present in the form of TAGs that are stored in oil bodies corresponding to small vesicles composed of an inner TAG core surrounded by lipid monolayer containing dedicated proteins such as oleosins,caleosins and steroleosins (Murphy, 1993; Murphy et al., 2001; Jolivet et al., 2009). Lipids were detected in both rice endosperm and embryo with local enrichments in the aleurone/subaleurone layers (Figure 2B) and scutellum epidermis (Figure 2C). This precise localization is very similar to that of barley (Neuberger et al., 2008). Yet, the composition of unsaturated fatty acids such as oleic acid (C18:1), linoleic acid (C18:2) and linolenic acid (C18:3) were more abundant in the endosperm (Figure 2D).

The nutritional and organoleptic quality of the rice seed is highly dependent on polyunsaturated fatty acid peroxidation caused by lipoxygenases (LOXs). The Aldo-Keto Reductases (AKRs) protein family detoxifies a wide variety of lipid peroxidation compounds. Correspondingly, overexpression of Aldo-ketoreductase-1 from *Pseudomonas* strain (PsAKR1) in rice improved seed viability and germination vigor (Narayana et al., 2017). Two AKRs were specifically present in the embryo (Os04g26910 and Os05g38230) while one AKR was common to both embryo and endosperm (Os01g43090). Plant AKRs were recently proposed as potential breeding targets for developing stress tolerant varieties (Sengupta et al., 2015). In our proteomic data, we detected two rice lipoxygenases (OsLOX2, Os03g52860; OsLOX3, Os03g49350) specifically in the dry embryo in accordance with ancient biochemical results (Table S3C; Ida et al., 1983). Functional analysis showed that OsLOX2 and OsLOX3 negatively affect the germination performance of seeds submitted to artificial or natural aging (Huang et al., 2014; Xu et al., 2015). The higher occurrence of lipid hydroperoxides in aged seeds has been linked with a decrease in seed longevity (Sattler et al., 2004). In addition, the overexpression or silencing of *OsLOX2* accelerates or slows germination *sensu stricto* (Huang et al., 2014). Rice *LOX2* gene expression is induced upon germination presumably to degrade TAGs present in oil bodies and fuel seedling establishment. Therefore, TAGs that are stored in the AL oil bodies would release free fatty acids, fuel carbohydrate synthesis and energy metabolism in the AL and ETC cells during rice seed germination. It has been shown that TAG degradation occurs very early in the embryo and AL cells during the germination process (Clarke et al., 1983; Leonova et al., 2010). Following that hypothesis, the present work reveals in an unexpected way that several mRNAs encoding for glyoxylate cycle enzymes i.e., glyoxysomal malate dehydrogenase (MDH, Os12g43630 and Os05g50940), malate synthase (MLS, Os04g40990) and isocitrate lyase (ICL, Os07g34520) were more abundant in the rice endosperm (Table S2). Notably, ICL was also found as endosperm-favored at the protein level (Table 2). These glyoxylate cycle enzymes preferentially found in the rice endosperm might be associated with anoxia and stressful conditions (Lu et al., 2005). This would also explain the endosperm-favored accumulated of the pyruvate phosphate dikinase 1 (PPDK1, Table 2). In contrast, at
<figure><img src="image_13.png" /></figure>

**FIGURE 6** | Translational activity and machinery in rice seed tissues. (A–C) Translational activities in isolated embryos and embryoless endosperms during germination sensu stricto. (A) Germination of isolated rice embryos (top panel) or complete seeds (lower panel) as evidenced by coleoptile emergence (arrows). Scale bars represent 2 mm. (B,C) Radiolabeled [³⁵S]-Methionine incorporation during 24 h of imbibition in 20 isolated embryos (E0) or five embryoless endosperm (A0). (B) [³⁵S]-Methionine incorporation per seed in isolated embryo (E0) and embryoless endosperm (A0). (C) [³⁵S]-Methionine incorporation per microgram of total protein in E0 and A0. Results are the mean (± SD) of three biological replicates and are expressed on a seed equivalent basis for comparison. Signal integration was performed during 10 min. Student's t-tests were applied to identify statistically significant differences (* means statistically significant as *P* < 0.05 and ** means statistically highly significant as *P* < 0.01). (D) Embryo and endosperm proteins involved in translation initiation identified in the dry mature rice seed. For each protein family, the number of proteins found in the rice seed proteome is indicated. Embryo-specific proteins (E0) are indicated along with their locus number. eIF, eukaryotic translation initiation factor; Met, methionine; PABP, polyadenylate-binding protein.
<figure><img src="image_14.png" /></figure>

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

**FIGURE 7** | Strong tissue-specificity of starch biosynthesis enzymes between endosperm and embryo. Specific endosperm (A) and embryo (B) proteins related to starch biosynthesis are displayed on their corresponding enzymatic reactions. AGPL, glucose-1-phosphate adenylyltransferase large subunit; AGPS, glucose-1-phosphate adenylyltransferase small subunit; BT1-1, brittle 1; PHO, Starch phosphorylase; SS, soluble starch synthase; BE, branching enzyme; ISA, isoamylase; GBSS, granule-bound starch synthase; PEP, phosphoenolpyruvate, PPT, phosphoenolpyruvate/phosphate translocator; NTT, ATP carrier protein.

the protein level, the rice embryo appears favorably equipped with enzymes involved in glycolysis, tricarboxylic acid cycle and ATP synthesis (Figure 8). The degradation of membrane phospholipids during seed storage is also detrimental to seed quality (Devaiah et al., 2007). In particular, phospholipase D (PLD) enzymes that cleave membrane phospholipids to phosphatidic acid (PA) are proposed to be one the earliest event of deterioration. In Arabidopsis, silencing of the most abundant PLD enzyme, phospholipase D alpha 1 (AtPLDα1), improves seed longevity (Devaiah et al., 2007). Among the 17 phospholipase D (PLD) rice genes (Li et al., 2007), the PLDα1 (OsPLDα1, Os01g07760) is the only one expressed and the corresponding enzyme is among the most abundant embryo-specific proteins (Table S3C). Among all rice PLD proteins, OsPLDα1 is closest relative of AtPLDα1 suggesting that their negative roles on seed longevity are probably conserved. In conclusion, the genetic manipulation of TAGs and phospholipid-related enzymes, in particular OsLOX2 and OsPLDα1, present in the rice embryo, have the potential to improve rice seed storability and organoleptic value.

## A Tissue-Specific Equipment Possibly Involved in Desiccation Tolerance and Seed Storability

### Glassy State Establishment

Historically, rice seed longevity has been strongly linked to desiccation tolerance (Ellis and Hong, 1994). Among the mechanisms involved, the accumulation of non-reducing sugars (sucrose, trehalose) and raffinose family oligosaccharides (RFO) at the end of seed development converts the cellular cytoplasm into a “glassy state” that restricts molecule mobility and halt enzymatic reactions (Buitink and Leprince, 2004; Rajjou and Debeaujon, 2008; Farrant and Moore, 2011; Hand et al., 2011). RFO are galactosyl-sucrose carbohydrates that are formed by the sequential addition of galactose moieties by galactinol synthase, raffinose synthase, and stachyose synthase. In the present metabolomic data, sucrose, fructose, glucose, and raffinose were the major simple carbohydrates detected in the dry mature rice seed (Table 1 and Table S1). Furthermore, raffinose is preferentially accumulated in the embryo while non-RFO carbohydrates such as sucrose and trehalose are more highly accumulated in the endosperm (Table 1). Interestingly, a colocalization between QTL of longevity and QTL controlling oligosaccharide contents (sucrose, raffinose, and stachyose) has been pointed out in Arabidopsis and rice (Bentsink et al., 2000; Zhu et al., 2007). Together with raffinose, trehalose can also protect proteins and membranes damages induced by desiccation (Fernandez et al., 2010). It was recently found that the trehalose-6-phosphate phosphatase 7 (Os09g20390, OsTPP7) allele from an Indica cultivar was likely to be the underlying QTL for enhanced seed longevity in two Nipponbare near-isogenic lines (Sasaki et al., 2015). The rice genome harbors at least nine trehalose-6-phosphate synthase (OsTPS) and nine OsTPP genes (Fernandez et al., 2010). One natural candidate for favored trehalose synthesis in the endosperm aleurone layer could be OsTPP10 (Os07g30160) whose transcript is only reliably detected in the endosperm (Table S2E) whereas OsTPP8 protein (Os05g50940) was only detected in the embryo (Table S3C). Thus, the present work provides novel knowledge on the spatial regulation of genes involved in trehalose accumulation in rice seed and possibly related to desiccation tolerance and seed longevity. Thus, in the rice endosperm and aleurone layer in particular, the glassy state seems to be dependent on trehalose, sucrose and raffinose while in the embryo, it depends mostly on raffinose (Table 1).

## Protein Folding Protection by Molecular Chaperones

From the mapping of preferentially accumulated proteins, seed categories related to heat stress, protein folding and LEA proteins were quite noticeable (Figure 8). These categories contain protein chaperone roles such as the LEA proteins, (Wang et al., 2007; Hincha and Thalhammer, 2012), the small HSP (Sarkar et al., 2009; Waters, 2013), annexins (Clark et al., 2012), lipocalins (Grzyb et al., 2006), and ClpB chaperones. Most of these proteins were described to be involved in the maintenance of protein folding, prevent membrane aggregation and can also have a synergic effect with non-reducing sugars
<figure><img src="image_16.png" /><figcaption>**FIGURE 8** | MapMan overview of differentially accumulated proteins with a preferential abundance in the endosperm or embryo. The log₂ ratios of the 335 proteins with a log₂ ratio superior to the median were mapped to a seed specific visualization (Joosen et al., 2011). A total of 241 proteins are visible. Red and blue colors represent endosperm and embryo-favored protein abundance respectively.</figcaption></figure>

and RFO to promote glassy state establishment (Boudet et al., 2006; Rajjou and Debeaujon, 2008; Hand et al., 2011). Commonly related with desiccation and abiotic stress tolerance, LEA proteins are members of intrinsically disordered proteins in aqueous solution. They undergo desiccation-induced folding during cell drying suggesting that these proteins could carry out distinct functions under different water states. Rice comprises 34 LEA proteins encoding genes (Wang et al., 2007). Our proteomic analysis identified 12 LEA proteins detected in both embryo and endosperm and 9 only detected in the embryo (Table 3). LEA proteins were previously associated to seed longevity (Chatelain et al., 2012). It was remarkable to note that several LEA proteins, including the dehydrin family were detected exclusively in the embryo, and could be involved in dry storage survival (Table 3). Indeed, a previous study showed that dehydrin RAB18 was very abundant in Arabidopsis dry mature seeds. The abundance of this protein progressively disappeared in aged seeds (Rajjou et al., 2008). Furthermore, it has been demonstrated that down-regulation of seed-specific dehydrins reduced Arabidopsis seed survival in the dry state (Hundertmark et al., 2011). Out

of the 23 predicted sHSP proteins (Sarkar et al., 2009), the HSP17.4 was exclusively found in the endosperm (Table 3). This contrasts with the eight sHSP proteins exclusively found in the embryo of which the OsHSP18.2 (Os01g08860) is capable of protecting the Arabidopsis embryo during artificial aging (Kaur et al., 2015). Furthermore, all three common sHSP proteins (HSP16.9, HSP17.9 and HSP26.7) are preferentially more accumulated in the embryo (Table 3). Remarkably, the HSP16.9 protein was shown to stabilize rice soluble proteins from heat denaturation under in vitro conditions (Yeh et al., 1995). Altogether, these results support the finding that such proteins would primarily serve to protect the embryo against desiccation injuries during late maturation program. The present proteomic analysis also revealed several other categories of chaperone proteins such as annexins, lipocalins, and Clp (caseinolytic protease) chaperones which have never been characterized in cereal dry seeds. First, annexins are probably essential for seed longevity since the overexpression of a sacred lotus (Nelumbo nucifera) isofom in Arabidopsis proved to enhance seed viability under heat stress (Chu et al., 2012). In the present proteomic
<table><caption>TABLE 3 | Proteins involved in folding and chaperone functions found in the embryo and endosperm.</caption><thead><tr><th>Protein</th><th>Family</th><th>Description</th><th>Specificity</th><th>Rank<sup>a</sup></th><th>Endosperm to embryo log<sub>2</sub> ratio<sup>b</sup></th></tr></thead><tbody><tr><td>Os08g36150.1</td><td>ASA1</td><td>Heat shock protein 90 co-chaperone</td><td>Common</td><td>-</td><td>-3.9</td></tr><tr><td>Os05g44340.1</td><td>clpA/clpB</td><td>Chaperone protein ClpB1</td><td>Common</td><td>-</td><td>-6.9</td></tr><tr><td>Os10g42439.1</td><td>DnaJ</td><td>DnaJ homolog</td><td>Embryo</td><td>1,312</td><td>-</td></tr><tr><td>Os01g04370.1</td><td>HSP20</td><td>16.9 kDa class I heat shock protein 1</td><td>Common</td><td>-</td><td>-5.4</td></tr><tr><td>Os01g08860.1</td><td>HSP20</td><td>18.0 kDa class II heat shock protein</td><td>Embryo</td><td>74</td><td>-</td></tr><tr><td>Os02g52150.1</td><td>HSP20</td><td>24.1 kDa heat shock protein</td><td>Embryo</td><td>362</td><td>-</td></tr><tr><td>Os02g54140.1</td><td>HSP20</td><td>18.6 kDa class III heat shock protein</td><td>Embryo</td><td>49</td><td>-</td></tr><tr><td>Os03g14180.1</td><td>HSP20</td><td>Chloroplastic 26.7 kDa heat shock protein</td><td>Common</td><td>-</td><td>-7.5</td></tr><tr><td>Os03g15960.1</td><td>HSP20</td><td>17.9 kDa class I heat shock protein</td><td>Common</td><td>-</td><td>-7.5</td></tr><tr><td>Os03g16020.1</td><td>HSP20</td><td>17.4 kDa class I heat shock protein</td><td>Endosperm</td><td>75</td><td>-</td></tr><tr><td>Os03g16030.1</td><td>HSP20</td><td>18.1 kDa class I heat shock protein</td><td>Embryo</td><td>394</td><td>-</td></tr><tr><td>Os03g16040.1</td><td>HSP20</td><td>17.7 kDa class I heat shock protein</td><td>Embryo</td><td>195</td><td>-</td></tr><tr><td>Os04g36750.1</td><td>HSP20</td><td>23.2 kDa heat shock protein</td><td>Embryo</td><td>93</td><td>-</td></tr><tr><td>Os06g11610.1</td><td>HSP20</td><td>Mitochondrial 26.2 kDa heat shock protein</td><td>Embryo</td><td>1,210</td><td>-</td></tr><tr><td>Os11g13980.1</td><td>HSP20</td><td>21.9 kDa heat shock protein</td><td>Embryo</td><td>1,350</td><td>-</td></tr><tr><td>Os01g08560.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-4.5</td></tr><tr><td>Os01g62290.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-3.6</td></tr><tr><td>Os02g48110.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-2.5</td></tr><tr><td>Os02g53420.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-2.1</td></tr><tr><td>Os03g02260.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Embryo</td><td>135</td><td>-</td></tr><tr><td>Os03g11910.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Embryo</td><td>218</td><td>-</td></tr><tr><td>Os03g16860.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Embryo</td><td>280</td><td>-</td></tr><tr><td>Os03g16920.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-3.5</td></tr><tr><td>Os03g60620.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>0.6</td></tr><tr><td>Os05g08840.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-4.0</td></tr><tr><td>Os05g23740.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>0.8</td></tr><tr><td>Os05g38530.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Embryo</td><td>206</td><td>-</td></tr><tr><td>Os11g47760.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>-2.3</td></tr><tr><td>Os12g14070.1</td><td>HSP70</td><td>70 kDa Heat shock protein</td><td>Common</td><td>-</td><td>7.9</td></tr><tr><td>Os04g01740.1</td><td>HSP90</td><td>90 kDa heat shock protein</td><td>Embryo</td><td>92</td><td>-</td></tr><tr><td>Os06g50300.1</td><td>HSP90</td><td>Endoplasmin homolog</td><td>Common</td><td>-</td><td>-2.9</td></tr><tr><td>Os08g38086.3</td><td>HSP90</td><td>90 kDa heat shock protein</td><td>Endosperm</td><td>109</td><td>-</td></tr><tr><td>Os08g39140.1</td><td>HSP90</td><td>Heat shock protein 81-1</td><td>Embryo</td><td>129</td><td>-</td></tr><tr><td>Os09g29840.1</td><td>HSP90</td><td>90 kDa heat shock protein</td><td>Common</td><td>-</td><td>0.5</td></tr><tr><td>Os09g30412.1</td><td>HSP90</td><td>Heat shock protein 81-2</td><td>Embryo</td><td>551</td><td>-</td></tr><tr><td>Os09g30418.1</td><td>HSP90</td><td>Heat shock protein 81-3</td><td>Common</td><td>-</td><td>0.3</td></tr><tr><td>Os12g32986.1</td><td>HSP90</td><td>90 kDa heat shock protein</td><td>Embryo</td><td>179</td><td>-</td></tr><tr><td>Os01g06630.1</td><td>LEA</td><td>Late embryogenesis abundant LEA_5</td><td>Embryo</td><td>21</td><td>-</td></tr><tr><td>Os01g12580.1</td><td>LEA</td><td>Late embryogenesis abundant LEA_2</td><td>Common</td><td>-</td><td>-5.5</td></tr><tr><td>Os01g50700.1</td><td>LEA</td><td>Dehydrin Rab25</td><td>Common</td><td>-</td><td>-4.4</td></tr><tr><td>Os01g50910.1</td><td>LEA</td><td>Late embryogenesis abundant LEA_4</td><td>Common</td><td>-</td><td>-6.9</td></tr><tr><td>Os02g15250.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Common</td><td>-</td><td>-4.5</td></tr><tr><td>Os02g44870.1</td><td>LEA</td><td>Dehydrin</td><td>Embryo</td><td>723</td><td>-</td></tr><tr><td>Os03g06360.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Common</td><td>-</td><td>-7.2</td></tr><tr><td>Os03g20680.1</td><td>LEA</td><td>Late embryogenesis abundant protein 1</td><td>Common</td><td>-</td><td>-4.7</td></tr><tr><td>Os03g53620.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Embryo</td><td>262</td><td>-</td></tr><tr><td>Os03g62620.2</td><td>LEA</td><td>Late embryogenesis abundant LEA_2</td><td>Common</td><td>-</td><td>-3.3</td></tr><tr><td>Os04g52110.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Common</td><td>-</td><td>-5.0</td></tr><tr><td>Os05g28210.1</td><td>LEA</td><td>Embryonic abundant protein 1</td><td>Embryo</td><td>14</td><td>-</td></tr><tr><td>Os05g46480.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Common</td><td>-</td><td>-5.9</td></tr></tbody></table>

(Continued)
<table><caption>**TABLE 3** | Continued</caption><thead><tr><th>Protein</th><th>Family</th><th>Description</th><th>Specificity</th><th>Rank<sup>a</sup></th><th>Endosperm to embryo log<sub>2</sub> ratio<sup>b</sup></th></tr></thead><tbody><tr><td>Os05g50710.1</td><td>LEA</td><td>Late embryogenesis abundant LEA_2</td><td>Embryo</td><td>321</td><td>-</td></tr><tr><td>Os06g23350.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Common</td><td>-</td><td>-4.5</td></tr><tr><td>Os08g23870.1</td><td>LEA</td><td>Late embryogenesis abundant LEA_1</td><td>Embryo</td><td>110</td><td>-</td></tr><tr><td>Os11g26570.1</td><td>LEA</td><td>Dehydrin</td><td>Common</td><td>-</td><td>-2.98</td></tr><tr><td>Os11g26750.1</td><td>LEA</td><td>Dehydrin Rab16D</td><td>Embryo</td><td>1,396</td><td>-</td></tr><tr><td>Os11g26780.1</td><td>LEA</td><td>Dehydrin Rab16B GN = RAB16B</td><td>Embryo</td><td>330</td><td>-</td></tr><tr><td>Os12g43140.1</td><td>LEA</td><td>Late embryogenesis abundant protein</td><td>Embryo</td><td>69</td><td>-</td></tr><tr><td>Os02g43020.1</td><td>TPR</td><td>Heat shock protein stress-induced TPLR_2</td><td>Common</td><td>-</td><td>-2.4</td></tr><tr><td>Os08g34150.1</td><td>Lipocalin</td><td>OsTIL-2 Temperature-induced lipocalin 2</td><td>Embryo</td><td>644</td><td>-</td></tr><tr><td>Os02g39930.1</td><td>Lipocalin</td><td>OsTIL-1 Temperature-induced lipocalin 1</td><td>Embryo</td><td>1,424</td><td>-</td></tr><tr><td>Os03g31300.1</td><td>clpA/clpB</td><td>Chloroplastic chaperone protein ClpB2</td><td>Common</td><td>-</td><td>-0.12</td></tr><tr><td>Os02g08490.1</td><td>clpA/clpB</td><td>Mitochondrial chaperone protein ClpB3</td><td>Embryo</td><td>165</td><td>-</td></tr></tbody></table>

data, three annexins with one exclusively present in the embryo (Os09g23160) and one significantly more accumulated in the embryo (Os02g51750, *p* < 0.05) were detected in the present study (Table 3). The last annexin (Os06g11800) was reported to be up-accumulated during germination suggesting a possible role on the embryo membrane dynamics (Yang et al., 2007). Secondly, lipocalins, a family of proteins that transport small hydrophobic molecules such as steroids, bilins, retinoids, and lipids, are classified in plants as temperature-induced lipocalins (TILs) and chloroplastic lipocalins (CHLs) (Charron et al., 2005). It has been demonstrated that both TILs and CHLs are involved in lipid protection, which is critical for stress adaptation. Two TIL proteins are predicted from the rice genome sequence (Charron et al., 2005) and they were detected in the mature rice embryo while the plastidial form OsCHL was undetectable in this tissue (Table 3). These results on the relative abundance of the OsTILs and OsCHL are consistent with those showed in Arabidopsis since the accumulation of AtCHL protein in the AtTIL KO mutant and *vice versa* suggests a functional overlap between these two lipocalin types (Boca et al., 2013). Interestingly, seed longevity is correlated with the accumulation of these proteins in Arabidopsis (Boca et al., 2013).

## Protein Repair Systems

Several enzymes involved in protein repair, were presently detected specifically in the rice embryo proteome. This was the case for three Methionine Sulfoxide Reductases (MSR) proteins namely MSRB5, MSRA2-1 and MSRA4 (Table S3C; Rouhier et al., 2006). MSRs are involved in the reversal of oxidized Met residues (Met sulfoxide, MetSO) in altered proteins thereby preventing aging-associated diseases in all organisms (Moskovitz, 2005). The MSRA4.1 is a plastidial enzyme potentially involved in oxidative stress resistance and that can repair free and protein-bound MetSO *in vitro* (Guo et al., 2009). MSR repair system in Medicago and Arabidopsis promote seed longevity (Châtelain et al., 2013). Secondly, the protein-L-isoaspartate O-methyltransferase can repair abnormal isoaspartyl occurring in damaged proteins (Thapar et al., 2001). In seeds, PIMT

are actively involved with the maintenance of seed viability in Arabidopsis (Ogé et al., 2008) and rice (Petla et al., 2016). In wheat, PIMT activity is very high in dry mature seeds, increase up to 4 h after imbibition and then decrease during subsequent germination (Mudgett and Clarke, 1994). Amongst the two rice PIMT genes, we found the OsPIMT2 (Os04g40540; Petla et al., 2016) among the most abundant embryo-specific proteins (rank #188). More precisely, this could be ΔOsPIMT2, a truncated yet functional version of OsPIMT2 (Petla et al., 2016). This rice PIMT protein is accumulated during the very late stages of seed development in relation due to the formation of aspartyl residues during desiccation (Petla et al., 2016).

## ROS Homeostasis

The control of Reactive Oxygen Species (ROS, e.g., H₂O₂) homeostasis during both desiccation and early germination is of paramount importance for seed vigor and longevity (Sattler et al., 2004; Bailly et al., 2008). Proteomic and Metabolomic results emphasized several mechanisms that could help the embryo to cope with desiccation-induced oxidative stress.

First, tocopherols and tocotrienols participate to seed longevity by limiting lipid peroxidation (Sattler et al., 2004). In our data, α- and γ-tocopherols were found to be preferentially, if not exclusively, accumulated in the embryo (Table 1). Several proteins involved in tocopherol biosynthesis pathway were specifically found in the embryo such as the 4-hydroxyphenylpyruvate dioxygenase (HPPD, Os02g07160; Table S3C), which is involved in the production of both plastoquinone and tocopherol essential for plant survival (Sano et al., 2016). In the same way, OsVTE1 protein (Os02g17650), that is responsible for γ-tocopherol synthesis, was specifically detected in the rice embryo (Table S3C). Finally, part of the same pathway, we found expression of Arabidopsis *VTE2* (Os06g44840) and *VTE3* (Os12g42090) homologs with a favored embryo gene expression (5 and 2.5 fold respectively, Table S2F).

Along with vitamin E, ascorbate is also a very important antioxidant molecule. Ascorbate (AsA) and dehydroascorbate (DHA) were specifically or favorably detected in the endosperm

<footnote><sup>a</sup>Rank of the protein respectively to the other tissue-specific proteins.</footnote>

<footnote><sup>b</sup>The average seed protein abundance in each tissue (endosperm or embryo) was used to calculate the log<sub>2</sub> ratio. A positive log<sub>2</sub> ratio indicates an endosperm-favored metabolite accumulation (negative for embryo-favored log<sub>2</sub> ratio).</footnote>
(Table 1 and Table S1). Ascorbate and DHA can be degraded to threonate upon non-enzymatic reaction with H₂O₂ or enzymatically. Interestingly, threonate is also present in high amounts in the endosperm compared to the embryo (Table 1). This suggests that a complete AsA to threonate pathway exist in the endosperm. In developing wheat kernels, ascorbate level decrease from mid- to final seed developmental stage and the ascorbate pool becomes progressively oxidized (Paradiso et al., 2012). From our data, it seems that ascorbate *de novo* synthesis could be restricted to the embryo since a putative mannose-1-phosphate guanyltransferase (Os03g11050) and the two GDP-mannose 3,5-epimerase 1 and 2 (GME1, GME2) are specifically detected in the embryo at similar abundances (Table S3C). In contrast, the ascorbate salvage pathway from monohydroascorbate (MDHA) by MDHA reductase (MDHAR) is present in both rice endosperm and embryo (Table S3C). In addition, the conversion of DHA back to AsA is also possible thanks to DHA reductases (OsDHAR1, Os05g02530) present in both rice embryo and endosperm proteomes at similar levels (Table S3C). Ascorbate could interact with ABA metabolism and/or signaling to modulate seed germination ability. Indeed, exogenous application of low concentrations of ascorbate is able to rescue rice seed germination from abscisic acid treatment (Ye et al., 2012).

The present proteome reveals a wide diversity of antioxidant enzymes that are already present in the dry seed with an embryo-favored accumulation (Figure 8). These enzymes include 4 superoxide dismutases (SODCC1, SODCC2, SODA, SODCP, Table S3C), two embryo-specific catalases (CATA and CATB, Table S3C) and 11 embryo-specific peroxidases including several ascorbate peroxidase i.e., the cytosolic OsAPX2, the peroxisomal OsAPX4 and the stromal OsAPX7 (Table S3C). These results seem to argue in favor of a more abundant ROS detoxification enzymatic apparatus in the rice embryo.

## The Mature Seed Is the Crossroad of Post-transcriptional and Translational Regulations Essential for Germination Success

Studies in various species demonstrated that the developmental transition from a maturing to a germinating seed is the place of strong post-transcriptional and translational regulations (Gallardo et al., 2007; Hajduch et al., 2010; Verdier et al., 2013; Galland et al., 2014b; Layat et al., 2014). In this study, we investigated that post-transcriptional and translational regulations occuring in both tissues of the rice seed at the end of its development (Figure 5). Since germination sensu stricto in both monocot and dicot seeds is only dependent on mRNA translation (Rajjou et al., 2004; Sano et al., 2012), we took a closer look at the translational machinery at the tissue-level.

### Stored mRNAs and the Translational Machinery

We were wondering whether these comparable translational activities in the endosperm and embryo relied on different translational machinery sets. Indeed, we could show that, in Arabidopsis germinating seeds, stored mRNAs were differentially

translated (Galland et al., 2014b) making selective mRNA translation a way to distinguish stored and neosynthesized mRNAs. For these reasons, proteins involved in translation and present in the embryo and endosperm were screened. Thus, 292 proteins related to mRNA translation processes (BIN 2.2.1-2.2.4) were identified in the embryo (Table S3F, Figure 6D). Specifically, 109 different ribosomal proteins could be identified. Initiation of translation is controlled by specific cap-dependent initiation factors (Roy and von Arnim, 2013). First, the 43S pre-initiation complex is formed through association between the 40S ribosomal subunit, a charged methionyl-tRNA and the eIF1, eIF2, eIF3 and eIF5 translation initiation factors. A complete set of the 43S pre-initiation complex was retrieved in the proteomic data (Table S3F). Out of the eight eIF3 protein subunits (B-C-D-E-F-H-K-H) monitored, the present data showed that two isoforms of eIF3E (Os07g12110 and Os07g07250) and one isoform of eIF3F (Os05g01450) are restricted to the embryo. In Arabidopsis, the eIF3F is a key regulator of embryo development particularly in actively developing tissues (Xia et al., 2010). In addition, it interacts with eIF3E suggesting that the observed rice eIF3E/3F proteins would play significant roles during embryonic or post-embryonic cellular processes. Upon formation, the multifactor complex (MFC) associates with the 40S ribosome thereby establishing the 43S pre-initiation complex (Figure 6D).

In parallel, the mRNA cap is recognized by the eIF4F complex composed of both a eIF4E and a eIF4G protein family (Roy and von Arnim, 2013). In plants, a very important feature is the presence of isoforms of eIF4E and eIF4G named eIF(iso)4E and eIF(iso)4G. These different isoforms participate to the mRNA translational selectivity (Mayberry et al., 2009; Martinez-Silva et al., 2012). In our data, we found that the eIF4E and eIF4(iso)4E subunits, responsible for recognition of the mRNA cap, were restricted to the rice embryo and present in relatively similar abundances (Table S3F, Figure 6D). The 43S pre-initiation complex subsequently associates to the eIF4 complex and binds to the mRNA. Translation initiation factors eIF4A with ATP-dependent helicase activity unwind mRNA 5'UTR secondary structures. The eIF4A-1 and eIF4A-3 were detected in our proteomic data (Table S3F, Figure 6D). The eIF4A-1 translation initiation factor is present in both endosperm and embryo, but, in contrast to almost all translation initiation factors, this protein is significantly and strongly up-accumulated in the endosperm (Table S3F, Figure 6D). The eIF4A-3 protein was strictly observed in the embryo suggesting that embryo and endosperm use different eIF4A helicases during mRNA translation initiation. The translation initiation factor eIF6-2 interacts with RACK1, a negative regulator of ABA response and positive regulator of GA signaling (Guo et al., 2011; Fennell et al., 2012). Especially relevant in seed biology, it was demonstrated that ABA inhibited RACK1 and eIF6 gene expressions (Guo et al., 2011). In Arabidopsis, three homologs of the mammalian RACK1, namely RACK1A, RACK1B and RACK1C were characterized. In our proteomic description of rice seeds, it was interesting to observe the presence of two RACK1 proteins with OsRACK1A detected in both tissues and OsRACK1B only detected in the embryo along with the embryo-specific eIF6-2 (Table S3C, Figure 6D). Recently, the OsRACK1A gene has been shown to positively
regulate rice seed germination through promotion of ABA catabolism and H₂O₂ synthesis (Zhang et al., 2014). Thus, in addition to RACK1A, RACK1B could also play a major role in the embryo during seed germination. RACK1A and B could also link ABA and GA signaling with mRNA translation. We only detected the eIF6-2 protein in the embryo suggesting that RACK1A/B regulations of mRNA translation only apply to the embryo and not to the endosperm. Together with the absence of the two cap recognition eIF4E proteins in the endosperm, it further confirms that embryo and endosperm may have contrasted qualitative mRNA translation regulations. Very recent evidences showed that eIFiso4G1 translation initiation factor has a role in the fatty acid profile of Arabidopsis developing seeds through the balance of plastic and nucleus-encoded mRNAs involved in fatty acid biosynthesis (Li et al., 2017). Moreover, this role of eIFisoG1 is not be compensated by eIFiso4G2 suggesting a very specific effect. Further work should distinguish the translational machinery at the tissue level and its consequences on seed metabolism.

## CONCLUDING REMARKS

Starch and SSPs were long associated with the endosperm storage function. Thus, it was remarkable to pinpoint, in our proteomic dataset, the presence of glutelins and starch biosynthesis enzymes at non-negligible level also in the embryo (Figures 4, 7). These results refine and expand previous proteomic results on whole developing rice seeds (Koller et al., 2002). Altogether, this would also support the supernumerary embryo hypothesis with both ancestral tissues being equipped with different molecular apparatus before divergence. Classically, the inner SE has been seen as a dead storage tissue since the central parts of the endosperm undergo PCD (Young and Gallie, 2000). Yet, we could show that, as expected, the embryoless endosperm (aleurone and SE) showed an important translational activity however, the embryo display a higher protein synthesis (Figures 6B,C). While the functional consequences of this mRNA translational activity remain to be established, it is clear that the seed endosperm has emerging new roles regarding the control of seed germination and environmental adaptation (Yan et al., 2014; Bassel, 2016). New determinants of agricultural seed quality both in monocots and dicots crops will undoubtedly benefit from tissue-specific combined “multi-omics.” The genetic and tissue heterogeneity of the mature seed is a considerable challenge to seed biologists. In addition, the seed definitely constitute a fascinating plant organ in which post-transcriptional regulations and translational selectivity fine-tune the biological processes that are spatially and temporally regulated within a few hours. A renewed vision of seed biology by integrative systems biology would certainly dig out meaningful new genetic determinants of seed quality.

## AUTHOR CONTRIBUTIONS

MG, DH, IL, and LR designed and performed the experimental work. GCI completed the metabolome analyses. SB and SH

realized the transcriptomic analysis while BV gathered the proteomic data. GC, BC were involved in the preparation of samples. FG, JT, and EA helped with R statistics. HM and BG provided help with cytological observations. MG and LR wrote the manuscript. All authors read and approved the final manuscript.

## FUNDING

This work was supported by the French Ministry of Industry (FUI, NUTRICE agreement # 092906334).

## ACKNOWLEDGMENTS

We are thankful to Loïc Fontaine and Jean-Benoît Morel (UMR BGPI, CIRAD, Montpellier, France) for supply of rice Nipponbare seeds. We acknowledge Dr Francisco Cubillos Riffo (Laboratorio de Microbiología y Biotecnológica Aplicada, Centro de Estudios en Ciencia y Tecnología de los Alimentos, Santiago, Chile) for help with statistical analysis. Thanks to Olivier Langella and Thierry Balliau, from PAPPSO (Plateforme d'Analyse Protéomique de Paris Sud-Ouest), for making the proteomic data available on PROTICdb. The author gratefully acknowledges the support of K. C. Wong Education Foundation, Hong Kong for the visit of DH at IJPB. Thanks to the Saclay Plant Sciences (SPS) LabEx supporting IJPB, IPS2 and GQE-Le Moulon (ANR-10-LABX-0040-SPS).

## SUPPLEMENTARY MATERIAL

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

**Figure S1** | Protein SDS-PAGE gels used for label-free shotgun proteomics. (A,B) Proteins were extracted, quantified by the Bradford assay and 25 µg of embryo (A) or endosperm (B) proteins were subjected to SDS-PAGE (10% acrylamide) and Coomassie staining. (C) Protein quantification in each lane according to the Multi Gauge 3.0 software. Areas quantified are shown on the gels.

**Figure S2** | Transcriptome (A) and proteome (B) analysis workflows applied in this study.

**Figure S3** | Gene Singular Enrichment Analysis of genes preferentially expressed in the endosperm. From the list of the differentially expressed genes, we extracted 787 probe sets with a preferred expression in the endosperm (log₂ ratio > +1.4, *p* < 0.01). The resulting list of 787 probe sets was then submitted to the AgriGO Gene Singular Enrichment Analysis (Hypergeometric test corrected by Yekutieli False Discovery Rate correction, Affymetrix Rice Genome Array set as background, *p* < 0.05) to detect enriched Biological Process (A), Cellular Compartment (B), and Molecular Function (C) GO categories. 722 probes were classified. Light yellow, dark yellow and red boxes indicates GO terms significant enrichment at *p* < 0.05, < 0.01, and < 0.001 respectively.

**Figure S4** | Gene Singular Enrichment Analysis of genes preferentially expressed in the embryo. From the list of the differentially expressed genes, we extracted 1,921 probe sets with a preferred expression in the embryo (log₂ ratio < −1, *p* < 0.01). The resulting list of 1,921 probe sets was then submitted to the AgriGO Gene Singular Enrichment Analysis (Hypergeometric test corrected by Yekutieli False Discovery Rate correction, Affymetrix Rice Genome Array set as background, *p* < 0.05) to detect enriched Biological Process (A), Cellular Compartment (B), and Molecular Function (C) GO categories. 1,712 probes were
classified. Light yellow, dark yellow and red boxes indicates GO terms significant enrichment at *p* < 0.05, < 0.01, and < 0.001 respectively.

**Figure S5** | Embryo and endosperm peptide and protein distributions. (A–D) Distributions of the gene-specific peptide and protein abundances from (A,C) embryo and (B,D) endosperm on a log₂-transformed axis. (E,F) Quantile-quantile comparison of embryo and endosperm (E) peptide and (F) protein abundances (log₂-transformed). (G) Distribution of the log₂ ratios (endosperm vs. embryo) for the 673 common proteins. The comparison with the theoretical normal law (red line with a mean equal to the estimated mean i.e., −1.89 and standard deviation 2.89) is shown.

**Figure S6** | Gene Singular Enrichment Analysis of endosperm-specific or endosperm-favored proteins. The list of the 113 endosperm-specific and 76 endosperm-favored (log₂ ratio > 1.7) was subjected to a Gene Singular Enrichment Analysis tool (hypergeometric test corrected by Yekutieli False Discovery rate, whole transcriptome set as background, *p* < 0.05) to detect enriched endosperm protein Biological Process (A), Cellular Compartment (B), and Molecular Function (C) GO categories. 164 proteins were classified. Light yellow, dark yellow and red boxes indicates GO terms significant enrichment at *p* < 0.05, < 0.01, and < 0.001 respectively.

**Figure S7** | Gene Singular Enrichment Analysis of embryo-specific or embryo-favored proteins. The list of the 1,426 embryo-specific and 267 embryo-favored (log₂ ratio < −2.9) was subjected to a Gene Singular Enrichment Analysis tool (hypergeometric test corrected by Yekutieli False Discovery rate, whole transcriptome set as background, *p* < 0.05) to detect enriched embryo protein Biological Process (A), Cellular Compartment (B), and Molecular Function (C) GO categories. 1,627 proteins were classified. Light yellow, dark yellow and red boxes indicates GO terms significant enrichment at *p* < 0.05, < 0.01, and < 0.001 respectively.

## REFERENCES

Alpert, P. (2005). The limits and frontiers of desiccation-tolerant life. *Integr. Comp. Biol.* 45, 685–695. doi: 10.1093/icb/45.5.685

Arc, E., Galland, M., Cueff, G., Godin, B., Lounifi, I., Job, D., et al. (2011). Reboot the system thanks to protein post-translational modifications and proteome diversity: how quiescent seeds restart their metabolism to prepare seedling establishment. *Proteomics* 11, 1606–1618. doi: 10.1002/pmic.201000641

Avila-Ospina, L., Clément, G., and Masclaux-Daubresse, C. (2017). Metabolite profiling for leaf senescence in barley reveals decreases in amino acids and glycolysis intermediates. *Agronomy* 7:15. doi: 10.3390/agronomy7010015

Bailly, C., El-Maarouf-Bouteau, H., and Corbineau, F. (2008). From intracellular signaling networks to cell death: the dual role of reactive oxygen species in seed physiology. *C. R. Biol.* 331, 806–814. doi: 10.1016/j.crvi.2008.07.022

Bao, J., Jin, L., Xiao, P., Shen, S., Sun, M., and Corke, H. (2008). Starch physicochemical properties and their associations with microsatellite alleles of starch-synthesizing genes in a rice RIL population. *J. Agric. Food Chem.* 56, 1589–1594. doi: 10.1021/jf073128+

Barthole, G., Lepiniec, L., Rogowsky, P. M., and Baud, S. (2012). Controlling lipid accumulation in cereal grains. *Plant Sci.* 185–186, 33–39. doi: 10.1016/j.plantsci.2011.09.002

Bassel, G. W. (2016). To Grow or not to Grow? *Trends Plant Sci.* 21, 498–505. doi: 10.1016/j.2016.02.001

Beckles, D. M., Smith, A. M., and ap Rees, T. (2001). A cytosolic ADP-glucose pyrophosphorylase is a feature of graminaceous endosperms, but not of other starch-storing organs. *Plant Physiol.* 125, 818–827. doi: 10.1104/pp.125.2.818

Belmonte, M. F., Kirkbride, R. C., Stone, S. L., Pelletier, J. M., Bui, A. Q., Yeung, E. C., et al. (2013). Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed. *Proc. Natl. Acad. Sci. U.S.A.* 110, E435–444. doi: 10.1073/pnas.122206110

Bentsink, L., Alonso-Blanco, C., Vreugdenhil, D., Tesnier, K., Groot, S. P., and Koornneef, M. (2000). Genetic analysis of seed-soluble oligosaccharides in relation to seed storability of Arabidopsis. *Plant Physiol.* 124, 1595–1604. doi: 10.1104/pp.124.4.1595

Bethke, P. C., Liboureel, I. G., Aoyama, N., Chung, Y. Y., Still, D. W., and Jones, R. L. (2007). The Arabidopsis aleurone layer responds to nitric oxide, gibberellin, and abscisic acid and is sufficient and necessary for seed dormancy. *Plant Physiol.* 143, 1173–1188. doi: 10.1104/pp.106.093435

Boca, S., Koestler, F., Ksas, B., Chevalier, A., Leymarie, J., Fekete, A., et al. (2013). Arabidopsis lipocalins AtCHL and AtTIL have distinct but overlapping functions essential for lipid protection and seed longevity. *Plant Cell Environ.* 37, 368–381. doi: 10.1111/pce.12159

Bonhomme, L., Valot, B., Tardieu, F., and Zivy, M. (2012). Phosphoproteome dynamics upon changes in plant water status reveal early events associated with rapid growth adjustment in maize leaves. *Mol. Cell. Proteomics* 11, 957–972. doi: 10.1074/mcp.M111.015867

Boudet, J., Buitink, J., Hoekstra, F. A., Rogniaux, H., Larre, C., Satour, P., et al. (2006). Comparative analysis of the heat stable proteome of radicles of *Medicago truncatula* seeds during germination identifies late embryogenesis

abundant proteins associated with desiccation tolerance. *Plant Physiol.* 140, 1418–1436. doi: 10.1104/pp.105.074039

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal. Biochem.* 72, 248–254. doi: 10.1016/0003-2697(76)90527-3

Buitink, J., and Leprince, O. (2004). Glass formation in plant anhydrobiotes: survival in the dry state. *Cryobiology* 48, 215–228. doi: 10.1016/j.cryobiol.2004.02.011

Buitink, J., and Leprince, O. (2008). Intracellular glasses and seed survival in the dry state. *C. R. Biol.* 331, 788–795. doi: 10.1016/j.crvi.2008.08.002

Charron, J. B., Ouellet, F., Pelletier, M., Danyluk, J., Chauve, C., and Sarhan, F. (2005). Identification, expression, and evolutionary analyses of plant lipocalins. *Plant Physiol.* 139, 2017–2028. doi: 10.1104/pp.105.070466

Chatelain, E., Hundertmark, M., Leprince, O., Le Gall, S., Satour, P., Deligny-Penninck, S., et al. (2012). Temporal profiling of the heat-stable proteome during late maturation of Medicago truncatula seeds identifies a restricted subset of late embryogenesis abundant proteins associated with longevity. *Plant Cell Environ.* 35, 1440–1455. doi: 10.1111/j.1365-3040.2012.02501.x

Chétgiorgi, J., Satour, P., Laugier, E., Ly Vu, B., Payet, N., Rey, P., et al. (2013). Evidence for participation of the methionine sulfoxide reductase repair system in plant seed longevity. *Proc. Natl. Acad. Sci. U.S.A.* 110, 3633–3638. doi: 10.1073/pnas.1220589110

Craig, R., and Beavis, R. C. (2004). TANDEM: matching proteins with tandem mass spectra. *Bioinformatics* 20, 1466–1467. doi: 10.1093/bioinformatics/bth092

Chu, P., Chen, H., Zhou, Y., Li, Y., Ding, Y., Jiang, L., et al. (2012). Proteomic and functional analyses of Nelumbo nucifera annexins involved in seed thermotolerance and germination vigor. *Planta* 235, 1271–1288. doi: 10.1007/s00425-011-1573-y

Clark, G. B., Morgan, R. O., Fernandez, M. P., and Roux, S. J. (2012). Evolutionary adaptation of plant annexins has diversified their molecular structures, interactions and functional roles. *New Phytol.* 196, 695–712. doi: 10.1111/j.1469-8137.2012.04308.x

Clarke, N. A., Wilkinson, M. C., and Laidman, D. L. (1983). "Lipids in cereal technology," in Lipids in Cereal Technology, ed P. J. Barnes (London: London Academic Press), 57–92.

De Giorgi, J., Piskurewicz, W., Loubery, S., Utz-Pugin, A., Bailly, C., Mène-Saffrané, L., et al. (2015). An endosperm-associated cuticle is required for Arabidopsis seed viability, dormancy and early control of germination. *PLoS Genet.* 11:e1005708. doi: 10.1371/journal.pgen.1005708

Devaiah, S. P., Pan, X., Hong, Y., Roth, M., Welti, R., and Wang, X. (2007). Enhancing seed quality and viability by suppressing phospholipase D in Arabidopsis. *Plant J.* 50, 950–957. doi: 10.1111/j.1365-313X.2007.03103.x

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

Ehrenshaft, M., and Brambl, R. (1990). Respiration and mitochondrial biogenesis in germinating embryos of maize. *Plant Physiol.* 93, 295–304. doi: 10.1104/pp.93.1.295
Ellis, R. H., and Hong, T. D. (1994). Desiccation tolerance and potential longevity of developing seeds of Rice (*Oryza sativa* L.). *Ann. Bot.* 73, 501–506. doi: 10.1006/anbo.1994.1062

Farrant, J. M., and Moore, J. P. (2011). Programming desiccation-tolerance: from plants to seeds to resumption plants. *Curr. Opin. Plant Biol.* 14, 340–345. doi: 10.1016/j.pbi.2011.03.018

Fennell, H., Olawin, A., Mizanur, R. M., Izumori, K., Chen, J. G., and Ullah, H. (2012). Arabidopsis scaffold protein RACK1A modulates rare sugar D-alloose regulated gibberellin signaling. *Plant Signal. Behav.* 7, 1407–1410. doi: 10.4161/psb.21995

Fernandez, O., Bethencourt, L., Quero, A., Sangwan, R. S., and Clement, C. (2010). Trehalose and plant stress responses: friend or foe? *Trends Plant Sci.* 15, 409–417. doi: 10.1016/j.tplants.2010.04.004

Fiehn, O., Wohlgemuth, G., Scholz, M., Kind, T., Lee, D. Y., Lu, Y., et al. (2008). Quality control for plant metabolomics: reporting MSI-compliant studies. *Plant J.* 53, 691–704. doi: 10.1111/j.1365-313X.2007.03387.x

Fincher, G. B. (1989). Molecular and cellular biology associated with endosperm mobilization in germinating cereal grains. *Annu. Rev. Plant Physiol. Plant Mol. Biol.* 40, 305–346. doi: 10.1146/annurev.pp.40.060189.001513

Folsom, J. J., Begcy, K., Hao, X., Wang, D., and Walia, H. (2014). Rice FIE1 regulates seed size under heat stress by controlling early endosperm development. *Plant Physiol.* 165, 238–248. doi: 10.1104/pp.113.232413

Friedman, W. (1998). The evolution of double fertilization and endosperm: an “historical” perspective. *Sex. Plant Reprod.* 11, 6–16. doi: 10.1007/s004970050114

Galland, M., Boutet-Mercey, S., Lounifi, I., Godin, B., Balzergue, S., Grandjean, O., et al. (2014a). Compartmentation and dynamics of flavone metabolism in dry and germinated rice seeds. *Plant Cell Physiol.* 55, 1646–1659. doi: 10.1093/pcp/pcu095

Galland, M., Huguet, R., Arc, E., Cueff, G., Job, D., and Rajjou, L. (2014b). Dynamic proteomics emphasizes the importance of selective mRNA translation and protein turnover during Arabidopsis seed germination. *Mol. Cell. Proteomics* 13, 252–268. doi: 10.1074/mcp.M113.032227

Galland, M., and Rajjou, L. (2015). Regulation of mRNA translation controls seed germination and is critical for seedling vigor. *Front. Plant Sci.* 6:284. doi: 10.3389/fpls.2015.00284

Gallardo, K., Firnhaber, C., Zuber, H., Hericher, D., Belghazi, M., Henry, C., et al. (2007). A combined proteome and transcriptome analysis of developing Medicago truncatula seeds: evidence for metabolic specialization of maternal and filial tissues. *Mol. Cell. Proteomics* 6, 2165–2179. doi: 10.1074/mcp.M700171-MCP200

Greenspan, P., Mayer, E. P., and Fowler, S. D. (1985). Nile red: a selective fluorescent stain for intracellular lipid droplets. *J. Cell Biol.* 100, 965–973. doi: 10.1083/jcb.100.3.965

Grzyb, J., Latowski, D., and Strzalka, K. (2006). Lipocalins-a family portrait. *J. Plant Physiol.* 163, 895–915. doi: 10.1016/j.jplph.2005.12.007

Guglielminetti, L., Yamaguchi, J., Perata, P., and Alpi, A. (1995). Amylolytic activities in cereal seeds under aerobic and anaerobic conditions. *Plant Physiol.* 109, 1069–1076. doi: 10.1104/pp.109.3.1069

Guo, J., Wang, S., Valerius, O., Hall, H., Zeng, Q., Li, J. F., et al. (2011). Involvement of Arabidopsis RACK1 in protein translation and its regulation by abscisic acid. *Plant Physiol.* 155, 370–383. doi: 10.1104/pp.110.160663

Guo, X., Wu, Y., Wang, Y., Chen, Y., and Chu, C. (2009). OsMSRA4.1 and OsMSRB1.1, two rice plastidial methionine sulfoxide reductases, are involved in abiotic stress responses. *Planta* 230, 227–238. doi: 10.1007/s00425-009-0934-2

Hajduch, M., Hearne, L. B., Miernyk, J. A., Casteel, J. E., Joshi, T., Agrawal, G. K., et al. (2010). Systems analysis of seed filling in Arabidopsis: using general linear modeling to assess concordance of transcript and protein expression. *Plant Physiol.* 152, 2078–2087. doi: 10.1104/pp.109.152413

Han, C., He, D., Li, M., and Yang, P. (2014). In-depth proteomic analysis of rice embryo reveals its important roles in seed germination. *Plant Cell Physiol.* 55, 1826–1847. doi: 10.1093/pcp/pcu114

Hand, S. C., Menze, M. A., Toner, M., Boswell, L., and Moore, D. (2011). LEA proteins during water stress: not just for plants anymore. *Annu. Rev. Physiol.* 73, 115–134. doi: 10.1146/annurev-physiol-012110-142203

Hay, F. R., and Whitehouse, K. J. (2017). Rethinking the approach to viability monitoring in seed genebanks. *Conserv. Physiol.* 5:cox009. doi: 10.1093/conphys/cox009

Hincha, D. K., and Thalhammer, A. (2012). LEA proteins: IDPs with versatile functions in cellular dehydration tolerance. *Biochem. Soc. Trans.* 40, 1000–1003. doi: 10.1042/BST20120109

Hoekstra, F. A., Golovina, E. A., and Buitink, J. (2001). Mechanisms of plant desiccation tolerance. *Trends Plant Sci.* 6, 431–438. doi: 10.1016/S1360-1385(01)02052-0

Howell, K. A., Millar, A. H., and Whelan, J. (2006). Ordered assembly of mitochondria during rice germination begins with pro-mitochondrial structures rich in components of the protein import apparatus. *Plant Mol. Biol.* 60, 201–223. doi: 10.1007/s11103-005-3688-7

Howell, K. A., Narsai, R., Carroll, A., Ivanova, A., Lohse, M., Usadel, B., et al. (2009). Mapping metabolic and transcript temporal switches during germination in rice highlights specific transcription factors and the role of RNA instability in the germination process. *Plant Physiol.* 149, 961–980. doi: 10.1104/pp.108.129874

Huang, J., Cai, M., Long, Q., Liu, L., Lin, Q., Jiang, L., et al. (2014). OsLOX2, a rice type I lipoxygenase, confers opposite effects on seed germination and longevity. *Transgenic Res.* 23, 643–655. doi: 10.1007/s11248-014-9803-2

Huber, W., Carey, V. J., Gentleman, R., Anders, S., Carlson, M., Carvalho, B. S., et al. (2015). Orchestrating high-throughput genomic analysis with Bioconductor. *Nat. Methods* 12, 115–121. doi: 10.1038/nmeth.3252

Hundertmark, M., Buitink, J., Leprince, O., and Hincha, D. K. (2011). The reduction of seed-specific dehydrins reduces seed longevity in Arabidopsis thaliana. *Seed Sci. Res.* 21, 165–173. doi: 10.1017/S0960258511000079

Ida, S., Masaki, Y., and Morita, Y. (1983). The isolation of multiple forms and product specificity of Rice Lipoxygenase. *Agric. Biol. Chem.* 47, 637–641.

International Rice Genome Sequencing Project (2005). The map-based sequence of the rice genome. *Nature* 436, 793–800. doi: 10.1038/nature03895

Jeon, J. S., Ryoo, N., Hahn, T. R., Walia, H., and Nakamura, Y. (2010). Starch biosynthesis in cereal endosperm. *Plant Physiol. Biochem.* 48, 383–392. doi: 10.1016/j.plaphy.2010.03.006

Jiao, Y., Tausta, S. L., Gandotra, N., Sun, N., Liu, T., Clay, N. K., et al. (2009). A transcriptome atlas of rice cell types uncovers cellular, functional and developmental hierarchies. *Nat. Genet.* 41, 258–263. doi: 10.1038/ng.282

Jolivet, P., Boulard, C., Bellamy, A., Larre, C., Barre, M., Rogniaux, H., et al. (2009). Protein composition of oil bodies from mature Brassica napus seeds. *Proteomics* 9, 3268–3284. doi: 10.1002/pmic.200800449

Joosen, R. V. L., Ligterink, W., Dekkers, B. J., and Hilhorst, H. W. M. (2011). Visualization of molecular processes associated with seed dormancy and germination using MapMan. *Seed Sci. Res.* 21, 143–152. doi: 10.1017/S09602585110000449

Jung, K. H., An, G., and Ronald, P. C. (2008). Towards a better bowl of rice: assigning function to tens of thousands of rice genes. *Nat. Rev. Genet.* 9, 91–101. doi: 10.1038/nrm2343

Kaur, H., Petla, B. P., Kamble, N. U., Singh, A., Rao, V., Salvi, P., et al. (2015). Differentially expressed seed aging responsive heat shock protein OsHSP18.2 implicates in seed vigor, longevity and improves germination and seedling establishment under abiotic stress. *Front. Plant Sci.* 6:713. doi: 10.3389/fpls.2015.00713

Kaur, H., Petla, B. P., and Majee, M. (2016). "Small heat shock proteins: roles in development, dessication tolerance and seed longevity," in *Heat Shock Proteins and Plants*, eds A. A. A. Asea, P. Kaur, and S. Calderwood (New York, NY: Springer International Publishing), 3–18.

Kawahara, Y., de la Bastide, M., Hamilton, J. P., Kanamori, H., McCombie, W. R., Ouyang, S., et al. (2013). Improvement of the *Oryza sativa* Nipponbare reference genome using next generation sequence and optical map data. *Rice* 6:4. doi: 10.1186/1939-8433-6-4

Kawakatsu, T., Yamamoto, M. P., Hirose, S., Yano, M., and Takaiwa, F. (2008). Characterization of a new rice glutelin gene GluD-1 expressed in the starchy endosperm. *J. Exp. Bot.* 59, 4233–4245. doi: 10.1093/jxb/ern265

Kim, S. T., Wang, Y., Kang, S. Y., Kim, S. G., Rakwal, R., Kim, Y. C., et al. (2009). Developing rice embryo proteomics reveals essential role for embryonic proteins in regulation of seed germination. *J. Proteome Res.* 8, 3598–3605. doi: 10.1021/pr900358s
Koller, A., Washburn, M. P., Lange, B. M., Andon, N. L., Deciu, C., Haynes, P. A., et al. (2002). Proteomic survey of metabolic pathways in rice. *Proc. Natl. Acad. Sci. U.S.A* 99, 11969–11974. doi: 10.1073/pnas.172183199

Kumamaru, T., Uemura, Y., Inoue, Y., Takemoto, Y., Siddiqui, S. U., Ogawa, M., et al. (2010). Vacuolar processing enzyme plays an essential role in the crystalline structure of glutelin in rice seed. *Plant Cell Physiol.* 51, 38–46. doi: 10.1093/pcp/pcp165

Lafon-Placette, C., and Kohler, C. (2014). Embryo and endosperm, partners in seed development. *Curr. Opin. Plant Biol.* 17, 64–69. doi: 10.1016/j.pbi.2013.11.008

Lawlor, K. E., and Vince, J. E. (2014). Ambiguities in NLRP3 inflammasome regulation: is there a role for mitochondria? *Biochim. Biophys. Acta* 1840, 1433–1440. doi: 10.1016/j.bbagen.2013.08.014

Layat, E., Leymarie, J., El-Maarouf-Bouteau, H., Caius, J., Langlade, N., and Bailly, C. (2014). Translatome profiling in dormant and nondormant sunflower (*Helianthus annuus*) seeds highlights post-transcriptional regulation of germination. *New Phytol.* 204, 864–872. doi: 10.1111/nph.13002

Lee, J., and Koh, H. J. (2011). A label-free quantitative shotgun proteomics analysis of rice grain development. *Proteome Sci.* 9:61. doi: 10.1186/1477-5956-9-61

Leonova, S., Grimberg, A., Marttila, S., Stymne, S., and Carlsson, A. S. (2010). Mobilization of lipid reserves during germination of oat (*Avena sativa* L.), a cereal rich in endosperm oil. *J. Exp. Bot.* 61, 3089–3099. doi: 10.1093/jxb/erq141

Li, G., Lin, F., and Xue, H. W. (2007). Genome-wide analysis of the phospholipase D family in *Oryza sativa* and functional characterization of PLD beta 1 in seed germination. *Cell Res.* 17, 881–894. doi: 10.1038/cr.2007.77

Li, Q., Shen, W., Zheng, Q., Tan, Y., Gao, J., Shen, J., et al. (2017). Effects of eIFiso4G1 mutation on seed oil biosynthesis. *Plant J.* 90, 966–978. doi: 10.1111/tpj.13522

Logan, D. C., Millar, A. H., Sweetlove, L. J., Hill, S. A., and Leaver, C. J. (2001). Mitochondrial biogenesis during germination in maize embryos. *Plant Physiol.* 125, 662–672. doi: 10.1104/pp.125.2.662

Lu, Y., Wu, Y. R., and Han, B. (2005). Anaerobic induction of isocitrate lyase and malate synthase in submerged rice seedlings indicates the important metabolic role of the glyoxylate cycle. *Acta Biochim. Biophys. Sin.* 37, 406–414. doi: 10.1111/j.1745-7270.2005.00060.x

Martinez-Silva, A. V., Aguirre-Martinez, C., Flores-Tinoco, C. E., Alejandri-Ramirez, N. D., and Dinkova, T. D. (2012). Translation initiation factor AteIF(iso)4E is involved in selective mRNA translation in *Arabidopsis thaliana* seedlings. *PLoS ONE* 7:e31606. doi: 10.1371/journal.pone.0031606

Mayberry, L. K., Allen, M. L., Dennis, M. D., and Browning, K. S. (2009). Evidence for variation in the optimal translation initiation complex: plant eIF4B, eIF4F, and eIF(iso)4F differentially promote translation of mRNAs. *Plant Physiol.* 150, 1844–1854. doi: 10.1104/pp.109.138438

Moskovitz, J. (2005). Methionine sulfoxide reductases: ubiquitous enzymes involved in antioxidant defense, protein regulation, and prevention of aging-associated diseases. *Biochim. Biophys. Acta* 1703, 213–219. doi: 10.1016/j.bbagap.2004.09.003

Mudgett, M. B., and Clarke, S. (1994). Hormonal and environmental responsiveness of a developmentally regulated protein repair L-isoaspartyl methyltransferase in wheat. *J. Biol. Chem.* 269, 25605–25612.

Murphy, D. J. (1993). Structure, function and biogenesis of storage lipid bodies and oleosins in plants. *Prog. Lipid Res.* 32, 247–280. doi: 10.1016/0163-7827(93)90009-L

Murphy, D. J., Hernandez-Pinzon, I., and Patel, K. (2001). Role of lipid bodies and lipid-body proteins in seeds and other tissues. *J. Plant Physiol.* 158, 471–478. doi: 10.1078/0176-1617-00359

Narayana, N. K., Vemanna, R. S., Kodekallu Chandrashekar, B., Rao, H., Vennapusa, A. R., Narasimaha, A., et al. (2017). Aldo-ketoreductase 1 (AKR1) improves seed longevity in tobacco and rice by detoxifying reactive cytotoxic compounds generated during ageing. *Rice* 10:11. doi: 10.1186/s12284-017-0148-3

Neuberger, T., Sreenivasulu, N., Rokitta, M., Rolletschek, H., Gobel, C., Rutten, T., et al. (2008). Quantitative imaging of oil storage in developing crop seeds. *Plant Biotechnol. J.* 6, 31–45. doi: 10.1111/j.1467-7652.2007.00294.x

Nguyen, H. C., Hoefgen, R., and Hesse, H. (2012). Improving the nutritive value of rice seeds: elevation of cysteine and methionine contents in rice plants by ectopic expression of a bacterial serine acetyltransferase. *J. Exp. Bot.* 63, 5991–6001. doi: 10.1093/jxb/ers253

Nguyen, T. P., Cueff, G., Hegedus, D. D., Rajjou, L., and Bentsink, L. (2015). A role for seed storage proteins in Arabidopsis seed longevity. *J. Exp. Bot.* 66, 6399–6413. doi: 10.1093/jxb/erv348

Nowack, M. K., Ungru, A., Bjerkan, K. N., Grini, P. E., and Schnittger, A. (2010). Reproductive cross-talk: seed development in flowering plants. *Biochem. Soc. Trans.* 38, 604–612. doi: 10.1042/BST0380604

Ogé, L., Bourdais, G., Bove, J., Collet, B., Godin, B., Granier, F., et al. (2008). Protein repair L-isoaspartyl methyltransferase 1 is involved in both seed longevity and germination vigor in Arabidopsis. *Plant Cell* 20, 3022–3037. doi: 10.1105/tpc.108.058479

Olsen, O. A. (2004). Nuclear endosperm development in cereals and *Arabidopsis thaliana*. *Plant Cell* 16, 214–227. doi: 10.1105/tpc.017111

Ouyang, S., Zhu, W., Hamilton, J., Lin, H., Campbell, M., Childs, K., et al. (2007). The TIGR Rice Genome Annotation Resource: improvements and new features. *Nucleic Acids Res.* 35(suppl_1):D883–D887. doi: 10.1093/nar/gkl976

Paradiso, A., de Pinto, M. C., Locato, V., and De Gara, L. (2012). Galactone-gamma-lactone-dependent ascorbate biosynthesis alters wheat kernel maturation. *Plant Biol.* 14, 652–658. doi: 10.1111/j.1438-8677.2011.00543.x

Penfield, S., Li, Y., Gilday, A. D., Graham, S., and Graham, I. A. (2006). Arabidopsis ABA INSENSITIVE4 regulates lipid mobilization in the embryo and reveals repression of seed germination by the endosperm. *Plant Cell* 18, 1887–1899. doi: 10.1105/tpc.106.041277

Petla, B. P., Kamble, N. U., Kumar, M., Verma, P., Ghosh, S., Singh, A., et al. (2016). Rice PROTEIN I-ISOASPARTYL METHYLTRANSFERASE isoforms differentially accumulate during seed maturation to restrict deleterious isoAsp and reactive oxygen species accumulation and are implicated in seed vigor and longevity. *New Phytol.* 211, 627–645. doi: 10.1111/nph.13923

Rajjou, L., Belghazi, M., Huguet, R., Robin, C., Moreau, A., Job, C., et al. (2006). Proteomic investigation of the effect of salicylic acid on Arabidopsis seed germination and establishment of early defense mechanisms. *Plant Physiol.* 141, 910–923. doi: 10.1104/pp.106.082057

Rajjou, L., and Debeaujon, I. (2008). Seed longevity: survival and maintenance of high germination ability of dry seeds. *C. R. Biol.* 331, 796–805. doi: 10.1016/j.crvi.2008.07.021

Rajjou, L., Duval, M., Gallardo, K., Catusse, J., Bally, J., Job, C., et al. (2012). Seed germination and vigor. *Annu. Rev. Plant Biol.* 63, 507–533. doi: 10.1146/annurev-arplant-042811-105550

Rajjou, L., Gallardo, K., Debeaujon, I., Vandekerckhove, J., Job, C., and Job, D. (2004). The effect of alpha-amanitin on the Arabidopsis seed proteome highlights the distinct roles of stored and neosynthesized mRNAs during germination. *Plant Physiol.* 134, 1598–1613. doi: 10.1104/pp.103.036293

Rajjou, L., Lovigny, Y., Groot, S. P., Belghazi, M., Job, C., and Job, D. (2008). Proteome-wide characterization of seed aging in the Arabidopsis: a comparison between artificial and natural aging protocols. *Plant Physiol.* 148, 620–641. doi: 10.1104/pp.108.123141

Rolletschek, H., Melkus, G., Grafahrend-Belau, E., Fuchs, J., Heinzel, N., Schreiber, F., et al. (2011). Combined noninvasive imaging and modeling approaches reveal metabolic compartmentation in the barley endosperm. *Plant Cell* 23, 3041–3054. doi: 10.1105/tpc.111.087015

Rouhier, N., Vieira Dos Santos, C., Tarrago, L., and Rey, P. (2006). Plant methionine sulfoxide reductase A and B multigenic families. *Photosyn. Res.* 89, 247–262. doi: 10.1007/s11120-006-9097-1

Roy, B., and von Arnim, A. G. (2013). Translational regulation of cytoplasmic mRNAs. *The Arabidopsis Book* 11:e0165. doi: 10.1199/tab.0165

Sano, N., Masaki, S., Tanabata, T., Yamada, T., Hirasawa, T., Kashiwagi, M., et al. (2013). RNA-binding proteins associated with desiccation during seed development in rice. *Biotechnol. Lett.* 35, 1945–1952. doi: 10.1007/s10529-013-1283-6

Sano, N., Permana, H., Kumada, R., Shinozaki, Y., Tanabata, T., Yamada, T., et al. (2012). Proteomic analysis of embryonic proteins synthesized from long-lived mRNAs during germination of rice seeds. *Plant Cell Physiol.* 53, 687–698. doi: 10.1093/pcp/pcs024

Sano, N., Rajjou, L., North, H. M., Debeaujon, I., Marion-Poll, A., and Seo, M. (2016). Staying alive: molecular aspects of seed longevity. *Plant Cell Physiol.* 57, 660–674. doi: 10.1093/pcp/pcv186

Sarkar, N. K., Kim, Y. K., and Grover, A. (2009). Rice sHsp genes: genomic organization and expression profiling under stress and development. *BMC Genomics* 10:393. doi: 10.1186/1471-2164-10-393
Sasaki, K., Takeuchi, Y., Miura, K., Yamaguchi, T., Ando, T., Ebitani, T., et al. (2015). Fine mapping of a major quantitative trait locus, qLG-9, that controls seed longevity in rice (*Oryza sativa* L.). *Theor. Appl. Genet.* 128, 769-778. doi: 10.1007/s00122-015-2471-7

Sattler, S. E., Gilliland, L. U., Magallanes-Lundback, M., Pollard, M., and DellaPenna, D. (2004). Vitamin E is essential for seed longevity and for preventing lipid peroxidation during germination. *Plant Cell* 16, 1419-1432. doi: 10.1105/tpc.021360

Sengupta, D., Naik, D., and Reddy, A. R. (2015). Plant aldo-keto reductases (AKRs) as multi-tasking soldiers involved in diverse plant metabolic processes and stress defense: a structure-function update. *J. Plant Physiol.* 179, 40-55. doi: 10.1016/j.jplph.2015.03.004

Shin, J. H., Kim, S. R., and An, G. (2009). Rice aldehyde dehydrogenase7 is needed for seed maturation and viability. *Plant Physiol.* 149, 905-915. doi: 10.1104/pp.108.130716

Sievers, F., Wilm, A., Dineen, D., Gibson, T. J., Karplus, K., Li, W., et al. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. *Mol. Syst. Biol.* 7:539. doi: 10.1038/msb.2011.75

Sreenivasulu, N., Usadel, B., Winter, A., Radchuk, V., Scholz, U., Stein, N., et al. (2008). Barley grain maturation and germination: metabolic pathway and regulatory network commonalities and differences highlighted by new MapMan/PageMan profiling tools. *Plant Physiol.* 146, 1738-1758. doi: 10.1104/pp.107.111781

Subbarao, K. V., Datta, R., and Sharma, R. (1998). Amylases synthesis in scutellum and aleurone layer of maize seeds. *Phytochemistry* 49, 657-666. doi: 10.1016/S0031-9422(97)00964-3

Takai, K., and Endo, Y. (2010). The cell-free protein synthesis system from wheat germ. *Methods Mol. Biol.* 607, 23-30. doi: 10.1007/978-1-60327-331-2_3

Tarpley, L., Duran, A. L., Kebrom, T. H., and Sumner, L. W. (2005). Biomarker metabolites capturing the metabolite variance present in a rice plant developmental period. *BMC Plant Biol.* 5:8. doi: 10.1186/1471-2229-5-8

Thapar, N., Kim, A. K., and Clarke, S. (2001). Distinct patterns of expression but similar biochemical properties of protein L-isoaspartyl methyltransferase in higher plants. *Plant Physiol.* 125, 1023-1035. doi: 10.1104/pp.125.2.1023

Valot, B., Langella, O., Nano, E., and Zivy, M. (2011). MassChroQ: a versatile tool for mass spectrometry quantification. *Proteomics* 11, 3572-3577. doi: 10.1002/pmic.201100120

Verdier, J., Lalanne, D., Pelletier, S., Torres-Jerez, I., Righetti, K., Bandyopadhyay, K., et al. (2013). A regulatory network-based approach dissects late maturation processes related to the acquisition of desiccation tolerance and longevity of *Medicago truncatula* seeds. *Plant Physiol.* 163, 757-774. doi: 10.1104/pp.113.222380

Walbot, V., and Evans, M. M. (2003). Unique features of the plant life cycle and their consequences. *Nat. Rev. Genet.* 4, 369-379. doi: 10.1038/nrg1064

Wang, L., Xie, W., Chen, Y., Tang, W., Yang, J., Ye, R., et al. (2010). A dynamic gene expression atlas covering the entire life cycle of rice. *Plant J.* 61, 752-766. doi: 10.1111/j.1365-313X.2009.04100.x

Wang, X.-S., Zhu, H.-B., Jin, G.-L., Liu, H.-L., Wu, W.-R., and Zhu, J. (2007). Genome-scale identification and analysis of LEA genes in rice (*Oryza sativa* L.). *Plant Sci.* 172, 414-420. doi: 10.1016/j.plantsci.2006.10.004

Wang, Y., Zhu, S., Liu, S., Jiang, L., Chen, L., Ren, Y., et al. (2009). The vacuolar processing enzyme OsVPE1 is required for efficient glutelin processing in rice. *Plant J.* 58, 606-617. doi: 10.1111/j.1365-313X.2009.03801.x

Waters, E. R. (2013). The evolution, function, structure, and expression of the plant sHSPs. *J. Exp. Bot.* 64, 391-403. doi: 10.1093/jxb/ers355

Westengen, O. T., Jeppson, S., and Guarino, L. (2013). Global ex-situ crop diversity conservation and the Svalbard Global Seed Vault: assessing the current status. *PLoS ONE* 8:e64146. doi: 10.1371/journal.pone.0064146

Wu, X., Gong, F., Yang, L., Hu, X., Tai, F., and Wang, W. (2014). Proteomic analysis reveals differential accumulation of small heat shock proteins and late embryogenesis abundant proteins between ABA-deficient mutant vp5 seeds and wild-type Vp5 seeds in maize. *Front. Plant Sci.* 5:801. doi: 10.3389/fpls.2014.00801

Wyse Jackson, P., and Kennedy, K. (2009). The global strategy for plant conservation: a challenge and opportunity for the international community. *Trends Plant Sci.* 14, 578-580. doi: 10.1016/j.tplants.2009.08.011

Xia, C., Wang, Y. J., Li, W. Q., Chen, Y. R., Deng, Y., Zhang, X. Q., et al. (2010). The Arabidopsis eukaryotic translation initiation factor 3, subunit F (AtelF3f), is required for pollen germination and embryogenesis. *Plant J.* 63, 189-202. doi: 10.1111/j.1365-313X.2010.04237.x

Xu, H., Wei, Y., Zhu, Y., Lian, L., Xie, H., Cai, Q., et al. (2015). Antisense suppression of LOX3 gene expression in rice endosperm enhances seed longevity. *Plant Biotechnol. J.* 13, 526-539. doi: 10.1111/pbi.12277

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

Xue, L. J., Zhang, J. J., and Xue, H. W. (2012). Genome-wide analysis of the complex transcriptional networks of rice developing seeds. *PLoS ONE* 7:e31081. doi: 10.1371/journal.pone.0031081

Yan, D., Duermeyer, L., Leoveanu, C., and Nambara, E. (2014). The functions of the endosperm during seed germination. *Plant Cell Physiol.* 55, 1521-1533. doi: 10.1093/pcp/pcu089

Yang, P., Li, X., Wang, X., Chen, H., Chen, F., and Shen, S. (2007). Proteomic analysis of rice (*Oryza sativa*) seeds during germination. *Proteomics* 7, 3358-3368. doi: 10.1002/pmic.200700207

Ye, N., Zhu, G., Liu, Y., Zhang, A., Li, Y., Liu, R., et al. (2012). Ascorbic acid and reactive oxygen species are involved in the inhibition of seed germination by abscisic acid in rice seeds. *J. Exp. Bot.* 63, 1809-1822. doi: 10.1093/jxb/err336

Yeh, C. H., Yeh, K. W., Wu, S. H., Chang, P. F., Chen, Y. M., and Lin, C. Y. (1995). A recombinant rice 16.9-kDa heat shock protein can provide thermoprotection in vitro. *Plant Cell Physiol.* 36, 1341-1348.

Young, T. E., and Gallie, D. R. (2000). Programmed cell death during endosperm development. *Plant Mol. Biol.* 44, 283-301. doi: 10.1023/A:1026588408152

Zhan, J., Thakare, D., Ma, C., Lloyd, A., Nixon, N. M., Arakaki, A. M., et al. (2015). RNA sequencing of laser-capture microdissected compartments of the maize kernel identifies regulatory modules associated with endosperm cell differentiation. *Plant Cell* 27, 513-531. doi: 10.1105/tpc.114.135657

Zhang, D., Chen, L., Li, D., Lv, B., Chen, Y., Chen, J., et al. (2014). OsRACK1 is involved in abscisic acid- and H₂O₂-mediated signaling to regulate seed germination in rice (*Oryza sativa* L.). *PLoS ONE* 9:e97120. doi: 10.1371/journal.pone.0097120

Zhang, Y. X., Xu, H. H., Liu, S. J., Li, N., Wang, W. Q., Moller, I. M., et al. (2016). Proteomic analysis reveals different involvement of embryo and endosperm proteins during aging of Yliangyou 2 hybrid rice seeds. *Front. Plant Sci.* 7:1394. doi: 10.3389/fpls.2016.01394

Zhu, C., Cheng, C., Liu, X., and He, J. Y. (2007). Accumulation of soluble sugars related to desiccation tolerance during rice (*Oryza sativa* L.) seed development and maturation. *Seed Sci. Technol.* 35, 649-659. doi: 10.15258/sst.2007.35.3.12

**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.

The reviewer ME and handling Editor declared their shared affiliation.

Received: 24 June 2017; Accepted: 03 November 2017; Published: 22 November 2017

Citation: Galland M, He D, Lounifi I, Arc E, Clément G, Balzergue S, Huguet S, Cueff G, Godin B, Collet B, Granier F, Morin H, Tran J, Valot B and Rajjou L (2017) *An Integrated "Multi-Omics" Comparison of Embryo and Endosperm Tissue-Specific Features and Their Impact on Rice Seed Quality*. *Front. Plant Sci.* 8:1984. doi: 10.3389/fpls.2017.01984

This article was submitted to Plant Systems and Synthetic Biology, a section of the journal *Frontiers in Plant Science*

Copyright © 2017 Galland, He, Lounifi, Arc, Clément, Balzergue, Huguet, Cueff, Godin, Collet, Granier, Morin, Tran, Valot and Rajjou. 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.