# Distinctive Gene Expression Patterns Define Endodormancy to Ecodormancy Transition in Apricot and Peach

Jiali Yu, Anna O Conrad, Véronique Decroocq, Tetyana Zhebentyayeva, Daniel E Williams, Dennis Bennett, Guillaume Roch, Jean-Marc Audergon, Christopher Dardick, Zongrang Liu, et al.

► To cite this version:

Jiali Yu, Anna O Conrad, Véronique Decroocq, Tetyana Zhebentyayeva, Daniel E Williams, et al. Distinctive Gene Expression Patterns Define Endodormancy to Ecodormancy Transition in Apricot and Peach. Frontiers in Plant Science, 2020, 11, pp.180. 10.3389/fpls.2020.00180. hal-02534855

HAL Id: hal-02534855

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

Submitted on 7 Apr 2020

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

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

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

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

Check for
updates

# Distinctive Gene Expression Patterns Define Endodormancy to Ecodormancy Transition in Apricot and Peach

Jiali Yu¹, Anna O. Conrad²,³, Véronique Decroocq⁴, Tetyana Zhebentyayeva⁵,
Daniel E. Williams⁶, Dennis Bennett⁷, Guillaume Roch⁸, Jean-Marc Audergon⁸,
Christopher Dardick⁷, Zongrang Liu⁷, Albert G. Abbott² and Margaret E. Staton¹,⁹*

OPEN ACCESS

**Edited by:**
Michael James Considine,
University of Western Australia,
Australia

**Reviewed by:**
Hisayo Yamane,
Kyoto University, Japan
Gabino Rios,
Instituto Valenciano de Investigaciones
Agrarias, Spain
Carmen Leida,
Fondazione Edmund Mach, Italy

*Correspondence:
Margaret E. Staton
mstaton1@utk.edu

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

**Received:** 05 November 2019

**Accepted:** 06 February 2020

**Published:** 28 February 2020

**Citation:**
Yu J, Conrad AO, Decroocq V, Zhebentyayeva T, Williams DE, Bennett D, Roch G, Audergon J-M, Dardick C, Liu Z, Abbott AG and Staton ME (2020) Distinctive Gene Expression Patterns Define Endodormancy to Ecodormancy Transition in Apricot and Peach. *Front. Plant Sci.* 11:180. doi: 10.3389/fpls.2020.00180

¹ Genome Science and Technology Program, University of Tennessee, Knoxville, TN, United States, ² Forest Health Research and Education Center, University of Kentucky, Lexington, KY, United States, ³ Department of Plant Pathology, The Ohio State University, Columbus, OH, United States, ⁴ UMR 1332 Biologie du Fruit et Pathologie, Equipe de Virologie, INRA, Université de Bordeaux, Villenave d'Ornon, France, ⁵ Department of Ecosystem Science and Management, Schatz Center for Tree Molecular Genetics, the Pennsylvania State University, University Park, PA, United States, ⁶ Center for Environmental Biotechnology, University of Tennessee, Knoxville, TN, United States, ⁷ Appalachian Fruit Research Station, United States Department of Agriculture—Agriculture Research Service, Kearneysville, WV, United States, ⁸ GAFL Fruit and Vegetable Genetics and Breeding, INRA Centre PACA, Montfavet, France, ⁹ Department of Entomology and Plant Pathology, Institute of Agriculture, University of Tennessee, Knoxville, TN, United States

Dormancy is a physiological state that plants enter for winter hardiness. Environmental-induced dormancy onset and release in temperate perennials coordinate growth cessation and resumption, but how the entire process, especially chilling-dependent dormancy release and flowering, is regulated remains largely unclear. We utilized the transcriptome profiles of floral buds from fall to spring in apricot (*Prunus armeniaca*) genotypes with contrasting bloom dates and peach (*Prunus persica*) genotypes with contrasting chilling requirements (CR) to explore the genetic regulation of bud dormancy. We identified distinct gene expression programming patterns in endodormancy and ecodormancy that reproducibly occur between different genotypes and species. During the transition from endo- to eco-dormancy, 1,367 and 2,102 genes changed in expression in apricot and peach, respectively. Over 600 differentially expressed genes were shared in peach and apricot, including three DORMANCY ASSOCIATED MADS-box (DAM) genes (*DAM4*, *DAM5*, and *DAM6*). Of the shared genes, 99 are located within peach CR quantitative trait loci, suggesting these genes as candidates for dormancy regulation. Co-expression and functional analyses revealed that distinctive metabolic processes distinguish dormancy stages, with genes expressed during endodormancy involved in chromatin remodeling and reproduction, while the genes induced at ecodormancy were mainly related to pollen development and cell wall biosynthesis. Gene expression analyses between two *Prunus* species highlighted the conserved
transcriptional control of physiological activities in endodormancy and ecodormancy and
revealed genes that may be involved in the transition between the two stages.

Keywords: dormancy, *Prunus*, floral buds, transcriptome, chill requirement, bloom date, co-expression network, RNASeq

INTRODUCTION

In temperate regions, many perennial plants protect their buds and reproductive tissues from winter cold by entering a state of dormancy. To transit from dormancy to bud-burst, these trees must be exposed to a particular period of chilling temperatures and warm temperatures, referred to as chilling requirements (CR) and heat requirements (HR) respectively. The CR for bud break is an intrinsic part of temperate tree phenology that varies between and within species based on adaptations to the climate of their native range (Luedeling, 2012). Various models exist for calculating CR, but they are usually measured as the number of accumulated hours in a specific range, around 0-7°C (Richardson, 1974; Erez, 1979; Shaltout, 1983; Hänninen, 1990; Cesaraccio et al., 2004). For agricultural tree crops, cultivars must be carefully selected for each growing region in order to provide a sufficient amount of cold for regular flowering and reliable fruit or nut production (Luedeling and Brown, 2011). With global climate change, understanding dormancy and its control are critical for the productivity of fruit trees, as increasing or irregular seasonal temperatures impact optimal timing for important phenological traits (Hatfield and Prueger, 2015). Lack of fulfillment of CR due to warm winters affects bud break, resulting in low flowering rates and thus low fruit yield as the fruit trees do not fulfill their CR (Viti and Monteleone, 1991). In addition, early blooming trees can suffer significant yield loss from frost damage to floral blooms during unusually late periods of freezing in the spring (Rieger, 1989).

Winter dormancy is divided into two stages, endodormancy and ecodormancy (Lang et al., 1987). Floral bud endodormancy is a state defined by physiological inhibition of flowering and is induced by cold temperatures and/or short photoperiods (Horvath et al., 2003; Foley et al., 2009). Fulfillment of CR marks the transition point from endodormancy to ecodormancy. In ecodormancy, buds regain competency to respond to external environmental factors. They remain ecodormant under unfavorable growth conditions (e.g. cold temperatures), but quickly progress toward budburst when favorable conditions are present. Plants that do not receive sufficient chilling fail to transition to ecodormancy, which leads to failing to bloom or blooming erratically in the spring (Luedeling et al., 2011). The transition from endodormancy to ecodormancy is irreversible, however, there is not a phenotype which uniquely identifies the physiological state of endodormancy other than measurements of chill accumulation and more recently, studies of starch accumulation in sweet cherry floral buds (Anderson et al., 1985; Citadin et al., 2001; Chavarria et al., 2009; Fadón et al., 2018a). More studies linking floral bud endodormancy with specific physiological networks

are necessary to determine what drives the transition from endo-
to eco-dormancy and how this transition is regulated at both the
genetic and physiological levels.

Prunus spp. have genotypes exhibiting CRs ranging from approximately 100 hours to over 1,000 hours (Audergon, 1993; Caruso et al., 1997; Valentini et al., 2002; Alburquerque et al., 2008; Dirlewanger et al., 2012; Quero-Garcia et al., 2016). Genotypes with a broad phenotypic range for this trait provide excellent materials for the study of genetic control of CR and molecular activities at different stages of dormancy. Previous genetic studies on dormancy associated traits suggest that CR is a complex quantitative trait and has a major effect on flowering time in *Prunus* spp. including almond (*P. dulcis*) (Sánchez-Pérez et al., 2012), apricot (*P. armeniaca*) (Olukolu et al., 2009), peach (*P. persica*) (Fan et al., 2010; Zhebentyayeva et al., 2014), and sweet cherry (*P. avium*) (Castède et al., 2014). Some major quantitative trait loci (QTL) for CR appear to overlap orthologous genomic regions across species, suggesting shared underlying molecular mechanisms (Dirlewanger et al., 2012). For example, the most significant QTL responsible for CR and flowering time was identified on LG4 in sweet cherry (Castède et al., 2014) and almond (Sánchez-Pérez et al., 2012). Despite this extensive research, the only QTL region for CR and blooming date (BD) with strongly supported candidate genes is a major QTL on LG1 from peach (Fan et al., 2010; Zhebentyayeva et al., 2014). The strong association between this QTL region on chromosome 1 and CR was also identified from a genome-wide association study on over 400 peach genotypes (Li et al., 2019). This region contains six tandemly repeated DORMANCY ASSOCIATED MADS-box (*DAM1-6*) transcription factors, previously mapped to the *evg* (*evergrowing*) locus (Wang et al., 2002; Bielenberg et al., 2004). The deletion of four of the six *DAM* genes is associated with the *evergrowing* phenotype in peach, which exhibits an inability to cease growth in winter (Bielenberg et al., 2008). Gene profiling reveals *DAM1-4* peak in expression during bud set and may influence dormancy onset (Li et al., 2009; Yamane et al., 2011) while levels of *DAM5* and *DAM6* are high at the beginning of endodormancy and decreased steadily during chilling period (Li et al., 2009; Jiménez et al., 2010b; Yamane et al., 2011; Leida et al., 2012a). This pattern suggests that these genes likely play a significant role in maintaining the endodormant state. Recent studies of *DAM* genes and their homologs *SHORT VEGETATIVE PHASE-like* (*SVL*) genes suggest they are cooperatively regulated by several transcription factors, hormones, and epigenetic factors (da Silveira Falavigna et al., 2018). Wang et al. (2019) propose peach *TEOSINTE BRANCHED1* transcription factor negatively regulates *DAM5* and *DAM6* expression, resulting in dormancy release (Wang et al., 2019). Leida et al. (2012b) reported that
histone modifications in *DAM6* promoter, the second exon and the second intron are involved in dormancy release (Leida et al., 2012b). Sequencing of the *DAM* gene region in selected members of a QTL population of peach revealed that large intronic insertions in *DAM5* and *DAM6* are associated with the low chill phenotype, and Zhebentyayeva et al. (2014) hypothesize these insertions change the epigenetic factors of the locus to influence expression (Zhebentyayeva et al., 2014).

In *Prunus*, winter buds were originally considered to be in a resting state and thus this state was termed “dormancy.” However, both endodormant and ecodormant *Prunus* buds are now known to have unique and changing transcriptional profiles, continued metabolic activities including starch accumulation and hormone fluctuations, and ongoing floral structure differentiation (Reinoso et al., 2002; Zhu et al., 2015; Chmielewski et al., 2017; Zhang et al., 2018). This differs from other species such as pear, which have a true resting state at endodormancy (Saito et al., 2015b). Previous studies have profiled expression for a subset of genes during bud dormancy in peach (*Prunus persica*) using suppression subtractive hybridization (Leida et al., 2010) and in apricot (*Prunus armeniaca*) using cDNA amplified fragment length polymorphism (Čechová et al., 2012). However, no high throughput global gene expression studies are available in either species. In Japanese apricot (*Prunus mume*), a more recent transcriptome study has profiled global gene expression patterns during bud dormancy (Zhang et al., 2018). Transcriptomic and metabolic profiles indicate that gibberellins (GA) and abscisic acid (ABA) likely play a role in bud dormancy (Zhang et al., 2018). Based on their data and previous results from others, Zhang et al. (2018) proposed a molecular model of dormancy control: cold temperatures induce c-repeat binding factor (*CBF*) transcripts, which promote the expression of *DAM* genes and inhibit the GA signaling pathway, resulting in endodormancy. The long-term cold period then reduces *CBF* and *DAM* gene expression, leading to dormancy release (Zhang et al., 2018). However, Japanese apricot floral timing, unlike other *Prunus* species, are more strongly correlated to HR rather than CR (Kitamura et al., 2017). Research from a wide variety of tree species has yielded other proposed mechanisms for dormancy regulation, including epigenetic regulation of key gene networks through DNA and histone methylation (Lloret et al., 2018). Studies in poplar also support that hormones such as ABA participate in the establishment and release of dormancy (Rinne et al., 2011; Singh et al., 2018; Tylewicz et al., 2018). Low temperature and short photoperiod induce poplar bud dormancy by upregulating ABA signaling, which increases poplar *SVL* expression (Singh et al., 2018; Singh et al., 2019). *SVL* is likely involved in dormancy initiation and bud break (Busov, 2019). ABA has been reported to induce dormancy by increasing the frequency of plasmodesmata closure, which inhibits growth signal transport to the apical meristem (Tylewicz et al., 2018). Overexpression of *Prunus mume* *DAM6* in apple delays bud break by increasing ABA and decreasing cytokinin levels (Yamane et al., 2019). While ethylene has mainly been characterized as a phytohormone controlling dormancy onset

in trees such as poplar and birch (Ruonala et al., 2006; Ruttink et al., 2007), it also may be important for later stages as exogenous application can delay budbreak (Liu and Sherif, 2019).

This study examined global gene expression profiles over dormancy states in the contrasting CR and BD of peach and apricot genotypes, allowing us to differentiate candidate genes responsible for endodormancy maintenance and/or endodormancy to ecodormancy transition. Species of *Prunus* including peach (*P. persica*), apricot (*P. armeniaca*), and plum (*P. domestica*) have high genomic synteny (Dirlewanger et al., 2004; Jung et al., 2009). Based on extensive genomic data, the current study proposes to characterize shared molecular mechanisms during the endodormancy to ecodormancy transition in floral bud tissues across two related *Prunus* species, peach and apricot. RNASeq analysis was performed to explore the transcriptomic changes as dormancy progressed in floral buds from four apricot genotypes with contrasting BD, along with floral buds from four peach genotypes with contrasting CR phenotypes. After integrating the peach and apricot RNASeq data, shared differential expression genes between endodormancy and ecodormancy stages were identified. Moreover, co-expression networks revealed similar biological pathways triggered by endodormancy or ecodormancy in both species. Therefore, examination of gene expression patterns in different *Prunus* species is expected to shed light on the physiological and molecular similarities of dormancy across species and potentially help to generate new candidate genes for previously identified QTLs controlling CR.

# MATERIALS AND METHODS

## Plant Materials

Five apricot genotypes grown at INRA-PACA Domaine des Pins de l'Amarine, France were used for this study: A1956 (Palsteyn, n=2), A2137 (Bakour, n=3), A1267 (Badami, n=3), A660 (Bergeron, n=3), and A2312 (Flamingold, n=1). Two late blooming genotypes (A660, A1267) and three early blooming genotypes (A1956, A2137, A2312) were selected, with the historical record of bloom dates reported in Conrad et al. (Conrad et al., 2019). While CR has not been directly measured for these genotypes, the bloom date (BD) phenotype is correlated with CR (Fan et al., 2010): early BD genotypes have low CR while late BD genotypes have high CR. Floral buds from clonally propagated (grafted) trees of each genotype were collected at 0, 100, 400, and 800 chill hours regardless of their developmental stages, starting from October 29th, 2015. Two additional collections were made at specific developmental stages, when sepals became visible (C2) and when petals became visible (D1) (Conrad et al., 2019).

Peach genotypes A209 (300 CR), A340 (300 CR), A318 (850 CR), and A323 (1,100 CR) were selected from an F₂ population derived from low CR cultivar "Fla.92-2C" by high CR cultivar "Contender" (Fan et al., 2010; Zhebentyayeva et al., 2014) located at Clemson University. Floral buds from four clones of each genotype were collected at 0, 100, 600, and 1,000 chill hours and the stage of pre-bloom, from October 14th, 2015 until
March 17th, 2016 (Table 1). No collection was made at the 1,000 chill hours time point for the low chill genotypes as the buds had already flowered.

## RNA Extraction and Sequencing

RNA was extracted from floral buds of five apricot genotypes (four as specified below for RNAseq and one genotype A2312 for qPCR) as previously described (Conrad et al., 2019). Not all samples yielded high quality RNA. Of the 24 genotype by time point collections, 14 had successful extractions of all three biological replicates, eight had successful extractions of only two biological replicates, and two had only a single biological replicate (A2137 at 0 chill hours and A1267 at 0 chill hours) (Supplementary Tables 1 and 2). The 60 RNA samples from four apricot genotypes (A1267, A660, A1956, A2137) were submitted to MOgene, LLC (St. Louis, MO, USA) for library preparation and sequencing by an Illumina NextSeq 500, yielding 50 bp single-end reads.

Peach RNA samples were prepared by grinding 1 g of flower buds to a fine powder in the presence of liquid nitrogen followed immediately by the addition of 5 ml cold Invitrogen Plant RNA Reagent (Thermo Fisher Scientific, Waltham, MA, Cat# 12322-012), aliquoted 1 ml into microcentrifuge tubes and following the manufacturer's instructions. RNA pellet was dissolved in 100 µl of RNAsecure (Thermo Fisher Scientific, Waltham, MA, Cat# AM7006). Additional removal of DNA was required by adding an equal volume Sigma Tri-Reagent (Sigma, St. Louis, MO, Cat# T9424) using manufacturer's instructions with two additional purification steps using Phenol: Chloroform and final elution in 100 µl of RNAsecure (Thermo Fisher Scientific, Waltham, MA, Cat# AM7006). One µg of RNA sample were loaded on a Low EEO/Multipurpose Agarose gel (Thermo Fisher Scientific, Waltham, MA, Cat# BP160) and normalized using a Typhoon 9600 FLS (GE Healthcare, Marlborough, MA) with ImageQuant TL Image Analysis Software version 8.1 (GE Healthcare, Marlborough, MA). All genotype by time point combinations had at least three biological replicates and most had four biological replicates (Supplementary Tables 1 and 2). Sixty-

nine RNA samples were submitted to GeneWiz (South Plainfield, NH, USA) for library preparation and sequencing on the HiSeq platform to generate 150 bp paired-end reads.

## Quality Control

RNA-seq yielded approximately 19–35 million reads per library. Raw reads were analyzed for quality by FastQC (Andrews, 2010). Adaptors and low-quality reads were trimmed by skewer (Jiang et al., 2014).

## Read Alignment and Gene Quantification

As an apricot reference genome is not available, and peach and apricot are closely related, the peach reference genome was used as a reference for both species. Initial read alignments revealed the peach RNAseq libraries contained significant portions of reads derived from rRNA. rRNA filtering was conducted before re-aligning to reference genome. Raw reads were mapped to small and large subunit rRNAs from SILVA ribosomal RNA database (Quast et al., 2013) and the mapped reads were removed. STAR (Dobin et al., 2013) was used to align the clean apricot and peach RNA-seq reads to both *Prunus persica* genomes: *Prunus persica* v1.0 (The International Peach Genome Initiative et al., 2013) and *Prunus persica* v2.0 (Verde et al., 2017). The aligned reads were analyzed by HTSeq-count (Anders, Pyl et al., 2015), which counts for the aligned reads mapped to the exons. The *Prunus persica* v2.0 was used for all analyses except for the DAM region, which is correctly annotated in the original version of the genome *Prunus persica* v1.0 (Supplementary Figure 1). For the HTSeq-count results, the DAM region in *Prunus persica* v2.0 (from Prupe.1G531100 to Prupe.1G531700, peach chromosome Pp01 from 43,417,246 to 43,480,648 bases) was removed and replaced by the read counts from the six *DAM* genes (ppa018667m, ppb017585m, ppa010758m, ppa011123m, ppa010822m, and ppa010714m) annotated in *Prunus persica* v1.0 assembly.

## Differential Expression Analysis and Functional Enrichment

Gene counts were analyzed for gene-level differential expression using DESeq2 (Love et al., 2014). Sample distances were calculated using the regularized log transformation of gene counts available from the DESeq2 package. The principal component (PC) analysis used the first and second largest sample distances as the first PC and second PC. Two Wald tests were performed. The first utilized one factor, stage of development. A second test was run with one factor, phenotype, coded as either early bloom or late bloom for each sample. Significantly differentially expressed genes (DEGs) were identified with fold change > 2 and adjusted p-value < 0.05. Genes of interest were annotated with gene ontology (GO) by AgriGO v2.0 (Tian et al., 2017) using both peach *Prunus persica* v2.0 and *Arabidopsis* TAIR10 genomes. The *Arabidopsis* genes corresponding to peach genes were obtained as the best hit of *Arabidopsis* as recorded in the annotation information file from Phytozome12 (Goodstein et al., 2012). GO enrichment analysis was conducted using the hypergeometric statistical model and p-values were adjusted by the Hochberg method. The GO terms with FDR < 0.05 were regarded as significantly enriched.

<table><caption>**TABLE 1** | Sampling time points and chilling requirements (CRs) for *Prunus persica* trees.</caption><thead><tr><th rowspan="2">Phenotype Genotype</th><th colspan="2">Low-CR</th><th colspan="2">High-CR</th></tr><tr><th>A209</th><th>A340</th><th>A318</th><th>A323</th></tr></thead><tbody><tr><td>Chill requirement<fnref n="1" /> (hours)</td><td>300</td><td>300</td><td>850</td><td>1,100</td></tr><tr><td>Time point</td><td colspan="4">Sampling date</td></tr><tr><td>0 CH<fnref n="2" /></td><td>14/16-Oct-15</td><td>14/16-Oct-15</td><td>14/16-Oct-15</td><td>14/16-Oct-15</td></tr><tr><td>100 CH</td><td>23-Nov-15</td><td>23-Nov-15</td><td>23-Nov-15</td><td>23-Nov-15</td></tr><tr><td>600 CH</td><td>28-Jan-16</td><td>28-Jan-16</td><td>28-Jan-16</td><td>28-Jan-16</td></tr><tr><td>1,000 CH<fnref n="3" /></td><td>NS</td><td>NS</td><td>27-Feb-16</td><td>27-Feb-16</td></tr><tr><td>Pre-bloom (PB)</td><td>10-Feb-16</td><td>10-Feb-16</td><td>12-Mar-16</td><td>17-Mar-16</td></tr></tbody></table>

<footnote n="1">¹CR reported in (Zhebentyayeva et al., 2014).</footnote>

<footnote n="2">²Samples were collected over a 2 d period, from October 14th to 16th, 2015.</footnote>

<footnote n="3">³NS, not sampled.</footnote>
## Differential Co-Expression Analysis and Network Enrichment Analysis

Weighted Gene Correlation Network Analysis (WGCNA) (Zhang and Horvath, 2005; Langfelder and Horvath, 2008) was used to identify modules of highly correlated genes based on the DESeq2 normalized gene expression data following the steps in the documentation for the WGCNA Package (Langfelder and Horvath, 2014). Genes with 0 normalized count were removed. The expression profiles of the remaining genes were hierarchically clustered into modules containing at least 100 genes. The expression profiles of each module were summarized using the first PC as the module eigengene (ME). Module-traits correlations were calculated by Pearson correlation coefficient (PCC) where each module was represented by its ME, and each developmental stage was represented with a numeric vector with “1” for the trait, and “0” for all the others. Genes from selected modules were annotated with the *Arabidopsis* gene as described above. The GO enrichment networks were analyzed by BiNGO in Cytoscape (Maere et al., 2005).

## Reverse Transcription Quantitative PCR (RT-qPCR)

Total RNA was extracted as previously described (Conrad et al., 2019). One RNA sample from each of three apricot genotypes (A2137, A1956, A2312) at time points 100, 400, 800, and sepal was subjected to real-time qPCR analysis using the Luna® Universal one step RT-qPCR kit (E3005S, New England Biolabs Inc, MA) with the primers listed in Table 2. Expression levels of mRNAs of interest were normalized for the 18S RNA levels. Three technical replicates were used for each sample.

## Floral Bud Dissection and Imaging

Ten flower buds were randomly selected from each of the four peach genotypes and dissected longitudinal, slightly off-center using razor blades (Electron Microscopy Sciences Stainless Steel “PTFE” coated Double Edge, PERSONNA brand.004” thick, Fisher Scientific, Waltham, MA, USA, Cat# 50-949-411). The dissected buds were imaged using Nikon SMZ1500 Stereoscopic Zoom Microscope (Nikon Instruments Inc, Melville, NY, USA) with a Pax-it! ARC Model PS-CM camera (Villa Park, IL, USA).

<table><caption>TABLE 2 | Primers used for qPCR.</caption><thead><tr><td>Target transcript</td><td>Primer sequence (5' to 3')</td></tr></thead><tbody><tr><td>18S_F</td><td>GTTACTTTTAGGACTCCGCC</td></tr><tr><td>18S_R</td><td>TTCCTTTAAGTTTCAGCCTTG</td></tr><tr><td>DAM6_F</td><td>TACTGGACCTGCGTTTGTGGAGCC</td></tr><tr><td>DAM6_R</td><td>TGTTGCAGCTGGTGGAGGTGGCAATT</td></tr><tr><td>Prupe.1G104900_F</td><td>TCATCTTCCGCTGCCTTTGTAGCCT</td></tr><tr><td>Prupe.1G104900_R</td><td>GACACTGCCAAGAACACCAAGGACA</td></tr><tr><td>Prupe.2G122600_F</td><td>GGAGAAATTGGAACGCCTGTGC</td></tr><tr><td>Prupe.2G122600_R</td><td>TGAGCCCTCAGTTGCTAGTTCAG</td></tr><tr><td>Prupe.5G014900_F</td><td>TCCCTTTGGACAGATTCCAGTGC</td></tr><tr><td>Prupe.5G014900_R</td><td>GCAGCCTCTTTCAGGTTGTTGTG</td></tr><tr><td>Prupe.7G084900_F</td><td>ACCATTTGCCGGATGGATGGAAG</td></tr><tr><td>Prupe.7G084900_R</td><td>CAACCATGTCAGCTGGAACCAC</td></tr></tbody></table>

## Multivariate Analysis of Variance

Permutational multivariate analysis of variance (PERMANOVA) (Anderson, 2001; Anderson, 2017) of peach samples was analyzed by “adonis” function from R package “vegan” (Oksanen et al., 2019). The regularized log-transformed gene counts were used for the PERMANOVA and the ratios of rRNA content, genotypes, and time points were used as factors.

## RESULTS

### Quantification of Gene Expression During Dormancy in Apricot

In *Prunus* spp., individual genotypes have specific CR accommodating various climatic zones, with trees with lower CR blooming earlier on average and trees with higher CR blooming later (Alburquerque et al., 2008; Fan et al., 2010; Castède et al., 2014; Benmoussa et al., 2017; Fadón et al., 2018b; Balogh et al., 2019). For this experiment, genotypes with varying average bloom dates when planted in the same location were selected. To determine the phenotypes of the four apricot genotypes, the bloom dates from 1999 to 2013 were recorded. Genotypes A1956 and A2137 reached 50% bloom after an average of 60 Julian calendar days and are hereafter referred to as “early blooming” (Audergon, 1993; Dirlewanger et al., 2012; Andreini et al., 2014; Quero-Garcia et al., 2016). Genotypes A660 and A1267 reached 50% bloom after more than 70 Julian days on average (“late blooming”) (Supplementary Figure 2A). A2137 was the earliest genotype entering sepal visible stage at 577 chill hours, while the latest one to enter the visible sepal stage at 1,002 chill hours was genotype A1267 (Supplementary Figure 2B).

Over the course of a single winter, bud samples were collected at 0, 100, 400, and 800 chill hours as well as at the sepal visible stage and the petal visible stage. Due to the variation in dormancy release time, early blooming trees had already bloomed at 800 chill hours while the late blooming trees were still dormant. Therefore, samples from A660 and A1267 at 800 chill hours are labeled as Bud-800, and samples from A2137 and A1956 are labeled as Flower-800 for clarity in developmental stage of tissues collected.

RNA sequencing was performed with individual libraries for each set of buds collected from each individual tree, yielding from 19.1 million to 35.2 million raw reads per library. Peach and apricot have structurally similar genomes as identified by the comparison of apricot and peach genetic maps (Dirlewanger et al., 2004; Arús et al., 2006). As a *Prunus armeniaca* reference genome and annotations are not publicly available, the peach genome was tested as a possible alternative. Over 90% of the raw reads were successfully mapped, and approximately 70%-80% of the reads map to annotated genes (Supplementary Figure 3). All additional apricot results are reported in reference to the *Prunus persica* v2.0 genome.

The expression of the *DAM* genes, a tandem duplication of six genes that are well-known regulators of dormancy in peach, were examined. The six *DAM* genes were found to be misannotated in *Prunus persica* v2.0 genome, with exons from
DAM 1-3 grouped together in a single gene model (Supplementary Figure 1). Therefore, we replaced the gene counts in this region using the *Prunus persica* whole genome assembly v1.0, where the *DAM* genes have been manually annotated (The International Peach Genome Initiative et al., 2013). The six *DAM* genes had lower expression levels in the early blooming genotypes than in the late blooming genotypes throughout dormancy to bloom (Figure 1). *DAM1* to *DAM3* genes had relatively high expression levels at the beginning of chilling, and were downregulated as chill hours increased. *DAM1* in the late blooming genotypes was increased at 400 chill hours, while *DAM1* expression was downregulated at 400 chill hours in early blooming genotypes (Figure 1). Notably, *DAM4* and *DAM5* in the late blooming genotypes were highly expressed at the beginning of endodormancy, and downregulated after 100 chill hours (i.e. *DAM5*) or 400 chill hours (i.e. *DAM4*). In the early blooming genotype A2137, the expression levels of *DAM4* had little change during dormancy. *DAM6* expression levels were decreased consistently in all genotypes as chill hours accumulated. However, the four genotypes have different expression levels at 0 chill hour: the earliest blooming genotype A2137 had the lowest *DAM6* expression while the late blooming genotype A660 had the highest expression level (Figure 1). Overall, the expression of the six *DAM* genes were downregulated at 400 chill hours in the early blooming genotypes, while in the late blooming genotypes, they reached the same level at 800 chill hours (Figure 1). This suggests the expression level is correlated with the timing of growth resumption. The *DAM* genes are highly expressed at the onset of dormancy followed by downregulation at bud break. This is consistent with previous studies that *DAM* genes are associated with dormancy regulation (Bielenberg et al., 2008; Li et al., 2009; Jiménez et al., 2010b). Similar gene expression profiles of *DAM5* and *DAM6* were observed in other apricot genotypes from the same year (Balogh et al., 2019). The difference of *DAM5*

expression levels in different apricot genotypes at 400 chill hours was also seen in different peach genotypes (Leida et al., 2012c). However, our results for *DAM5* and *DAM6* expression levels differ from the findings of Li et al. (2009) who reported that *DAM5* and *DAM6* were upregulated until CR was met in peach (Li et al., 2009).

Transcriptome Profiles Clustered by Bud
Developmental Stages, Regardless of Chill
Hours

To evaluate the overall transcriptomes of four apricot genotypes, the sample distances were represented using a principal component analysis (PCA). The first two PCs represented the two largest variances in the gene expression profiles, which explain 78% and 10% of the variance for the first and second PCs, respectively. The 60 transcriptome profiles were clearly separated into five groups (Figure 2A) correlated with developmental stages. As CR has not been measured for the four apricot genotypes, the endodormant and ecodormant developmental stages were unknown at the time of collection. However, the PCA plot reveals one main cluster of most dormant bud samples (cluster 1) and a smaller set of dormant bud samples in a separate cluster (cluster 2). Cluster 2 is formed of samples of the earliest blooming genotype (A2137) at 400 hours and one of the later blooming genotypes (A660) at 800 hours. Both the timing and comparison to bloom date data suggest this cluster represents ecodormancy, i.e., as these two genotypes have different bloom dates, the early blooming genotype moved into ecodormancy earlier than the late blooming genotype. Likely, the ecodormant stage was not captured in the other two genotypes. Although it is difficult to distinguish endodormancy and ecodormancy physically from the buds, the gene expression profiles indicated that the trees were progressing towards the same developmental stages but at different rates. From here on the five PCA clusters will be referred to by these development

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | The expression of six *DAM* genes during dormancy through bud break in the four apricot genotypes. Top: *DAM1* to *DAM3*, Bottom: *DAM4* to *DAM6*. Two red lines represent early blooming genotypes, two blue lines represent late blooming genotypes (mean ± SE). Biological replicates range from 1 to 3 for each genotype and time point (**Supplementary Tables 1** and **2**).</figcaption></figure>
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Transcriptome variances distinguish endodormancy and ecodormancy. (A) A principal component analysis (PCA) of samples by transcriptome profile in apricot. PC1 and PC2 represent the first two largest sample variances from overall gene expression. The samples were clustered into five groups. 1: endodormancy, 2: ecodormancy, 3: sepal, 4: petal, 5: flower. Biological replicates ranged from 1 to 3 for each genotype and time point (Supplementary Tables 1 and 2). (B) The number of DEGs between early blooming and late blooming genotypes at 0, 100, 400, and 800 chill hours, with red and blue bars representing the number of upregulated and downregulated genes, respectively. Numbers of DEGs upregulated and downregulated are displayed above the bars. (C) Volcano plot of DEGs (red) and non-DEGs (black) (log2 fold change| > 1, adjusted p-value < 0.05) with 1,203 genes upregulated and 164 genes downregulated at ecodormancy compared to endodormancy.

time points: 1- endodormancy, 2- ecodormancy, 3- sepal visible, 4- petal visible, and 5- flower.

To understand the phenotypic effect on dormancy transition between early and late blooming phenotypes, we compared the gene expression of early blooming trees and late blooming trees at each of three time points. After experiencing 0, 100, and 400 chill hours, we found 25, 114, and 553 genes, respectively, were differentially expressed between the two phenotypes (Figure 2B and Supplementary Tables 3–5), suggesting relatively small differences in the transcriptome profiles at the earliest sampling times. However, at 800 chill hours, early blooming genotypes had passed the bud break stage, while the late blooming genotypes were still in dormancy. At this point, the expression of 2,073 (1,354 upregulated and 719 downregulated) genes significantly changed (Figure 2B). However, since different genotypes progressed developmentally at different rates, comparing genotypes at the same time point does not yield

information on the developmental stages of dormancy. To investigate further, we compared the expression of genes between endodormancy and ecodormancy stages based on the PCA clusters (regardless of the date of sampling) by contrasting the 28 samples in endodormancy and five samples in ecodormancy. In total, 1,367 out of 26,872 genes were significantly differentially expressed with the cutoff for false discovery rate < 0.05 and an absolute fold change > 2 (log2FC| > 1). This analysis yielded 1,203 genes that were upregulated and 164 genes that were downregulated at the ecodormancy stage (Figure 2C).

To identify the functions and biological pathways that the DEGs are involved in, GO enrichment analysis was applied. The 1,367 genes differentially expressed between endodormancy and ecodormancy were analyzed by AgriGO v2.0 (Tian et al., 2017), yielding significant GO terms (FDR < 0.05) involved in pathways such as oxidation reduction (GO:0055114, FDR = 0.00075),
carbohydrate metabolic process (GO:0005975, FDR = 0.0057), and hydrolase activity (GO:0016787, FDR = 0.014) (Table 3). The GO enrichment analysis of 640 DEGs between early and late blooming genotypes at the 0, 100, and 400 chill hours indicated enrichment of GO terms related to cell wall metabolism (GO:0044036, FDR= 0.00089) and chitin metabolic process (GO:0006030, FDR = 0.00089) (Table 4). Our results indicate both oxidative stress and cell wall modification are important during the period of CR acquisition during endodormancy. Changes in cell wall structure and components are critical for the processes of growth and stress resistance (Zwiazek, 1991), and previous studies have shown the increase of cell wall formation and reactive oxygen species (ROS) signaling (Jian et al., 1997; Considine and Foyer, 2014) during bud dormancy. Chitinases were reported to be associated with protection against freezing and promote cell wall biosynthesis during dormancy in white spruce (González et al., 2015).

## Co-Expression Analysis Identified Endodormancy and Ecodormancy Related Genes

Co-expression clustering was performed to provide better insight into coordinated gene expression modules that contribute to the endodormancy to ecodormancy transition. 25,960 genes detected by RNA-seq were clustered into 23 co-expression modules based on the TOMsimilarity algorithm in the WGCNA R package (Langfelder and Horvath, 2008). Module size ranged from 126

genes (ME22) to 5,967 genes (ME1) (Figure 3B). Module-factor relationship analysis correlated the expression level of the modules to the five developmental stages. Higher correlation p-values of a module at a particular stage indicate the genes in that module have higher expression at that stage. According to the module-factor relationship, we identified two modules, module 2 (ME2) and module 15 (ME15) most correlated to dormancy stages (Figure 3A). Module 2 is the most positively related to endodormancy and contains 2,350 genes that are expressed highly at the endodormancy stage (Figures 3A, C and Supplementary Table 6). GO network by BiNGO (Maere et al., 2005) associates the genes in module 2 with biological regulation, stress response, reproduction including embryonic development and cell cycle, metabolic process, and chromatin modification (Figure 3D and Supplementary Table 7). Module 15 is the most related to ecodormancy (Figures 3A, C) and the GO network of the 422 genes in the module indicates that they are associated with hormone biosynthesis including ABA catabolism and secondary metabolism such as terpenoid and isoprenoid catabolic process, as well as flower development such as pollen exine formation, sporopollenin biosynthesis, and anatomical structure morphogenesis (Figure 3E and Supplementary Tables 8 and 9). The co-expression network patterns indicated that the genes related to pollen and flower development are activated at ecodormancy, suggesting an increase in development of reproductive tissues in the bud at ecodormancy after the CRs are fulfilled.

<table><caption>TABLE 3 | Significantly enriched (FDR > 0.05) GO terms of DEGs comparing endodormancy and ecodormancy.</caption><thead><tr><th>GO term</th><th>Category</th><th>Functional description</th><th>Number in DEGs</th><th>Number in genome</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0055114</td><td>P</td><td>Oxidation reduction</td><td>107</td><td>1195</td><td>0.00075</td></tr><tr><td>GO:0005975</td><td>P</td><td>Carbohydrate metabolic process</td><td>58</td><td>591</td><td>0.0057</td></tr><tr><td>GO:0004180</td><td>F</td><td>Carboxypeptidase activity</td><td>13</td><td>44</td><td>0.00012</td></tr><tr><td>GO:0004185</td><td>F</td><td>Serine-type carboxypeptidase activity</td><td>13</td><td>42</td><td>0.00012</td></tr><tr><td>GO:0070008</td><td>F</td><td>Serine-type exopeptidase activity</td><td>13</td><td>48</td><td>0.00018</td></tr><tr><td>GO:0016491</td><td>F</td><td>Oxidoreductase activity</td><td>118</td><td>1344</td><td>0.00018</td></tr><tr><td>GO:0008236</td><td>F</td><td>Serine-type peptidase activity</td><td>24</td><td>161</td><td>0.0007</td></tr><tr><td>GO:0017171</td><td>F</td><td>Serine hydrolase activity</td><td>24</td><td>161</td><td>0.0007</td></tr><tr><td>GO:0004553</td><td>F</td><td>Hydrolase activity, hydrolyzing O-glycosyl compounds</td><td>43</td><td>377</td><td>0.0007</td></tr><tr><td>GO:0008238</td><td>F</td><td>Exopeptidase activity</td><td>13</td><td>59</td><td>0.0011</td></tr><tr><td>GO:0016798</td><td>F</td><td>Hydrolase activity, acting on glycosyl bonds</td><td>43</td><td>396</td><td>0.0013</td></tr><tr><td>GO:0005506</td><td>F</td><td>Iron ion binding</td><td>44</td><td>423</td><td>0.0027</td></tr><tr><td>GO:0020037</td><td>F</td><td>Heme binding</td><td>40</td><td>385</td><td>0.0051</td></tr><tr><td>GO:0046906</td><td>F</td><td>Tetrapyrrole binding</td><td>40</td><td>386</td><td>0.0051</td></tr><tr><td>GO:0009055</td><td>F</td><td>Electron carrier activity</td><td>18</td><td>129</td><td>0.011</td></tr><tr><td>GO:0016787</td><td>F</td><td>Hydrolase activity</td><td>144</td><td>1945</td><td>0.014</td></tr><tr><td>GO:0016651</td><td>F</td><td>Oxidoreductase activity, acting on NADH or NADPH</td><td>8</td><td>34</td><td>0.016</td></tr><tr><td>GO:0016705</td><td>F</td><td>Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen</td><td>33</td><td>323</td><td>0.018</td></tr><tr><td>GO:0016209</td><td>F</td><td>Antioxidant activity</td><td>15</td><td>106</td><td>0.023</td></tr><tr><td>GO:0070011</td><td>F</td><td>Peptidase activity, acting on L-amino acid peptides</td><td>43</td><td>467</td><td>0.024</td></tr><tr><td>GO:0003824</td><td>F</td><td>Catalytic activity</td><td>412</td><td>6346</td><td>0.04</td></tr><tr><td>GO:0008233</td><td>F</td><td>Peptidase activity</td><td>43</td><td>485</td><td>0.045</td></tr></tbody></table>

P. biological processes, F. molecular functions.
<table><caption>TABLE 4 | Significantly enriched GO terms of DEGs comparing early blooming and late blooming genotypes in endodormancy time points (0, 100 and 400 chill hours).</caption><thead><tr><th>GO term</th><th>Category</th><th>Functional description</th><th>Number in DEGs</th><th>Number in genome</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0044036</td><td>P</td><td>Cell wall macromolecule metabolic process</td><td>6</td><td>20</td><td>0.00018</td></tr><tr><td>GO:0016998</td><td>P</td><td>Cell wall macromolecule catabolic process</td><td>6</td><td>19</td><td>0.00018</td></tr><tr><td>GO:0006030</td><td>P</td><td>Chitin metabolic process</td><td>5</td><td>13</td><td>0.00018</td></tr><tr><td>GO:0006032</td><td>P</td><td>Chitin catabolic process</td><td>5</td><td>13</td><td>0.00018</td></tr><tr><td>GO:0006026</td><td>P</td><td>Aminoglycan catabolic process</td><td>5</td><td>13</td><td>0.00018</td></tr><tr><td>GO:0006022</td><td>P</td><td>Aminoglycan metabolic process</td><td>5</td><td>17</td><td>0.00066</td></tr><tr><td>GO:0000272</td><td>P</td><td>Polysaccharide catabolic process</td><td>5</td><td>22</td><td>0.0022</td></tr><tr><td>GO:0071554</td><td>P</td><td>Cell wall organization or biogenesis</td><td>10</td><td>104</td><td>0.0032</td></tr><tr><td>GO:0005976</td><td>P</td><td>Polysaccharide metabolic process</td><td>8</td><td>91</td><td>0.02</td></tr><tr><td>GO:0005975</td><td>P</td><td>Carbohydrate metabolic process</td><td>27</td><td>591</td><td>0.024</td></tr><tr><td>GO:0004568</td><td>F</td><td>Chitinase activity</td><td>5</td><td>13</td><td>0.00082</td></tr><tr><td>GO:0004553</td><td>F</td><td>Hydrolase activity, hydrolyzing O-glycosyl compounds</td><td>23</td><td>377</td><td>0.0033</td></tr><tr><td>GO:0016798</td><td>F</td><td>Hydrolase activity, acting on glycosyl bonds</td><td>24</td><td>396</td><td>0.0033</td></tr><tr><td>GO:0016787</td><td>F</td><td>Hydrolase activity</td><td>75</td><td>1945</td><td>0.011</td></tr><tr><td>GO:0008236</td><td>F</td><td>Serine-type peptidase activity</td><td>12</td><td>161</td><td>0.011</td></tr><tr><td>GO:0017171</td><td>F</td><td>Serine hydrolase activity</td><td>12</td><td>161</td><td>0.011</td></tr><tr><td>GO:0003824</td><td>F</td><td>Catalytic activity</td><td>207</td><td>6346</td><td>0.014</td></tr><tr><td>GO:0016684</td><td>F</td><td>Oxidoreductase activity, acting on peroxide as acceptor</td><td>8</td><td>95</td><td>0.023</td></tr><tr><td>GO:0004601</td><td>F</td><td>Peroxidase activity</td><td>8</td><td>95</td><td>0.023</td></tr><tr><td>GO:0070001</td><td>F</td><td>Aspartic-type peptidase activity</td><td>7</td><td>84</td><td>0.027</td></tr><tr><td>GO:0016747</td><td>F</td><td>Transferase activity, transferring acyl groups other than amino-acyl groups</td><td>12</td><td>198</td><td>0.027</td></tr><tr><td>GO:0016209</td><td>F</td><td>Antioxidant activity</td><td>8</td><td>106</td><td>0.027</td></tr><tr><td>GO:0004190</td><td>F</td><td>Aspartic-type endopeptidase activity</td><td>7</td><td>84</td><td>0.027</td></tr><tr><td>GO:0004180</td><td>F</td><td>Carboxypeptidase activity</td><td>5</td><td>44</td><td>0.027</td></tr><tr><td>GO:0004185</td><td>F</td><td>Serine-type carboxypeptidase activity</td><td>5</td><td>42</td><td>0.027</td></tr><tr><td>GO:0004175</td><td>F</td><td>Endopeptidase activity</td><td>14</td><td>258</td><td>0.03</td></tr><tr><td>GO:0070008</td><td>F</td><td>Serine-type exopeptidase activity</td><td>5</td><td>48</td><td>0.031</td></tr><tr><td>GO:0016746</td><td>F</td><td>Transferase activity, transferring acyl groups</td><td>13</td><td>237</td><td>0.032</td></tr><tr><td>GO:0004252</td><td>F</td><td>Serine-type endopeptidase activity</td><td>7</td><td>97</td><td>0.045</td></tr></tbody></table>

P. biological processes, F. molecular functions.

## Validation of RNAseq by qPCR

To validate our transcriptomic analysis, we performed a reverse transcription qPCR (RT-qPCR) experiment using floral bud RNAs from three early blooming apricot genotypes: A2137, A1956, and A2312. We aimed to validate the RNASeq using the same RNA samples for two of the genotypes (A2137 and A1956), and a previously unexamined early blooming genotype (A2312). A2312 had not previously been analyzed by RNASeq, however, bud samples were collected at the same time points and from the same orchard as the other genotypes. As samples at the 0 chill hour and 100 chill hours were very similar in expression pattern in the previous RNASeq analysis, we chose samples from 100, 400, and 800 chill hours as well as sepal stage for the RT-qPCR validation (Supplementary Table 2). The five DEGs with the highest fold change between endodormancy and ecodormancy stage were selected. Prupe.7G084900 is a fatty acyl-CoA reductase, which is associated with male sterility and is involved in oxidoreductase activity and pollen exine formation (Chen et al., 2011). Prupe.2G122600 encodes a member of the chalcone and stilbene synthase family that most closely matches the *Arabidopsis thaliana* gene less adhesive pollen 6 (LAP6, AT1G02050), involved in tapetosome development and

pollen viability [The Arabidopsis Information Resource (TAIR); Dobritsa et al., 2010; Wang et al., 2018]. Prupe.5G014900 is a glutathione S-transferase, many of which provide cellular detoxification and have been profiled in stress response in many plant species (Seppänen et al., 2000; Anderson and Davis, 2004; Jain et al., 2010). Prupe.1G104900 is a late embryogenesis abundant protein (LEA), a family known to be involved in cold and drought stress response (Pedrosa et al., 2015) and previously profiled in dormant apricot (Yamane et al., 2006) and oak buds (Ueno et al., 2013). The final RT-qPCR gene was *DAM6* (ppa010714m) discussed above.

The RT-qPCR analysis of the five genes all matched the RNASeq results in the two previously sequenced genotypes. The expression pattern for A2312, the newly examined genotype, was similar to the other two early blooming genotypes, confirming the replicability of these results (Figure 4). For all three genotypes, the *DAM6* gene was reduced in expression as chill hours accumulated until the release of dormancy at the sepal stage. The expression level of *DAM6* in the genotype A2137, the earliest blooming tree, was already lower at 100 chill hours than that in the other two genotypes. Two of the genes, jojoba acyl CoA reductase-related male sterility protein (Prupe.7G084900) and chalcone and
A Module-factors relationships

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

B Number of genes in each module

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

C Module expression profiles

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

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

D Response to stimulus

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

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

ME15

**FIGURE 3** | Co-expression modules detected genes induced at endodormancy and ecodormancy. (A) The module-factor relationships between modules and developmental stages. (B) The number of differentially expressed genes (DEGs) and non-DEGs in each module. (C) The expression profiles of ME2 (top) and ME15 (bottom) eigengenes in four apricot genotypes (mean ± SE). The eigengene represents all genes in the module as a single normalized expression pattern. (D) The enriched GO network for ME2 genes, and (E) the enriched network for ME15 genes. Nodes represent individual GO terms, and node color indicates the *p*-value. General GO categories for the nodes are labeled here with boxes. The individually labeled node networks are available in **Supplementary Figures 4** and **5**. Detailed descriptions of nodes and their *p*-values are in **Supplementary Tables 7** and **9**.

stilbene synthase (Prupe.2G122600), were only induced at the ecodormancy stage in a single genotype, A2137.

## Expression Patterns in Peach Buds Were Similar to Apricot During Dormancy Transition

To confirm the transferability of these results to peach, a similar transcriptome profiling experiment was performed using a set of F₂ peach trees with different CRs. Four genotypes were selected with low CR and high CR: two genotypes with 300 hours CR (A209 and A340; low CR) and two genotypes with 850 hours and 1,100 hours of CR (A318 and A323, respectively; high CR) (Fan et al., 2010; Zhebentyayeva et al., 2014). Floral buds from the low CR genotypes were collected at 0, 100, 600 chill hours, and pre-

bloom, while the high CR genotypes were collected with an additional time point, 1,000 chill hours. Based on imaging of dissected buds, (Supplementary Figure 6), the carpels from low CR buds were more mature than those from high CR buds at 600 chill hours while the high CR buds developed similar carpal sizes at 1,000 chill hours (Supplementary Figure 6). This result suggests tissue differentiation and development proceeds at different rates based on CR with significant floral development activated at ecodormancy. This agrees with previous reports of cell division and differentiation in anthers during winter dormancy, with microsporogenesis beginning after reaching CR in apricot (Julian et al., 2011).

To examine the changes in transcriptome profiles at sample time points in these peach floral buds, RNASeq analysis was conducted with 69 peach libraries using a similar experimental approach to
<figure><img src="image_11.png" /><figcaption>**FIGURE 4** | Validation of differentially expressed genes (DEGs) with RT-qPCR. The expression pattern of DEGs in the RNA-seq (top) matched the qPCR (bottom) results. Pink: A2137 (Bakour); red: A1956 (Palsteyn); grey: A2312 (Flamingold). *DAM6*, DORMANCY-ASSOCIATED MADS-box gene 6; *LEA*, Late embryogenesis abundant protein gene; *MS2*, Jojoba acyl CoA reductase-related male sterility gene; *LAP6*, Chalcone and stilbene synthase; *GSTF13*, Glutathione S-transferase family gene. One RNA sample from each genotype at each time point was assayed with three technical replicates.</figcaption></figure>

that of the apricot. Unlike the apricot data, significant ribosomal RNA (rRNA) contamination was identified in more than 50% of the libraries (Supplementary Figure 7). To inspect whether the rRNA contamination significantly affected the transcriptome profiles, a PERMANOVA was conducted with transcriptome distance matrices. The result indicated that the percentage of rRNA in the library is a significant factor in the transcriptome variation, but it does not show a significant interaction with the genotype or time point factors (Supplementary Table 10). The

rRNA content was likely caused by the failure of rRNA depletion during library preparation before sequencing and, based on the statistical results, can be adequately factored out of the experimental analysis. This is further supported by PCA plotting, showing the gene expression profiles clustering into developmental stages regardless of rRNA content (Figure 5). In a strikingly similar pattern to the apricot data, the clustering reveals that the low CR buds were in an ecodormant state at 600 chill hours while the high CR buds were still at endodormancy. One genotype may indicate

<figure><img src="image_12.png" /><figcaption>**FIGURE 5** | Transcriptome variances of peach floral buds in dormancy during chill accumulation. PCA analysis of peach bud transcriptome profiles. PC1 and PC2 represented the largest two components of variation identified among samples. Solid circles indicate samples at endodormancy (bottom left), ecodormancy (bottom right) and pre-bloom (top) stages. The dashed circle indicates the four A318 samples taken 1,000 chill hours. Despite this genotype having a previously measured CR of 850 hours, these samples cluster more closely with endodormant buds. All genotypes and time points had three or four biological replicates (Supplementary Tables 1 and 2).</figcaption></figure>
this pattern does not perfectly correspond to endo- and
ecodormancy. Buds at 1,000 chill hours from genotype A318 with
850 CR (i.e. putatively ecodormant) were clustered more closely to
other endodormant buds (Figure 5). Whether the buds were
actually endodormant or ecodormant can only be predicted from
chill hours and was not directly assayed. However, the overall
transcriptome profile variations in peach and apricot both indicate
two distinct transcriptional programs that largely correspond to
endodormancy and ecodormancy. To ensure that the differential
gene expression analysis of peach samples clearly represents these
two stages, the differential expression analysis was performed on
samples from the 600 hour time point, in which both the PCA-
based clustering and predicted chill hours for each cultivar both
agree on the endodormant and ecodormant stage of the samples.

The expression profiles of the *DAM* genes were comparable in peach and apricot. Similar to apricot *DAM5* and *DAM6*, peach *DAM5* and *DAM6* were expressed at a high level at the beginning of endodormancy then downregulated after 100 chill hours. However, a different expression profile for *DAM3* was found between the species. In peach, *DAM3* shared the same expression pattern with *DAM4*, generally upregulated in endodormancy, then downregulated at ecodormancy (**Figure 1** and **Supplementary Figure 8**). This is different from *DAM3* in apricot, which is consistently downregulated from the beginning of dormancy to the end.

Transcriptome profiles for low chill and high chill genotypes
were compared to identify the DEGs at each time point. As
expected, at 0 and 100 chill hours, when low chill and high chill
genotypes were both in endodormancy, few gene expression
differences were found: 49 DEGs at 0 chill hour and 77 DEGs at
100 chill hours (**Supplementary Figure 9A**). At the time point of
600 chill hours, when the two low chill genotypes moved into
ecodormancy and the high chill genotypes were still in
endodormancy, 2,102 genes were differentially expressed, with
1,603 upregulated and 499 downregulated (**Supplementary
Figures 9A, B**). Four of the six *DAM* genes (*DAM3*, *DAM4*,
*DAM5*, and *DAM6*) were included in the DEGs (**Supplementary
Table 11**). The DEGs were enriched in the pathways of oxidation
reduction process (GO:0055114, FDR = 1.01e-08), carbohydrate
metabolic process (GO: 0005975, FDR = 1.7e-07), and cell wall
metabolic process (GO: 0044036, FDR = 4.89e-04), which were
also enriched in the apricot DEGs (**Supplementary Figure 9C**).
These results confirm that the genes involved in stress response,
sugar metabolism, and cell wall assembly contribute to the
endodormancy to ecodormancy transition, and that this pattern
of expression is replicable across at least two *Prunus* species.

Comparison of Peach and Apricot
Transcriptome Profiles

To compare peach and apricot, we considered only the DEGs
found to differentiate endodormancy and ecodormancy stages.
The DEGs in both apricot and peach were distributed across all
eight chromosomes (Figure 6A). Of the 1,367 DEGs identified in
apricot and the 2,102 DEGs in peach, 608 genes were statistically
significant in both species (Figure 6B). When comparing with the
fold change of DEGs in peach and apricot, over 99% of the DEGs

were either upregulated or downregulated consistently in both
species, with only six genes expressed inconsistently (Figure 6C,
Supplementary Table 11). Prupe.8G238000, Prupe.6G047900,
and Prupe.5G057900 were downregulated at ecodormancy in
apricot but were upregulated in peach. Prupe.3G259300,
Prupe.4G069300, and Prupe.3G127600 were upregulated in
apricot but downregulated in peach. These very consistent
patterns between peach and apricot gene expression suggest
there is a highly conserved gene regulatory system governing
dormancy stages and transitions across *Prunus* spp. as well as a
few potentially important variations.

CR and bud break QTLs have previously been identified in
both apricot and peach linkage maps (Olukolu et al., 2009; Fan
et al., 2010; Zhebentyayeva et al., 2014). Using these genetic
mapped positions, we sought to filter our DEGs for possible
candidate genes in these QTL regions. For apricot, there are
1,059 markers [43 simple sequence repeats (SSR) and 1,016
amplified fragment length polymorphism (AFLP)] located on
the Perfection and A.1740 linkage map (Olukolu et al., 2009).
The detected CR QTLs regions were delimited by AFLP markers,
which were unable to be located in the peach reference genome
due to a lack of marker-associated sequence data. For peach, 15
QTLs associated with CR from a cross of peach cultivars "Fla.92-
2C" and "Contender" (Fan et al., 2010; Zhebentyayeva et al.,
2014) encompass a total of 3,571 genes. These genes were
extracted to compare with apricot and peach DEGs. Of 1,367
apricot DEGs, 213 were detected in the QTLs while 298 out of
2,102 peach DEGs overlapped with these QTLs (Figure 6A). A
final set of 99 genes was identified as being located within the
QTL regions and statistically significant in dormancy-associated
expression changes in both apricot and peach (Figure 6A).
*DAM4* (ppa011123m), *DAM5* (ppa010822m), and *DAM6*
(ppa010714m) were part of this final set and were the only
genes located on QTL qCR1a-2009.

In order to broadly understand the functions of DEGs in
peach and apricot individually as well as their shared
functions, a GO enrichment analysis was conducted for
three sets of genes: the shared DEGs (608), the DEGs
unique to apricot (759), and the DEGs unique to peach
(1,494). When examining GO term enrichment of the DEGs
shared in both peach and apricot, significantly enriched terms
include lipid localization (GO:0010876, FDR=1.5e-13),
response to stimulus (GO:0050895, FDR=5.3e-9), pollen
wall assembly (GO:0010208, FDR=2.6e-11), and cellular
component assembly involved in morphogenesis processes
(GO:0010927, FDR=2.6e-11) (Figure 6D). The 2,253 DEGs
identified in either peach or apricot but not both, were
enriched in the same biological pathways as the DEGs that
were only significant in one species. Despite the two lists not
sharing any genes, when examining the unique peach DEGs
and the unique apricot DEGs, four of the top 10 enriched GO
terms based on the lowest FDR values were shared
(Figures 6E, F). This shared set of GO terms suggests that
peach and apricot share many of the same functions and
pathways during dormancy stages and transitions, even if
different individual genes are being identified as
<figure><img src="image_13.png" /></figure>

**FIGURE 6** | Comparison of differentially expressed genes (DEGs) between endodormancy and ecodormancy in apricot and peach. (A) The distribution of DEGs and quantitative trait loci (QTLs) on the reference genome. Red and black dots represent up and down fold changes, respectively. (B) Venn diagram of apricot and peach DEGs overlapping with peach CR QTLs. (C) A comparison between DEGs in apricot and peach based on the log2-fold change (logFC). Colors indicate genes that were statistically significant in apricot only (green), peach only (purple), or both (orange). (D–F) Top 10 GO terms enriched in DEGs found in both species (D), peach only (E), and apricot only (F).

differentially expressed. For example, 232 DEGs involved in stimulus response were in peach but not in apricot, while 119 stimulus response genes were differentially expressed in apricot but not in peach. A similar situation was found for phenylpropanoid biosynthetic pathway, response to abiotic stimulus, and response to chemical stimulus, i.e. they are significant in both species but mostly consisted of different underlying genes from peach and apricot. In contrast, pollen and floral development genes tended to be found as DEGs in both species, and thus GO terms relevant to those processes appear on the shared list.

There are pathways that were enriched only for one species but not for the other. For example, GO terms for photosynthesis and carbon fixation were significant in apricot. Fourteen and five DEGs, respectively, were found to be associated with those GO terms (Figure 6E), and all are found in the chloroplast genome. A different set of five DEGs involved in photosynthesis was found in the peach transcripts, but was not significantly differentially expressed between stages. While chloroplasts occur in green floral tissues such as sepals, they have also been found to be photosynthetically active in other floral organs such as petals and corollas in other plants (Weiss et al., 1989; Pyke and Page, 1998).
Chloroplasts have also been proposed as essential for redox processes and developmental signaling during flower development (Muñoz and Munné-Bosch, 2018).

## Comparison of Co-Expression Networks Between Apricot and Peach

Similar to apricot co-expression networks, peach bud transcriptome profiles distinguished co-expressed genes that were upregulated at endodormancy or ecodormancy, respectively. ME6, the cluster of genes with the highest expression level at endodormancy was identified as positively related to endodormancy stage (correlation = 0.63). (Supplementary Figure 10A). The co-expression networks by Cytoscape indicated that the ME6 genes were enriched in phenylpropanoid biosynthetic process (GO:0009699, FDR=0.024), regionalization (GO:0003002, FDR=0.0319), and positive regulation of flavonoid biosynthetic process (GO:0009963, FDR=0.0319) pathways (Supplementary Figure 10C and Supplementary Table 12). ME4 and ME10, were both identified as positively related with the ecodormancy stage (ME4 correlation = 0.94, ME10 correlation = 0.75) (Supplementary Figures 10A, B). ME4 genes were enriched in organic acid catabolic process (GO:0016054, FDR=7.59e-7), oxoacid metabolic process (GO:0043436, FDR=9.68e-7), and response to stimulus (GO:0050896, FDR=9.68e-7). ME10 genes were enriched in pollen wall assembly (GO:0010208, FDR=1.19e-11), cellular component assembly involved in morphogenesis (GO:0010927, FDR=1.19e-11), and pollen exine formation (GO:0010584, FDR=1.86e-11) pathways (Supplementary Figures 10C, D). Both significant and insignificant GO descriptions are provided in the Supplementary Tables 12–14.

Co-expression networks developed independently from the peach and apricot data identified genes involved in many shared pathways, such as pollen wall assembly and stress response. To further investigate the common expression patterns in both apricot and peach, a co-expression analysis integrating the peach and apricot datasets was performed. Twenty-two thousand three hundred seventeen genes expressed in more than half of the samples in the dataset were clustered into 16 modules spanning from 215 genes (ME15) to 6,089 genes (ME1). DEGs were distributed unevenly across the 16 modules, with the largest number as 174 DEGs in ME11, and no DEGs in ME7 or ME15 (Figure 7B). To understand the major biological processes that genes are involved in for each module, genes were then enriched by AgriGO v2.0 using the *Prunus persica* background. Three hundred thirty-nine significant GO terms from the biological process category were enriched across the 16 modules. ME1 had the most with 84 enriched GO functions, while ME0, ME6, ME13, and ME14 had no significant GO terms (Figure 7D). Most of the modules do not share the same GO terms, suggesting multiple genes in important dormancy-related biological functions are co-expressed in similar patterns across our experimental time points. For example, cell wall macromolecule catabolic process (GO:0016998) was only enriched in ME4, while cell cycle process (GO:0022402) and mitotic cell cycle process (GO:1903047) were only enriched in ME10.

Two modules, ME4 and ME11, are highly correlated to ecodormancy in peach (correlation = 0.9) and apricot (correlation = 0.77) in the co-expression analysis (Figure 7A). A total of 35.5% of the genes in the ME11 were DEGs shared between apricot and peach, and of 1,000 genes in ME4, 129 were shared DEGs (Figure 7B). These two modules contained half of the shared DEGs (303 out of 608), suggesting that genes in these two modules are essential to defining transcriptional patterns in ecodormancy. The expression profile of the ME4 eigengene showed an upregulation at ecodormancy in peach low chill genotypes but up-regulation in only one apricot genotype, A660 (Figure 8C). The expression of ME11 eigengene was upregulated at ecodormancy in both peach and apricot genotypes, however, the eigengene in apricot A660 continued increasing at the sepal stage (Figure 8D). The GO enrichment network of ME4 indicates the genes are involved in response to abiotic and biotic stimuli and response to hormones (ethylene, abscisic acid, gibberellin, and jasmonic acid). A group of genes in ME4 were involved in organic acid transport and secondary metabolic processes, specifically sesquiterpenoid metabolism (Figure 8A and Supplementary Table 15). These pathways indicate the ME4 genes are most likely responding to the abiotic stress imposed during winter (Supplementary Table 15). The other module, ME11, is enriched in pollen wall development, flower development, lignan metabolism, and lipid transport (Figure 8B and Supplementary Table 16). These results suggest that both peach and apricot activate new floral development activities including pollen wall and flower morphogenesis after CRs are fulfilled.

## DISCUSSION

*Prunus* flowers are sensitive to frost damage, and maintenance of winter dormancy for a sufficient time period is essential to prevent frost damage to delicate reproductive tissues, especially during the spring months when temperatures may fluctuate. This is one of the most important adaptation mechanisms that ensure perennial plants survive low winter temperatures and coordinates the timing of reproductive activities (Campoy et al., 2011; Luedeling et al., 2011; Leida, 2012). There have been many efforts to find the possible regulation factors of winter dormancy in perennial trees, including examinations of internal factors such as phase change, circadian clock, and hormones, and external factors including photoperiod, temperature, light radiation, and water (Perry, 1971; Rohde and Bhalerao, 2007; Allona et al., 2008; Cooke et al., 2012; Maurya and Bhalerao, 2017). Our time-series transcriptome profiles in apricot and peach buds allowed us to investigate the gene expression changes as floral bud dormancy progressed towards bloom and to compare complementary data from two closely related species and among genotypes varying in CR/BD.

## Buds Stay Active During Winter Dormancy

Meristems are the functional units for plant development, including root and shoot apical meristems, as well as cambium
<figure><img src="image_14.png" /><figcaption>**FIGURE 7** | Co-expression modules identified genes induced commonly at endodormancy and ecodormancy in peach and apricot. (A) The module-factors relationships with the correlation of module eigengenes (MEs) with time points and developmental stages. LC, low chill; HC, high chill. (B) The number of genes in each module, including differentially expressed genes (DEGs) (in red) shared in both species and non-DEGs (blue). (C) A heatmap with the top 5 enriched gene ontology (GO) terms for each module. (D) The number of significant GO terms shared in both apricot and peach in each module.</figcaption></figure>

and intercalary meristems (Esau 1965). Dormancy of buds at the cellular level is defined as eco-developmental quiescence of a shoot meristematic organ, whereby tissues fail to respond to the environment (Considine and Considine, 2016). However, based on our transcriptome profiles and imaging of internal tissues, buds continued development during endodormancy. With the co-expression analysis in this study, results showed that cell cycle genes were upregulated during endodormancy (Figure 3 and Supplementary Figure 4). For example, the co-expression module with the highest correlated score to endodormancy stage in apricot (ME2) included cyclin D genes (Prupe.1G430500, Prupe.6G229700, and Prupe.8G146700) and cell cycle checkpoint control protein (Prupe.6G321300), indicating cells were actively dividing during endodormancy. Cell division can be induced by cytokinins through increasing cyclin expression (Dewitte et al., 2007). A cytokinin mediated gene, cytokinin response factor 2 (CRF2, Prupe.5G114100) in apricot buds was co-

expressed with cyclin D genes, indicating potentially high levels of cytokinins during endodormancy. Treatment of a synthetic cytokinin has been shown to terminate potato tuber dormancy through downregulation of the AGAMOUS-like MADS-box transcription factor, a homolog to the peach *DAM* genes (Campbell et al., 2014). CRF mutants reduced pistil and ovule growth (Cucinotta et al., 2016) and *CRF2* is transcriptionally induced by cold stress in *Arabidopsis* (Jeon et al., 2016), implying cytokinin signaling may play a critical role in modulating cell division in cold. High expression levels of cell cycle genes and cytokinin response genes during endodormancy suggest that cells in the dormant floral buds were actively dividing and/or differentiating. A previous study of anatomical development of floral structures in peach buds reported slow but continuous differentiation throughout the winter, especially of androecium and gynoecium tissues (Reinoso et al., 2002). Our bud images show increasing definition of floral structures (Supplementary Figure 6).
<figure><img src="image_15.png" /><figcaption>**FIGURE 8** | The enriched gene ontology (GO) networks for ME4 (A) and ME11 (B). (C) The expression profiles of the ME4 eigengene (mean ± SE) in peach (top) and apricot (bottom). (D) The expression profiles of the ME11 eigengene (mean ± SE) in peach (top) and apricot (bottom).</figcaption></figure>

## Epigenetic Modification Plays an Important Role in Dormancy

One set of major dormancy regulators known in fruit trees are *DAMs*, with previous studies reporting that increased trimethylation of H3K27 and decreased H3K4 trimethylation represses the expression of *DAMs* during dormancy (Leida, 2012; Leida et al., 2012b; Saito et al., 2015a). In our study, all six copies of *DAM* genes decreased at the transition from endodormancy to ecodormancy, and were maintained at a low level after ecodormancy (**Figure 1** and **Supplementary Figure 8**). The co-expression networks identified a group of genes co-expressed with *DAM1* (ppa018667m) and *DAM3* (ppa010758m), and that were involved in chromatin modification and organization (**Figure 3D** and **Supplementary Figure 4**), including Prupe.2G042400 and Prupe.1G177800 related to histone H3K4 methylation, Prupe.4G091400, and Prupe.6G153800 and Prupe.1G533300 related to production of small RNA involved in gene silencing. However, none of these genes were differentially expressed between endodormancy and ecodormancy stages. Studies reported H3K4 trimethylation was enriched during bud break in pear (Anh Tuan et al., 2016) and a transcription factor EARLY BUD-BREAK 1 (*EBB1*) associated with H3K4 positively regulated bud break in poplar (Yordanov et al., 2014), confirming that histone methylation

is involved in winter dormancy and bud break of woody trees. Small RNA profiles in Japanese pear have also been shown to be associated with bud dormancy transition (Bai et al., 2016), however, whether the epigenetic effects are causative or consequential remains unclear.

## Genes Responsive to Cold Stress Play a Key Role in Maintaining Endodormancy

One major function of endodormancy is resistance to cold temperatures and the protection of new tissues developing inside of buds. Cold-response genes were identified as DEGs during the endodormancy to ecodormancy transition in apricot and peach. For example, Prupe.2G294400 is a catabolite activator protein (CAP160), which has previously been associated with cold stress response (Kaye et al., 1998). The expression of Prupe.2G294400 initially increased as chilling accumulated and then declined once ecodormancy was reached (**Supplementary Figure 11**), and this expression pattern was consistently observed in all four genotypes of both species. This mirrors experiments in other plants, for example, expression of *CAP160* transcripts in spinach increased when plants were exposed to low temperatures, and a transgenic tobacco overexpressing *CAP160* displayed a higher tolerance to freezing stress (Kaye et al., 1998).
Similar to the CAP160, one LEA gene (Prupe.1G104900) was significantly downregulated during ecodormancy but increased during chill hour accumulation in endodormancy (Figure 4 and Supplementary Figure 11). The LEA protein family has been shown to respond to abiotic stress such as cold, drought and salt (Hara et al., 1999; Chandra Babu et al., 2004; Checker et al., 2012; Duan and Cai, 2012; Pedrosa et al., 2015). The LEA expression pattern in our experiments is consistent with previously observed LEA gene expression, which was induced by short photoperiod and downregulated as dormancy progressed (Jiménez et al., 2010a; Leida et al., 2012a).

CBFs are known to respond to cold and drought in plants and are responsible for 10%-20% of transcriptional changes induced by cold (Stockinger et al., 1997; Vogel et al., 2005; Medina et al., 2011; Shi et al., 2018). Two CBFs were significantly downregulated during the ecodormancy stage in peach but not in apricot (Supplementary Table 11). Overexpression of CBF in peach delayed bud break and induced LEA genes and other cold-induced genes (Wisniewski et al., 2015; Artlip et al., 2019). CBFs bind upstream of DAM genes and induce their expression at early time points (Niu et al., 2016; Zhao et al., 2018). While we did see upregulation of DAM genes 4 and 5 in apricot and DAM genes 3, 4, 5, and 6 in peach at early time points especially in late blooming genotypes, all DAM genes were relatively highly expressed at the beginning of the experiment, prior to chilling that would induce CBFs. These patterns match previous Prunus DAM gene expression profiles (Li et al., 2009; Yamane et al., 2011) and may indicate regulatory mechanisms other than CBF that operate prior to chill. Cold-responsive (COR) genes are also responsive to CBF signals and are induced by low temperatures, in order to increase plant cold tolerance, leading to cold acclimation (Nylander et al., 2001; Chinnusamy et al., 2007; Chinnusamy et al., 2010). Peach Prupe.7G161100, an AtCOR ortholog, was differentially expressed between dormancy phases in peach but not apricot. It also shared similar expression patterns with the LEA gene (Prupe.1G104900) (Supplementary Figure 11). The induction of cold responsive genes including COR, LEA, and CAP160 during endodormancy suggests that these genes may contribute to cold tolerance of flower buds during endodormancy. The up-regulation of cold tolerance genes lasted longer in the high chill genotypes than in the low chill genotypes in peach, suggesting cold tolerance is consistently maintained throughout endodormancy. In a process analogous to vernalization, it has been proposed that these genes may also contribute to sensing the amount of chill accumulation during endodormancy (Shi et al., 2018). However, after trees fulfilled their CR and moved into ecodormancy, the genes were decreased to a low level of expression, despite the continuing vulnerability to frost or freeze damage during ecodormancy. Therefore, these particular cold responsive genes may only respond to cold and induce tree cold tolerance during endodormancy. Other mechanisms may be involved in combating cold stress during ecodormancy. As these trees are grown outside their climatically adapted ranges, this could also be due to a misalignment of the phenology of the trees with the locale and would be interesting to investigate further in native populations.

## ROS Responsive Genes Are Differentially Expressed Between Endodormancy and Ecodormancy

Oxidative stress response has been found during dormancy and is hypothesized to be not only a stress response but also an important signaling mechanism (Considine and Foyer, 2014). Our results showed that 63 DEGs found in both Prunus species were enriched in response to stress, which included oxidation stress (Figure 6). Co-expression networks of apricot and peach also indicated ME4, which was induced during ecodormancy, was enriched in genes involved in the oxidation reduction process (Figure 7C). The production of ROS, especially hydrogen peroxide ($H_2O_2$) during endodormancy release was identified as the major signal that triggered antioxidation pathways, such as an increase of peroxidase and superoxide dismutase (SOD) (Leida et al., 2010; Tan et al., 2010; Viti et al., 2010; Prassinos et al., 2011; Vergara et al., 2012; Bai et al., 2013). In our analyses, six oxidoreductase genes (Prupe.2G288400, Prupe.3G168200, Prupe.5G089100, Purpe.6G262200, Prupe.8G156700, and Prupe.8G195100) and five peroxidase family proteins (Prupe.1G081600, Prupe.6G191900, Prupe.6G192000, Prupe.6G239400, and Prupe.7G035800) were significantly upregulated during ecodormancy in both species; two of them (Prupe.1G081600 and Prupe.5G089100) are within the peach CR QTL regions (qCR1d-2008 and qCR5-2009, respectively). Previous studies have shown that exogenous application of $H_2O_2$ can induce endodormancy break in Japanese pear (Kuroda et al., 2005); an increase of $H_2O_2$ was also observed in sweet cherry and peach after applying hydrogen cyanamide, a chemical that induces dormancy-break (Pérez et al., 2008; Ionescu et al., 2017). Our results using field-grown apricot and peach trees found an upregulation of antioxidant defense genes, consistent with the findings from hydrogen cyanamide induced dormancy break (Tang et al., 2019). This supports the hypothesis that ROS induces endodormancy release, which consequently induces bud break. The upregulation of antioxidants and other defense genes likely function as protection from the increase of ROS during dormancy phase transition, which was observed in both apricot and peach buds.

## Flower Development-Related Genes Have Specific Patterns of Expression Corresponding to Endodormancy and Ecodormancy Stages

DEGs and co-expression clusters included a number of genes involved in flower development that are upregulated during ecodormancy. Jojoba acyl CoA reductase-related male sterility proteins (Prupe.7G084900), chalcone and stilbene synthase (Prupe.2G122600 and Prupe.8G159600), Pollen Ole e 1 allergen and extensin family proteins (Prupe.4G061100 and Prupe.3G204300), and gamete expressed protein 1 (Prupe.1G133700) were identified as DEGs in both species, and were clustered into ME11 where gene expression was induced at ecodormancy (Figure 7D, Supplementary Tables
11 and 16). Previous studies reported stamen and microsore development after dormancy break in apricot flower buds, and sweet cherry anthers showed developmental activation after CR fulfillment (Julian et al., 2011; Fadón et al., 2019). Pollen development-related genes have also been shown to be upregulated after dormancy break in peach cultivars (Ríos et al., 2013).

Besides the upregulated DEGs at ecodormancy, genes involved in flower development were also found to be downregulated during ecodormancy. Prupe.1G388300, an ortholog of VERDANDI (VDD), which is a transcription factor that regulates female gametophyte differentiation in *Arabidopsis* (Matias-Hernandez et al., 2010), was differentially expressed in both apricot and peach. VDD was upregulated in endodormancy during chill accumulation and then downregulated at the ecodormancy stage (**Supplementary Table 12**). In *Arabidopsis*, this gene is induced by ovule identity MADS-box transcription factor complexes consisting of SEEDSTICK (STK), SEPALLATA3 (SEP3), and SHATTERPROOF during ovule development (Favaro et al., 2003). However, the expression of VDD in peach and apricot was opposite to the expression of STK, as STK (Prupe.1G549600) was upregulated from dormancy transition to flowering (**Supplementary Table 11** and **Supplementary Figure 12**). This suggests that VDD and female gametogenesis may be regulated by other factors rather than STK complexes during dormancy in *Prunus* spp. Another jojoba acyl CoA reductase-related male sterility protein (Prupe.6G126800) involved in flower wax biosynthesis (Busta and Jetter, 2017) and ovate family protein 2 (OFP2, Prupe.6G290900) involved in vascular formation (Schmitz et al., 2015; Wang et al., 2016) were also downregulated during ecodormancy, indicating that some floral development activities may be completed prior to entering ecodormancy. Julian et al. (2011) reported a rapid development of stamen and vascular differentiation upon CR fulfillment in apricot (Julian et al., 2011). These transcriptional patterns, coupled with our images of developing buds and previous work, demonstrate that floral development progresses in stages closely tied to dormancy progression.

## The DEGs Located in QTL Regions May Be Involved in Dormancy Regulation

Of the shared DEGs for apricot and peach, 99 were located within the CR QTL regions, including three of the *DAM* genes (*DAM4*, *DAM5*, and *DAM6*) in the QTL qCR1a-2009. However, these three *DAM* genes were clustered into two different co-expression modules. *DAM4* fell into ME1, which did not show obvious up or down-regulation during CR accumulation in either peach or apricot, although the *DAM4* gene was downregulated during dormancy release (**Figure 1** and **Supplementary Figure 13**). *DAM5* and *DAM6* were clustered into ME5, where genes were expressed highly at 0 chill hour and decreased during dormancy release (**Supplementary Figure 13**).

Other DEGs overlapping with CR QTLs include genes involved in transcriptional regulation, sugar transportation, stress response, and cell wall development. As a response to abiotic stress, the heat shock transcription factor A-2 (*HSFA2*, Prupe.7G206900) found

on LG 7 QTL qCR7-2008 was upregulated at ecodormancy. *HSFA2* accumulates under heat stress, similar to other heat shock proteins (Schramm et al., 2006). The targets of HSFA2 include the HSP70 (70-kDa heat shock) protein family (Prupe.7G108000 and Prupe.7G108400), which were also upregulated during ecodormancy in apricot. This agrees with a previous study reporting an induction of HSP70 before floral bud break in Japanese pear (*Pyrus pyrifolia*) (Takemura et al., 2015). However, some other heat shock proteins such as Hsp21, Hsp83, and Hsp40 were downregulated during ecodormancy in both apricot and peach, in contrast to the expression patterns of HSP70 (**Supplementary Table 11**).

Three orthologs of sugar transporter protein 1 (Prupe.5G083900, Prupe.5G090900, and Prupe.5G091100 in LG5 QTL qCR5-2008 and qCR5-2009), which are regulated by sugar levels in plants, were significantly upregulated during ecodormancy in both apricot and peach, indicating that the sugar concentrations were likely increased during ecodormancy. This is consistent with a study in Japanese apricot (*Prunus mume*), which found that soluble sugars such as glucose and sucrose were at low levels at the beginning of dormancy, and gradually increased until dormancy release (Zhang et al., 2018). Sugar transportation and synthesis are important to budbreak in Japanese pear, Japanese apricot, and cherry (Marafon et al., 2011; Zhuang et al., 2015; Chmielewski et al., 2017; Fadón et al., 2018a). Dormant cherry buds accumulate starch throughout endodormancy and chilling induces the increase of starch and sucrose concentrations in Japanese pear buds (Marafon et al., 2011; Fadón et al., 2018a). Sugar transport may also have an active role in regulating dormancy transitions, as transgenic poplar expressing *Arabidopsis* sucrose phosphate synthase exhibits early bud flush (Park et al., 2009). The activity of α-amylase which hydrolases starch into sucrose is reported to increase after floral and vegetative buds release from endodormancy (Hussain et al., 2015). To resume growth from dormancy, buds not only use local sugars, but also transport carbohydrates from long distances (Tixier et al., 2019). Tixier et al. (2017) propose long-distance sugar transport in walnut sustains the fast growth resumption after dormancy release (Tixier et al., 2017). In this study two beta-galactosidase genes found in the QTLs qCR7-2008 and qCR-2009 (Fan et al., 2010), *BGAL7* (Prupe.7G194500), and *BGAL1* (Prupe.7G210000), were downregulated during ecodormancy in both apricot and peach. BGAL enzymes produce free galactose, which is the most dynamic sugar residue of the cell wall during fruit development in tomato (Smith and Gross, 2000). Down-regulation of *BGAL* during ecodormancy also supports the idea that the bud is increasing cell wall formation activities during ecodormancy. In *Arabidopsis* thaliana, *BGAL7* and *BGAL8* have been profiled as increasing expression during bud development, particularly in stamens and mature pollen for *BGAL7* (Schmid et al., 2005; Berardini et al., 2015).

## Plant Hormones Are Involved in the Endodormancy to Ecodormancy Transition

Co-expression networks identified genes involved in signaling pathways of growth-related hormones during ecodormancy, including ABA responsive genes and GA responsive genes (**Figure 8A** and **Supplementary Tables 11** and **15**). Two GA
2-oxidase genes (Prupe.4G080700 and Prupe.4G204600) were significantly downregulated during ecodormancy in both apricot and peach (Supplementary Table 11). GA 2-oxidase is the enzyme controlling the bioactive GA levels during plant development. It deactivates GA₁ and GA₄, two major bioactive GAs in plants (Yamaguchi, 2008). Downregulation of GA 2-oxidase genes may increase the bioactive GA levels after the transition from endodormancy to ecodormancy, reactivating plant growth. In *Prunus mume*, GA₃ concentrations were found to increase during the dormancy phase transition, and the low chill genotype had higher GA₃ levels than the high chill genotype during ecodormancy (Wen et al., 2016). The GA biosynthesis gene GA 20-oxidase was also expressed at a high level during ecodormancy (Wen et al., 2016). Unlike *Prunus mume*, there were no significant changes in GA 20-oxidase gene in our study in either peach or apricot, suggesting that the increase of GA levels in peach and apricot may be due to the suppression of GA deactivation rather than an increase of GA biosynthesis. In *Populus*, chilling was found to increase GA biosynthesis and signaling genes, leading to a suggested model where GA is required for dormancy release and acts by inducing removal of callose thus reopening pores and plasmodesmata for signaling to resume growth (Rinne et al., 2011).

ABA has an opposite effect as GA on dormancy. Our co-expression networks among peach and apricot identified three genes involved in ABA metabolic process (GO:0009687, FDR=0.0258), and 15 genes responsive to ABA stimulus (GO:0009737, FDR=1.71E-04); all were upregulated during ecodormancy. Among these 18 genes, two Myb-containing domain genes *Myb108* (Prupe.1G111700) and *Myb102* (Prupe.4G192000) and one homeobox 7 (*HB7*, Prupe.3G316600) were differentially expressed in both species. Past studies reported that *HB7* is a stress response protein positively regulated by the ABA pathway (Söderman et al., 1996; Georgii et al., 2019). The high expression level of *HB7* transcripts during ecodormancy in apricot and peach suggests high concentrations of ABA, which may be acting to repress bud break during ecodormancy. However, the two Myb proteins are shown to inhibit ABA accumulation and down-regulate ABA signaling (Cui et al., 2013; Piao et al., 2019). These genes may indicate counteracting or balancing influences on the ABA pathway during ecodormancy, to promote dormancy release and bud growth.

The small auxin up RNA (SAUR)-like auxin responsive gene family was also significantly upregulated during ecodormancy in both species. The upregulation of six SAUR-like genes (Prupe.2G194600, Prupe.7G167000, Prupe.8G080300, Prupe.8G081100, Prupe.8G081700, and Prupe.8G157800) indicated that auxin levels may be increasing during ecodormancy compared to during endodormancy, consistent with the results found in grapevine that auxin concentrations rise from ecodormancy until budburst (Aloni et al., 1990). Moreover, the downregulation of auxin efflux carrier genes (Prupe.1G071800 and Prupe.5G004300) and upregulation of auxin influx carrier gene (Prupe.1G503100) indicated that increased auxin transporters contribute to increasing auxin levels. Auxin transport has been proposed as a regulatory gate

for dormancy control in apple (Porto et al., 2015). Our analyses identified auxin efflux carrier Prupe.1G071800 as a DEG located on a peach CR QTL qCR1d-2008, making it a strong candidate as a dormancy regulator in peach as well.

## Potential Implications

The seasonal transition from endodormancy to ecodormancy is a complex and tightly coordinated network of many processes, including environmental sensing, cold acclimation, abiotic stress response, hormone fluctuations, chill hour accumulation, and floral development. Examining the transcriptional patterns at multiple time points during dormancy in two species from two locations each with phenotypically varying genotypes enabled us to elucidate strong, reproducible patterns governing dormancy. Our co-expression networks and DEGs begin to untangle this process at the transcript level into two major gene expression patterns shown in **Figure 9**: 1) gene expression was gradually downregulated as dormancy progresses; 2) low gene expression during endodormancy that then peaks at ecodormancy. Functional annotation and gene ontology enrichment analysis defined the genes in the first pattern as involved in pathways including reproduction and chromatin modification. Genes expressed in the second pattern are involved in pollen development, cell wall formation, defense systems, oxidation-reduction, and hormone metabolic and transport pathways. Genes following the first pattern are likely responsible for chill accumulation and release of endodormancy, while genes following the second pattern may be responsible for preparing the bud for the final steps of floral maturation, such as development of pollen. These patterns form a framework of conserved biological pathways that determine winter dormancy and spring flower timing. This transcriptome knowledge contributes to the overall

<figure><img src="image_16.png" /><figcaption>**FIGURE 9** | Diagram of two distinctive expression patterns during dormancy phase transition.</figcaption></figure>
understanding molecular controls of phenological traits, which is particularly critical for tree crops that suffer major crop losses due to warm winters followed by late frosts in spring when trees have already bloomed and are damaged.

## Future Directions

Our analyses using time-series transcriptomic profiles from dormant buds over winter in peach and apricot identified 99 genes located on the CR QTL region and differentially expressed in both species. These candidate genes need to be examined in mapping populations for sequence or structural variants that might yield phenotypic differences and their transcriptional profiles examined against more dormancy time points. Further, the same material needs to be experimentally verified for the specific endo-vs-ecodormancy transition point. The biological pathways identified by the co-expression networks also point to chromatin remodeling and cell cycle pathways during chilling in the endodormancy stage as likely factors regulating dormancy transition and warrant additional research. The finding of increasing defense response at ecodormancy only in peach also needs further validation in other *Prunus* species.

## DATA AVAILABILITY STATEMENT

All reported raw read data is available at the National Center for Biotechnology Information (NCBI) under the study of project accession PRJNA567655. Raw reads for apricot samples are available from the Sequence Read Archive (SRA) with accession IDs SAMN12791244 to SAMN12791303. Peach samples are available under the accession IDs SAMN12791304 to SAMN12791372. The raw counts of reads per gene derived after read mapping were deposited to NCBI Gene Expression Omnibus (GEO) with accession GSE138792.

## REFERENCES

* Alburquerque, N., García-Montiel, F., Carrillo, A., and Burgos, L. (2008). Chilling and heat requirements of sweet cherry cultivars and the relationship between altitude and the probability of satisfying the chill requirements. *Environ. Exp. Bot.* 64, 162–170. doi: 10.1016/j.envexpbot. 2008.01.003

* Allona, I., Ramos, A., Ibáñez, C., Contreras, A., Casado, R., and Aragoncillo, C. (2008). Molecular control of winter dormancy establishment in trees: a review. *Span. J. Agric. Res.* 6, 201–210. doi: 10.5424/sjar/200806S1-389

* Aloni, R., Baum, S. F., and Peterson, C. A. (1990). The role of cytokinin in sieve tube regeneration and callose production in wounded coleus internodes. *Plant Physiol.* 93, 982–989. doi: 10.1104/pp.93.3.982

* Anderson, J. V., and Davis, D. G. (2004). Abiotic stress alters transcript profiles and activity of glutathione S-transferase, glutathione peroxidase, and glutathione reductase in Euphorbia esula. *Physiol. Plant* 120, 421–433. doi: 10.1111/j.0031-9317.2004.00249.x

* Anderson, J. L., Richardson, E. A., and Kesner, C. D. (1985). 'Validation of chill unit and flower bud phenology models for "Montmorency" sour cherry,' in *I International Symposium on Computer Modelling in Fruit Research and Orchard Management. Acta Hortic.* 184, 71–78. doi: 10.17660/ActaHortic.1986.184.7

* Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. *Austral Ecol.* 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x

## AUTHOR CONTRIBUTIONS

AA, CD, ZL, MS, TZ, and VD conceived the project. VD, GR, and J-MA collected apricot samples and VD and AA extracted RNA. AC and TZ collected peach samples and CD, ZL, and DB extracted RNA. ZL and DB provided the peach bud imaging. DW performed RT-qPCR. JY and MS performed the data analyses and wrote the manuscript with contributions from all co-authors. All authors approved the manuscript.

## FUNDING

This project is supported by the Agriculture and Food Research Initiative Competitive Grant No. 2016-67014-24577 from the USDA National Institute of Food and Agriculture (Albert G Abbott, Christopher Dardick, Zongrang Liu, Margaret E Staton) and the French ANR CHEX ABRIWG No. ANR-11-CHEX-0002 (Albert G Abbott, Véronique Decroocq, and Jean-Marc Audergon).

## ACKNOWLEDGMENTS

We thank Dr. Shenghua Fan (Forest Health Research and Education Center, University of Kentucky, USA) for the assistance with Prunus QTL maps and markers.

## SUPPLEMENTARY MATERIAL

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

* Anderson, M. J. (2017). Permutational Multivariate Analysis of Variance (PERMANOVA). In *Wiley StatsRef: Statistics Reference Online* (eds N. Balakrishnan, T. Colton, B. Everitt, W. Piegorsch, F. Ruggeri and J. L. Teugels) 1–15. doi: 10.1002/9781118445112.stat07841

* Andreini, L., de Cortázar-Atauri, I. G., Chuine, I., Viti, R., Bartolini, S., Ruiz, D., et al. (2014). Understanding dormancy release in apricot flower buds (Prunus armeniaca L.) using several process-based phenological models. *Agric. For. Meteorol.* 184, 210–219. doi: 10.1016/j.agrformet.2013.10.005

* Andrews, S. (2010). "Babraham bioinformatics-FastQC a quality control tool for high throughput sequence data." https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

* Anh Tuan, P., Bai, S., Saito, T., Imai, T., Ito, A., and Moriguchi, T. (2016). Involvement of EARLY BUD-BREAK, an AP2/ERF transcription factor gene, in bud break in Japanese pear (Pyrus pyrifolia Nakai) lateral flower buds: expression, histone modifications and possible target genes. *Plant Cell Physiol.* 57, 1038–1047. doi: 10.1093/pcp/pcw041

* Arús, P., Yamamoto, T., Dirlewanger, E., and Abbott, A. G. (2006). Synteny in the rosaceae. *Plant Breed. Rev.* 27, 175.

* Artlip, T., McDermaid, A., Ma, Q., and Wisniewski, M. (2019). Differential gene expression in non-transgenic and transgenic "M.26" apple overexpressing a peach CBF gene during the transition from eco-dormancy to bud break. *Hortic. Res.* 6, 1–16. doi: 10.1038/s41438-019-0168-9

* Audergon, J. M. (1993). 'Variety and breeding' in *X International Symposium on*
Apricot Culture. Acta Hortic. 384, 35-46. doi: 10.17660/ActaHortic.1995.384.2

Bai, S., Saito, T., Sakamoto, D., Ito, A., Fujii, H., and Moriguchi, T. (2013). Transcriptome analysis of Japanese pear (Pyrus pyrifolia Nakai) flower buds transitioning through endodormancy. *Plant Cell Physiol.* 54, 1132-1151. doi: 10.1093/pcp/pct067

Bai, S., Saito, T., Ito, A., Tuan, P. A., Xu, Y., Teng, Y., et al. (2016). Small RNA and PARE sequencing in flower bud reveal the involvement of sRNAs in endodormancy release of Japanese pear (Pyrus pyrifolia "Kosui"). *BMC Genomics* 17, 12. doi: 10.1186/s12864-016-2514-8

Balogh, E., Halász, J., Soltész, A., Erős-Honti, Z., Gutermuth, Á., Szalay, L., et al. (2019). Identification, structural and functional characterization of dormancy regulator genes in apricot (Prunus armeniaca L.). *Front. Plant Sci.* 10, 402. doi: 10.3389/fpls.2019.00402

Benmoussa, H., Ghrab, M., Ben Mimoun, M., and Luedeling, E. (2017). Chilling and heat requirements for local and foreign almond (Prunus dulcis Mill.) cultivars in a warm Mediterranean location based on 30 years of phenology records. *Agric. For. Meteorol.* 239, 34-46. doi: 10.1016/j.agrformet.2017.02.030

Berardini, T. Z., Reiser, L., Li, D., Mezheritsky, Y., Muller, R., Strait, E., et al. (2015). The arabidopsis information resource: Making and mining the "gold standard" annotated reference plant genome. *genesis* 53, 474-485. doi: 10.1002/dvg.22877

Bielenberg, D. G., Wang, Y., Fan, S., Reighard, G. L., Scorza, R., and Abbott, A. G. (2004). A deletion affecting several gene candidates is present in the Evergrowing peach mutant. *J. Hered.* 95, 436-444. doi: 10.1093/jhered/esh057

Bielenberg, D. G., Wang, Y. E., Li, Z., Zhebentyayeva, T., Fan, S., Reighard, G. L., et al. (2008). Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. *Tree Genet. Genomes* 4, 495-507. doi: 10.1007/s11295-007-0126-9

Busov, V. B. (2019). Plant development: dual roles of poplar svl in vegetative bud dormancy. *Curr. Biol.* 29, R68-R70. doi: 10.1016/j.cub.2018.11.061

Busta, L., and Jetter, R. (2017). Structure and biosynthesis of branched wax compounds on wild type and wax biosynthesis mutants of arabidopsis thaliana. *Plant Cell Physiol.* 58, 1059-1074. doi: 10.1093/pcp/pcx051

Campbell, M., Suttle, J., Douches, D. S., and Buell, C. R. (2014). Treatment of potato tubers with the synthetic cytokinin 1-(*α*-ethylbenzyl)-3-nitroguanidine results in rapid termination of endodormancy and induction of transcripts associated with cell proliferation and growth. *Funct. Integr. Genomics* 14, 789-799. doi: 10.1007/s10142-014-0404-1

Campoy, J. A., Ruiz, D., and Egea, J. (2011). Dormancy in temperate fruit trees in a global warming context: a review. *Sci. Hortic.* 130, 357-372. doi: 10.1016/j.scienta.2011.07.011

Caruso, T., Giovannini, D., Marra, F. P., and Sottile, F. (1997). Two new planting systems for early ripening peaches (Prunus persica L. Batsch): yield and fruit quality in four low-chill cultivars. *J. Hortic. Sci.* 72, 873-883. doi: 10.1080/14620316.1997.11515579

Castède, S., Campoy, J. A., Garcia, J. Q., Le Dantec, L., Lafargue, M., Barreneche, T., et al. (2014). Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements. *New Phytol.* 202, 703-715. doi: 10.1111/nph.12658

Čechová, J., Baránek, M., Krška, B., and Pidra, M. (2012). Screening of differentially expressed genes during the end of endogenous dormancy of flower buds in Prunus armeniaca L. *Plant Growth Regul.* 67, 141-150. doi: 10.1007/s10725-012-9671-9

Cesaraccio, C., Spano, D., Snyder, R. L., and Duce, P. (2004). Chilling and forcing model to predict bud-burst of crop and forest species. *Agric. For. Meteorol.* 126, 1-13. doi: 10.1016/j.agrformet.2004.03.002

Chandra Babu, R., Zhang, J., Blum, A., David Ho, T.-H., Wu, R., and Nguyen, H. T. (2004). HVA1, a LEA gene from barley confers dehydration tolerance in transgenic rice (Oryza sativa L.) *via* cell membrane protection. *Plant Sci.* 166, 855-862. doi: 10.1016/j.plantsci.2003.11.023

Chavarria, G., Herter, F. G., do Carmo Bassols Raseira, M., Rodrigues, A. C., Reisser, C., and da Silva, J. B. (2009). Mild temperatures on bud breaking dormancy in peaches. *Ciência Rural* 39, 2016-2021. doi: 10.1590/s0103-84782009000700010

Checker, V. G., Chhibbar, A. K., and Khurana, P. (2012). Stress-inducible expression of barley Hval gene in transgenic mulberry displays enhanced tolerance against drought, salinity and cold stress. *Transgenic Res.* 21, 939-957. doi: 10.1007/s11248-011-9577-8

Chen, W., Yu, X.-H., Zhang, K., Shi, J., De Oliveira, S., Schreiber, L., et al. (2011). Male sterile2 encodes a plastid-localized fatty acyl carrier protein reductase required for pollen exine development in arabidopsis. *Plant Physiol.* 157, 842-853. doi: 10.1104/pp.111.181693

Chinnusamy, V., Zhu, J., and Zhu, J.-K. (2007). Cold stress regulation of gene expression in plants. *Trends Plant Sci.* 12, 444-451. doi: 10.1016/j.tplants.2007.07.002

Chinnusamy, V., Zhu, J.-K., and Sunkar, R. (2010). "Gene Regulation During Cold Stress Acclimation in Plants," in *Plant Stress Tolerance: Methods and Protocols*. Ed. R. Sunkar (Totowa, NJ: Humana Press), 39-55.

Chmielewski, F. M., Gotz, K. P., Homann, T., Huschek, G., and Rawel, H. M. (2017). Identification of endodormancy release for cherries (Prunus Avium L.) by abscisic acid and sugars. *J. Hortic.* 4, 3. doi: 10.4172/2376-0354.1000210

Citadin, I., Raseira, M. C. B., Herter, F. G., and da Silva, J. B. (2001). Heat requirement for blooming and leafing in peach. *HortScience* 36, 305-307. doi: 10.21273/HORTSCI.36.2.305

Conrad, A. O., Yu, J., Staton, M. E., Audergon, J.-M., Roch, G., Decroocq, V., et al. (2019). Association of the phenylpropanoid pathway with dormancy and adaptive trait variation in apricot (Prunus armeniaca). *Tree Physiol.* 39 (7), 1136-1148. doi: 10.1093/treephys/tpz053

Considine, M. J., and Considine, J. A. (2016). On the language and physiology of dormancy and quiescence in plants. *J. Exp. Bot.* 67, 3189-3203. doi: 10.1093/jxb/erw138

Considine, M. J., and Foyer, C. H. (2014). Redox regulation of plant development. *Antioxid. Redox Signal.* 21, 1305-1326. doi: 10.1089/ars.2013.5665

Cooke, J. E. K., Eriksson, M. E., and Junttila, O. (2012). The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms. *Plant Cell Environ.* 35, 1707-1728. doi: 10.1111/j.1365-3040.2012.02552.x

Cucinotta, M., Manrique, S., Guazzotti, A., Quadrelli, N. E., Mendes, M. A., Benkova, E., et al. (2016). Cytokinin response factors integrate auxin and cytokinin pathways for female reproductive organ development. *Development* 143, 4419-4424. doi: 10.1242/dev.143545

Cui, F., Brosché, M., Sipari, N., Tang, S., and Overmyer, K. (2013). Regulation of ABA dependent wound induced spreading cell death by MYB108. *New Phytol.* 200, 634-640. doi: 10.1111/nph.12456

da Silveira Falavigna, V., Guitton, B., Costes, E., and Andrés, F. (2018). I want to (bud) break free: the potential role of DAM and SVP2 like genes in regulating dormancy cycle in temperate fruit trees. *Front. Plant Sci.* 9, 1990. doi: 10.3389/fpls.2018.01990

Dewitte, W., Scofield, S., Alcasabas, A. A., Maughan, S. C., Menges, M., Braun, N., et al. (2007). Arabidopsis CYCD3 D-type cyclins link cell proliferation and endocycles and are rate-limiting for cytokinin responses. *Proc. Natl. Acad. Sci. U. S. A.* 104, 14537-14542. doi: 10.1073/pnas.0704166104

Dirlewanger, E., Graziano, E., Joobeur, T., Garriga-Calderé, F., Cosson, P., Howad, W., et al. (2004). Comparative mapping and marker-assisted selection in Rosaceae fruit crops. *Proc. Natl. Acad. Sci. U. S. A.* 101, 9891-9896. doi: 10.1073/pnas.0307937101

Dirlewanger, E., Quero-García, J., Le Dantec, L., Lambert, P., Ruiz, D., Dondini, L., et al. (2012). Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry. *Heredity* 109, 280-292. doi: 10.1038/hdy.2012.38

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: ultrafast universal RNA-seq aligner. *Bioinformatics* 29, 15-21. doi: 10.1093/bioinformatics/bts635

Dobritsa, A. A., Lei, Z., Nishikawa, S.-I., Urbanczyk-Wochniak, E., Huhman, D. V., Preuss, D., et al. (2010). LAP5 and LAP6 encode anther-specific proteins with similarity to chalcone synthase essential for pollen exine development in Arabidopsis. *Plant Physiol.* 153, 937-955. doi: 10.1104/pp.110.157446

Duan, J., and Cai, W. (2012). OsLEA3-2, an abiotic stress induced gene of rice plays a key role in salt and drought tolerance. *PloS One* 7, e45117. doi: 10.1371/journal.pone.0045117

Erez, A. (1979). Quantitative chilling enhancement and negation in peach bud by high temperatures in a daily cycle. *J. Am. Soc Hortic. Sci.* 104, 536-540.

Fadón, E., Herrero, M., and Rodrigo, J. (2018a). Dormant flower buds actively accumulate starch over winter in sweet cherry. *Front. Plant Sci.* 9, 171. doi: 10.3389/fpls.2018.00171

Fadón, E., Rodrigo, J., and Herrero, M. (2018b). Is there a specific stage to rest? Morphological changes in flower primordia in relation to endodormancy in sweet cherry (Prunus avium L.). *Trees* 32, 1583-1594. doi: 10.1007/s00468-018-1735-7
Fadón, E., Herrero, M., and Rodrigo, J. (2019). Anther and pollen development in sweet cherry (Prunus avium L.) in relation to winter dormancy. *Protoplasma* 256, 733-744. doi: 10.1007/s00709-018-01332-4

Fan, S., Bielenberg, D. G., Zhebentayeva, T. N., Reighard, G. L., Okie, W. R., Holland, D., et al. (2010). Mapping quantitative trait loci associated with chilling requirement, heat requirement and bloom date in peach (Prunus persica). *New Phytol.* 185, 917-930. doi: 10.1111/j.1469-8137.2009.03119.x

Favaro, R., Pinyopich, A., Battaglia, R., Kooiker, M., Borghi, L., Ditta, G., et al. (2003). MADS-box protein complexes control carpel and ovule development in Arabidopsis. *Plant Cell* 15, 2603-2611. doi: 10.1105/tpc.015123

Foley, M. E., Anderson, J. V., and Horvath, D. P. (2009). The effects of temperature, photoperiod, and vernalization on regrowth and flowering competence in Euphorbia esula (Euphorbiaceae) crown buds. *Botany* 87, 986-992. doi: 10.1139/B09-055

Georgii, E., Kugler, K., Pfeifer, M., Vanzo, E., Block, K., Domagalska, M. A., et al. (2019). The systems architecture of molecular memory in poplar after abiotic stress. *Plant Cell* 31, 346-367. doi: 10.1105/tpc.18.00431

González, L. M. G., El Kayal, W., Morris, J. S., and Cooke, J. E. K. (2015). Diverse chitinases are invoked during the activity-dormancy transition in spruce. *Tree Genet. Genomes* 11, 41. doi: 10.1007/s11295-015-0871-0

Goodstein, D. M., Shu, S., Howson, R., Neupane, R., Hayes, R. D., Fazo, J., et al. (2012). Phytozome: a comparative platform for green plant genomics. *Nucleic Acids Res.* 40, D1178–D1186. doi: 10.1093/nar/gkr944

Hänninen, H. (1990). Modelling bud dormancy release in trees from cool and temperate regions. *Acta Forestalia Fennica* 47. doi: 10.14214/aff.7660

Hara, M., Wakasugi, Y., Ikoma, Y., Yano, M., Ogawa, K., and Kuboi, T. (1999). cDNA sequence and expression of a cold-responsive gene in Citrus unshiu. *Biosci. Biotechnol. Biochem.* 63, 433-437. doi: 10.1271/bbb.63.433

Hatfield, J. L., and Prueger, J. H. (2015). Temperature extremes: effect on plant growth and development. *Weather Climate Extremes* 10, 4–10. doi: 10.1016/j.wace.2015.08.001

Horvath, D. P., Anderson, J. V., Chao, W. S., and Foley, M. E. (2003). Knowing when to grow: signals regulating bud dormancy. *Trends Plant Sci.* 8, 534–540. doi: 10.1016/j.tplants.2003.09.013

Hussain, S., Niu, Q., Yang, F., Hussain, N., and Teng, Y. (2015). The possible role of chilling in floral and vegetative bud dormancy release in Pyrus pyrifolia. *Biol. Plant.* 59, 726-734. doi: 10.1007/s10535-015-0547-5

Ionescu, I. A., López-Ortega, G., Burow, M., Bayo-Canha, A., Junge, A., Gericke, O., et al. (2017). Transcriptome and metabolite changes during hydrogen cyanamide-induced floral bud break in sweet cherry. *Front. Plant Sci.* 8, 1233. doi: 10.3389/fpls.2017.01233

Jain, M., Ghanashyam, C., and Bhattacharjee, A. (2010). Comprehensive expression analysis suggests overlapping and specific roles of rice glutathione S-transferase genes during development and stress responses. *BMC Genomics* 11, 73. doi: 10.1186/1471-2164-11-73

Jeon, J., Cho, C., Lee, M. R., Van Binh, N., and Kim, J. (2016). Cytokinin response factor2 (CRF2) and CRF3 regulate lateral root development in response to cold stress in arabidopsis. *Plant Cell* 28, 1828-1843. doi: 10.1105/tpc.15.00909

Jian, L.-C., Li, P. H., Sun, L.-H., and Chen, T. H. H. (1997). Alterations in ultrastructure and subcellular localization of Ca2+ in poplar apical bud cells during the induction of dormancy. *J. Exp. Bot.* 48, 1195-1207. doi: 10.1093/jxb/48.6.1195

Jiang, H., Lei, R., Ding, S.-W., and Zhu, S. (2014). Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. *BMC Bioinf.* 15, 182. doi: 10.1186/1471-2105-15-182

Jiménez, S., Li, Z., Reighard, G. L., and Bielenberg, D. G. (2010a). Identification of genes associated with growth cessation and bud dormancy entrance using a dormancy-incapable tree mutant. *BMC Plant Biol.* 10, 25. doi: 10.1186/1471-2229-10-25

Jiménez, S., Reighard, G. L., and Bielenberg, D. G. (2010b). Gene expression of DAM5 and DAM6 is suppressed by chilling temperatures and inversely correlated with bud break rate. *Plant Mol. Biol.* 73, 157-167. doi: 10.1007/s11103-010-9608-5

Julian, C., Rodrigo, J., and Herrero, M. (2011). Stamen development and winter dormancy in apricot (Prunus armeniaca). *Ann. Bot.* 108, 617–625. doi: 10.1093/aob/mcr056

Jung, S., Jiwan, D., Cho, I., Lee, T., Abbott, A., Sosinski, B., et al. (2009). Synteny of Prunus and other model plant species. *BMC Genomics* 10, 76. doi: 10.1186/1471-2164-10-76

Kaye, C., Neven, L., Hofig, A., Li, Q. B., Haskell, D., and Guy, C. (1998). Characterization of a gene for spinach CAP160 and expression of two

spinach cold-acclimation proteins in tobacco. *Plant Physiol.* 116, 1367-1377. doi: 10.1104/pp.116.4.1367

Kitamura, Y., Yamane, H., Yukimori, A., Shimo, H., Numaguchi, K., and Tao, R. (2017). Blooming date predictions based on japanese apricot "nanko" flower bud responses to temperatures during dormancy. *HortScience* 52, 366-370. doi: 10.21273/hortsci11253-16

Kuroda, H., Sugiura, T., and Sugiura, H. (2005). Effect of hydrogen peroxide on breaking endodormancy in flower buds of japanese pear (Pyrus pyrifolia Nakai). *J. Jpn. Soc. Hortic. Sci.* 74, 255-257. doi: 10.2503/jjshs.74.255

Lang, A. G., Early, J. D., Martin, G. C., and Darnell, R. L. (1987). Endo-, para- and ecodormancy: Physiological terminology and classification for dormancy research. *HortScience* 22, 371-377.

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

Langfelder, P., and Horvath, S. (2014). *Tutorials for the WGCNA Package*, Available at: https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/Tutorials/.

Leida, C., Terol, J., Martí, G., Agustí, M., Llácer, G., Badenes, M. L., et al. (2010). Identification of genes associated with bud dormancy release in Prunus persica by suppression subtractive hybridization. *Tree Physiol.* 30, 655-666. doi: 10.1093/treephys/tpq008

Leida, C., Conejero, A., Arbona, V., Gómez-Cadenas, A., Llácer, G., Badenes, M. L., et al. (2012a). Chilling-dependent release of seed and bud dormancy in peach associates to common changes in gene expression. *PloS One* 7, e35777. doi: 10.1371/journal.pone.0035777

Leida, C., Conesa, A., Llácer, G., Badenes, M. L., and Ríos, G. (2012b). Histone modifications and expression of DAM6 gene in peach are modulated during bud dormancy release in a cultivar-dependent manner. *New Phytol.* 193, 67–80. doi: 10.1111/j.1469-8137.2011.03863.x

Leida, C., Romeu, J. F., García-Brunton, J., Ríos, G., and Badenes, M. L. (2012c). Gene expression analysis of chilling requirements for flower bud break in peach. *Plant Breed.* 131, 329-334. doi: 10.1111/j.1439-0523.2011.01946.x

Leida, C. A. (2012). Molecular aspects of dormancy in peach (Prunus persica [L.] Batsch.) [Tesis doctoral no publicada]. Universitat Politècnica de València. doi: 10.4995/Thesis/10251/15864

Li, Z., Reighard, G. L., Abbott, A. G., and Bielenberg, D. G. (2009). Dormancy-associated MADS genes from the EVG locus of peach [Prunus persica (L.) Batsch] have distinct seasonal and photoperiodic expression patterns. *J. Exp. Bot.* 60, 3521-3530. doi: 10.1093/jxb/erp195

Li, Y., Cao, K., Zhu, G., Fang, W., Chen, C., Wang, X., et al. (2019). Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. *Genome Biol.* 20, 36. doi: 10.1186/s13059-019-1648-9

Liu, J., and Sherif, S. M. (2019). Combating spring frost with ethylene. *Front. Plant Sci.* 10, 1408. doi: 10.3389/fpls.2019.01408

Lloret, A., Badenes, M. L., and Ríos, G. (2018). Modulation of dormancy and growth responses in reproductive buds of temperate trees. *Front. Plant Sci.* 9, 1368. doi: 10.3389/fpls.2018.01368

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* 15, 550. doi: 10.1186/s13059-014-0550-8

Luedeling, E., and Brown, P. H. (2011). A global analysis of the comparability of winter chill models for fruit and nut trees. *Int. J. Biometeorol.* 55, 411-421. doi: 10.1007/s00484-010-0352-y

Luedeling, E., Girvetz, E. H., Semenov, M. A., and Brown, P. H. (2011). Climate change affects winter chill for temperate fruit and nut trees. *PloS One* 6, e20155. doi: 10.1371/journal.pone.0020155

Luedeling, E. (2012). Climate change impacts on winter chill for temperate fruit and nut production: a review. *Sci. Hortic.* 144, 218-229. doi: 10.1016/j.scienta.2012.07.011

Maere, S., Heymans, K., and Kuiper, M. (2005). BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. *Bioinformatics* 21, 3448-3449. doi: 10.1093/bioinformatics/bti551

Marafon, A. C., Citadin, I., do Amarante, L., Herter, F. G., and Hawerroth, F. J. (2011). Chilling privation during dormancy period and carbohydrate mobilization in Japanese pear trees. *Sci. Agricola* 68, 462-468. doi: 10.1590/s0103-90162011000400011

Matias-Hernandez, L., Battaglia, R., Galbiati, F., Rubes, M., Eichenberger, C., Grossniklaus, U., et al. (2010). VERDANDI is a direct target of the MADS
domain ovule identity complex and affects embryo sac differentiation in
Arabidopsis. *Plant Cell* 22, 1702–1715. doi: 10.1105/tpc.109.068627

Maurya, J. P., and Bhalerao, R. P. (2017). Photoperiod- and temperature-mediated control of growth cessation and dormancy in trees: a molecular perspective. *Ann. Bot.* 120, 351–360. doi: 10.1093/aob/mcx061

Medina, J., Catalá, R., and Salinas, J. (2011). The CBFs: three arabidopsis transcription factors to cold acclimate. *Plant Sci.* 180, 3–11. doi: 10.1016/j.plantsci2010.06.019

Muñoz, P., and Munné-Bosch, S. (2018). Photo-Oxidative Stress during Leaf, Flower and Fruit Development. *Plant Physiol.* 176, 1004–1014. doi: 10.1104/pp.17.01127

Niu, Q., Li, J., Cai, D., Qian, M., Jia, H., Bai, S., et al. (2016). Dormancy-associated MADS-box genes and microRNAs jointly control dormancy transition in pear (Pyrus pyrifolia white pear group) flower bud. *J. Exp. Bot.* 67, 239–257. doi: 10.1093/jxb/erv454

Nylander, M., Svensson, J., Palva, E. T., and Welin, B. V. (2001). Stress-induced accumulation and tissue-specific localization of dehydrins in Arabidopsis thaliana. *Plant Mol. Biol.* 45, 263–279. doi: 10.1023/A:1006469128280

Oksanen, J., Guillaume Blanchet, F., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). *Vegan: community ecology package*. R package version 2.5-6, Available at: https://CRAN.R-project.org/package=vegan.

Olukolu, B. A., Trainin, T., Fan, S., Kole, C., Bielenberg, D. G., Reighard, G. L., et al. (2009). Genetic linkage mapping for molecular dissection of chilling requirement and budbreak in apricot (Prunus armeniaca L.). *Genome* 52, 819–828. doi: 10.1139/G09-050

Pérez, F. J., Vergara, R., and Rubio, S. (2008). H2O2 is involved in the dormancy-breaking effect of hydrogen cyanamide in grapevine buds. *Plant Growth Regul.* 55, 149–155. doi: 10.1007/s10725-008-9269-4

Park, J.-Y., Canam, T., Kang, K.-Y., Unda, F., and Mansfield, S. D. (2009). Sucrose phosphate synthase expression influences poplar phenology. *Tree Physiol.* 29, 937–946. doi: 10.1093/treephys/tpp028

Pedrosa, A. M., Martins, C., de, P. S., Gonçalves, L. P., and Costa, M. G. C. (2015). Late embryogenesis abundant (lea) constitutes a large and diverse family of proteins involved in development and abiotic stress responses in sweet orange (Citrus sinensis L. Osb.). *PloS One* 10, e0145785. doi: 10.1371/journal.pone.0145785

Perry, T. O. (1971). Dormancy of trees in winter. *Science* 171, 29–36. doi: 10.1126/science.171.3966.29

Piao, W., Kim, S.-H., Lee, B.-D., An, G., Sakuraba, Y., and Paek, N.-C. (2019). Rice transcription factor OsMYB102 delays leaf senescence by down-regulating abscisic acid accumulation and signaling. *J. Exp. Bot.* 70, 2699–2715. doi: 10.1093/jxb/erz095

Porto, D. D., Bruneau, M., Perini, P., Anzanello, R., Renou, J.-P., dos Santos, H. P., et al. (2015). Transcription profiling of the chilling requirement for bud break in apples: a putative role for FLC-like genes. *J. Exp. Bot.* 66, 2659–2672. doi: 10.1093/jxb/erv061

Prassinos, C., Rigas, S., Kizis, D., Vlahou, A., and Hatzopoulos, P. (2011). Subtle proteome differences identified between post-dormant vegetative and floral peach buds. *J. Proteomics* 74, 607–619. doi: 10.1016/j.jprot.2011.01.018

Pyke, K. A., and Page, A. M. (1998). Plastid ontogeny during petal development in Arabidopsis. *Plant Physiol.* 116, 797–803. doi: 10.1104/pp.116.2.797

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. *Nucleic Acids Res.* 41, D590–D596. doi: 10.1093/nar/gks1219

Quero-Garcia, J., Roch, G., Charlot, G., Hilaire, C., Barreneche, T., Campoy-Corbalan, J. A., et al. (2016). *Etude des déterminants génétiques et moléculaires de caractères influencés par le changement climatique chez l'abricotier et le cerisier.* Available at: https://hal.archives-ouvertes.fr/hal-01578240/.

Ríos, G., Tadeo, F. R., Leida, C., and Badenes, M. L. (2013). Prediction of components of the sporopollenin synthesis pathway in peach by genomic and expression analyses. *BMC Genomics* 14, 40. doi: 10.1186/1471-2164-14-40

Reinoso, H., Luna, V., Pharis, R. P., and Bottini, R. (2002). Dormancy in peach (Prunus persica) flower buds. V. Anatomy of bud development in relation to phenological stage. *Can. J. Bot.* 80, 656–663. doi: 10.1139/b02-052

Richardson, E. A. (1974). A model for estimating the completion of rest for "Redhaven" and "Elberta" peach trees. *HortScience* 9, 331–332.

Rieger, M. (1989). Freeze protection for horticultural crops. *Hortic. Rev.* 11, 45–109. doi: 10.1002/9781118060841.ch3

Rinne, P. L. H., Welling, A., Vahala, J., Ripel, L., Ruonala, R., Kangasjärvi, J., et al. (2011). Chilling of dormant buds hyperinduces flowering locus t and recruits ga-inducible 1,3-β-glucanases to reopen signal conduits and release dormancy in populus. *Plant Cell* 23, 130–146. doi: 10.1105/tpc.110.081307

Rohde, A., and Bhalerao, R. P. (2007). Plant dormancy in the perennial context. *Trends Plant Sci.* 12, 217–223. doi: 10.1016/j.tplants.2007.03.012

Ruonala, R., Rinne, P. L. H., Baghour, M., Moritz, T., Tuominen, H., and Kangasjärvi, J. (2006). Transitions in the functioning of the shoot apical meristem in birch (Betula pendula) involve ethylene. *Plant J.* 46, 628–640. doi: 10.1111/j.1365-313X.2006.02722.x

Ruttink, T., Arend, M., Morreel, K., Storme, V., Rombauts, S., Fromm, J., et al. (2007). A molecular timetable for apical bud formation and dormancy induction in poplar. *Plant Cell* 19, 2370–2390. doi: 10.1105/tpc.107.052811

Sánchez-Pérez, R., Dicenta, F., and Martinez-Gómez, P. (2012). Inheritance of chilling and heat requirements for flowering in almond and QTL analysis. *Tree Genet. Genomes* 8, 379–389. doi: 10.1007/s11295-011-0448-5

Söderman, E., Mattsson, J., and Engström, P. (1996). The Arabidopsis homeobox gene ATHB-7 is induced by water deficit and by abscisic acid. *Plant J.* 10, 375–381. doi: 10.1046/j.1365-313X.1996.10020375.x

Saito, T., Bai, S., Imai, T., Ito, A., Nakajima, I., and Moriguchi, T. (2015a). Histone modification and signalling cascade of the dormancy-associated MADS-box gene, PpMADS13-1, in Japanese pear (Pyrus pyrifolia) during endodormancy. *Plant Cell Environ.* 38, 1157–1166. doi: 10.1111/pc.12469

Saito, T., Tuan, P. A., Katsumi-Horigane, A., Bai, S., Ito, A., Sekiyama, Y., et al. (2015b). Development of flower buds in the Japanese pear (Pyrus pyrifolia) from late autumn to early spring. *Tree Physiol.* 35, 653–662. doi: 10.1093/treephys/tpv043

Schmid, M., Davison, T. S., Henz, S. R., Pape, U. J., Demar, M., Vingron, M., et al. (2005). A gene expression map of Arabidopsis thaliana development. *Nat. Genet.* 37, 501–506. doi: 10.1038/ng1543

Schmitz, A. J., Begcy, K., Sarath, G., and Walia, H. (2015). Rice Ovate Family Protein 2 (OFP2) alters hormonal homeostasis and vasculature development. *Plant Sci.* 241, 177–188. doi: 10.1016/j.plantsci.2015.10.011

Schramm, F., Ganguli, A., Kiehlmann, E., Englich, G., Walch, D., and von Koskull-Döring, P. (2006). The heat stress transcription factor HsfA2 serves as a regulatory amplifier of a subset of genes in the heat stress response in Arabidopsis. *Plant Mol. Biol.* 60, 759–772. doi: 10.1007/s11103-005-5750-x

Seppänen, M. M., Cardi, T., Borg Hyökki, M., and Pehu, E. (2000). Characterization and expression of cold-induced glutathione S-transferase in freezing tolerant Solanum commersonii, sensitive S. tuberosum and their interspecific somatic hybrids. *Plant Sci.* 153, 125–133. doi: 10.1016/S0168-9452(99)00252-6

Shaltout, A. D. (1983). Rest completion prediction model for "Starkrimson Delicious apples." *J. Am. Soc Hortic. Sci.* 108, 957–961.

Shi, Y., Ding, Y., and Yang, S. (2018). Molecular Regulation of CBF Signaling in Cold Acclimation. *Trends Plant Sci.* 23, 623–637. doi: 10.1016/j.tplants.2018.04.002

Singh, R. K., Maurya, J. P., Azeez, A., Miskolczi, P., Tylewicz, S., Stojković, K., et al. (2018). A genetic network mediating the control of bud break in hybrid aspen. *Nat. Commun.* 9, 4173. doi: 10.1038/s41467-018-06696-y

Singh, R. K., Miskolczi, P., Maurya, J. P., and Bhalerao, R. P. (2019). A Tree Ortholog of short vegetative phase floral repressor mediates photoperiodic control of bud dormancy. *Curr. Biol.* 29, 128–133.e2. doi: 10.1016/j.cub.2018.11.006

Smith, D. L., and Gross, K. C. (2000). A family of at least seven beta-galactosidase genes is expressed during tomato fruit development. *Plant Physiol.* 123, 1173–1183. doi: 10.1104/pp.123.3.1173

Stockinger, E. J., Gilmour, S. J., and Thomashow, M. F. (1997). Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. *Proc. Natl. Acad. Sci.* 94, 1035–1040. doi: 10.1073/pnas.94.3.1035

Takemura, Y., Kuroki, K., Jiang, M., Matsumoto, K., and Tamura, F. (2015). Identification of the expressed protein and the impact of change in ascorbate peroxidase activity related to endodormancy breaking in Pyrus pyrifolia. *Plant Physiol. Biochem.* 86, 121–129. doi: 10.1016/j.plaphy.2014.11.016

Tan, Y., Yue, T. A. N., Dong-sheng, G. A. O., Ling, L. I., Chen, X., and Ai-hong, X. U. (2010). Respiratory response of dormant nectarine floral buds on chilling deficiency. *Agric. Sci. In China* 9, 1448–1454. doi: 10.1016/s1671-2927(09)60236-4

Tang, L., Chhajed, S., Vashist, T., Olmstead, M. A., Olmstead, J. W., and Colquhoun, T. A. (2019). Transcriptomic study of early responses to the bud dormancy-breaking agent hydrogen cyanamide in "tropicbeauty" peach. *J. Am. Soc. Hortic. Sci.* 144, 244–256. doi: 10.21273/jashs04686-19
The Arabidopsis Information Resource (TAIR) www.arabidopsis.org. Available at: https://www.arabidopsis.org/servers/TairObject?id=137756&type=locus [Accessed September 13, 2019].

The International Peach Genome Initiative, Verde, I., Abbott, A. G., Scalabrin, S., Jung, S., Shu, S., et al. (2013). The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. *Nat. Genet.* 45, 487. doi: 10.1038/ng.2586

Tian, T., Liu, Y., Yan, H., You, Q., Yi, X., Du, Z., et al. (2017). agriGO v2.0: a GO analysis toolkit for the agricultural community 2017 update. *Nucleic Acids Res.* 45, W122–W129. doi: 10.1093/nar/gkx382

Tixier, A., Roxas, A. A., Godfrey, J., Saa, S., Lightle, D., Maillard, P., et al. (2017). Role of bark color on stem temperature and carbohydrate management during dormancy break in persian walnut. *J. Am. Soc. Hortic. Sci.* 142, 454–463. doi: 10.21273/jashs04216-17

Tixier, A., Gambetta, G. A., Godfrey, J., Orozco, J., and Zwieniecki, M. A. (2019). Non-structural carbohydrates in dormant woody perennials; the tale of winter survival and spring arrival. *Front. For. Global Change* 2, 18. doi: 10.3389/f fig.2019.00018

Tylewicz, S., Petterle, A., Marttila, S., Miskolczi, P., Azeez, A., Singh, R. K., et al. (2018). Photoperiodic control of seasonal growth is mediated by ABA acting on cell-cell communication. *Science* 360, 212–215. doi: 10.1126/science.aan8576

Ueno, S., Klopp, C., Leplé, J. C., Derory, J., Noirot, C., Léger, V., et al. (2013). Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing. *BMC Genomics* 14, 236. doi: 10.1186/1471-2164-14-236

Valentini, N., Me, G., Spanna, F., and Lovisetto, M. (2002). Chilling and heat requirement in apricot and peach varieties, in: *XXVI International Horticultural Congress: Key Processes in the Growth and Cropping of Deciduous Fruit and Nut Trees*, Vol. 636. pp. 199–203.

Verde, I., Jenkins, J., Dondini, L., Micali, S., Pagliarani, G., Vendramin, E., et al. (2017). The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. *BMC Genomics* 18, 225. doi: 10.1186/s12864-017-3606-9

Vergara, R., Rubio, S., and Pérez, F. J. (2012). Hypoxia and hydrogen cyanamide induce bud-break and up-regulate hypoxic responsive genes (HRG) and VvFT in grapevine-buds. *Plant Mol. Biol.* 79, 171–178. doi: 10.1007/s11103-012-9904-3

Viti, R., and Monteleone, P. (1991). Observations on flower bud growth in some low yield varieties of apricot. *Acta Hortic.* 293, 319–326. doi: 10.17660/ actahort.1991.293.37

Viti, R., Bartolini, S., and Zanol, G. C. (2010). Biological changes and active oxygen-scavenging enzymes activities in apricot (Prunus armeniaca L.) flower buds during dormancy transitions, in: *in XXVIII International Horticultural Congress on Science and Horticulture for People (IHC2010): International Symposium on the*, Vol. 940. pp. 331–339.

Vogel, J. T., Zarka, D. G., Van Buskirk, H. A., Fowler, S. G., and Thomashow, M. F. (2005). Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. *Plant J.* 41, 195–211. doi: 10.1111/j.1365-313X.2004.02288.x

Wang, Y., Georgi, L. L., Reighard, G. L., Scorza, R., and Abbott, A. G. (2002). Genetic mapping of the evergrowing gene in peach [Prunus persica (L.) Batsch]. *J. Hered.* 93, 352–358. doi: 10.1093/jhered/93.5.352

Wang, S., Chang, Y., and Ellis, B. (2016). Overview of Ovate family proteins, a novel class of plant-specific growth regulators. *Front. Plant Sci.* 7, 417. doi: 10.3389/fpls.2016.00417

Wang, K., Guo, Z.-L., Zhou, W.-T., Zhang, C., Zhang, Z.-Y., Lou, Y., et al. (2018). The Regulation of sporopollenin biosynthesis genes for rapid pollen wall formation. *Plant Physiol.* 178, 283–294. doi: 10.1104/pp.18.00219

Wang, Q., Xu, G., Zhao, X., Zhang, Z., Wang, X., Liu, X., et al. (2019). TCP transcription factor PpTCP20 is involved in peach bud endodormancy by inhibiting PpDAM5/PpDAM6 and interacting with PpABF2. *J. Exp. Bot.* 71 (4), 1585–1597. doi: 10.1093/jxb/erz516

Weiss, D., Schonfeld, M., and Halevy, A. H. (1989). Photosynthetically active chloroplasts in petunia corollas. *Tech. New Dev. In Photosynth. Res.* 168, 489–492. doi: 10.1007/978-1-4684-8571-4_58

Wen, L. H., Zhong, W. J., Huo, X. M., Zhuang, W. B., Ni, Z. J., and Gao, Z. H. (2016). Expression analysis of ABA- and GA-related genes during four stages of bud dormancy in Japanese apricot (Prunus mume Sieb. et Zucc.). *J. Hortic. Sci. Biotechnol.* 91, 362–369. doi: 10.1080/14620316.2016.1160546

Wisniewski, M., Norelli, J., and Artlip, T. (2015). Overexpression of a peach CBF gene in apple: a model for understanding the integration of growth, dormancy, and cold hardiness in woody plants. *Front. Plant Sci.* 6, 85. doi: 10.3389/fpls.2015.00085

Yamaguchi, S. (2008). Gibberellin metabolism and its regulation. *Annu. Rev. Plant Biol.* 59, 225–251. doi: 10.1146/annurev.arplant.59.032607.092804

Yamane, H., Kashiwa, Y., Kakehi, E., Yonemori, K., Mori, H., Hayashi, K., et al. (2006). Differential expression of dehydrin in flower buds of two Japanese apricot cultivars requiring different chilling requirements for bud break. *Tree Physiol.* 26, 1559–1563. doi: 10.1093/treephys/26.12.1559

Yamane, H., Ooka, T., Jotatsu, H., Hosaka, Y., Sasaki, R., and Tao, R. (2011). Expressional regulation of PpDAM5 and PpDAM6, peach (Prunus persica) dormancy-associated MADS-box genes, by low temperature and dormancy-breaking reagent treatment. *J. Exp. Bot.* 62, 3481–3488. doi: 10.1093/jxb/err028

Yamane, H., Wada, M., Honda, C., Matsuura, T., Ikeda, Y., Hirayama, T., et al. (2019). Overexpression of Prunus DAM6 inhibits growth, represses bud break competency of dormant buds and delays bud outgrowth in apple plants. *PloS One* 14, e0214788. doi: 10.1371/journal.pone.0214788

Yordanov, Y. S., Ma, C., Strauss, S. H., and Busov, V. B. (2014). EARLY BUDBREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees. *Proc. Natl. Acad. Sci.* 111, 10001–10006. doi: 10.1073/pnas.1405621111

Zhang, B., and Horvath, S. (2005). A general framework for weighted gene co-expression network analysis. *Stat. Appl. Genet. Mol. Biol.* 4 (1), 1544–6115. doi: 10.2202/1544-6115.1128

Zhang, Z., Zhuo, X., Zhao, K., Zheng, T., Han, Y., Yuan, C., et al. (2018). Transcriptome profiles reveal the crucial roles of hormone and sugar in the bud dormancy of prunus mume. *Sci. Rep.* 8, 5090. doi: 10.1038/s41598-018-23108-9

Zhao, K., Zhou, Y., Ahmad, S., Yong, X., Xie, X., Han, Y., et al. (2018). PmCBFs synthetically affect PmDAM6 by alternative promoter binding and protein complexes towards the dormancy of bud for Prunus mume. *Sci. Rep.* 8, 4527. doi: 10.1038/s41598-018-22537-w

Zhebentyayeva, T. N., Fan, S., Chandra, A., Bielenberg, D. G., Reighard, G. L., Okie, W. R., et al. (2014). Dissection of chilling requirement and bloom date QTLs in peach using a whole genome sequencing of sibling trees from an F2 mapping population. *Tree Genet. Genomes* 10, 35–51. doi: 10.1007/s11295-013-0660-6

Zhu, Y., Li, Y., Xin, D., Chen, W., Shao, X., Wang, Y., et al. (2015). RNA-Seq-based transcriptome analysis of dormant flower buds of Chinese cherry (Prunus pseudocerasus). *Gene* 555, 362–376. doi: 10.1016/j.gene.2014.11.032

Zhuang, W., Gao, Z., Wen, L., Huo, X., Cai, B., and Zhang, Z. (2015). Metabolic changes upon flower bud break in Japanese apricot are enhanced by exogenous GA4. *Hortic. Res.* 2, 15046. doi: 10.1038/ hortres.2015.46

Zwiatzek, J. J. (1991). Cell wall changes in white spruce (Picea glauca) needles subjected to repeated drought stress. *Physiol. Plant.* 82, 513–518. doi: 10.1034/j.1399-3054.1991.820405.x

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

Copyright © 2020 Yu, Conrad, Decroocq, Zhebentyayeva, Williams, Bennett, Roch, Audergon, Dardick, Liu, Abbott and Staton. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.