# Integrative analysis of the late maturation programme and desiccation tolerance mechanisms in intermediate coffee seeds.

Stéphane Dussert, Julien Serret, Aldecinei Bastos-Siqueira, Fabienne Morcillo, Eveline Déchamp, Valerie Rofidal, Philippe Lasermes, Hervé Etienne, Thierry Joët

► To cite this version:

Stéphane Dussert, Julien Serret, Aldecinei Bastos-Siqueira, Fabienne Morcillo, Eveline Déchamp, et al. Integrative analysis of the late maturation programme and desiccation tolerance mechanisms in intermediate coffee seeds. Journal of Experimental Botany, 2018, 69 (7), pp.1583-1597. 10.1093/jxb/erx492. hal-01720116

HAL Id: hal-01720116

https://hal.science/hal-01720116

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
RESEARCH PAPER

# Integrative analysis of the late maturation programme and desiccation tolerance mechanisms in intermediate coffee seeds

Stéphane Dussert<fnref n="1" />, Julien Serret<fnref n="1" />, Aldecinei Bastos-Siqueira<fnref n="1" />, Fabienne Morcillo<fnref n="2" />, Eveline Déchamp<fnref n="3" />, Valérie Rofidal<fnref n="4" />, Philippe Lashermes<fnref n="1" />, Hervé Etienne<fnref n="3" /> and Thierry Joët<fnref n="1" />,*

Received 30 August 2017; Editorial decision 20 December 2017; Accepted 20 December 2017

Editor: Steve Penfield, John Innes Centre, UK

## Abstract

The ‘intermediate seed’ category was defined in the early 1990s using coffee (Coffea arabica) as a model. In contrast to orthodox seeds, intermediate seeds cannot survive complete drying, which is a major constraint for seed storage and has implications for both biodiversity conservation and agricultural purposes. However, intermediate seeds are considerably more tolerant to drying than recalcitrant seeds, which are highly sensitive to desiccation. To gain insight into the mechanisms governing such differences, changes in desiccation tolerance (DT), hormone contents, and the transcriptome were analysed in developing coffee seeds. Acquisition of DT coincided with a dramatic transcriptional switch characterised by the repression of primary metabolism, photosynthesis, and respiration, and the up-regulation of genes coding for late-embryogenesis abundant (LEA) proteins, heat-shock proteins (HSPs), and antioxidant enzymes. Analysis of the heat-stable proteome in mature coffee seeds confirmed the accumulation of LEA proteins identified at the transcript level. Transcriptome analysis also suggested a major role for ABA and for the transcription factors CaHSFA9, CaDREB2G, CaANAC029, CaPLATZ, and CaDOG-like in DT acquisition. The ability of CaHSFA9 and CaDREB2G to trigger HSP gene transcription was validated by Agrobacterium-mediated transformation of coffee somatic embryos.

**Keywords:** Agrobacterium transformation, Coffea arabica, desiccation tolerance, heat-stable proteome, intermediate seed, late-embryogenesis abundant protein, late maturation, heat-shock factor, heat-shock protein, transcriptome.

## Introduction

The seed is a key structure in the life cycle of higher plants; it facilitates dispersal in space and over time, and consequently ensures the survival of the species. It has been proposed that

the majority of flowering plants produce seeds that are able to withstand almost complete loss of cellular water. Together with certain desiccation-tolerant invertebrate animals, such

<footnote n="1">¹ IRD, Université Montpellier, UMR DIADE, BP 64501, F-34394 Montpellier, France</footnote>

<footnote n="2">² CIRAD, Université Montpellier, UMR DIADE, BP 64501, F-34398 Montpellier, France</footnote>

<footnote n="3">³ CIRAD, Université Montpellier, UMR IPME, BP 64501, F-34394 Montpellier, France</footnote>

<footnote n="4">⁴ Biochimie et physiologie moléculaire des plantes, CNRS, INRA, Montpellier Supagro, Université Montpellier, 2 place Viala, IBIP, Bâtiment 7,34060 Montpellier, France</footnote>

<footnote>* Correspondence: thierry.joet@ird.fr</footnote>
as bdelloid rotifers, tardigrades, and nematodes, resurrection plants, many yeast cells, and fungal spores, most plant seeds are anhydrobiotic organisms (Crowe et al., 1992). In the dry quiescent state, these seeds are endowed with an exceptional capacity to endure extreme conditions after dispersal, such as intense frost/heat events or long drought periods (Buitink and Leprince, 2008). In addition, the longevity of such seeds increases in a predictable way by reducing seed moisture content and temperature (Roberts, 1973), and these so-called ‘orthodox’ seeds can therefore survive *ex situ* storage for very long periods under conventional gene-bank conditions (Li and Pritchard, 2009).

However, about 8% of the world's flowering plants produce seeds that do not tolerate desiccation (Tweddle et al., 2003; Wyse and Dickie, 2017). Two other categories of seeds – intermediate and recalcitrant – have thus been defined with respect to their storability in gene banks (Roberts, 1973; Ellis et al., 1990). These terms are now commonly used to describe the level of sensitivity to drying of non-orthodox seeds. Although the level of desiccation tolerance (DT) may vary considerably within each of these two seed categories (Dussert et al., 1999; Berjak and Pammenter, 2008), it is commonly acknowledged that recalcitrant seeds do not survive if dried to water potentials of about –15 MPa, equivalent to about 90% relative humidity (RH), while intermediate seeds are able to withstand drying to about 30% RH (Black and Pritchard, 2002). Plants that produce non-orthodox seeds prevail in specific habitats, such as tropical rain forests where they represent >45% of tree species (Tweddle et al., 2003; Hamilton et al., 2013). Seeds of many major tropical crops are recalcitrant (e.g. cocoa, coconut, rubber tree) or intermediate (e.g. tea, coffee, oil palm, citrus), which represents a major constraint for growers, seed companies, and germplasm repositories (Li and Pritchard, 2009; Walters et al., 2013).

At least four types of mechanisms, which act synergically and appear to be shared by all anhydrobiotic organisms, mitigate the deleterious effects of desiccation (Crowe et al., 1992; Hoekstra et al., 2001; Leprince and Buitink, 2010; Gaff and Oliver, 2013; Tapia and Koshland, 2014): (i) prevention of oxidative damage through accumulation of antioxidant compounds, such as tocols and glutathione, together with co-ordinated repression of basal metabolism during drying; (ii) avoidance of structural stress through cell wall modification, reorganisation of intracellular membranes and cytoskeleton, and chromatin condensation; (iii) stabilisation of membranes and proteins by non-reducing sugars, late-embryogenesis abundant proteins (LEA) and heat-shock proteins (HSP) through the formation of hydrogen bonds with polar residues of proteins and membrane phospholipids, which maintain dry molecules in a physical state similar to that seen in the presence of water (a mechanism known as the ‘water replacement hypothesis’, reviewed in Crowe, 2007); and (iv) efficient detoxication and repair systems for damaged DNA and proteins upon rehydration. Several recent comparative transcriptionic, proteomic, or metabolomic studies have determined how these mechanisms are orchestrated during dehydration and rehydration in anhydrobiotes, such as resurrection plants, bryophytes, tardigrades, and midges (Oliver et al., 2011; Yobi

et al., 2012; Dinakar and Bartels, 2013; Gusev et al., 2014; Wang et al., 2014).

In orthodox seeds, DT is acquired *in planta* during the late maturation stage of seed development. Again, proteome, transcriptome, and gene co-expression network analyses have greatly improved our understanding of DT acquisition in orthodox seeds (Buitink et al., 2006; Chatelain et al., 2012; Verdier et al., 2013). In maturing orthodox seeds, the down-regulation of genes involved in the cell cycle, DNA processing, and primary metabolism, and the up-regulation of genes coding for small HSPs (sHSPs) and LEAs, and of genes involved in stress defence (antioxidant system, secondary metabolites) coincide with DT acquisition (Buitink et al., 2006). The specific role of certain LEA proteins in DT has been further characterised through proteome analysis of developing *Medicago truncatula* seeds (Chatelain et al., 2012) and of pre-germinated *M. truncatula* and *Arabidopsis thaliana* seeds treated with polyethylene glycol to re-induce DT, which was lost during germination (Boudet et al., 2006; Maia et al., 2011).

The co-ordinated induction of these protective mechanisms involves complex regulatory networks and signalling pathways (Verdier et al., 2013). In *A. thaliana*, seed maturation is controlled by several master regulators, which interact in a complex manner linked with ABA signalling (Karssen et al., 1983; Gutierrez et al., 2007). They include the CCAAT-box binding factor LEC1 and the three B3 domain-containing proteins ABI3, FUS3, and LEC2 (To et al., 2006). Their downstream targets include other transcription factors (TFs), hormonal pathways, and LEA genes (To et al., 2006; Santos-Mendoza et al., 2008). Part of the seed maturation process is therefore indirectly controlled by master regulators via secondary TFs, which activate their own transcriptional programme (Mönke et al., 2012). Such a cascade has been described between ABI3 and Heat Shock Factor A9 (HSFA9), a secondary TF responsible for HSP transcriptional activation (Kotak et al., 2007). Analysis of gene co-expression networks during *M. truncatula* seed development also revealed close associations between ABI3, secondary TFs such as ABI4, ABI5, AP2/EREBP, and HSFs, and many LEA genes (Verdier et al., 2013; Zinsmeister et al., 2016).

In contrast to the considerable progress made in understanding DT acquisition during the development of orthodox seeds, our knowledge of the factors involved in desiccation sensitivity of intermediate and recalcitrant seeds remains very limited. Reactive oxygen species (ROS) are directly implemented in recalcitrant seed death following partial dehydration (Leprince et al., 1999; Bailly, 2004; Berjak and Pammenter, 2008). The lack of co-ordinated repression of metabolism during drying, e.g. unabated respiratory activity and subsequent metabolic dysfunction, is thought to play a predominant role in the oxidative burst that triggers lipid oxidation, membrane disruption, and ultimately the death of recalcitrant seeds. Comparative analyses with orthodox seeds has also revealed quantitative differences in key protective components such as LEAs (Delahaie et al., 2013) and the lack of structural adaptation to prevent mechanical damage (Berjak and Pammenter, 2008).
The basis of intermediate seed sensitivity to desiccation is even less well documented than that of recalcitrant seeds. The features that intermediate seeds would need to be as tolerant to desiccation as orthodox seeds and those that make them considerably less sensitive to drying than recalcitrant seeds are totally unknown. The intermediate seed category has been defined using coffee (*Coffea arabica*) as a model (Ellis et al., 1990). Coffee seeds die when dried below 20% RH (Dussert et al., 2006). The coffee seed is composed of a copious living endosperm surrounding a tiny spatulated embryo (1% of the mature seed mass) (Eira et al., 2006). The embryo is more tolerant to drying than the endosperm and the level of desiccation sensitivity of whole mature seeds actually corresponds to that of the endosperm (Dussert et al., 2006). Species of the genus *Coffea* vary considerably in the level of desiccation sensitivity of their seeds (Dussert et al., 2000), but neither the total lipid and fatty acid content of the endosperm (Dussert et al., 2001) nor the content of non-reducing sugars (Chabrillange et al., 2000) explain the interspecific variation in seed DT. Since these early studies, no progress has been made in identifying the reasons for coffee seed sensitivity to drying, whereas their developmental features and the transcriptional programme for the accumulation of reserves are now better understood (Joët et al., 2009, 2010, 2014). In the present study, changes in transcript abundance and hormone contents of developing coffee seeds were investigated to better understand DT acquisition in intermediate seeds. The roles of two TFs in sHSP accumulation were validated by *Agrobacterium*-mediated transformation of coffee somatic embryos, and changes in LEA transcript accumulation were further characterised by analysing the heat-stable proteome of mature seeds.

## Material and methods

### Plant material, desiccation tolerance, and viability assays

Experiments were performed with seeds of *Coffea arabica* cv. 'Laurina'. Developing seeds were harvested from plants grown at Grand Tampon, Reunion Island (1015 m above sea level) with an average temperature of 18.5 °C and mean annual rainfall of 1300 mm. For transcriptome and hormone analyses, three independent biological samples (pools of ~200 seeds) were collected from 20 trees randomly selected in the plot for each of the seven developmental stages studied (Salmona et al., 2008). After being cross-sectioned, the seed was separated from the pericarp and immediately frozen in liquid nitrogen and stored at -80 °C. For germination capacity and DT assays, one pool of 300 seeds was collected at 10-d intervals from 140 to 260 d after flowering (DAF), as described above. For DT measurement, batches of 50 seeds were desiccated for 20 d at 25 °C in the dark over saturated solutions of KOH (9% RH), K acetate (23% RH), MgCl₂ (32% RH), K₂CO₃ (45% RH), NH₄NO₃ (62% RH), or (NH₄)₂SO₄ (81% RH). Seed moisture content, seed germination and viability, zygotic embryo extraction, *in vitro* culture and viability were then assessed as described in Dussert et al. (2006). For

DT measurement, somatic embryos (SE) were desiccated for 6 h at 25 °C in the dark over saturated solutions of (NH₄)₂SO₄ (81% RH), KCl (85% RH), KNO₃ (92% RH), or K₂SO₄ (97% RH). SE viability was assessed using the criterion of normal autotrophic seedling conversion after 8 weeks of *in vitro* culture at 25 °C in the dark on the maturation medium M described by Etienne (2005).

### Hormone analysis

Seed tissues were freeze-dried and ground to fine powder in an analytical grinder (IKA A10, Staufen, Germany) before being sent to the National Research Council, Canada for hormone analysis (50 mg of dry powder per sample). ABA and ABA metabolites, cytokinins, auxins, and GAs were quantified by ultra-performance liquid chromatography-ESI-tandem mass spectrometry as described in detail by Chiwocha et al. (2003). The values presented are means of triplicate determinations.

### Transcriptome analysis

Gene identifiers are from the Coffee Genome Hub (coffee-genome.org). For microarray analyses, total RNA was extracted from whole seeds as described in Joët et al. (2014). However, since the embryo represents less than 1% of total seed mass, it was assumed that its contribution to the total RNA population was negligible. Microarrays were performed at the MGX transcriptomic platform (Montpellier-GenomiX, Institut de Génomique Fonctionnelle, Montpellier) using the 15k PUCE CAFE and the protocols described in Privat et al. (2011). The microarray data are publicly available at http://www.ncbi.nlm.nih.gov/geo/ (GEO accession number GSE107949). For 454 pyrosequencing, total RNA was extracted from separated endosperm and embryos using the RNeasy™ Lipid Tissue Kit (Qiagen). cDNAs were tagged independently (Titanium Kit, Roche Diagnostics) then mixed for sequencing by GATC Biotech AG (Germany) using GS FLX Titanium (Roche). Raw data have been deposited at the European Nucleotide Archive (ENA) under the project number PRJEB23959. Trimmed reads were mapped on CDS using the BWA-MEM package with default parameters (Li, 2013). Samtools (Li et al., 2009) was used to count mapped reads and the number of reads per kilobase and million reads (rpkm) were then calculated. Statistical analyses of microarray and 454 data were performed as described in Dussert et al. (2013) and Joët et al. (2014). Briefly, differences in transcript accumulation between two consecutive stages were tested using the Limma Package (Wettenhall and Smyth, 2004) and Audic-Claverie statistics (Audic and Claverie, 1997), respectively. In both cases, P-values were Bonferroni-corrected to control for false discovery rates. Hierarchical cluster analysis (HCA) was used to group genes according to their transcription profile using the tool developed by Eisen et al. (1998). Gene ontology (GO) annotation was performed using Blast2GO software (Götz et al., 2008) with default parameters. Fisher's exact test was then used to evaluate the significance of GO term enrichment. Biological enrichment analysis was also performed using Arabidopsis homologues of differentially expressed
coffee genes and the Classification Supervisor tool (http://bar.utoronto.ca; Provart et al., 2003) for Mapman categorisation. The method used for qPCR analysis is described in Joët et al. (2014). The set of primers that enabled amplification of target genes is detailed in Supplementary Table S1 at JXB online. The level of expression of each gene was normalised to the geometric mean of expression levels of three validated reference coffee genes (Cc08_g05690 Ubiquitin UBQ10, Cc00_g15790 40S ribosomal protein S24, and Cc00_g17460 14-3-3 protein; Cruz et al., 2009; Joët et al., 2014).

## Production of transgenic coffee somatic embryos overexpressing CaHSFA9 and CaDREB2 genes

Coffea arabica HSFA9 and DREB2 cDNAs (Genbank accession numbers JQ687374 and JQ687375, respectively) were used as a template for PCR amplification using the primers described in Supplementary Table S1 and were cloned into the plant overexpression vector pMDC32. Coffee embryogenic calli (C. arabica cv. 'Caturra') were genetically transformed using recombinant Agrobacterium tumefaciens strain LB1119 containing the recombinant plasmid, as previously described (Ribas et al., 2011). About 10 somatic embryos regenerated from each hygromycin-resistant callus (independent transgenic lines) and four regenerated plantlets were tested for transformation by PCR amplification of the selection gene, HPTII. Transgenic lines were then tested for both gene expression and desiccation tolerance.

## Heat-stable proteome analysis

Proteins were extracted from 500 mg of ground coffee endosperm obtained from 25 mature seeds. Heat-stable soluble proteins were extracted according to the method described by Boudet et al. (2006) and optimised as described by Chatelain et al. (2012). The only modification made in the present work was the addition of dithiothreitol (DTT, 5 mM) and polyvinylpolypyrrolidone (PVPP, 4% w/v) to the extraction buffer to prevent tissue oxidation. Proteome analysis included protein fractionation and tryptic digestion, as described in Decourcelle et al. (2015), and LC-MS/MS analysis using the mass spectrometer QExactive Plus (ThermoFisher Scientific Inc) (see Supplementary Table S2). Raw files were analysed using Mascot (2.4.0) with a threshold P-value of 0.05 and the predicted peptide database from the Coffee Genome Hub.

## Results

### Acquisition of desiccation tolerance, hormonal changes, and transcriptional switches in the developing and maturing coffee seed

Under the field conditions used in this study, seeds completed their development in approximately 260 d. At 150 DAF, most seeds were already able to germinate and to develop into normal seedlings (Fig. 1A) but they could not withstand drying in any of the four RH tested (Fig. 1B). In contrast, at this stage, 75%

of embryos extracted from seeds dried in 23% RH germinated normally (Fig. 1C), demonstrating that desiccation-induced loss of viability of immature seeds was only due to damage to the endosperm. Between 150 and 190 DAF, almost all seeds acquired the capacity to be dried in 62% RH and half of them also became tolerant to drying in 23% RH (Fig. 1B). The level of DT of mature coffee seeds was thus mostly gained between 150 and 190 DAF, which corresponds to the transition between stages 5 and 6 according to the anatomic and metabolic criteria defined in previous studies, in which seven stages are distinguished during seed development (Salmona et al., 2008; Joët et al., 2009). These seven developmental stages can be briefly described as follows (Fig. 1D): from stage 1 to stage 2, the perisperm undergoes significant growth, which determines the final size of the seed; at stage 3, the endosperm develops rapidly and replaces the perisperm in the locule; endosperm growth ends by stage 4 and oil starts to accumulate; stage 5 is characterised by endosperm hardening due to the massive deposition of galactomannans in cell walls; accumulation of reserves ends by stage 6, as illustrated by the transcription pattern of three genes representative of oil (OLE-1), protein (SSP1), and galactomannan (ManS1) storage (Fig. 1D), when the pericarp of the fruit turns yellow; finally, fruit and seed maturity is completed at stage 7, when the pericarp becomes red.

In addition, DT acquisition was concomitant with a slight decrease in seed water content (Fig. 1A), most likely caused by the end of reserve accumulation (Fig. 1D). It did not correspond to a desiccation stage sensu stricto since seed water content remained high until maturity, i.e. around 50% on a fresh weight basis (Fig. 1A), neither was DT acquisition coupled with a peak in ABA (Fig. 1E, Supplementary Table S3). The highest ABA level was observed at the onset of the accumulation of storage compounds (stage 4) and endosperm ABA content then decreased gradually until maturity. The IAA level showed a similar pattern during endosperm development and maturation. By contrast, a spectacular drop in active cytokinins (trans-zeatin and 2iP) occurred concomitantly with the acquisition of DT between stages 5 and 6, while the level of the inactive glycosylated form of trans-zeatin-O-glucoside remained high until maturity (see Supplementary Table S3). No active gibberellins were detected at any time during endosperm development (Supplementary Table S3). Among the 15 522 genes represented on the PUCE CAFE microarray (Privat et al., 2011), 9573 (61.7%) displayed significant variation in expression during seed development (FDR<0.05; Supplementary Table S4). Among these genes, 5857 were found to be significantly up- or down-regulated between two consecutive stages. This transcriptome survey revealed two major transcriptional switches during seed development (Fig. 1F). The first switch occurred at the transition between stages 3 and 4 (with 1396 and 1312 genes up- and down-regulated, respectively) before the maturation stage (accumulation of reserves). The second major transcriptional switch, of even greater magnitude, occurred between stages 5 and 6 and thus coincided with endosperm DT acquisition and the end of the reserve accumulation programme. About 2000 genes were up-regulated during this stage and approximately the same number were down-regulated (Fig. 1F).
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Acquisition of desiccation tolerance (DT) during coffee seed development. (A) Acquisition of germinative capacity and changes in moisture content (on a fresh weight basis) during coffee seed maturation. (B) Changes in viability (%) after equilibration drying at 81%, 62%, 45%, and 23% relative humidity (RH) during seed development. (C) Viability of zygotic embryos at 150 d after flowering (DAF, stage 5) and 200 DAF (stage 6) after equilibration drying at 9%, 23%, 32%, 45% and 62% RH. (D) Schematic representation of the seven developmental stages studied and transcript profile of three genes chosen to illustrate accumulation of storage proteins, triacylglycerols, and cell wall polysaccharides: 11S globulin (SSP1, Cc03_g05570), oleosin (OLE-1, Cc02_g04750), and mannan synthase (ManS1, Cc06_g04240), respectively. (E) Changes in ABA, cytokinins (isopentenyladenine, trans-Zeatin), and auxin (IAA) during seed maturation. The DT acquisition timeframe, i.e. the stage 5–6 transition, is indicated by the dashed vertical lines. (F) Numbers of up- and down-regulated genes at each of the six developmental transitions.</figcaption></figure>

Acquisition of desiccation tolerance coincided
with transcriptional activation of cellular protection
mechanisms and co-ordinated repression of metabolic
activity

Hierarchical clustering analysis (HCA) resulted in 26 major
clusters of temporally co-expressed genes (Supplementary
Fig. S1). Of these, nine clusters (C18 to C26) grouped 1521

genes specifically up-regulated during DT acquisition (stage
5-6 transition), while cluster C5 contained 613 genes that
were specifically down-regulated during this developmental
stage (see Supplementary Table S5). All the top 100 up-reg-
ulated genes (Fig. 2) displayed a considerable fold-change in
gene expression during the stage 5-6 transition: i.e. from 5- to
400-fold changes (Supplementary Table S5). About half of
them code for proteins directly involved in cell protection and rescue, including sHSPs, LEAs, defence-, cold-, and drought-induced proteins, ROS-scavenging enzymes, as well as the metabolism of trehalose, a major sugar-sensing and stress-signalling molecule (Fig. 2A). The top 100 up-regulated genes also contained several genes coding for cell wall-modifying enzymes or involved in secondary metabolism, such as genes encoding cytochrome P450 enzymes, and phenylpropanoid and flavonoid biosynthetic enzymes. Finally, transcription regulators represented a substantial proportion of both the top 100 up-regulated genes (10%) and of the full set of genes specifically activated during the stage 5–6 transition (16%; Supplementary Table S6).

The composition of the top 100 down-regulated genes (fold-change between 0.1 and 0.5, Supplementary Table S5) suggests a transcriptional co-ordination of metabolic slow-down during the stage 5–6 transition, with a decrease in the mRNA abundance of genes involved in the fate of proteins (biosynthesis and degradation), primary metabolism, photosynthesis, and respiration (Fig. 2B). The 138 protein-fate genes account for the largest proportion of the transcript repertoire down-regulated during late maturation (32%; Supplementary Table S6). Genes involved in basal cellular activity, such as those encoding components of cytoskeleton

or membrane transporters, were also found to be down-regulated. In addition, functional category enrichment analysis performed on the 424 annotated genes of cluster C5 highlighted global down-regulation of genes involved in energy processes, such as photosynthesis, the tricarboxylic acid cycle, respiratory electron transport, and ATP synthase, and glycolysis (Fig. 2D). Similarly, the most significantly over-represented GO terms in cluster C5 were related to translation, proteasome, photosynthesis, cytoskeleton organisation, and vesicle-mediated transport, while translational activity and organisation of ribosomes were the most significantly under-represented GO terms in up-regulated genes (see Supplementary Table S6).

### Re-routing of energy metabolism and induction of a maturation-specific ROS-scavenging programme

The analysis of up- and down-regulated genes during DT acquisition revealed complex re-routing of energy metabolism, with down-regulation of photosynthesis- and respiration-related transcripts, including key mitorespiration genes, such as those encoding subunits of ATP synthase and electrogenic transporters of the electron transfer chain, and, in counterpart, large transcript accumulation of genes involved

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Functional analysis of genes differentially expressed during acquisition of desiccation tolerance (DT) in developing coffee seeds. Pie-charts of functions assigned to the top 100 up-regulated (A) and top 100 down-regulated (B) genes. Functional enrichment analysis (normalised frequencies of Mapman categories and bootstrap standard deviations) of the 1521 up-regulated (C) and the 677 down-regulated genes (D).</figcaption></figure>
in gluconeogenesis, fermentation, glycolysis-related methylglyoxal detoxication, and the glyoxylate cycle (Fig. 2C, Supplementary Table S6). The very few respiratory or photosynthetic genes up-regulated during DT acquisition encode proteins that act as energy-dissipating systems and are involved in the fine-tuning of the redox poise, such as mitochondrial and plastidial alternative terminal oxidases (Cc08_g05330 and Cc09_g01070, respectively), or non-electrogenic (type II) external NADH dehydrogenases (i.e. Cc01_g14360 and Cc02_g12400), which oxidise extra-mitochondrial NADH accumulated in the case of intense glycolytic activity.

Transcriptome analysis also revealed a major switch in redox homeostasis during the stage 5-6 transition, with the transcriptional activation of many key ROS-scavenging enzymes, such as superoxide dismutases, catalases, glutathione reductases, glutaredoxins, and glutathione peroxidases (Fig. 3). Moreover, many genes encoding the antioxidant and detoxifying enzyme glutathione-S-transferase, which is involved in the reversibility of S-glutathionylation, degradation of xenobiotics, and protein protection against glycation, were also up-regulated during the stage 5-6 transition (Fig. 3). Other genes involved in protection against protein oxidation and glycation, or lipid oxidation, were also highly up-regulated during DT acquisition. These included genes encoding methionine sulfoxide reductase (MSR, Cc02_g04840), glyoxalases, and an aldehyde dehydrogenase (Cc04_g11660) with high homology with Osaldh7, which has been shown to play an important role in seed longevity by detoxifying aldehydes generated by lipid peroxidation (Shin et al., 2009). The induction of a lipocalin gene (Cc06_g15300) homologous to the gene encoding the TIL protein in Arabidopsis, which is involved in seed longevity through lipid protection against oxidation, is also worth noting (Boca et al., 2014). By contrast, an ascorbate peroxidase gene (Cc01_g11800) was down-regulated, suggesting that ascorbate does not play a pivotal role at this stage of seed development.

## Hormone homeostasis and signalling pathways during late maturation

DT acquisition at the stage 5-6 transition coincided with a strong transcriptional activation of several genes primarily involved in the biosynthesis of ABA (zeaxanthin epoxidase, Cc07_g06010, 9-cis-epoxycarotenoid dioxygenase NCED1, Cc05_g10210, xanthoxin dehydrogenase ABA2, Cc01_g08700, and abscisic-aldehyde oxidase, Cc04_g10210), jasmonate (linoleate 9S-lipoxygenase 5, Cc02_g33800), and ethylene (aminocyclopropane-1-carboxylate oxidase ACO1, Cc05_g02890 and Cc09_g01110) (see Supplementary Table S5). However, several genes that control the inactivation and homeostatic regulation of endogenous phytohormone levels were also induced at this stage. Among them, we identified a gene coding for ABA 8'-hydroxylase (Cc08_g05920). The other hormones likely subjected to catabolic or homeostatic processes during DT acquisition were brassinosteroids (cytochrome P450 734A1, Cc04_g03930), gibberellins (gibberellin 2-beta-dioxygenase GA2OX1, Cc05_g10350), jasmonate (amidohydrolase ILL6, Cc10_g12030), cytokinins

(zeatin O-glucosyltransferase, ZOG1, Cc06_g17790), and auxin [indole-3-acetic acid (IAA)-amido synthetase GH3.3, Cc02_g19470; UDP-glycosyltransferase 74E2, Cc01_g01490]. In Arabidopsis, the latter two enzymes modulate auxin levels in the presence of excess auxin and under abiotic stress, respectively, while ZOG1 regulates active versus storage forms of cytokinins, and ILL6 contributes to jasmonoyl-isoleucine turnover and jasmonate homeostasis (Staswick et al., 2005, Tognetti et al., 2010, Widemann et al., 2013). Finally, transcript profiles of the cytokinin biosynthetic enzymes adenylate isopentenyltransferase (IPT1, Cc11_g10030) and cytokinin riboside 5'-monophosphate phosphoribohydrolase (LOG1, Cc06_g03960) mirrored seed cytokinin contents, with high transcript levels at stages 4 and 5 followed by a significant decline in expression at stages 6 and 7.

The repertoire of genes up-regulated during DT acquisition also included proteins involved in intracellular signalling and protein kinase cascades (e.g. MAPK). Among them, we detected several type-2C protein phosphatase genes (PP2Cs; Cc01_g12800, Cc02_g34520, Cc06_g03340), which play important roles in stress signal transduction in plants by negatively regulating the ABA response and MAPK cascade pathways. SRK2I (sucrose non-fermenting 1-related subfamily 2, SnRK2.3, Cc07_g05710), a direct component of the ABA signalling pathway, was also activated during DT acquisition. The induction of the ABA receptor gene PYR1 (Cc08_g02750) and of the cytokinin receptor AHK3 (Cc11_g17030), which is a negative regulator of drought and salt stress responses and of ABA signalling (Tran et al., 2007; Kumar and Verslues, 2015), is also worth noting. Several key players of energy homeostasis were also transcriptionally activated, such as trehalose phosphate synthase genes (TPS; Cc11_g06690), KINB1 (Cc01_g06750), and KING1 (Cc04_g10020, Cc11_g13370), which encode β and γ regulatory subunits of SnRK1, the plant orthologue of the evolutionarily conserved SNF1/AMPK/SnRK1 protein kinase family.

## LEA accumulation in developing coffee seeds

The amino acid sequences of the 51 Arabidopsis LEAs (Hundertmark and Hincha, 2008) were used to conduct BLAST searches against the coffee genome, which led to the identification of 29 LEA proteins (see Supplementary Table S7). Phylogenetic analysis of all coffee and Arabidopsis LEA proteins (Supplementary Fig. S2) showed that coffee LEAs are distributed in the seven distinct groups defined according to PFAM domains (LEA_1 to LEA_5, dehydrins and seed maturation proteins, SMP; Hundertmark and Hincha, 2008). Using publicly available coffee transcriptomes of various vegetative and reproductive tissues (Denoeud et al., 2014), HCA revealed that only five LEA genes are ubiquitously expressed (LEA_2-1, LEA_2-2, LEA_3-1, LEA_Dh3, LEA_Dh4; Fig. 4). Most LEAs displayed marked tissue-specificity, with maximum expression in either the stamen (Cluster II-a, seven LEA genes) or the endosperm (Cluster I-a, nine LEA genes). Since not all LEA genes were represented on the PUCE CAFE microarray, we then conducted a 454 pyrosequencing-based transcriptome analysis of the coffee embryo and endosperm at the stage
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Transcript profiles of desiccation tolerance (DT)-associated reactive oxygen species (ROS) scavenging and detoxification enzymes. Gene expression levels at the seven developmental stages are indicated by the different levels of shading, as indicated in the key, and are presented as percentages of the maximum normalised transcript abundance of the gene. ALDH, aldehyde dehydrogenase; APX, ascorbate peroxidase; AOX, alternative oxidase; CAT, catalase; DHAR, dehydroascorbate reductase; GPX, glutathione peroxidase; GR, glutathione reductase; GRX, glutaredoxin; GST, glutathione-S-transferase; LAC, laccase (diphenol oxidase); LOX, lipoxygenase; MDAR, monodehydroascorbate reductase; MSR, methionine sulfoxide reductase; NDH2, type-2 NAD(P)H dehydrogenase; OXO, oxalate oxidase; PDI, protein disulphide isomerise; PPO, polyphenol oxidase; PTOX, plastid terminal oxidase; SOD, superoxide dismutase; THRX, thioredoxin.</figcaption></figure>

5-6 transition (Fig. 4, Supplementary Table S8). This analysis showed that 19 LEA genes were very highly expressed in the endosperm at stage 6. Significant amounts of these 19 transcripts were already present at stage 5 in the embryo (150 DAF), in which two additional LEA genes were significantly expressed (LEA_5-3, Cc02_g06270 and LEA_4-9, Cc07_g04790). LEA genes expressed in the endosperm at stage 6 belonged to all seven groups of LEAs (i.e. three genes in group LEA_1, two in LEA_2, one in LEA_3, six in LEA_4, two in LEA_5, one SMP, and four dehydrins). Within this subset, 13 LEA genes

were significantly up-regulated during DT acquisition at the stage 5-6 transition according to Audic-Claverie (454 data) and Limma (microarray) statistics (Supplementary Table S7). To ascertain whether the lower DT of intermediate seeds compared to orthodox seeds was caused by incomplete translation of LEA transcripts during late maturation, we analysed the heat-stable soluble proteome of mature coffee seeds (Supplementary Table S7). Proteins were detected for 17 of the 19 transcripts accumulated at stage 6 (Fig. 4). Except for group LEA_3, all the LEA groups were represented in the mature coffee seed proteome.
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Genome-wide analysis of late-embryogenesis abundant (LEA) proteins in coffee. LEA genes were classified according to hierarchical clustering analysis of transcript abundance in *C. canephora* roots, stamens, pistils, leaves, and perisperm at 180 d after flowering (DAF) and in endosperm at 260 DAF. Expression levels in *C. canephora* tissues and in *C. arabica* zygotic embryos and endosperm at 150 and 200 DAF are indicated by the blue dots, the size and intensity of which represent the number of reads per kilobase and million reads (rpkm; see key). LEA peptides detected in the heat-stable proteome of mature seeds are labelled with an orange star. LEA genes differentially expressed during seed development as shown by microarray data analysis are labelled with a diamond (closed diamonds indicate LEA genes differentially expressed during the stage 5–6 transition).</figcaption></figure>

## TFs that putatively govern the late maturation programme in coffee seed

Transcriptome profiling revealed intense transcriptional activity during the establishment of DT and, among the up-regulated genes, 128 displayed putative functions related to

transcriptional regulation (representing 16% of genes with an assigned putative function; Supplementary Table S5). The TF showing the highest induction was a NAC domain-containing protein 29 (ANAC029, Cc11_g12740) with a 45-fold change during DT acquisition. In addition, the most
represented TF classes were AP2/ERF, MYB, WRKY, and HSF families (with more than six genes in each class). Some examples of very highly up-regulated TFs during DT acquisition are homologues of DREB2G (Cc09_g03140, 17-fold change), MYB113 (Cc07_g08030, 15-fold change), HSFA9 (Cc06_g11570, 4-fold change), and four WRKY genes with fold-changes higher than 3 (Cc04_g05080, Cc11_g12480, Cc09_g01430, and Cc01_g14950). Among these candidates, it is worth noting that the *DREB2G* and *HSFA9* genes had large sequence identity with their homologues in Arabidopsis and sunflower (see Supplementary Fig. S3), which have been shown to be involved in drought response and desiccation tolerance (Sakuma et al., 2006; Prieto-Dapena et al., 2006, 2008). Two other TFs with high homology to PLATZ (Cc08_g02880) and Delay of Germination-like (DOG-like 4, Cc02_g15140), which are involved in DT acquisition in Arabidopsis seeds (Gonzalez-Morales et al., 2016), were also strongly induced at the stage 5-6 transition (8-fold change for *DOGL4*). Finally, with regards to the master regulators of seed development and maturation, it is worth noting that *ABI3* (Cc01_g17380) and *FUS3* (Cc07_g01190) peaked at stages 4 and 5, respectively, i.e. during seed filling and before DT acquisition (see Supplementary Table S4). Unfortunately, the PUCE CAFE microarray lacks representative probes for coffee genes homologous to *AtABI4* (Cc02_g06230) and *AtABI5* (Cc06_g05680). *CaABI4* was not detected at either stage 5 or 6 in the 454 pyrosequencing-based transcriptome data, and very low transcript amounts were observed for *CaABI5* (Supplementary Table S8).

## Functional validation of HSFA9 and DREB2G

Transgenic coffee somatic embryos overexpressing the HSFA9 or the DREB2G cDNA sequences under the control of the doubled constitutive 35S promoter of the Cauliflower mosaic virus were regenerated by Agrobacterium-mediated transformation. All the lines tested (four independent lines for each gene) were positively transformed, as shown by the integration of the *hptII* gene (Supplementary Fig. S4). Real-time RT-PCR analysis demonstrated strong expression of the transgenes in the transformed lines. HSFA9 expression showed a 23-fold increase in HSFA9⁺ overexpressing somatic embryos, while an 8-fold increase was observed for DREB2G expression in DREB2G⁺ somatic embryos (Fig. 5A). Transgene functionality was then investigated through the analysis of putative downstream target expression. Six sHSP genes belonging to the top 100 up-regulated genes during DT acquisition were chosen as putative targets. All six sHSP genes were significantly up-regulated, with a 2- to 8-fold increase in somatic embryos overexpressing HSFA9 (Fig. 5B). A significant increase in gene expression level was also observed for two sHSP in somatic embryos overexpressing DREB2G. However, overexpression of HSFA9 or DREB2G and downstream sHSP genes was not sufficient to increase the DT of coffee somatic embryos, since there was no statistical difference in survival after equilibrium drying between the HSFA9⁺ or DREB2G⁺ somatic embryos and the controls (Fig. 5C).

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Functional validation of HSFA9 and DREB2G. Small heat-shock protein (sHSP) gene expression and desiccation tolerance. (A) Real-time RT-PCR analysis of HSFA9 and DREB2G ectopic expression in coffee somatic embryos overexpressing HSFA9 or DREB2G (HSFA9⁺ and DREB2⁺, respectively). Values are means (±SD) of four independent transgenic lines. (B) Real-time RT-PCR analysis of sHSP gene expression in HSFA9⁺ and DREB2⁺ transgenic somatic embryos. (C) Viability of HSFA9⁺ (circles), DREB2⁺ (squares), and control (diamonds) somatic coffee embryos after drying to equilibrium at 81%, 85%, 92%, and 97% relative humidity (RH).</figcaption></figure>

## Discussion

Using coffee seeds as a model, we searched for developmental processes that confer desiccation tolerance (DT) in orthodox seeds but that are lacking in intermediate seeds. We showed
that DT acquisition occurs late in coffee seed development, i.e. at the end of seed filling, while these events may overlap in seeds whose development is very short, such as Arabidopsis seeds (Leprince et al., 2017 for a review). As in all fleshy fruits, coffee seeds are not subjected to dry ambient air during maturation. DT is therefore acquired without dehydration in planta in developing coffee seeds. The programme that governs DT acquisition therefore does not appear to require an environmental cue to start, but appears to be strictly under developmental control. Transcriptome analysis of developing coffee seeds revealed the existence of a major transcriptional switch coinciding with DT acquisition. The period that we define as the late maturation stage in coffee seeds starts with this transcriptional switch at the stage 5-6 transition. It is characterised by the transcriptional activation of several stress-defence mechanisms (antioxidants, detoxication and repair systems, LEAs, HSPs, cell-wall reorganisation) and down-regulation of genes involved in energy processes, primary metabolism, and protein fate. Transcriptome analyses of developing *M. truncatula* and Arabidopsis seeds showed that the same activation and repression processes occur during the late maturation phase of orthodox seeds (see Leprince and Buitink, 2010 for a review). Moreover, none of the DT mechanisms inferred from orthodox seed transcriptomes appeared to be missing in maturing coffee seeds. This observation led us to scrutinise the role played by hormones, LEAs, sHSPs, and TFs thought to play a role in orthodox seed DT.

As in orthodox seeds, ABA appears to play a major role in controlling the late maturation programme in the coffee seed. First, almost all enzymatic steps in ABA biosynthesis were significantly upregulated during DT acquisition. Second, the TF (Cc11_g12740) that displayed the highest up-regulation during DT acquisition (45-fold change) showed very high sequence homology with the Arabidopsis gene *AtANAC029* (NAC III subfamily), which is required for normal Arabidopsis seed development and morphology (Kunieda et al., 2008) and was recently shown to control ABA biosynthesis, notably through trans-activation of the promoter of the abscisic aldehyde oxidase *AAO3* gene (Yang et al., 2014). Our results also suggest the existence of crosstalk between ABA and cytokinins during late maturation. DT acquisition coincided with a marked decrease in cytokinin levels and down-regulation of the cytokinin biosynthetic genes *IPT1* and *LOG1*. Antagonist actions of cytokinins and ABA have already been demonstrated in several biological processes, including seed germination (Guan et al., 2014) and plant survival under abiotic salt and drought stresses (Nishiyama et al., 2011).

Several other transcripts that play a key role in ABA homeostasis and signalling pathways strongly support a role for ABA in coffee seed DT acquisition. ABA 8'-hydroxylase (Cc08_g05920), which is involved in the main ABA catabolic route and regulates ABA levels in mature Arabidopsis seeds (Okamoto et al., 2006), was significantly up-regulated during the stage 5-6 transition. Key components of the ABA sensing (PYR1) and signalling pathways (PP2Cs and SRK12) also showed very high transcript accumulation during DT acquisition. They included the E3 ubiquitin-protein ligase RGLG1 (Cc07_g18890), which functions as a positive regulator of

ABA signalling in Arabidopsis through ABA-dependent degradation of PP2CA (Wu et al., 2016). PP2Cs have been shown to co-ordinate both ABA signalling and the SnRK1-mediated energy-sensing response (Rodrigues et al., 2013). SnRK1, a sensor of energy status that mediates extensive metabolic and developmental reprogramming in response to environmental stresses (Baena-González and Hanson, 2017), is also involved in the co-ordination of metabolic, hormonal, and developmental processes during seed development (Radchuk et al., 2010). SnRK1 is a heterotrimeric protein kinase complex composed of a catalytic subunit, α, and two regulatory subunits, β and γ. In our coffee seed transcriptome, gene expression for the SnRK1 subunit α (*KIN10*, Cc06_g20520) peaked at stage 5, together with *FUS3* with which it has been reported to interact physically to regulate transitions between developmental stages in Arabidopsis (Tsai and Gazzarrini, 2012), while genes for the regulatory subunits KINB1 and KING1 were specifically up-regulated during coffee seed DT acquisition. It is worth noting that silencing of SnRK1β (SNF4) impairs pollen hydration mechanisms in Arabidopsis through modulation of ROS signalling and mitochondrial biogenesis (Gao et al., 2016) and alters seed longevity in *M. truncatula* (Rosnoblet et al., 2007).

It should, however, be noted that the changes in ABA content measured during seed maturation did not reflect transcript accumulation of ABA biosynthetic genes at the stage 5-6 transition: the maximum ABA content was in fact observed at stage 4. It is therefore very likely that most ABA detected during seed filling (stages 4 and 5) was translocated from maternal tissues and then masked local biosynthesis during DT acquisition. This mechanism has been observed in many other plants (Frey et al., 2004; Kanno et al., 2010). The question as to whether translocated ABA and ABA synthesised locally play the same roles in late maturation induction requires further investigation. For instance, in many orthodox seed species, ABA produced by maternal tissues is not sufficient to induce dormancy, which strictly depends on ABA synthesised in the endosperm and the embryo (Kucera et al., 2005).

Several *in vitro* studies have shown that LEAs are intrinsically disordered proteins that can fold reversibly upon desiccation (Tolleter et al., 2007; Hundertmark et al., 2012), and either stabilise membranes (Tolleter et al., 2010; Popova et al., 2011) or prevent protein misfolding and aggregation (Pouchkina-Stantcheva et al., 2007; Boucher et al., 2010; Dang et al., 2014). We identified 29 LEA proteins in the coffee genome. Our transcriptome and proteome analyses revealed that 19 of them were massively transcribed and 17 were translated during coffee seed maturation, and 15 of them were seed-specific. The closest Arabidopsis homologue of each of this subset of 19 LEA genes is also highly transcribed in the maturing Arabidopsis seed (Hundertmark and Hincha, 2008; also visualised with the eFP browser by Winter et al., 2007; Supplementary Table S7). Within this maturation subset, 13 LEA genes were significantly up-regulated during DT acquisition. The closest Arabidopsis homologues of nine LEAs of this coffee seed DT subset are associated with DT re-induction in pre-germinated Arabidopsis seeds (Supplementary Table
S7; Maia et al., 2011) and three of them also display high homology with one of the LEAs involved in DT re-induction in *M. truncatula* pre-germinated seeds (Supplementary Table S7; Boudet et al., 2006). Also within this DT subset, we found a LEA homologous to EM6 (73% identity at the amino acid level), which plays a critical role in water binding during Arabidopsis seed maturation (Manfre et al., 2009). In conclusion, our focus on LEA proteins did not reveal a specific lack in coffee seeds in comparison with orthodox seeds. Although accumulation of LEAs is a process shared by most anhydروbiotes, and ectopic expression of plant LEA genes may enhance DT in yeast (Dang et al., 2014), the specific role of LEAs in seed DT acquisition remains to be revealed. The timing of accumulation of most LEAs coincides with longevity rather than with DT acquisition during orthodox seed development (Chatelain et al., 2012; Leprince et al., 2017). Several Arabidopsis mutants that lack one or two LEA proteins are still desiccation tolerant (Carles et al., 2002; Hundertmark et al., 2011). Silencing of dehydrins alters Arabidopsis seed longevity under controlled deterioration without affecting DT (Hundertmark et al., 2011). Similarly, *Castanospermum australe* seeds accumulate large amounts of dehydrins but are recalcitrant (Delahaie et al., 2013). How LEAs contribute to seed DT thus requires further clarification.

sHSPs are ATP-independent chaperones that prevent protein aggregation by maintaining them in a folding-competent state and that refold misfolded proteins alone or in concert with other ATP-dependent chaperones. Several seed-specific heat-shock factors, including HSFA9 (Kotak et al., 2007), are responsible for HSP transcriptional activation and accumulation during late maturation in Arabidopsis (Wehmeyer et al., 1996). The hypothesis that sHSPs contribute to seed DT has been supported by enhanced DT in desiccation-sensitive vegetative organs ectopically expressing HSFs (Prieto-Dapena et al., 2008; Personat et al., 2014). In the present study, we also demonstrated increased sHSP transcript accumulation in coffee somatic embryos overexpressing HSFA9 or DREB2. By contrast, HSFA9 or DREB2 overexpression did not confer desiccation tolerance to transformed embryos. In common with LEA proteins, while several studies suggest that sHSPs play a major role in seed longevity (Prieto-Dapena et al., 2006; Almoguera et al., 2009; Kaur et al., 2015), their contribution to seed DT requires further investigation.

In conclusion, the present work revealed that the cellular and regulation processes that occur during late maturation in the coffee seed are remarkably similar to those known or thought to be involved in orthodox seed DT. Two hypotheses, which are not mutually exclusive, thus deserve to be explored to understand the impaired tolerance to drying of intermediate seeds. First, the differences are quantitative. For instance, the metabolic shutdown may be only partial in intermediate seeds whereas it is hypothesised to be complete in orthodox seeds, or the amount of a certain type of LEA accumulated in intermediate seeds may be smaller than in orthodox seeds. Second, intermediate seeds may have certain particularities that are not modified by transcriptional changes during DT acquisition. For instance, a peculiar membrane lipid composition of these tropical seeds could result in specific lethal

phospholipid phase transitions at low RH that would not be prevented by membrane stabilisers up-regulated during DT acquisition. The interspecific variability of seed DT that exists within the genus *Coffea* (Dussert et al., 1999) could be very useful to test these two hypotheses.

## Supplementary data

Supplementary data are available at JXB online.

Table S1. List of specific primers used for real-time RT-PCR and PCR.

Table S2. Proteins detected in the heat-stable proteome of the mature coffee seed.

Table S3. Phytohormone profiling in developing coffee seed: ESI-MS-MS analysis of ABA, auxin, cytokinin, and gibberellin families.

Table S4. Relative expression of the 15 522 genes represented on the oligonucleotide PUCE CAFE microarray.

Table S5. Gene composition of differentially down- (cluster C5) and up-regulated genes (clusters C18–C26) during DT acquisition.

Table S6. GO and MapMan enrichment analyses of clusters of down-regulated genes (cluster C5) and up-regulated genes (clusters C18–C26) during DT acquisition.

Table S7. Genome-wide analysis of the coffee LEA gene family.

Table S8. Endosperm and embryo transcriptomes (454) during the stage 5–6 transition.

Fig. S1. Hierarchical cluster analysis of the 5857 genes displaying significant up- or down-regulation at one of the six developmental transitions.

Fig. S2. Phylogenetic tree of the 29 coffee and 51 *Arabidopsis thaliana* LEA proteins.

Fig. S3. Comparisons of the amino acid sequence of HSFA9 proteins from *Helianthus annuus* (AAM43804) and coffee (Cc06_g11570), and of the amino acid sequence of DREB2G proteins from *Arabidopsis thaliana* (At5g18450) and coffee (Cc09_g03140).

Fig. S4. PCR detection of the HPTII hygromycin-resistance gene in transgenic coffee plants.

## Acknowledgments

Aldecinei Bastos-Siqueira benefited from a CNPq post-doctoral grant from the Science Without Borders Program, Brazil.

## References

**Almoguera C**, Prieto-Dapena P, Díaz-Martín J, Espinosa JM, Carranco R, Jordano J. 2009. The HaDREB2 transcription factor enhances basal thermotolerance and longevity of seeds through functional interaction with HaHSFA9. BMC Plant Biology **9**, 75.

**Audic S**, Claverie JM. 1997. The significance of digital gene expression profiles. Genome Research **7**, 986–995.

**Baena-González E**, Hanson J. 2017. Shaping plant development through the SnRK1-TOR metabolic regulators. Current Opinion in Plant Biology **35**, 152–157.

**Bailly C.** 2004. Active oxygen species and antioxidants in seed biology. Seed Science Research **14**, 93–107.
Berjak P, Pammenter NW. 2008. From *Avicennia* to *Zizania*: seed recalcitrance in perspective. *Annals of Botany* **101**, 213–228.

Black M, Pritchard H. eds. 2002. Desiccation and survival in plants: drying without dying. Wallingford: CABI Publishing.

Boca S, Koestler F, Ksas B, Chevalier A, Leymarie J, Fekete A, Mueller MJ, Havaux M. 2014. Arabidopsis lipocalins AtCHL and AtTIL have distinct but overlapping functions essential for lipid protection and seed longevity. *Plant, Cell & Environment* **37**, 368–381.

Boucher V, Buitink J, Lin X, Boudet J, Hoekstra FA, Hundertmark M, Renard D, Leprince O. 2010. MtPM25 is an atypical hydrophobic late embryogenesis-abundant protein that dissociates cold and desiccation-aggregated proteins. *Plant, Cell & Environment* **33**, 418–430.

Boudet J, Buitink J, Hoekstra FA, Rogniaux H, Larré C, Satour P, Leprince O. 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 Physiology* **140**, 1418–1436.

Buitink J, Leger JJ, Guisle I, et al. 2006. Transcriptome profiling uncovers metabolic and regulatory processes occurring during the transition from desiccation-sensitive to desiccation-tolerant stages in *Medicago truncatula* seeds. *The Plant Journal* **47**, 735–750.

Buitink J, Leprince O. 2008. Intracellular glasses and seed survival in the dry state. *Comptes Rendus Biologies* **331**, 788–795.

Carles C, Bies-Etheve N, Aspart L, Léon-Kloosterziel KM, Koornneef M, Echeverria M, Delseny M. 2002. Regulation of *Arabidopsis thaliana* Em genes: role of ABI5. *The Plant Journal* **30**, 373–383.

Chabrillange N, Dussert S, Engelmann F, Doulbeau S, Hamon S. 2000. Desiccation tolerance in relation to soluble sugar contents in seeds of ten coffee (*Coffea* L.) species. *Seed Science Research* **10**, 393–396.

Chatelain E, Hundertmark M, Leprince O, Le Gall S, Satour P, Deligny-Penninck S, Rogniaux H, Buitink J. 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 & Environment* **35**, 1440–1455.

Chiwocha SD, Abrams SR, Ambrose SJ, Cutler AJ, Loewen M, Ross AR, Kermode AR. 2003. A method for profiling classes of plant hormones and their metabolites using liquid chromatography-electrospray ionization tandem mass spectrometry: an analysis of hormone regulation of thermodormancy of lettuce (*Lactuca sativa* L.) seeds. *The Plant Journal* **35**, 405–417.

Crowe JH. 2007. Trehalose as a “chemical chaperone”. Fact and fantasy. In: Csermely P, Vigh L, eds. Molecular aspects of the stress response: chaperones, membranes and networks, vol. 594. New York: Springer, 143–158.

Crowe JH, Hoekstra FA, Crowe LM. 1992. Anhydrobiosis. Annual Review of Physiology **54**, 579–599.

Cruz F, Kalaoun S, Nobile P, Colombo C, Almeida J, Barros LMG, Romano E, Grossi-de-Sá MF, Vaslin M, Alves-Ferreira M. 2009. Evaluation of coffee reference genes for relative expression studies by quantitative real-time RT-PCR. *Molecular Breeding* **23**, 607–616.

Dang NX, Popova AV, Hundertmark M, Hincha DK. 2014. Functional characterization of selected LEA proteins from *Arabidopsis thaliana* in yeast and *in vitro*. *Planta* **240**, 325–336.

Decourcelle M, Perez-Fons L, Baulande S, et al. 2015. Combined transcript, proteome, and metabolite analysis of transgenic maize seeds engineered for enhanced carotenoid synthesis reveals pleotropic effects in core metabolism. *Journal of Experimental Botany* **66**, 3141–3150.

Delahaie J, Hundertmark M, Bove J, Leprince O, Rogniaux H, Buitink J. 2013. LEA polypeptide profiling of recalcitrant and orthodox legume seeds reveals ABI3-regulated LEA protein abundance linked to desiccation tolerance. *Journal of Experimental Botany* **64**, 4559–4573.

Denoeud F, Carretero-Paulet L, Dereeper A, et al. 2014. The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. *Science* **345**, 1181–1184.

Dinakar C, Bartels D. 2013. Desiccation tolerance in resurrection plants: new insights from transcriptome, proteome and metabolome analysis. *Frontiers in Plant Science* **4**, 482.

Dussert S, Chabrillange N, Engelmann F, Anthony F, Louarn J, Hamon S. 2000. Relationship between seed desiccation sensitivity, seed water content at maturity and climatic characteristics of native environments of nine *Coffea* L. species. *Seed Science Research* **10**, 293–300.

Dussert S, Chabrillange N, Engelmann F, Hamon S. 1999. Quantitative estimation of seed desiccation sensitivity using a quantal response model: application to nine species of the genus *Coffea* L. *Seed Science Research* **9**, 135–144.

Dussert S, Chabrillange N, Rocquelin G, Engelmann F, Lopez M, Hamon S. 2001. Tolerance of coffee (*Coffea* spp.) seeds to ultra-low temperature exposure in relation to calorimetric properties of tissue water, lipid composition, and cooling procedure. *Physiologia Plantarum* **112**, 495–504.

Dussert S, Davey MW, Laffargue A, Doulbeau S, Swennen R, Etienne H. 2006. Oxidative stress, phospholipid loss and lipid hydrolysis during drying and storage of intermediate seeds. *Physiologia Plantarum* **127**, 192–204.

Dussert S, Guerin C, Andersson M, et al. 2013. Comparative transcriptome analysis of three oil palm fruit and seed tissues that differ in oil content and fatty acid composition. *Plant Physiology* **162**, 1337–1358.

Eira MTS, Amaral da Silva EA, de Castro RD, Dussert S, Walters C, Bewley JD, Hilhorst WWM. 2006. Coffee seed physiology. *Brazilian Journal of Plant Physiology* **18**, 149–163.

Eisen MB, Spellman PT, Brown PO, Botstein D. 1998. Cluster analysis and display of genome-wide expression patterns. *Proceedings of the National Academy of Sciences, USA* **95**, 14863–14868.

Ellis RH, Hong TD, Roberts EH. 1990. An intermediate category of seed storage behaviour? 1. Coffee. *Journal of Experimental Botany* **41**, 1167–1174.

Etienne H. 2005. Somatic embryogenesis protocol: Coffee (*Coffea arabica* L. and *C. canephora* P.). In: Jain SM, Gupta PK, eds. *Protocol for somatic embryogenesis in woody plants*, vol. **77**. Springer Netherlands, 167–179.

Frey A, Godin B, Bonnet M, Sotta B, Marion-Poll A. 2004. Maternal synthesis of abscisic acid controls seed development and yield in *Nicotiana plumbaginifolia*. *Planta* **218**, 958–964.

Gaff DF, Oliver M. 2013. The evolution of desiccation tolerance in angiosperm plants: a rare yet common phenomenon. *Functional Plant Biology* **40**, 315–328.

Gao XQ, Liu CZ, Li DD, Zhao TT, Li F, Jia XN, Zhao XY, Zhang XS. 2016. The Arabidopsis KINβγ subunit of the SnRK1 complex regulates pollen hydration on the stigma by mediating the level of reactive oxygen species in pollen. *PLoS Genetics* **12**, e1006228.

Gonzalez-Morales SI, Chavez-Montes RA, Hayano-Kanashiro C, Alejo-Jacuinde G, Rico-Cambron TY, de Folter S, Herrera-Estrella L. 2016. Regulatory network analysis reveals novel regulators of seed desiccation tolerance in *Arabidopsis thaliana*. *Proceedings of the National Academy of Sciences, USA* **113**, E5232–E5241.

Götz S, García-Gómez JM, Terol J, et al. 2008. High-throughput functional annotation and data mining with the Blast2GO suite. *Nucleic Acids Research* **36**, 3420–3435.

Guan C, Wang X, Feng J, Hong S, Liang Y, Ren B, Zuo J. 2014. Cytokinin antagonizes abscisic acid-mediated inhibition of cotyledon greening by promoting the degradation of abscisic acid insensitive5 protein in Arabidopsis. *Plant Physiology* **164**, 1515–1526.

Gusev O, Suetsugu Y, Cornette R, et al. 2014. Comparative genome sequencing reveals genomic signature of extreme desiccation tolerance in the anhydrobiotic midge. *Nature Communications* **5**, 4784.

Gutierrez L, Van Wuytswinkel O, Castelain M, Bellini C. 2007. Combined networks regulating seed maturation. *Trends in Plant Science* **12**, 294–300.

Hamilton KN, Offord CA, Cuneo P, Deseo MA. 2013. A comparative study of seed morphology in relation to desiccation tolerance and other physiological responses in 71 Eastern Australian rainforest species. *Plant Species Biology* **28**, 51–62.

Hoekstra FA, Golovina EA, Buitink J. 2001. Mechanisms of plant desiccation tolerance. *Trends in Plant Science* **6**, 431–438.

Hundertmark M, Buitink J, Leprince O, Hincha DK. 2011. The reduction of seed-specific dehydrins reduces seed longevity in *Arabidopsis thaliana*. *Seed Science Research* **21**, 165–173.
**Hundertmark M, Hincha DK.** 2008. LEA (late embryogenesis abundant) proteins and their encoding genes in *Arabidopsis thaliana*. BMC Genomics **9**, 118.

**Hundertmark M, Popova AV, Rausch S, Seckler R, Hincha DK.** 2012. Influence of drying on the secondary structure of intrinsically disordered and globular proteins. Biochemical and Biophysical Research Communications **417**, 122–128.

**Joët T, Salmona J, Laffargue A, Descroix F, Dussert S.** 2010. Use of the growing environment as a source of variation to identify the quantitative trait transcripts and modules of co-expressed genes that determine chlorogenic acid accumulation. Plant, Cell & Environment **33**, 1220–1233.

**Joët T, Laffargue A, Salmona J, Doulbeau S, Descroix F, Bertrand B, de Kochko A, Dussert S.** 2009. Metabolic pathways in tropical dicotyledonous albuminous seeds: *Coffea arabica* as a case study. New Phytologist **182**, 146–162.

**Joët T, Laffargue A, Salmona J, Doulbeau S, Descroix F, Bertrand B, Lashermes P, Dussert S.** 2014. Regulation of galactomannan biosynthesis in coffee seeds. Journal of Experimental Botany **65**, 323–337.

**Kanno Y, Jikumaru Y, Hanada A, Nambara E, Abrams SR, Kamiya Y, Seo M.** 2010. Comprehensive hormone profiling in developing Arabidopsis seeds: examination of the site of ABA biosynthesis, ABA transport and hormone interactions. Plant & Cell Physiology **51**, 1988–2001.

**Karssen CM, Brinkhorst-van der Swan DL, Creekland AE, Koornneef M.** 1983. Induction of dormancy during seed development by endogenous abscisic acid: studies on abscisic acid deficient genotypes of *Arabidopsis thaliana* (L.) Heynh. Planta **157**, 158–165.

**Kaur H, Petla BP, Kamble NU, Singh A, Rao V, Salvi P, Ghosh S, Majee M.** 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. Frontiers in Plant Science **6**, 713.

**Kotak S, Vierling E, Bäumlein H, von Koskull-Döring P.** 2007. A novel transcriptional cascade regulating expression of heat stress proteins during seed development of Arabidopsis. The Plant Cell **19**, 182–195.

**Kucera B, Cohn MA, Leubner-Metzger G.** 2005. Plant hormone interactions during seed dormancy release and germination. Seed Science Research **15**, 281–307.

**Kumar MN, Verslues PE.** 2015. Stress physiology functions of the Arabidopsis histidine kinase cytokinin receptors. Physiologia Plantarum **154**, 369–380.

**Kunieda T, Mitsuda N, Ohme-Takagi M, Takeda S, Aida M, Tasaka M, Kondo M, Nishimura M, Hara-Nishimura I.** 2008. NAC family proteins NARS1/NAC2 and NARS2/NAM in the outer integument regulate embryogenesis in Arabidopsis. The Plant Cell **20**, 2631–2642.

**Leprince O, Buitink J.** 2010. Desiccation tolerance: from genomics to the field. Plant Science **179**, 554–564.

**Leprince O, Buitink J, Hoekstra FA.** 1999. Axes and cotyledons of recalcitrant seeds of *Castanea sativa* Mill. exhibit contrasting responses of respiration to drying in relation to desiccation sensitivity. Journal of Experimental Botany **50**, 1515–1524.

**Leprince O, Pellizzaro A, Berriri S, Buitink J.** 2017. Late seed maturation: drying without dying. Journal of Experimental Botany **68**, 827–841.

**Li DZ, Pritchard HW.** 2009. The science and economics of *ex situ* plant conservation. Trends in Plant Science **14**, 614–621.

**Li H.** 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997v2.

**Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R.** 2009. Subgroup 1000 genome project data processing: the sequence alignment/map format and SAMtools. Bioinformatics **25**, 2078–2079.

**Maia J, Dekkers BJ, Provart NJ, Ligterink W, Hilhorst HW.** 2011. The re-establishment of desiccation tolerance in germinated *Arabidopsis thaliana* seeds and its associated transcriptome. Plos ONE **6**, e29123.

**Manfre AJ, LaHatte GA, Climer CR, Marcotte WR Jr.** 2009. Seed dehydration and the establishment of desiccation tolerance during seed maturation is altered in the *Arabidopsis thaliana* mutant *atem6-1*. Plant & Cell Physiology **50**, 243–253.

**Mönke G, Seifert M, Keilwagen J, et al.** 2012. Toward the identification and regulation of the *Arabidopsis thaliana* ABI3 regulon. Nucleic Acids Research **40**, 8240–8254.

**Nishiyama R, Watanabe Y, Fujita Y, et al.** 2011. Analysis of cytokinin mutants and regulation of cytokinin metabolic genes reveals important regulatory roles of cytokinins in drought, salt and abscisic acid responses, and abscisic acid biosynthesis. The Plant Cell **23**, 2169–2183.

**Okamoto M, Kuwahara A, Seo M, Kushiro T, Asami T, Hirai N, Kamiya Y, Koshiba T, Nambara E.** 2006. CYP707A1 and CYP707A2, which encode abscisic acid 8'-hydroxylases, are indispensable for proper control of seed dormancy and germination in Arabidopsis. Plant Physiology **141**, 97–107.

**Oliver MJ, Guo L, Alexander DC, Ryals JA, Wone BW, Cushman JC.** 2011. A sister group contrast using untargeted global metabolomic analysis delineates the biochemical regulation underlying desiccation tolerance in *Sporobolus stapfianus*. The Plant Cell **23**, 1231–1248.

**Personat JM, Tejedor-Cano J, Prieto-Dapena P, Almoguera C, Jordano J.** 2014. Co-overexpression of two Heat Shock Factors results in enhanced seed longevity and in synergistic effects on seedling tolerance to severe dehydration and oxidative stress. BMC Plant Biology **14**, 56.

**Popova AV, Hundertmark M, Seckler R, Hincha DK.** 2011. Structural transitions in the intrinsically disordered plant dehydration stress protein LEA7 upon drying are modulated by the presence of membranes. Biochimica et Biophysica Acta-Biomembranes **1808**, 1879–1887.

**Pouchkina-Stantcheva NN, McGee BM, Boschetti C, et al.** 2007. Functional divergence of former alleles in an ancient asexual invertebrate. Science **318**, 268–271.

**Prieto-Dapena P, Castaño R, Almoguera C, Jordano J.** 2006. Improved resistance to controlled deterioration in transgenic seeds. Plant Physiology **142**, 1102–1112.

**Prieto-Dapena P, Castano R, Almoguera C, Jordano J.** 2008. The ectopic overexpression of a seed-specific transcription factor, HaHSFA9, confers tolerance to severe dehydration in vegetative organs. The Plant Journal **54**, 1004–1014.

**Privat I, Bardil A, Gomez AB, et al.** 2011. The 'PUCE CAFE' Project: the first 15K coffee microarray, a new tool for discovering candidate genes correlated to agronomic and quality traits. BMC Genomics **12**, 5.

**Provart NJ, Gil P, Chen W, Han B, Chang HS, Wang X, Zhu T.** 2003. Gene expression phenotypes of Arabidopsis associated with sensitivity to low temperatures. Plant Physiology **132**, 893–906.

**Radchuk R, Emery RJ, Weier D, Vigeolas H, Geigenberger P, Lunn JE, Feil R, Weschke W, Weber H.** 2010. Sucrose non-fermenting kinase 1 (SnRK1) coordinates metabolic and hormonal signals during pea cotyledon growth and differentiation. The Plant Journal **61**, 324–338.

**Ribas AF, Dechamp E, Champion A, Bertrand B, Combes MC, Verdeil JL, Lapeyre F, Lashermes P, Etienne H.** 2011. *Agrobacterium*-mediated genetic transformation of *Coffea arabica* (L.) is greatly enhanced by using established embryogenic callus cultures. BMC Plant Biology **11**, 92.

**Roberts EH.** 1973. Predicting the storage life of seeds. Seed Science and Technology **1**, 499–514.

**Rodrigues A, Adamo M, Crozet P, et al.** 2013. ABI1 and PP2CA phosphatases are negative regulators of Snf1-related protein kinase1 signaling in Arabidopsis. The Plant Cell **25**, 3871–3884.

**Rosnoblet C, Aubry C, Leprince O, Vu BL, Rogniaux H, Buitink J.** 2007. The regulatory gamma subunit SNF4b of the sucrose non-fermenting-related kinase complex is involved in longevity and stachyose accumulation during maturation of *Medicago truncatula* seeds. The Plant Journal **51**, 47–59.

**Sakuma Y, Maruyama K, Qin F, Osakabe Y, Shinozaki K, Yamaguchi-Shinozaki K.** 2006. Dual function of an Arabidopsis transcription factor DREB2A in water-stress-responsive and heat-stress-responsive gene expression. Proceedings of the National Academy of Sciences, USA **103**, 18822–18827.

**Salmona J, Dussert S, Descroix F, de Kochko A, Bertrand B, Joët T.** 2008. Deciphering transcriptional networks that govern *Coffea arabica* seed development using combined cDNA array and real-time RT-PCR approaches. Plant Molecular Biology **66**, 105–124.

**Santos-Mendoza M, Dubreucq B, Baud S, Parcy F, Caboche M, Lepiniec L.** 2008. Deciphering gene regulatory networks that control seed development and maturation in Arabidopsis. The Plant Journal **54**, 608–620.

**Shin JH, Kim SR, An G.** 2009. Rice aldehyde dehydrogenase7 is needed for seed maturation and viability. Plant Physiology **149**, 905–915.
Staswick PE, Serban B, Rowe M, Tiryaki I, Maldonado MT, Maldonado MC, Suza W. 2005. Characterization of an Arabidopsis enzyme family that conjugates amino acids to indole-3-acetic acid. The Plant Cell **17**, 616–627.

Tapia H, Koshland DE. 2014. Trehalose is a versatile and long-lived chaperone for desiccation tolerance. Current Biology **24**, 2758–2766.

To A, Valon C, Savino G, Guilleminot J, Devic M, Giraudat J, Parcy F. 2006. A network of local and redundant gene regulation governs Arabidopsis seed maturation. The Plant Cell **18**, 1642–1651.

Tognetti VB, Van Aken O, Morreel K, et al. 2010. Perturbation of indole-3-butyric acid homeostasis by the UDP-glucosyltransferase UGT74E2 modulates Arabidopsis architecture and water stress tolerance. The Plant Cell **22**, 2660–2679.

Tolleter D, Hincha DK, Macherel D. 2010. A mitochondrial late embryogenesis abundant protein stabilizes model membranes in the dry state. Biochimica et Biophysica Acta-Biomembranes **1798**, 1926–1933.

Tolleter D, Jaquinod M, Mangavel C, et al. 2007. Structure and function of a mitochondrial late embryogenesis abundant protein are revealed by desiccation. The Plant Cell **19**, 1580–1589.

Tran LS, Urao T, Qin F, Maruyama K, Kakimoto T, Shinozaki K, Yamaguchi-Shinozaki K. 2007. Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in Arabidopsis. Proceedings of the National Academy of Sciences, USA **104**, 20623–20628.

Tsai AY, Gazzarrini S. 2012. AKIN10 and FUSCA3 interact to control lateral organ development and phase transitions in Arabidopsis. The Plant Journal **69**, 809–821.

Tweddle JC, Dickie JB, Baskin CC, Baskin JM. 2003. Ecological aspects of seed desiccation sensitivity. Journal of Ecology **91**, 294–304.

Verdier J, Lalanne D, Pelletier S, 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 Physiology **163**, 757–774.

Walters C, Berjak P, Pammenter N, Kennedy K, Raven P. 2013. Preservation of recalcitrant seeds. Science **339**, 915–916.

Wang C, Grohme MA, Mali B, Schill RO, Frohme M. 2014. Towards decrypting cryptobiosis–analyzing anhydrobiosis in the tardigrade *Milnesium tardigradum* using transcriptome sequencing. Plos ONE **9**, e92663.

Wehmeyer N, Hernandez LD, Finkelstein RR, Vierling E. 1996. Synthesis of small heat-shock proteins is part of the developmental program of late seed maturation. Plant Physiology **112**, 747–757.

Wettenhall JM, Smyth GK. 2004. limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics **20**, 3705–3706.

Widemann E, Miesch L, Lugan R, Holder E, Heinrich C, Aubert Y, Miesch M, Pinot F, Heitz T. 2013. The amidohydrolases IAR3 and ILL6 contribute to jasmonoyl-isoleucine hormone turnover and generate 12-hydroxyjasmonic acid upon wounding in Arabidopsis leaves. The Journal of Biological Chemistry **288**, 31701–31714.

Winter D, Vinegar B, Nahal H, Ammar R, Wilson GV, Provart NJ. 2007. An “Electronic Fluorescent Pictograph” browser for exploring and analyzing large-scale biological data sets. Plos ONE **2**, e718.

Wu Q, Zhang X, Peirats-Llobet M, Belda-Palazon B, Wang X, Cui S, Yu X, Rodriguez PL, An C. 2016. Ubiquitin ligases RGLG1 and RGLG5 regulate abscisic acid signaling by controlling the turnover of phosphatase PP2CA. The Plant Cell **28**, 2178–2196.

Wyse SB, Dickie JB. 2017. Predicting the global incidence of seed desiccation sensitivity. Journal of Ecology **105**, 1082–1093.

Yang J, Worley E, Udvardi M. 2014. A NAP-AAO3 regulatory module promotes chlorophyll degradation via ABA biosynthesis in Arabidopsis leaves. The Plant Cell **26**, 4862–4874.

Yobi A, Wone BWM, Xu WX, Alexander DC, Guo LN, Ryals JA, Oliver MJ, Cushman JC. 2012. Comparative metabolic profiling between desiccation-sensitive and desiccation-tolerant species of *Selaginella* reveals insights into the resurrection trait. The Plant Journal **72**, 983–999.

Zinsmeister J, Lalanne D, Terrasson E, et al. 2016. ABI5 is a regulator of seed maturation and longevity in legumes. The Plant Cell **28**, 2735–2754.