A Conserved Carbon Starvation Response Underlies Bud Dormancy in Woody and Herbaceous Species

Plant shoot systems give rise to characteristic above-ground plant architectures. Shoots are formed from axillary meristems and buds, whose growth and development is modulated by systemic and local signals. These cues convey information about nutrient and water availability, light quality, sink/source organ activity and other variables that determine the timeliness and competence to maintain development of new shoots. This information is translated into a local response, in meristems and buds, of growth or quiescence. Although some key genes involved in the onset of bud latency have been identified, the gene regulatory networks (GRNs) controlled by these genes are not well defined. Moreover, it has not been determined whether bud dormancy induced by environmental cues, such as a low red-to-far-red light ratio, shares genetic mechanisms with bud latency induced by other causes, such as apical dominance or a short-day photoperiod. Furthermore, the evolution and conservation of these GRNs throughout angiosperms is not well established. We have reanalyzed public transcriptomic datasets that compare quiescent and active axillary buds of Arabidopsis, with datasets of axillary buds of the woody species Vitis vinifera (grapevine) and apical buds of Populus tremula x Populus alba (poplar) during the bud growth-to-dormancy transition. Our aim was to identify potentially common GRNs induced during the process that leads to bud para-, eco- and endodormancy. In Arabidopsis buds that are entering eco- or paradormancy, we have identified four induced interrelated GRNs that correspond to a carbon (C) starvation syndrome, typical of tissues undergoing low C supply. This response is also detectable in poplar and grapevine buds before and during the transition to dormancy. In all eukaryotes, C-limiting conditions are coupled to growth arrest and latency like that observed in dormant axillary buds. Bud dormancy might thus be partly a consequence of the underlying C starvation syndrome triggered by environmental and endogenous cues that anticipate or signal conditions unfavorable for sustained shoot growth.

Plant shoot systems give rise to characteristic above-ground plant architectures. Shoots are formed from axillary meristems and buds, whose growth and development is modulated by systemic and local signals. These cues convey information about nutrient and water availability, light quality, sink/source organ activity and other variables that determine the timeliness and competence to maintain development of new shoots. This information is translated into a local response, in meristems and buds, of growth or quiescence. Although some key genes involved in the onset of bud latency have been identified, the gene regulatory networks (GRNs) controlled by these genes are not well defined. Moreover, it has not been determined whether bud dormancy induced by environmental cues, such as a low red-to-far-red light ratio, shares genetic mechanisms with bud latency induced by other causes, such as apical dominance or a short-day photoperiod. Furthermore, the evolution and conservation of these GRNs throughout angiosperms is not well established. We have reanalyzed public transcriptomic datasets that compare quiescent and active axillary buds of Arabidopsis, with datasets of axillary buds of the woody species Vitis vinifera (grapevine) and apical buds of Populus tremula x Populus alba (poplar) during the bud growth-to-dormancy transition. Our aim was to identify potentially common GRNs induced during the process that leads to bud para-, eco-and endodormancy. In Arabidopsis buds that are entering eco-or paradormancy, we have identified four induced interrelated GRNs that correspond to a carbon (C) starvation syndrome, typical of tissues undergoing low C supply. This response is also detectable in poplar and grapevine buds before and during the transition to dormancy. In all eukaryotes, C-limiting conditions are coupled to growth arrest and latency like that observed in dormant axillary buds. Bud dormancy might thus be partly a consequence of the underlying C starvation syndrome triggered by environmental and endogenous cues that anticipate or signal conditions unfavorable for sustained shoot growth.

INTRODUCTION
Shoot branching patterns define overall above-ground plant architecture. In angiosperms, shoots are formed from axillary meristems initiated at the base of leaves. These meristems grow and develop into axillary buds that contain, preformed, most of the elements of adult branches (shoot meristems, leaf primordia, reproductive meristems). Axillary buds can enter a quiescent state, rather than growing out immediately to give a branch. In this latent or dormant state their metabolic activity and cell division are very limited (Shimizu and Mori, 1998;Ruttink et al., 2007). Bud dormancy and bud activation are influenced by environmental signals such as nutrient and water availability, light quality, day-length and temperature, and by endogenous signals such as sink/source organ activity and hormone signaling. Once dormant, buds require changes in specific developmental and/or environmental cues to resume growth and generate an elongated branch. These cues are monitored in different organs, and inform the plant as to when develop new shoots. This information is transduced to the bud and translated into a gene response that leads to quiescence or growth activation (Rameau et al., 2015).
Bud dormancy is therefore an adaptive trait that allows plants to endure adverse situations until conditions are favorable for development of new shoots. It has great impact on plant reproductive success and productivity, and on survival in temperate woody species. Evolution of this trait might have allowed plants to colonize habitats with fluctuating conditions not always suitable for sustained, uninterrupted growth. Depending on the type of stimulus that promotes growth arrest, Lang et al. (1985Lang et al. ( , 1987 distinguished three types of bud dormancy. When dormancy is induced by environmental factors, it is termed ecodormancy; when promoted by other plant organs it is paradormancy or correlative inhibition, and when it is maintained by signals internal to the bud and can only be reversed under certain conditions it is defined as endodormancy. In woody plants, axillary buds undergo transitions between different dormant states throughout the year. Paradormant buds enter endodormancy in response to changes in daylength and temperature. Chilling promotes transition from endo-to ecodormancy, after which the buds are susceptible to grow in response to mild temperatures (Rohde and Bhalerao, 2007).
Transcriptomic studies have been carried out in several herbaceous and woody species to define expression changes in buds during the transitions into and out of different types of dormancy, in response to changes in daylength, light quality, and apical dominance, and in mutant genotypes in which bud growth is affected (e.g., Tatematsu et al., 2005;Ruttink et al., 2007;González-Grandío et al., 2013;Reddy et al., 2013;Ueno et al., 2013;Porto et al., 2015). The GRNs that act inside the bud to control the stages leading to dormancy nonetheless remain little known. It is also largely unknown whether different types of dormancy share common underlying genetic mechanisms. Even less is known about the degree of conservation and evolution of the genetic control of this process in different plant species. Comparative analyses to identify common themes among different types of dormancy, or across species, are scarce (González-Grandío and Cubas, 2014;Fennell et al., 2015;Howe et al., 2015;Hao et al., 2017). Such comparisons could help us determine whether eco-, paraand endodormancy are variations of a single ancestral genetic program or whether each type is controlled by unrelated GRNs. It also will help elucidate whether GRNs that cause bud growth arrest are conserved in different herbaceous and woody plant species.
The master regulators that locally control the dormancy onset are also largely unknown. The best characterized are the genes that encode the TCP transcription factors (TF) teosinte branched1 (Tb1, Doebley et al., 1997), BRANCHED1 (BRC1, Aguilar-Martínez et al., 2007;Finlayson, 2007) and their orthologs in mono-and dicotyledonous species, respectively. These widely conserved factors play a very important role in the regulation of para-and ecodormancy in herbaceous plants. These genes are expressed in axillary buds and promote bud dormancy in response to fluctuating environmental cues such as light quality and quantity, and endogenous signals such as apical dominance, sugar availability and hormone signaling (reviewed in . In Arabidopsis thaliana, BRC1 controls transcription of several GRNs in buds; one, positively controlled by BRC1, leads to abscisic acid (ABA) accumulation and signaling, whereas another two that are downregulated by BRC1 are enriched in ribosomal protein genes in one case, and in cell division and DNA replication genes in the other (González-Grandío et al., 2013. Additional GRNs controlled by BRC1 remain to be characterized. In this study our aim was to identify potentially common GRNs induced during the process that leads to bud para-, eco-and endodormancy. For that we compared publicly available transcriptomic data from active paraand ecodormant axillary buds of Arabidopsis, and found, induced in dormant buds, a shared transcriptomic response typical of tissue undergoing C starvation. We then detected this response also in Populus tremula × Populus alba (poplar) apical buds undergoing endodormancy and in Vitis vinifera (grapevine) axillary buds entering para-, endo-and ecodormancy. This C starvation transcriptional response, activated shortly after exposure to conditions leading to bud dormancy, anticipates and underlies the growth-to-dormancy transition in the three species. The C starvation syndrome entails a suite of interconnected transcriptional responses that include sugar signaling, sugar metabolism reprogramming, senescence, autophagy, catabolism, and ABA and ethylene signaling. It also involves downregulation of cytokinin (CK) signaling, inhibition of anabolism, and repression of protein/DNA synthesis and cell division, conditions typical of cells in dormant buds. This conserved starvation response, genetically connected to cell growth arrest, may be one of the underlying forces driving the growth-to-dormancy transition of axillary buds in response to suboptimal conditions in herbaceous and woody species.

Arabidopsis Bud Dormancy Is Associated With the Induction of Four GRN
Three independent transcriptomic analyses have compared active and dormant buds in Arabidopsis. One study compared the transcriptional profiling of (dormant) buds of intact plants and of (active) buds of decapitated plants at 24 h post-treatment (Tatematsu et al., 2005). Two additional experiments compared the transcripts of active vs. dormant buds of plants exposed to high red-to-far-red light ratio (R:FR, active buds) or low R:FR (dormant buds) (González-Grandío et al., 2013;Reddy et al., 2013). Here we define dormancy as a state in which bud growth is reversibly interrupted, regardless of the requirements to resume development. A search for genes upregulated in dormant buds relative to active buds in the three experiments identified 78 genes termed bud dormancy genes (Supplementary Figure S1 and Dataset S1; González-Grandío and Cubas, 2014). These genes correspond to the least common denominator of the three studies and were induced in para-and ecodormant buds by either correlative inhibition or low R:FR, respectively. They were also differentially expressed at 3 h (Reddy et al., 2013), 8 h (González-Grandío et al., 2013 and 24 h (Tatematsu et al., 2005) after treatment onset.
We evaluated the degree of coregulation of these genes using the most updated co-expression database of ATTED-II (15,275 microarray experiments; Obayashi et al., 2007). Hierarchical clustering analysis revealed four clusters of coregulated genes (14, 20, 13, and 31 genes; Figure 1A and Supplementary Dataset S1). We then searched for additional genes coregulated with each cluster using CoExSearch (ATTED-II, Obayashi et al., 2007) and obtained four lists of highly coregulated genes (Supplementary Dataset S1). Analysis of their fold change (FC) induction in the three active-vs.-dormant bud experiments confirmed that a significant proportion of the genes in each list were induced (FC ≥ 1.2) in dormant buds in at least one experiment ( Figure 1B and Supplementary Figure S2). We termed the gene lists that comprised the bud dormancy genes of the original clusters plus their coregulated genes (induced in dormant buds in at least one experiment, red dots in Figure 1B) bud dormancy GRNI-IV, with 297, 283, 271, and 295 genes respectively (Supplementary Dataset S1).

Bud Dormancy GRNs Are Related to Hormone Signaling, Stress, Catabolism and Starvation Response
To elucidate the biological processes in which these GRNs were involved, we searched for enrichment in gene ontology (GO) terms using the Panther Classification System (Mi et al., 2017;Supplementary Dataset S2), complemented with a MapMan bin analysis (Thimm et al., 2004; Supplementary Dataset S1). GRNI was significantly enriched in terms related to ethylene, auxin and gibberellin signaling and response; GRNII in terms related to ABA, catabolism and response to abiotic stress; GRNIII in terms related to lipid and amino acid catabolism, senescence, response to starvation and biotic stress; and GRNIV in terms related to protein ubiquitination and response to sucrose starvation.
We evaluated the degree of overlap between these GRNs by seeking common genes. GRNIII and GRNIV shared one-third of their genes; GRNII and GRNIII shared 30%, and GRNI and GRNIV had 26% genes in common ( Figure 1C, Supplementary Figure S3 and Dataset S3). This suggested that these GRNs are not strictly independent, but correspond to related aspects of the same syndrome, probably coordinated or maintained by ethylene, auxin and ABA signaling ( Figure 1D).

Bud Dormancy GRNs Are Enriched in Genes Typical of a C Starvation Response
We observed that three robust sugar starvation gene markers, GIBBERELLIN-STIMULATED ARABIDOPSIS 6 (GASA6), DORMANCY-ASSOCIATED PROTEIN-LIKE 1 (DRM1/DYL1) and DARK INDUCIBLE 6 (DIN6) (Contento et al., 2004;Price et al., 2004;Gonzali et al., 2006;Zhong et al., 2015) were members of one or several GRNs (GASA6, GRNI; DRM1, GRNI, III and IV; DIN6, GRNIII and IV; Supplementary Dataset S3). As sugar has a prominent role in the control of shoot outgrowth (Mason et al., 2014;reviewed in Barbier F.F. et al., 2015) and GRNIV was significantly enriched in terms related to sucrose starvation, we studied this response further. The C starvation syndrome, triggered under C-limiting conditions (e.g., an extended night), helps to obtain an alternative energy source and C skeletons. In Arabidopsis, it comprises a suite of interconnected events that result in changes in C balance and growth. They include reprogramming of sugar sensing, transport, signaling and metabolism, increased protein ubiquitination and degradation, amino acid and lipid catabolism, induction of ABA and ethylene signaling, and recycling of cell components via autophagy and senescence. In addition, CK signaling, ribosomal gene expression, DNA synthesis and cell division are inhibited (Contento et al., 2004;Lin and Wu, 2004;Thimm et al., 2004;Gonzali et al., 2006;Rolland et al., 2006;Rose et al., 2006). Remarkably, many of the GO terms and/or MapMan bins enriched in the four GRNs matched categories induced by C-limiting conditions ( Table 1 and Supplementary Datasets S1, S2).
To test the possibility that these GRNs correspond to a C starvation response, we compared the GRN genes with four lists of genes induced in C-limiting conditions: (i) 26 genes of a robust core of C-signaling response shared by 21 Arabidopsis accessions (Supplementary Dataset S4; Sulpice et al., 2009), (ii) 57 sugar-responsive genes, proposed upstream components of the transcriptional response to sucrose (Supplementary Dataset S4; Osuna et al., 2007), (iii) 429 dark-induced, sugarrepressed genes (Supplementary Dataset S4; Gonzali et al., 2006) and (iv) 507 genes responsive to AKIN10, a catalytic subunit of the SUCROSE-NON-FERMENTING-1-RELATED PROTEIN KINASE (SnRK1), which integrates stress and C signals to coordinate energy balance, metabolism and growth (Supplementary Dataset S4;Baena-González et al., 2007).
Frontiers in Plant Science | www.frontiersin.org The GRN to which each gene belongs to is indicated.
Genes from these sets appeared in the GRNs at a much higher frequency than expected in a random list (pval 4.5E-11 to 7.5E-215; Table 2 and Supplementary Figure S4), indicating that the bud dormancy GRNs were very highly enriched in genes typical of a C starvation response.

GSEA Analyses Confirm a C Starvation Response in Dormant Buds
We assessed this potential C starvation response by performing a Gene Set Enrichment Analysis (GSEA) using all transcribed  genes from each experiment, rather than focusing on the bud dormancy GRNs. GSEA is a statistical approach that allows identification of overrepresented gene sets among differentially up-or downregulated genes of a transcriptomic experiment (Subramanian et al., 2005). For each "active-vs.-dormant bud" experiment, we generated a ranked gene list using relative gene expression levels and False Discovery Rate (FDR) values. We then tested whether gene sets related to a potential C starvation syndrome (sugar-, darkness-and AKIN10-responsive genes, ABA, ethylene and CK markers, ribosomal genes, cell cycle and cell division genes; Supplementary Dataset S4) were found toward the top (upregulated) or the bottom (downregulated) of the ranked gene lists. Analyses 1 confirmed significant overrepresentation of C-signaling, sugar-repressed, AKIN10induced and ABA and ethylene marker genes among those upregulated in the three experiments (Figures 2A, 3). In contrast, CK markers, ribosomal genes and S-phase genes were overrepresented among the downregulated genes (Figures 2B, 3).
Other cell division markers such as M-phase genes, histones and kinesins were overrepresented only in the 8 and 24 h experiments, which suggests they are downregulated at later stages of the process (Figures 2B, 3). All these results suggest that a C starvation syndrome is induced early in the growth-to-dormancy transition in para-and ecodormant axillary buds in Arabidopsis.

Regulation of the Bud Dormancy GRNs
To find potential master regulators of the C starvation-related bud dormancy GRNs, we searched for overrepresented motifs in the gene promoters (1 kilobase upstream of the transcription start site) of each GRN using Oligo-analysis and Pattern 1 http://bioinfogp.cnb.csic.es/files/projects/tarancon_et_al_2017_supp/ assembly (Rsat; Medina-Rivera et al., 2015). In GRNI and GRNIV, tcTTATCCAc was the most-overrepresented motif (Supplementary Figure S5 and Dataset S5); it contains the sucrose-repressible element TATCCA, bound in rice by the MYB factors OsMYBS1, 2, and 3, which mediate sugar-regulated gene expression (Lu et al., 1998(Lu et al., , 2002. We looked for TFs within GRNI and IV that could bind this motif, based on DNA affinity purification sequencing (DAP-Seq) data (O'Malley et al., 2016) or chromatin immunoprecipitation sequencing (ChIP-Seq) data , and that could act as master regulators of the GRNs. The Arabidopsis OsMYBS2 ortholog MYBS2 (At5g08520), which has a role in sugar and ABA signaling , pertains to GRNI and might bind this motif (Supplementary Figure S6A and Dataset S5). In GRNIV, MYBS2 and three other MYB-related proteins, MYBH/KUA1 (At5g47390), At1g19000 and At1g74840 could bind this motif (Supplementary Figure S6A and Dataset S5). MYBH/KUA1 has a critical role in dark-induced leaf senescence (Huang et al., 2015). These sugar-regulated genes could be instrumental in coordinating gene expression in GRNI and GRNIV.
We confirmed significant enrichment of the GRNs in the target genes of these TFs by using DAP-Seq and ChIP-Seq data (O'Malley et al., 2016;Song et al., 2016); their numbers in the GRNs were significantly higher than expected in a random gene list (pval < 0.01). For instance, the number of gene targets for NAC102, RD26 and ABF4 was 6, 3.2, and 3.4 times higher, respectively; for the remainder, this value was between 1.7 and 2.7 times higher than predicted (Supplementary Figure S6B).
All these results indicate that four interrelated GRNs associated to a C starvation response are induced in para-and ecodormant Arabidopsis buds. MYB-related, bZIP and NAC TFs could have a key role in the regulation of these GRNs. A large proportion of the genes in the GRNs are rapidly repressed by sugar and upregulated by AKIN10. They are tightly coregulated with or directly involved in sugar signaling and metabolism, autophagy, senescence, catabolism of lipids and proteins, and ABA and ethylene signaling. This response is also associated with downregulation of CK signaling, protein synthesis and cell division, all conditions that lead to the cell and tissue growth arrest typical of dormant buds.

Conservation of Bud Dormancy GRNs in Arabidopsis, Poplar and Grapevine
We investigated whether the GRNs related to a C starvation syndrome identified in Arabidopsis were also induced during the growth-to-dormancy transition in buds of the woody plant species, poplar and grapevine. We studied two public transcriptomic experiments in which apical buds of poplar (Ruttink et al., 2007) or axillary buds of grapevine (Díaz-Riquelme et al., 2012) underwent dormancy. To induce dormancy, shoot apices of poplar plants grown in long days (LD, 16 h light-8 h darkness) were exposed to 1-6 weeks of short days (w SD, 8 h light-16 h darkness) (Ruttink et al., 2007). During treatment, the shoot apices developed into buds (1-3 w SD), grew adapted to dehydration and cold (3-6 w SD), and became dormant (5-6 w SD). Samples were collected weekly. Díaz-Riquelme et al. (2012) collected monthly samples of axillary buds of grapevine plants grown in natural conditions in the northern hemisphere. Grapevine axillary buds are formed between April and May; in July and August they grow, undergo flowering and develop inflorescence meristems, enter endodormancy at the end of September, and exit dormancy by the end of November. They remain ecodormant throughout December, until environmental conditions become benign around March, when they sprout (Martìnez de Toda Fernaàndez, 1991;Díaz-Riquelme et al., 2012).
We analyzed the expression patterns of the poplar and grapevine orthologs of the GRNI-IV genes. Of 838 Arabidopsis genes in these GRNs, we identified 390 poplar and 421 grapevine orthologs (Supplementary Dataset S6). In both species, we studied gene expression relative to levels in the "active bud" sample (LD in poplar, April in grapevine). In general, a large proportion of the bud dormancy gene orthologs were significantly induced at most time points in poplar and grapevine buds (Supplementary Figure S7), which supports a conservation, during the growth-to-dormancy transition in these woody species, of the responses found in Arabidopsis. In poplar, the global induction appeared to increase over the weeks in SD, especially for GRNII and III genes. In contrast, grapevine gene induction was detectable throughout the year (Supplementary Figure S7).
In summary, a large proportion of the genes orthologous to Arabidopsis bud dormancy genes are also induced, either early and transiently or early and constantly during the growth-todormancy transition in poplar and grapevine, which supports their functional conservation in these woody species.

The C Starvation Response Is Conserved in Poplar and Grapevine Buds
To obtain a general view of the transcriptomic responses in these experiments, we performed GSEA similar to that for Arabidopsis, using all genes with proposed Arabidopsis orthologs (8023 genes in poplar and 8390 in grapevine) (Ruttink et al., 2007;Díaz-Riquelme et al., 2012).
The sugar-and AKIN10-responsive gene sets were overrepresented among upregulated genes from 1 w SD in poplar, and July in grapevine, and were also induced throughout the treatment/year (Figure 3). This finding confirms that the C starvation response begins early, long before endodormancy onset, and underlies the entire process. The ribosomal gene set was constitutively overrepresented among downregulated genes in all three species, which confirmed that inhibition of protein synthesis is an early and sustained response in buds entering dormancy. General downregulation of cell cycle and cell division genes was also observed in grapevine, whereas in poplar, cell division gene sets were repressed more gradually and reached maximum repression at 5 w SD. In contrast to Arabidopsis, histones were not significantly downregulated in the woody species. Nevertheless, C starvation response gene sets (upregulated) and cell growth-related gene sets (downregulated) clustered together in the three species.
Unlike the gene sets discussed above, hormone responses did not appear to be strongly conserved among species, suggesting more relaxed evolution of these pathways. Whereas ABArelated genes were induced constitutively in grapevine, ABA and ethylene responses were induced from 2 w SD onward in poplar, in accordance with previous observations (Ruttink et al., 2007). CK signaling is repressed in Arabidopsis, but not notably in poplar or grapevine. An early, extended response to the senescence-associated hormone jasmonate (MJ) was repressed in two Arabidopsis experiments, and was induced in most poplar and grapevine samples (Figure 3).
This results indicate that an early and sustained sugarstarvation response associated with downregulation of ribosomal and cell cycle proteins is conserved in buds of Arabidopsis, poplar and grapevine, and might constitute a core response of buds entering dormancy in the angiosperms.

Cell Type-Specific Gene Expression of Bud Dormancy Genes in the Shoot Apex
To further analyze the function of the genes induced during the C starvation response in buds, we selected those most highly expressed in Arabidopsis, poplar and grapevine, to determine the cell types in which they are expressed. We used a highresolution gene expression database of the Arabidopsis shoot apex, which contains the same tissues as axillary buds: meristem and leaf primordia. This database comprises gene expression profiles of different cell populations obtained by fluorescenceactivated cell sorting (L1, L2 and L3 layers, central (CZ) and peripheral zone (PZ), leaf primordia, xylem and phloem) (Yadav et al., 2014). As we cannot rule out that the expression levels of these genes change in dormant axillary buds, we used this database for qualitative rather than quantitative analysis, to identify the cell types in which these genes were expressed most abundantly.
In summary, whereas ABA signaling occurs mostly in the xylem, sugar signaling in the phloem, and ethylene in the meristem proper (Yadav et al., 2014), autophagy and arrest of cell growth take place throughout the meristem. This suggests that cell-to-cell communication and movement of signaling molecules, hormones and proteins must take place across different cell types in buds entering dormancy.

Bud Dormancy Early Markers
It is of great interest to identify robust, universal markers that allow diagnosis of axillary bud status. These markers should be induced early, and their expression be sustained in para-, eco-, and endodormant buds throughout angiosperms. Based on our analysis, several genes met these criteria. We tested further four of them: DRM1, HIS1-3, GID1C, and NAP. DRM1 and HIS1-3 were expressed at high levels in the Arabidopsis, poplar and grapevine experiments. DRM1 is a well-known dormancy marker in herbaceous and woody plants species (Stafstrom et al., 1998;Park and Han, 2003;Tatematsu et al., 2005;Aguilar-Martínez et al., 2007;Kebrom et al., 2010;Wood et al., 2013). It is repressed by sugar, which supports its strong association to low sugar levels and dormancy. ABA-responsive HIS1-3 is upregulated before ABA signaling in poplar and grapevine buds. GID1C, which encodes a gibberellin receptor, is expressed at high levels in the three Arabidopsis experiments and in poplar ( Figure 4A). The senescence-promoting gene NAP is expressed at very high levels in Arabidopsis and throughout the year in grapevine ( Figure 4B). Both GID1C and NAP belong to the four bud dormancy GRNs (Supplementary Dataset S1).
We tested whether expression of these genes also correlated with bud dormancy in axillary buds of potato (Solanum tuberosum, Solanaceae, Asteridae), a species only distantly related to Arabidopsis, poplar (both Rosidae) and grapevine (basal angiosperms). We identified the potato ortholog genes for the four candidates. In the case of NAP, we found two potato paralogs (NAPa and NAPb). We studied their mRNA levels in buds of plants treated for 10 h with white light (W) or with W supplemented with far red light (W+FR), a treatment that promotes axillary bud dormancy in potato. We also compared mRNA levels in buds of intact and decapitated plants. Whereas DRM1 and HIS1-3 were confirmed as reliable markers of bud dormancy in potato, GID1C, NAPa and NAPb did not respond as anticipated in decapitated plants; GID1C was not upregulated in low R:FR and NAPa/NAPb were highly induced after decapitation (Figure 6).

Carbon Availability, a Key Signal for Growth
In all eukaryotes, cell proliferation and growth demand high carbohydrate levels for energy generation and macromolecule synthesis. Correspondingly, low C availability promotes a reduction in growth rate in order to retain sufficient C to support essential maintenance functions (Rolland et al., 2006). In Arabidopsis, sugar availability has a great influence on growth and development, both in seedlings and adult plants. C-limiting conditions (e.g., sucrose depletion, night extensions, short-day photoperiods, starchless mutants) trigger a suite of transcriptional responses that lead to growth cessation, and which include repression of genes involved in anabolism, protein synthesis, cell division, cell cycle, and DNA synthesis and repair (Moore et al., 2003;Smith and Stitt, 2007;Wiese et al., 2007).
Regarding shoot branching, it has been shown that sugar availability to buds plays a major role in its control in pea and rose (Mason et al., 2014;. In agreement, we have found that induction of GRNs typical of tissues undergoing C starvation precedes and underlies the bud growthto-dormancy transition in Arabidopsis, poplar and grapevine. This is concomitant with transcriptional repression of ribosomal and cell-cycle genes, responses typical of buds entering dormancy as well as of tissues undergoing C limitations (e.g., Thimm et al., 2004;Smith and Stitt, 2007;González-Grandío et al., 2013). Indeed it is possible that bud dormancy is a manifestation and a consequence of the observed C starvation syndrome.

Dormancy-Promoting Stimuli and the C Starvation Response
How is this C starvation response induced? In apical dominance it has been proposed that the growing shoot apex acting as a sugar sink might limit sugar availability to axillary buds so this can be the direct trigger of the response (Mason et al., 2014). Tre6P, a metabolite that acts as a proxy for C status, may also promote signaling in addition to, or instead of direct sugar sensing (Paul   Lunn et al., 2014). This would be in agreement with the observation that plants that express microbial trehalosephosphate synthase (TPS) genes show increased shoot branching (Goddijn et al., 1997) and maize plants with a mutation in the trehalose-phosphate phosphatase gene RAMOSA3 have altered inflorescence branching (Satoh-Nagasawa et al., 2006). Nevertheless, it is likely that the syndrome is not only induced by an actual sugar shortfall, but also by cues that inform of current or future suboptimal conditions which may affect energy availability and/or interfere with respiration and C assimilation (Baena-González et al., 2007;Baena-González and Sheen, 2008). Seasonal environmental changes that perturb these processes (e.g., daylength shortening, light levels, temperature, water availability) may trigger acclimatory signaling pathways that anticipate C limitations (Smith and Stitt, 2007). Those and other stimuli could feed into regulatory networks that economize resources locally, to result in a moderation of growth rate in axillary meristems and buds.
In two of the Arabidopsis experiments examined, the dormancy-inducing stimuli was an exposure to low R:FR light ratio. Low R:FR light is interpreted by plants as a situation with limited light available for photosynthesis. It severely reduces the expression of photosynthesis-related genes (Cagnola et al., 2012) and induces cell-wall remodeling in stem and petioles, which may divert carbohydrates away from axillary buds (Sasidharan et al., 2010). Furthermore, low R:FR promotes ethylene and ABA signaling and CK degradation (Carabelli et al., 2007;Cagnola et al., 2012), hormonal responses tightly linked to the C starvation response (see below).
In the poplar and grapevine studies, the sugar-repressed networks are induced in buds soon after beginning of daylength shortening: in poplar at 1 w SD; in grapevine, in July, when daylength shortening has just begun (June 21), even though buds are still growing. Short-day photoperiod leads to localized flower and seed abortion associated with low levels of C in Arabidopsis (Lauxmann et al., 2016). Likewise, in poplar a measurable shortage of sugar availability is detectable after 1 w SD (Ruttink et al., 2007). Under the natural conditions in which grapevine plants are grown, daylength shortening and C limitations are progressive, but relatively small changes in C balance may trigger the response. Indeed in Arabidopsis even minor alterations in C status, well before C starvation, lead to notable changes in C-related signaling and response (Usadel et al., 2008). In addition, genetic pathways that sense photoperiod might help anticipate and adapt to impending C-limiting conditions in short days. These pathways, controlled by phytochromes, circadian clock, and genes controlling flowering time (Horvath, 2009), may regulate and establish crosstalk with the C starvation response. Indeed, sugars affect the expression of clock genes (Haydon et al., 2013) and conversely, the clock regulates carbohydrate metabolism (Smith and Stitt, 2007). Phytochromes, which monitor changes in R:FR and in day-length, also regulate SD-induced endodormancy in woody species (Johnson et al., 1994;Reed et al., 1994;Olsen et al., 1997;Neff and Chory, 1998;Monte et al., 2003;Ruonala et al., 2008;Franklin and Quail, 2010). Changes in low R:FR light ratio or photoperiod might therefore trigger partially overlapping responses, including potential anticipation of a C-limiting situation.
Although it has not been analyzed in this work, coordination between C and N metabolic pathways probably affect this process as well, as sugar responses depend significantly on the N status of the plant.

The C Starvation Syndrome in Axillary Buds: Sugar Signaling
The C starvation syndrome comprises a cascade of transcriptomic events that culminate in changes in growth and C balance (Figure 7). These events include induction of genes involved in transcriptional regulation, sugar sensing, transport and signaling, catabolism (i.e., amino acid and lipid degradation), protein ubiquitination and degradation, hormone signaling, autophagy and senescence. Genes required for growth, such as ribosomal, cell cycle and anabolism-related genes become downregulated (Thimm et al., 2004). Buds entering dormancy in Arabidopsis, poplar and grapevine show induction of genes of the former categories and repression of genes of the latter categories.
EXORDIUM-like (EXL)2 and EXL4 are bud dormancy genes potentially involved in sugar sensing. They are induced in extended night treatments in Arabidopsis seedlings, in accordance with a role under C-limiting conditions (Schröder et al., 2012). Their close paralogs, EXORDIUM (EXO) and EXL1, are proposed to integrate apoplastic C status with intracellular responses (EXO) (Lisso et al., 2013) and to control primary and long-term adaptation to C starvation (EXL1) (Schröder et al., 2011(Schröder et al., , 2012. Several sugar transporters are also induced in dormant buds. These are STP1, one of the most rapidly and prominently downregulated genes in response to sugars (Price et al., 2004;Cordoba et al., 2015), STP14, which is strongly repressed by sugars (Büttner, 2010) and the sucrose efflux transporters SWEET11 and SWEET12, which act with the sucrose/proton symporter SUT1/SUC2 for phloem loading and long-distance transport (Chen et al., 2012).
Sugar signaling involves Tre6P (Paul et al., 2008;Lunn et al., 2014) and four class-II TPS, TPS8, TPS9, TPS10, and TPS11, are induced in buds entering dormancy. They belong to a core C-signaling response, are usually strongly upregulated in C starvation, and are AKIN10-responsive (Contento et al., 2004;Price et al., 2004;Thimm et al., 2004;Baena-González et al., 2007). Although their proteins may be catalytically inactive, they might modulate other TPSs or act as Tre6-P sensors (Lunn, 2007). In addition, the SnRK1 protein-kinase, a central regulator of growth in response to C availability (Baena-González et al., 2007), is likely to have a key role in the induction of bud dormancy (see below).
One role of ABA and ethylene in the C starvation syndrome is induction of senescence, a genetically programmed process that promotes degradation of cell components and macromolecules, remobilizes nutrients, and optimizes resources to supply energy and C skeletons. Ethylene and ABA activate senescence-related genes and senescence induces ABA signaling (Abeles et al., 1988;Zeevaart and Creelman, 1988;Zacarias and Reid, 1990;Reid and Wu, 1992;Weaver et al., 1998;Seo et al., 2000;Yang et al., 2003;Buchanan-Wollaston et al., 2005;Lim et al., 2007). It is noteworthy that the potential master regulators of GRNII and III, ATAF1, ORE1/NAC6 and NAP, are ABA-induced factors that control senescence. ATAF1 induces a C starvation transcriptome and ABA biosynthesis (Jensen et al., 2013;Garapati et al., 2015). NAP activates SAG113/HAI1 and controls expression of ABSCISIC ALDEHYDE OXIDASE3 (AAO3), encoding an enzyme that catalyzes the final steps of ABA synthesis (Guo and Gan, 2006;Yang et al., 2014). ORE1 controls the expression of at least 78 SAGs and might also promote DNA degradation (Balazadeh et al., 2010;Matallana-Ramirez et al., 2013;Kim et al., 2014). HAT22, GBF2, GBF3, ABF3 and ABF4 are additional bud dormancy genes related both to senescence and ABA (Lin and Wu, 2004;Rivero et al., 2007;Song et al., 2016). Remarkably, MAX2/ORE9, which encodes an F-box involved directly in strigolactone perception and signaling and has a critical role in the control of shoot branching, also promotes senescence (Woo et al., 2001;Stirnberg et al., 2002). Finally, the MYB genes At1g19000 and At1g74840, proposed to be master regulators of GRNI and IV, are responsive to dark-induced senescence (Lin and Wu, 2004).
In contrast, CK signaling is antagonistic to senescence, and a reduction in CK levels is a key signal for senescence initiation in Arabidopsis (Gan and Amasino, 1995;Kim et al., 2006). Consistent with this, in Arabidopsis dormant buds CK signaling is reduced, and in other species CK levels have also been reported to be reduced relative to active buds (Turnbull et al., 1997;Dun et al., 2012;Roman et al., 2016). Consistently, four genes encoding F-box proteins that promote the ubiquitination and degradation of ARR factors [KISS ME DEADLY (KMD)1-4] are bud dormancy genes.
Autophagy is another process induced by C starvation (Izumi et al., 2013) and whose markers (ATG genes) are upregulated in buds entering dormancy. This is a process by which cytoplasmic components and organelles are transported to the vacuole, where they are broken down and recycled. Under C-limiting conditions it contributes to plant energy availability (Aubert et al., 1996;Rose et al., 2006;Izumi et al., 2013). Autophagy is associated with induction of lipid degradation and upregulation of E2-and E3-ubiquitin ligase components, which promote proteasomal-dependent protein degradation (Thompson and Vierstra, 2005). We have found a remarkable number of bud dormancy genes related to autophagy, ubiquitination, protein degradation and lipid catabolism, many of them controlled by SnRK1 (see below).

SnRK1
Could Have a Pivotal Role during the Bud Growth-to-Dormancy Transition SnRK1, a protein-kinase active in low energy conditions, promotes catabolism and represses anabolism, cell division and growth. Our transcriptomic data indicates that it may play an important role during the bud growth-to-dormancy transition. SnRK1 affects expression of robust dormancy markers such HIS1.3 and DRM1, and the potential master regulator of GRNI and GRNIV, MYBH/KUA1. In buds entering dormancy, the SnRK1 β subunit AKINBETA1, whose mRNA levels correlate directly with night duration (Pokhilko et al., 2014), is induced. Most importantly, our GSEA analysis indicates that the transcriptional network downstream of the catalytic SnRK1 α subunit, AKIN10, is significantly induced from the earliest stages of growth-to-dormancy transition in Arabidopsis, poplar and grapevine buds, and is maintained in para-, ecoand endodormant buds. Many of the abovementioned genes involved in sugar sensing, signaling, autophagy and repression of CK signaling are AKIN10-dependent, including EXL4, STP1/14, SWEET11/12, TPS8/9/19/11, AKINBETA, ATG8E/F/G/H, ATG18F/G, and F-box genes KMD1, 3, 4. SnRK1 also causes downregulation of a large number of ribosomal genes, another conserved significant effect detected by our GSEA analysis. SnRK1 could also be responsible for at least part of the observed induction of the ubiquitination machinery and lipid degradation.

A Conserved Core C Starvation Response Underlies Bud Dormancy in Angiosperms
Bud dormancy is an adaptive response present in all angiosperms. It prevents shoot development when endogenous or environmental conditions are unfavorable for sustained growth. It has great impact on reproductive success, productivity and survival, and must have been influential in the colonization of habitats with fluctuating conditions.
We have found induction of a conserved C starvation syndrome that precedes and underlies the growth-to-dormancy transition in buds of three distantly-related species, one herbaceous (Arabidopsis) and two woody (poplar and grapevine). This transcriptional response, composed by several interconnected GRNs, comprises ortholog genes in Arabidopsis, poplar and grapevine, as gene sets generated in Arabidopsis were used to detect the response in the woody species. Furthermore, this syndrome has been observed is several unrelated experiments, regardless the stimulus that promoted dormancy, either environmental (low R:FR, shortday photoperiods) or endogenous (apical dominance). This remarkable conservation suggests that a syndrome aimed at adapting to C-limiting situations is deeply rooted in the control of shoot meristem and bud development across angiosperms. Bud dormancy might thus be an ancestral response directly resulting from this C starvation syndrome, coordinated by different pathways that sense and/or anticipate situations on low C availability and feed into this core response to prevent untimely growth and development.

Functional Annotation of Bud Dormancy GRNs
Automated function prediction for the GRNs was carried out using GO analyses. The PANTHER classification system (Mi et al., 2017) was used to identify overrepresented biological process ontologies using a statistical overrepresentation test followed by Bonferroni correction for multiple testing. TAIR10 version of Arabidopsis thaliana genome was used as reference. We selected ontologies with a pval < 0.05. In addition, Mapman bins (Thimm et al., 2004) were added to all the genes in Supplementary Dataset S1.

Gene Set Enrichment Analysis
Gene Set Enrichment Analysis (Subramanian et al., 2005) was used to identify gene sets whose genes are overrepresented in different conditions. The GSEA method evaluates whether these genes occur preferentially toward the top or bottom of a ranked list. Enrichment scores are calculated using "weighted" statistics. For each sample, we calculated the log2 ratios of normalized gene intensities vs. normalized gene intensities of the "active bud" sample: white light-treated buds for González-Grandío et al. (2013); High R:FR-treated n-2 buds for Reddy et al. (2013); Buds of decapitated plants for Tatematsu et al. (2005); Buds of LD-grown poplars for Ruttink et al. (2007), April grapevine buds for Díaz-Riquelme et al. (2012). Genes were ranked by their log2 ratios calculated as the difference between normalized log intensity in the "dormant bud" condition minus normalized log2 intensity in the "active bud" condition. Intensity expression values were obtained from the references above. Gene sets for hormone markers were obtained from Nemhauser et al. (2006) and Mashiguchi et al. (2009). Gene sets related to sugar and AKIN10 responses were obtained from Sulpice et al. (2009) (Core C-signaling), Osuna et al. (2007) (Sugar-responsive), Gonzali et al. (2006) (Dark-induced, sugar-repressed), and Baena-González et al. (2007) (AKIN10responsive). The other gene sets are from González-Grandío et al. (2013). The GSEA Normalized Enrichment Score for all gene sets in all comparisons were clustered with TM4 Multi Experiment Viewer (MeV, Saeed et al., 2003). Tree was generated by the hierarchical clustering method (HCL) using Euclidean distance and average linkage options. Complete results are in http:// bioinfogp.cnb.csic.es/files/projects/tarancon_et_al_2017_supp/.

Promoter Motif Analysis
Sequences (1 kb) 5 of the transcription start site of the bud dormancy GRN genes were retrieved with Sequence Bulk Download 2 . Overrepresented 6-8mer motifs were identified with Motif discovery (RSAT, Medina-Rivera et al., 2015). The oligo-analysis tool was used to find significantly overrepresented motifs, which were assembled into frequency matrices with pattern-assembly and default parameters. Matrices were converted into consensus motifs with convert-matrix and represented using WebLogo (Crooks et al., 2004).

Generation and Visualization of Poplar and Grapevine Expression Datasets
For each time point we calculated the log2 ratios of normalized gene intensities vs. normalized gene intensities on LD (active buds). Expression data was visualized and clustered with MeV. Tree was generated by HCL, using Euclidean distance and average linkage options.

Cell-Type Specific Shoot Apex Expression of Bud Dormancy GRN Genes
For each sample, we calculated the log2 ratios of normalized gene intensities vs. normalized gene intensities of the "active bud" sample: LD for poplar, April for grapevine. Expression data for selected bud dormancy genes obtained from Yadav et al. (2014) was visualized and clustered with MeV. Trees were generated by HCL using Euclidean distance and average linkage options.

Identification of Solanum tuberosum Orthologs
The putative orthologs of Arabidopsis genes were identified by a tblastn search with protein sequences as query in the Spud DB Potato Genomics Resource website 3 . cDNAs showing a high similarity e-value with the query were selected. Proteins were aligned with those of Arabidopsis and phylogenetic trees (BioNeighbor joining method, 500 replicates; Gascuel, 1997) were built to identify the most likely orthologs, which were selected for expression studies (Supplementary Figure S9).

Quantitative-PCR Expression Analyses in Solanum tuberosum
Plant growth conditions, experimental design, light treatments, techniques and expression level normalization were as described in . For each biological replicate, 8 axillary buds from node 2, and 8 from node 3 were dissected (node 1 = lowest plant node); 4-5 biological replicates were collected for each condition. Primers used are listed in Suplemmentary  Table S1.

AUTHOR CONTRIBUTIONS
CT, EG-G, JO, and MN, performed experiments. PC, performed experiments, and wrote the manuscript.

ACKNOWLEDGMENTS
We thank Desmond Bradley and Elena Baena for constructive criticisms of the manuscript and Catherine Mark for editorial assistance. PC is supported by a MINECO grant (BIO2014-57011-R). CT is a La Caixa predoctoral fellow. EG-G was a predoctoral fellow of Fundación Ramón Areces and a CSIC JAE-Predoc fellow. MN is a Excelence Severo Ochoa (MINECO) postdoctoral researcher.