# The Identification of Small RNAs Differentially Expressed in Apple Buds Reveals a Potential Role of the Mir159-MYB Regulatory Module during Dormancy

Julio Garighan, Etienne Dvorak, Joan Estevan, Karine Loridon, Bruno Huettel, Gautier Sarah, Isabelle Farrera, Julie Leclercq, Priscila Grynberg,

Roberto Coiti Togawa, et al.

► To cite this version:

Julio Garighan, Etienne Dvorak, Joan Estevan, Karine Loridon, Bruno Huettel, et al. The Identification of Small RNAs Differentially Expressed in Apple Buds Reveals a Potential Role of the Mir159-MYB Regulatory Module during Dormancy. Plants, 2021, 10 (12), 10.3390/plants10122665. [https://doi.org/10.3390/plants10122665]

HAL Id: hal-03602721

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

Submitted on 9 Mar 2022

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>

Article

# The Identification of Small RNAs Differentially Expressed in Apple Buds Reveals a Potential Role of the Mir159-MYB Regulatory Module during Dormancy

Julio Garighan<fnref n="1" />,†, Etienne Dvorak <fnref n="1" />,†, Joan Estevan<fnref n="1" />, Karine Loridon<fnref n="1" />, Bruno Huettel<fnref n="2" />,
Gautier Sarah <fnref n="1" />, Isabelle Farrera <fnref n="1" />, Julie Leclercq <fnref n="1" />,<fnref n="3" />, Priscila Grynberg <fnref n="4" />, Roberto Coiti Togawa <fnref n="4" />,
Marcos Mota do Carmo Costa <fnref n="4" />, Evelyne Costes <fnref n="1" /> and Fernando Andrés <fnref n="1" />,*

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

check for
updates

**Citation:** Garighan, J.; Dvorak, E.; Estevan, J.; Loridon, K.; Huettel, B.; Sarah, G.; Farrera, I.; Leclercq, J.; Grynberg, P.; Coiti Togawa, R.; et al. The Identification of Small RNAs Differentially Expressed in Apple Buds Reveals a Potential Role of the Mir159-MYB Regulatory Module during Dormancy. *Plants* **2021**, *10*, 2665. https://doi.org/10.3390/plants10122665

Academic Editors: Mingjun Li and
Ting Wu

Received: 29 October 2021

Accepted: 2 December 2021

Published: 3 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Winter dormancy is an adaptive mechanism that temperate and boreal trees have developed to protect their meristems against low temperatures. In apple trees (*Malus domestica*), cold temperatures induce bud dormancy at the end of summer/beginning of the fall. Apple buds stay dormant during winter until they are exposed to a period of cold, after which they can resume growth (budbreak) and initiate flowering in response to warmer temperatures in spring. It is well-known that small RNAs modulate temperature responses in many plant species, but however, how small RNAs are involved in genetic networks of temperature-mediated dormancy control in fruit tree species remains unclear. Here, we have made use of a recently developed ARGONAUTE (AGO)-purification technique to isolate small RNAs from apple buds. A small RNA-seq experiment resulted in the identification of 17 micro RNAs (miRNAs) that change their pattern of expression in apple buds during dormancy. Furthermore, the functional analysis of their predicted target genes suggests a main role of the 17 miRNAs in phenylpropanoid biosynthesis, gene regulation, plant development and growth, and response to stimulus. Finally, we studied the conservation of the *Arabidopsis thaliana* regulatory miR159-MYB module in apple in the context of the plant hormone abscisic acid homeostasis.

**Keywords:** dormancy; apple tree; small RNAs; miR159

## 1. Introduction

Temperate trees spread over geographical regions presenting wide seasonal environmental fluctuations all over the year. In order to optimize their reproductive success, they adjust their growth and flowering cycles to these recurrent conditions. This is possible thanks to environment sensing mechanisms and signaling pathways that reprogram their meristems in response to changes of photoperiod and temperature. In winter, short photoperiods and low temperatures induce growth cessation, the formation of winter buds to protect the meristematic tissues, and a phase of meristem growth inhibition named endodormancy (or winter dormancy) [1]. Endodormant buds regain their competence to grow after exposure to a certain period of low temperatures that is known as chilling requirement (CR). Once CR is fulfilled, the shoot apical meristem (SAM) undergoes an

<footnote n="1">¹ UMR AGAP Institute, Institut Agro, CIRAD, INRAE, University of Montpellier, F-34398 Montpellier, France; julio.garighan@inrae.fr (J.G.); etienne.dvorak@inrae.fr (E.D.); joan.estevan@inrae.fr (J.E.); karine.loridon@inrae.fr (K.L.); gautier.sarah@inrae.fr (G.S.); isabelle.farrera@supagro.fr (I.F.); julie.leclercq@cirad.fr (J.L.); evelyne costes@inrae.fr (E.C.)</footnote>

<footnote n="2">² Genome Centre, Max Planck Institute for Plant Breeding Research, D-50829 Cologne, Germany; huettel@mpipz.mpg.de</footnote>

<footnote n="3">³ UMR AGAP Institute, CIRAD, F-34398 Montpellier, France</footnote>

<footnote n="4">⁴ Bioinformatica Laboratory, Embrapa Recursos Genéticos e Biotecnologia—Cenargen, Brasília 02372, Brazil; priscila.grynberg@embrapa.br (P.G.); roberto.togawa@embrapa.br (R.C.T.); marcos.costa@embrapa.br (M.M.d.C.C.)</footnote>

<footnote>* Correspondence: fernando.andres-lalaguna@inrae.fr</footnote>

<footnote>† Equal contribution.</footnote>
ecodormant phase. During ecodormancy, warm temperatures typical of springtime lead to growth resumption, budbreak and flowering [2,3].

The Rosaceae is the third most economically important plant family in temperate regions and responsible for the major part of the total worldwide-consumed fruits [4]. Iconic Rosaceous species such as apple (Malus domestica), peach (Prunus persica) and sweet cherry (Prunus avium) are characterized by high CRs. This feature limits the range of latitudes at which they can be productively grown and makes them highly susceptible to global warming. However, despite the importance of a well-adjusted dormancy cycle for flowering timing and fruit production, the understanding on how this process is controlled in Rosaceous species and other significant crops is still in its infancy. A group of genes encoding MADS transcription factors (TF) are believed to be major regulators of the dormancy cycle in many Rosaceous fruit trees [3,5–7]. Some of these MADS TF are named as DORMANCY-ASSOCIATED MADS-BOX (DAM) because of their genetic association with the non-dormant phenotype of the *evg* mutant of peach [6]. DAM genes belong to the MIKC<sup>C</sup> type of MADS TF genes and are similar to SHORT VEGETATIVE PHASE (SVP) from Arabidopsis thaliana [8]. They associate in transcriptional complexes that act in an intricate gene regulatory network (GRN) that control the transition between dormancy phases [3,5,9]. This GRN seems to be regulated, at least in part, by the hormonal and environment-mediated transcriptional modulation of its key components [5,7,10]. For example, the plant hormone abscisic acid (ABA) is a positive regulator of dormancy induction in several tree species [11–13]. Although some recent studies have shed light on mechanisms of transcriptional regulation of bud dormancy in trees, how it is controlled at the post-transcriptional level remains unclear.

Most of the plant small RNAs consist in 21 to 24-nucleotide RNA molecules pro-
duced by DICER-LIKE (DCL) proteins [14]. Small RNAs are loaded onto ARGONAUTE
(AGO)-like proteins to induce mRNA cleavage, translational repression, chromatin com-
paction or DNA elimination [15–17]. They are widely involved in the regulation of plant
development, reproduction, responses to the environment, biotic and abiotic stress and
genome reprogramming [18,19]. The micro RNAs (miRNAs) are a class of small RNAs
that in plants mediate post-transcriptional gene silencing (PTGS) by transcript cleavage
or translational repression [18]. In *A. thaliana*, the miR156 and miR172 are key regulators
of flowering through the post-transcriptional regulation of *SQUAMOSA PROMOTER-
BINDING PROTEIN-LIKE* (SPL) and *APETALA2*-like (AP2-like) genes, respectively [20–22].
Interestingly, miR156 and miR172 were found to be also involved in seed dormancy in
lettuce (*Lactuca sativa*) and *A. thaliana* [23]. However, reports describing the role of small
RNAs in bud dormancy are scarce. Only a few recent studies using genome-wide small
RNA sequencing (small RNA-seq) identified potential miRNAs that regulate dormancy in
trees [24–27]. In tree peony, several miRNAs were found to be differentially accumulated
in buds between endodormancy and ecodormancy after a period of chilling. Among the
identified miRNAs, five of them (miR156k, miR159a, miR167a, miR169a and miR172a)
showed an inverse expression pattern as compared to their target genes [27], suggesting
their role during endodormancy to ecodormancy transition in tree peony. Another recent
study reported the existence of miRNA-mediated regulating the photoperiod-dependent
dormancy induction in Vitis species [28]. Most of the studies concerning Rosaceous species
have been done on pear and peach [24,26,29,30]. More than hundred miRNAs were ex-
pressed in buds during dormancy of Japanese pear (*Pyrus pyrifolia* 'Kosui') [26]. In white
pear (*Pyrus pyrifolia*), it was found that miRNA6390 targeted several *DAM* genes as a part of
a potential mechanism of endodormancy release [24]. A *Prunus*-specific miRNA (miR6285)
was also found to be differentially expressed between endodormancy and ecodormancy in
peach. MiR6285 targets an asparagine-rich protein (NRP) that is involved in the regulation
of ABA signaling in *A. thaliana* [29].

Here, we have made use of a recently developed AGO-purification technique [31]
to isolate small RNAs from apple buds. Combining small RNA-seq and public RNA-seq
data allowed us to identify miRNA and their target genes potentially involved in the
control of dormancy cycle of fruit trees. Furthermore, we studied the conservation between *A. thaliana* and apple of a regulatory miRNA/target gene module in the context of the ABA homeostasis.

## 2. Results

### 2.1. Different Types of Known and Unknown Small RNAs Are Expressed in Apple Buds during the Dormancy Cycle

In order to identify small RNAs expressed during dormancy cycle in apple, bud samples were collected at different dates covering the transition between endodormancy to ecodormancy. Four dates were selected according to the physiological stage of the buds (Table 1 and Figure 1A,B). The physiological stage of the buds was defined using a forcing test (Tabuenca test [32,33]) that allowed us to distinguish between endodormant and ecodormant buds (Figure 1B). Based on the Tabuenca test that we periodically performed between November 2018 and March 2019, we determined that the transition from endodormancy to ecodormancy occurred around 04/02/2019 (Stage 2, Table 1). Before this date, buds remained fully endodormant (Stage 1), whereas after this date, buds became ecodormant (Stages 3 and 4), i.e., competent to resume growth upon exposure to warm T° (Figure 1A and Table 1).

<table><caption>Table 1. Samples used for the small RNA and sequencing results.</caption><thead><tr><th>Sample</th><th>Date</th><th>Avg T (°C)</th><th>CH (h)</th><th>Phase</th><th>Library</th><th>Reads</th></tr></thead><tbody><tr><td rowspan="4">Stage 1</td><td rowspan="4">7 January 2019</td><td rowspan="4">6.00</td><td rowspan="4">665</td><td rowspan="4">Early endo</td><td>A</td><td>22,353,097</td></tr><tr><td>B</td><td>23,372,379</td></tr><tr><td>C</td><td>22,980,048</td></tr><tr><td>D</td><td>20,623,140</td></tr><tr><td rowspan="3">Stage 2</td><td rowspan="3">4 February 2019</td><td rowspan="3">5.84</td><td rowspan="3">1098</td><td rowspan="3">Late endo</td><td>E</td><td>23,353,672</td></tr><tr><td>F</td><td>18,094,050</td></tr><tr><td>G</td><td>23,560,798</td></tr><tr><td rowspan="3">Stage 3</td><td rowspan="3">19 February 2019</td><td rowspan="3">7.98</td><td rowspan="3">1280</td><td rowspan="3">Early eco</td><td>H</td><td>16,446,157</td></tr><tr><td>I</td><td>20,745,045</td></tr><tr><td>J</td><td>20,396,876</td></tr><tr><td rowspan="2">Stage 4</td><td rowspan="2">4 March 2019</td><td rowspan="2">11.28</td><td rowspan="2">1395</td><td rowspan="2">Late eco</td><td>K</td><td>22,812,563</td></tr><tr><td>L</td><td>21,813,675</td></tr></tbody></table>

Sample: name of the sample, Date: date of sampling, Avg T: average of the day-temperature at the sampling date, CH: chilling hours at the sampling date, Phase: physiological stage of the trees based on a forcing test (endo: endodormant, eco: ecodormant), Library: name of the library for sequencing, Reads: number of reads obtained per library.

To isolate small RNAs from these samples, we have made use of a recently developed method called "TraPR", which purification protocol is based on AGO's ability to associate with small RNA. The sequencing results after small RNA isolation using this method are shown in Table 1. A total of 372 different small RNAs were detected across samples, among which, 361 corresponded to biologically relevant small RNAs (i.e., length between 20 to 24 nucleotides) (Table S1). The size distribution shown in Figure 1C indicates a predominant abundance of 21-nucleotides small RNAs. Notably, around one third (28%) of the detected unique small RNAs corresponded to 24-nucleotides small RNAs, which could mediate transcriptional silencing of transposons and pericentromeric repeats via the RdDM pathway [18]. A total of 233 small RNAs were annotated in public databases as known miRNAs, whereas the other 128 were annotated as unknown small RNAs (Table S1). Particularly abundant were the miRNAs belonging to the miR166, miR159, miR482, miR1511, miR319 and miR171 families. Remarkably, only the miR166 accounted for around 80% of the total recovered reads (Table S1).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Identification of small RNAs expressed during dormancy of apple buds. (A) Average temperature and chilling hours accumulation during the experiment done between 2018 and 2019. (B) Result of the Tabuenca test showing the around which endodormancy is released (red arrow) (C) Percentage of reads corresponding to the number of distinct types of small RNAs (Unique sRNAs), to the nucleotide lengths distribution (Reads for each sRNA) with or without mdm-miR166 (Reads for each sRNA w/o mdm-miR166).</figcaption></figure>

## 2.2. Several miRNAs Are Differentially Expressed between Endodormancy and Ecodormancy in Apple

We tested for small RNA differential expression using two R packages (DESeq2 and Edge2). An analysis performed with DESeq2 resulted in the detection of 26 differentially expressed small RNAs (DE-small RNAs), while EdgeR found 33 DE-small RNAs (Supplementary Table S2; *p*-value ≤ 0.05). A group of 20 DE-small RNAs was common between the two analyses (Figure 2A). This group was considered as a high confidence list of DE-small RNAs. The high confident list contained five small RNAs (t11329296_x5768, t05229383_x4081, t08515395_x8579, t00205877_x43290 and t00969484_x10829) that were not included in any public miRNA database. Therefore, we performed a further analysis to evaluate their stem-loop structures. In silico RNA folding analyses predicted that only t08515395_x8579 and t00205877_x43290 have an optimal stem-loop structure compatible with biologically functional miRNAs [34] (Figure 2B). Therefore, the other three small RNAs (t11329296_x5768, t05229383_x4081 and t00969484_x10829) were discarded for downstream
studies. Distinct patterns of expression were observed within the 17 DE-miRNAs of the high confidence list (Table 2 and Figure 2C). Four miRNAs were upregulated at Stage 4 (grey bar in Figure 2C). Two of them, mdm-miR159a and mdm-miR5225c increased their expression at the transition between endodormancy to ecodormancy. Similarly, a group of seven miRNAs (orange bar) were upregulated at Stage 3 but their expression levels dropped at Stage 4. Three miRNAs displayed reduced levels of expression at Stage 4 compared to Stages 1 and 2 (blue bar) and other three ones showed two peaks of expression at Stages 1 and 3 (green bar).

<table><caption>**Table 2.** Top DE-miRNAs.</caption><thead><tr><th rowspan="2">Name (ID-Accession)</th><th colspan="5">EDGER</th><th colspan="4">DESEQ2</th></tr><tr><th>S2 vs. S1</th><th>S3 vs. S1</th><th>S4 vs. S1</th><th>p Value</th><th>S2 vs. S1</th><th>S3 vs. S1</th><th>S4 vs. S1</th><th>p Value</th></tr></thead><tbody><tr><td>ath-miR858b</td><td>-0.70</td><td>1.12</td><td>-0.16</td><td>3.45 × 10<sup>-3</sup></td><td>-0.56</td><td>1.17</td><td>0.01</td><td>1.04 × 10<sup>-3</sup></td></tr><tr><td>mdm-miR159aMIMAT0025898</td><td>0.72</td><td>1.94</td><td>1.64</td><td>4.37 × 10<sup>-3</sup></td><td>0.89</td><td>1.99</td><td>1.93</td><td>1.24 × 10<sup>-3</sup></td></tr><tr><td>mdm-miR164cMIMAT0025909</td><td>-0.90</td><td>-0.22</td><td>-1.98</td><td>3.26 × 10<sup>-3</sup></td><td>-0.76</td><td>-0.17</td><td>-1.77</td><td>3.71 × 10<sup>-3</sup></td></tr><tr><td>mdm-miR171bMIMAT0025939</td><td>-0.83</td><td>0.74</td><td>-0.94</td><td>2.02 × 10<sup>-3</sup></td><td>-0.70</td><td>0.78</td><td>-0.78</td><td>6.28 × 10<sup>-4</sup></td></tr><tr><td>mdm-miR390eMIMAT0025973</td><td>-0.51</td><td>0.08</td><td>1.51</td><td>9.33 × 10<sup>-4</sup></td><td>-0.44</td><td>0.16</td><td>1.70</td><td>3.99 × 10<sup>-6</sup></td></tr><tr><td>mdm-miR482bMIMAT0026022</td><td>-0.75</td><td>0.30</td><td>-1.19</td><td>1.48 × 10<sup>-2</sup></td><td>-0.63</td><td>0.33</td><td>-1.03</td><td>1.19 × 10<sup>-2</sup></td></tr><tr><td>mdm-miR482cMIMAT0026023</td><td>-0.93</td><td>-0.04</td><td>-1.23</td><td>2.89 × 10<sup>-2</sup></td><td>-0.81</td><td>-0.01</td><td>-1.07</td><td>2.47 × 10<sup>-2</sup></td></tr><tr><td>mdm-miR5225cMIMAT0026052</td><td>-0.83</td><td>1.54</td><td>1.40</td><td>5.60 × 10<sup>-6</sup></td><td>-0.70</td><td>1.56</td><td>1.60</td><td>4.76 × 10<sup>-8</sup></td></tr><tr><td>mdm-miR7121cMIMAT0026042</td><td>-0.31</td><td>0.94</td><td>-0.59</td><td>1.17 × 10<sup>-2</sup></td><td>-0.16</td><td>1.01</td><td>-0.35</td><td>4.83 × 10<sup>-3</sup></td></tr><tr><td>mdm-miR7121eMIMAT0026044</td><td>0.30</td><td>1.59</td><td>-0.72</td><td>1.37 × 10<sup>-3</sup></td><td>0.44</td><td>1.62</td><td>-0.52</td><td>1.02 × 10<sup>-3</sup></td></tr><tr><td>mdm-miR858MIMAT0026070</td><td>0.59</td><td>2.51</td><td>-0.26</td><td>3.81 × 10<sup>-4</sup></td><td>0.77</td><td>2.57</td><td>-0.05</td><td>3.12 × 10<sup>-3</sup></td></tr><tr><td>mtr-miR4414a-5p</td><td>0.59</td><td>1.88</td><td>0.29</td><td>2.22 × 10<sup>-3</sup></td><td>0.75</td><td>1.94</td><td>0.47</td><td>9.98 × 10<sup>-4</sup></td></tr><tr><td>peu-miR2910</td><td>-0.01</td><td>1.33</td><td>0.17</td><td>2.05 × 10<sup>-2</sup></td><td>0.11</td><td>1.37</td><td>0.39</td><td>6.48 × 10<sup>-3</sup></td></tr><tr><td>ptc-miR159d</td><td>-0.11</td><td>0.70</td><td>1.92</td><td>2.53 × 10<sup>-4</sup></td><td>0.02</td><td>0.77</td><td>2.16</td><td>1.78 × 10<sup>-7</sup></td></tr><tr><td>sit-miR05-npr</td><td>0.44</td><td>-0.99</td><td>-0.54</td><td>4.77 × 10<sup>-2</sup></td><td>0.56</td><td>-0.94</td><td>-0.26</td><td>2.67 × 10<sup>-2</sup></td></tr><tr><td>t00205877_x43290</td><td>0.45</td><td>0.31</td><td>-1.00</td><td>1.33 × 10<sup>-2</sup></td><td>0.59</td><td>0.38</td><td>-0.79</td><td>3.48 × 10<sup>-3</sup></td></tr><tr><td>t08515395_x8579</td><td>0.20</td><td>-0.01</td><td>-1.67</td><td>6.17 × 10<sup>-3</sup></td><td>0.31</td><td>0.03</td><td>-1.43</td><td>7.52 × 10<sup>-3</sup></td></tr></tbody></table>

Name (ID-accession): ID and accession according to the miRBase nomenclature. S1: Stage 1; S2: Stage 2; S3: Stage 3; S4: Stage 4. Numeric values correspond to Log2 of the Fold Change (Log2FC) between stages.

## 2.3. Predicted Target Genes of DE-miRNAs Are Involved in Key Regulatory Biological Processes

We made use of the psRNATarget web server for the identification of target genes of the DE-miRNAs. This analysis resulted in the identification of 420 potential miRNA target genes (Table S3). One quarter of these genes (101 genes) were common targets among the 17 DE-miRNAs. Next, we made use of public RNA-seq data [35] to represent the expression profile of the 319 unique target genes in a heat map. Several target genes displayed changing expression patterns in apple buds collected at distinct dates during winter dormancy (Figures 3A and S1), suggesting their possible role in this biological process. To investigate the biological function of these target genes, we performed a gene ontology (GO) enrichment analysis. This analysis resulted in the classification of the 319 target genes into five main categories related to plant development, regulation of gene transcription, phenylpropanoid biosynthesis, regulation of growth and response to stimulus (Figure 3B and Table S4). Categories related to plant development included genes involved in meristem and flower development, as well as xylem formation and specification of symmetry. Remarkably, genes with roles related to response to hormones, i.e., ABA, jasmonic acid, gibberellin and auxin, were overrepresented in the category response to stimulus. Moreover, the category regulation of growth included genes related to cell development and proliferation among others. Finally, several genes encoding transcription factors were found within the category of regulation of gene expression.

We found that 18% of the genes showing a GO enrichment (57 genes) were annotated as MYB transcription factors. These MYB-encoding genes were targeted by many of the 17 DE-miRNAs, suggesting their role in a post-transcriptional network regulating bud dormancy. Corroborating this idea, several of these MYB genes were differentially expressed during the dormancy cycle (Figure 3C).
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Several small RNAs are differentially expressed between endodormancy to ecodormancy in apple buds. (A) Venn Diagram comparing the results between the EdgeR and Deseq2 analyses. (B) Stem-loop structure prediction of unknown DE-small RNAs. (C) Heat map Expression profile of the 17 DE-miRNAs.</figcaption></figure>

## 2.4. *The ABA Regulatory Module Mir159a-MYB33/MYB65 Is Conserved in Apple*

Our GO analysis suggested a role of ABA in apple bud dormancy control mediated by miRNAs. ABA is a well-known hormonal regulator of dormancy in tree species [13,36,37]. Interestingly, one of the identified 17 apple DE-miRNAs is annotated as mdm-miRNA159a, which homolog in *A. thaliana* modulates ABA signalling by targeting genes encoding R2R3 MYB transcription factors [38]. According to our small RNA-seq, mdm-miRNA159a expression stayed low during Stages 1 and 2 (endodormancy) and increased during Stages 3 and 4 (ecodormancy) (Figure 2B). This pattern of expression was confirmed by quantifying the pre-miRNA level of mdm-miRNA159a during the dormancy cycle in apple buds (Figure 4A). The mRNA level of its putative target genes *MdMYB33* and *MdMYB65* was also differentially expressed during the dormancy cycle (Figures 3C and 4B,C). As shown in the Figure 4B, *MdMYB33* mRNA was upregulated during endodormancy (between December and January) and strongly downregulated in the ecodormancy (from February to March), coinciding with the maximum accumulation of mdm-miRNA159a pre-miRNA levels (Figure 4A). Although the expression of *MdMYB65* mRNA was similar to that
shown for *MYB33* in the years 2016/2017, this pattern was not reproduced in 2018/2019 (Figure 4C). The downregulation of *MdMYB33* and *MdMYB65* expression could be caused by the activity of mdm-miR159a, since these loci contain a canonical miR159 binding site [38] (Figure 4D). Moreover, the miR159a binding site on *MdMYB33* was empirically demonstrated by degradome sequencing [39].

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

**Figure 3.** GO enrichment of DE-miRNA targets. (A) Heat map of DE-miRNA target genes patterns of expression. (B) Schematic map showing the GO enrichment of the DE-miRNA target genes. A high-resolution image is available in the Figure S1. (C) Heat map of target genes of DE-small RNA encoding MYB transcription factors. In (A,C), blue and yellow colours mean upregulation and downregulation, respectively.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Expression of mdm-miR159a and their targets during dormancy and in response to ABA. (A) Annual expression profile of pre-mdm-miR159a. (B) Annual expression profile of MdMYB33. (C) Annual expression profile of MdMYB65. (D) Comparison between *A. thaliana* and apple of target sites of miR159 on *MYBs*. (E) Alignment of the miR159 target site (red underline) and the sequences recovered from 28 independent clones resulted from the 5' RLM-RACE assay. In (D,E), arrows indicate the 5' position of the cleaved mRNA fragment identified by RLM-5' RACE, and the numbers refer to the number of independent clones analyzed. In (A–C), blue and orange shadows mean endodormancy and ecodormancy phases, respectively. Statistical analysis was done using a *t*-test. Letters shared in common between the dates indicate no significant differences (for *p* ≤ 0.01).</figcaption></figure>

Furthermore, we determined the precise mdm-miRNA159a cleavage site on the *Md-MYB65* gene by 5' RNA ligase-mediated rapid amplification of cDNA ends (5' RLM-RACE). Notably, we detected two potential cleavage sites within the miR159a binding sequence on *MdMYB65* (which is identical to the one in *MdMYB33*). The most frequent cleaved mRNA fragment, supported by 25 out of 28 independent clones, coincides with the same nucleotide sequence that was shown for *Arabidopsis* (Figure 4C,D).

The above-mentioned results suggested a role of the mdm-miR159a and its targets *MdMYB33* and *MdMYB65* during dormancy cycle. In *A. thaliana*, miR159a is induced
by ABA to reduce the sensitivity to this hormone in a feedback loop that involves the repression of these *MYB* genes [38]. In order to correlate the expression of miRNA159 and the ABA metabolism during dormancy, we monitored the level of expression of ABA biosynthetic and catabolic genes in buds in the same time-course experiment shown in Figure 4. Remarkably, the expression of *MdNCDE1*, a gene that encodes the ABA biosynthetic enzyme 9-cis-epoxycarotenoid dioxygenase [40], was induced at certain dates during endodormancy (winter period) and highly reduced during ecodormancy. This result suggested an ABA biosynthetic activity during the endodormancy phase. Furthermore, the mRNA level of *MdABA8'H*, which encodes an ABA 8'-hydroxylase enzyme involved in ABA catabolism [40], was clearly upregulated during endodormancy and strongly repressed in ecodormancy (Figure 5B). As *MdABA8'H* mRNA expression is induced by ABA [40], together with the *MdNCED1* expression profile (Figure 5A), this result corroborates previous reports proposing that highest ABA levels are found during endodormancy and negatively correlates with budbreak in different fruit tree species, including apple, sour cherry, sweet cherry, pear and peach [12,37,41,42]. This result suggested that miR159a is induced by ABA during apple bud dormancy.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Expression of ABA-responsive genes during dormancy. (A) Relative expression of *MdNCED1*. (B) Relative expression of *MdABA8'H*. Blue and orange shadows mean endodormancy and ecodormancy phases, respectively. Statistical analysis was done using a *t*-test. Letters shared in common between the dates indicate no significant differences (for *p* ≤ 0.01).</figcaption></figure>

**3. Discussion**

**3.1. Seventeen DE-miRNAs Could Be Part of GRNs Controlling Bud Dormancy in Apple**

The transition from endodormancy to ecodormancy and the activation of budbreak is believed to be controlled by complex GRNs that integrate endogenous and environmental signals [3,5,43]. TFs are key players in these GRNs, and a timed and coordinated control of their gene expression seems essential to orchestrate dormancy phase transitions [3,5]. This tight regulation is likely to occur not only at the transcriptional level but also post-transcriptionally. MiRNAs are regulatory molecules that exert a large part of their function through the post-transcriptional inhibition of major transcription factors [44]. Furthermore, miRNA can diffuse across tissues and act as long-distance signals [45] making them central players in developmental GRNs, as those regulating dormancy in trees. The sampling
of apical dormant buds from a homogeneous shoot type (spurs) may have favoured the detection of miRNAs, as such buds are synchronized in their developmental stages up to outgrowth. Our small-RNAseq approach and data analysis resulted in the identification of 17 DE-miRNAs possibly involved in GRN controlling bud dormancy of apple. Within these 17 miRNAs, some of them were already found DE during bud dormancy in other tree species. This is the case of members of the miRNA families miR390, miR858, miR164 and miR159 that were also DE during dormancy in tree peony [27] and Japanese pear [26]. However, the observed miRNA expression patterns are not always consistent among species, which suggest a partial evolutionary divergence on the mechanisms controlling dormancy across taxa. Supporting this idea, two recently described miRNAs, miR6390 and miR6285, in white pear [24] and peach [29], respectively, were not detected in apple buds during dormancy. Nevertheless, the 17 DE-miRNAs we identified here are likely to participate in dormancy-related GRNs. Confirming this hypothesis will require a further molecular and functional validation of these miRNAs and their target genes.

## 3.2. Potential Roles of Target Genes of DE-miRNAs during Bud Dormancy of Apple

In silico prediction identified more than 300 target genes for the 17 DE-miRNAs. The functional classification of these potential target genes based on GO enrichment, allow us to hypothesize on the role of the DE-miRNAs in apple dormancy. Corroborating the expected role of the DE-miRNAs in dormancy-associated GRNs, a large part of the predicted target genes encoded TFs, and in particular MYB TFs, that were included in the GO term of regulation of gene expression. The DE-miRNAs' potential targets were also associated with GO terms related to phenylpropanoid biosynthesis, plant development and growth and response to stimulus. Interestingly, phenylpropanoids are secondary metabolites that include flavonoids, monolignols, phenolic acids, stilbenes and coumarin, and have essential roles in plant development. Their function has been already related to bud dormancy in trees. For instance, changes in the expression of genes involved in phenylpropanoid pathway during bud dormancy were reported in apricot (*Prunus armeniaca*) grapevine (*Vitis* species) and pear (*Pyrus pyrifolia*) [46–49]. Moreover, some flavonoids (e.g., kaempferol and quercetin) have functions as auxin transport and redox scavengers [50–52] and could be involved thus in bud dormancy and growth resumption [53,54]. In *A. thaliana*, miR858 silences the expression of several MYB TFs involved in the phenylpropanoid pathway, including MYB11, 12, 111 and 123 [55]. We found mdm-miR858 and ath-miR858b among the identified 17 DE-miRNAs and both of them were specifically upregulated at the transition between endodormancy and ecodormancy and downregulated in the ecodormancy (Figure 2C). Most of the predicted targets of these two miRNAs were upregulated in ecodormancy and encoded MYB TF related to phenylpropanoid pathway (Figure 3C). This is the case of MD02G1087900, MD08G1070700, MD15G1215400 and MD15G1215500, which are predicted orthologues of the proanthocyanidins biosynthesis activator of *A. thaliana* MYB123 TF [56], and MD15G1215100, MD01G1084400 and MD15G1051000, encoding putative MYB3 TF orthologues, which in *A. thaliana* represses phenylpropanoid biosynthesis gene expression [57]. Thus, the miRNAs mdm-miR858 and ath-miR858b could function as modulators of phenylpropanoid during dormancy. This modulation seems rather complex, since potential targets of mdm-miR858 and ath-miR858b act as either activators or repressors in the phenylpropanoid pathway.

Several DE-miRNA target genes were included in groups of genes enriched for GO terms related to growth and plant development (Figure 3B). A group of these genes encoding NAC (NAM, ATAF and CUC) domain transcription factors are known to play functions related to meristem development in *A. thaliana* and are post-transcriptionally regulated by miR164 [58]. In particular, we found that the apple orthologues of *CUP SHAPED COTYLEDON2* (*CUC2*) (MD11G1253500), *NAC1* (MD10G1198400) and *NAC100* (MD06G1196100, MD09G1053700, MD14G1203000 and MD17G1051600) are potential targets of mdm-miR164c. MiR164 inhibits axillary formation via downregulation of targeted *CUC* genes and it is required for normal branching in *A. thaliana* [59]. The miR164-CUC
module seems conserved in other plant species and in cotton (Gossypium hirsutum) *gh*-miR164 and *gh-CUC2* form a regulatory gene network with *BRANCHED 1* (BRC1) that is believed to control branching. BRC1 encodes a TF that controls branching in several plant species [60,61] and it is part of a GRN controlling dormancy and budbreak in apple [5]. Therefore, the downregulation of mdm-miR164c that we observed at ecodormancy (stage 4) (Table 2 and Figure 2C) could be involved in a BRC1-mediated mechanism of budbreak control in apple. These data suggest a role of miRNAs in axillary meristem development, bud dormancy and budbreak. This notion is supported by the high abundance of the miR166 we observed in apple buds during endodormancy and ecodormancy (Table S1). In *A. thaliana*, miR166 targets genes encoding HD-ZIP III transcription factors involved in the formation of axillary meristems, root lateral meristems, and determination of lateral organ polarity required for laminar outgrowth [62]. Thus, miRNAs acting in meristem development seem to play a significant role in bud developmental switches of apple.

### 3.3. *A Possible Mechanism of ABA Hyposensibilization Mediated by the Mdm-miR159a-MYB33 Module during Ecodormancy*

The plant hormone ABA is a central regulator of dormancy by mechanisms involving repression of cell cycle and intercellular communication via plasmodesmata [52]. In many fruit tree and vine species such as apple, peach, pear, sweet cherry and grapevine (Vitis vinifera), ABA levels increase during the endodormancy phase and decline towards budbreak [12,36,41,42,53,63]. We could confirm this pattern of ABA accumulation in two consecutive years through the mRNA level quantification of the ABA-responding gene *MdABA8'H* [40] (Figure 5B). Additionally, the pattern of expression of *MdNCDE1* suggests that the observed ABA increase might be triggered by cold temperatures and/or drought typical of winter [40]. In *A. thaliana*, ABA induces the accumulation of miR159, which in turn reduces the sensitivity to this hormone by repressing the expression of *MYB33* and *MYB105*. This feedback loop is essential for ABA homeostasis in plant cells [40]. MiR159 also represses the expression of *MYB65*, which together with *MYB33*, promotes programmed cell death in endosperm and anthers [64]. Notably, we found that mdm-miR159a is upregulated from endodormancy to ecodormancy in apple buds, and that at least one of its potential post-transcriptional targets (*MdMYB33*) is downregulated at ecodormancy. In this context, we speculate that hyposensibilization to ABA by *MdMYB33* downregulation might be required to prevent ABA-mediated growth inhibition at the initiation of budbreak and flowering. This could be especially important to attenuate the effect of the rise of ABA levels in growing buds caused by punctual drought and/or cold stress during the springtime. The miR159-mediated repression of *MdMYB33* mRNA is also related to juvenile-to-adult phase transition in *A. thaliana* through the transcriptional control of miR156 and its target gene *SPL9* [65]. However, although we could find small RNAs annotated as miR156 within our sequencing data (Table S1), they were not differentially expressed during the endodormancy to ecodormancy transition, which makes their participation in this biological process unlikely.

## **4. Materials and Methods**

### 4.1. Plant Material and Experimental Conditions

Apical floral dormant buds from short shoots (spurs) were collected in 2018–2019 winter from *Malus domestica* cv. “Golden Delicious” trees in an orchard located at the experimental station SudExpé in Marsillargues, south of France. All trees were the same age and their phenology appeared tightly synchronized. Four time points were chosen to cover endodormancy to ecodormancy as well as the transition in between (stages 1–4, see Results). For each stage, three biological replicates were collected, each of them being a mix of nine buds from three trees. The dormant stage of the buds was evaluated under forcing conditions (16 h/8 h, light: dark photoperiod at 22 °C) using Tabuenca’s test [32,33].
## 4.2. RNA Extraction

Total RNA was extracted from ~100 mg bud powder with 'Spectrum Plant Total RNA Kit' (Sigma, Darmstadt, Germany) and quantified on a NanoQuant Plate instrument (TTecan, Männedorf, Switzerland). The quality of the RNA was verified in a TapeStation 4200 (Agilent Genomics, Santa Clara, CA, USA) and in an agarose gel. Small RNAs were extracted using Lexogen's TraPR Small RNA Isolation kit (Lexogen, Wien, Austria) [31].

## 4.3. Small RNA-Seq and Data Analysis

Small RNA libraries were prepared using the Small RNA-Seq Library Prep Kit for Illumina with TraPR (Lexogen) and sequenced on an Illumina HiSeq3000 instrument (1 × 150 bp single read) at the Max Planck Genome Centre in Cologne, Germany. Raw sRNA-seq data was analysed using Mir-Island 2.0 [66]. Reads were mapped to the apple "Golden Delicious" reference genome GDDH13 v1.1 [67]. Plant mature miRNAs from the literature were obtained by combining two miRNA databases: mirBase [68] and PNRD [69], which contain respectively 308 and 208 entries for mature apple miRNAs. Identification of miRNAs within the small RNA sequencing data was performed as follows: the Mir-Island program first identifies already known small RNAs based on data provided by the user (miRbase and PNRD). Then, it searches for new miRNAs by mapping small RNA reads to loci in the genome that can form a proper RNA hairpin structure that is coherent with a miRNA-producing locus. It also provides the reads of mature miRNA and miRNA* found in the sequencing for each small RNA. An independent miRNA annotation by sequence homology with known MiRNA described in miRBase and PNRD database was also done (Table S1). Mir-Island read counts were then used for differential expression analysis. This was conducted in R 4.0.3, using two common packages for RNA-seq data analysis: EdgeR 3.32.1 [70] and DESeq2 1.30.1 [71]. Counts were normalized using the trimmed mean of M values (TMM) method in EdgeR and the median of ratios method in DESeq2. Low-count tags were discarded according to each package pipeline. In both packages, counts are fitted to Negative Binomial generalized linear models. Fold-changes are estimated and miRNAs are tested for differential expression compared to the first stage (stage 1) with a likelihood ratio test (LRT). False discovery rates (FDR) are calculated using the Benjamini-Hochberg method. Heat maps of DE-miRNAs were constructed using Morpheus (https://software.broadinstitute.org/morpheus, 31 March 2021) and default parameters. RNA secondary structure was calculated using RNAfold [72].

## 4.4. In Silico Prediction of miRNA Target Genes

MiRNA targets were predicted using the psRNAtarget web interface (2017 release) [73]. We chose a cut-off of E = 3.0 for miRNA-target sequence matching. This only allows to find potential degradation by cleavage, as translation inhibition by miRNAs usually acts when there are more mismatches between the sequence. The "Golden Delicious" reference transcriptome GDDH13 v1.1 was used [67]. Public RNA-seq data [35] were used to build expression heat maps using Morpheus (https://software.broadinstitute.org/morpheus, 31 March 2021) and default parameters.

## 4.5. Gene Ontology Analyses

The GO enrichment test was done using the Bingo tool (v3.0.3) [74] implemented in Cytoscape (v3.8.2, cytoscape.org, April 2021). The GO term enrichment in the list of 319 potential DE-miRNA target genes was assessed using the hypergeometric test. Raw *p*-values were adjusted for multiple testing using the Benjamini-Hochberg method [75]. A high resolution GO map can be found in Figure S1.

## 4.6. Gene Expression Quantification

RT-qPCRs were performed on a LightCycler 480 using SYBR Green fluorescence detection (Roche, Basel, Switzerland). The parameters used for large were 50 cycles and 60 °C for annealing. PCR efficiencies were estimated using the LinRegPCR program [76].
Primers were designed using Primer3 online interface [77]. Relative quantification was calculated using the ddCT method [78]. WD40 and MDH were used as reference genes [79]. A list of all primers used in this study can be found in Table S5. All statistical analyses were performed using R stats package for R 4.0.3.

## 4.7. RLM-5'RACE

To validate the predicted miRNA-induced cleavage site, a modified RNA ligase-mediated 5' Rapid Amplification of cDNA ends was used as described in [80]. From 5 µg of total RNA, the 5'RACE protocol was performed using the GeneRacer kit (Invitrogen Life Technologies, Waltham, MA, USA), according to the supplier's instructions except for the first two steps (decapping and phosphorylation), which are specific of intact transcripts.

In order to maximise amplification specificity, two successive rounds of PCR amplification (nested PCR) were performed using the Phusion High-Fidelity DNA Polymerase with HF buffer (ThermoScientific), primers provided in the GeneRacer kit (GeneRacer 5' primer and GeneRacer nested 5' primer in the first and second round, respectively) as forward primers and gene-specific primers (MYB65-GSP and MYB65-nested-GSP) targeting MYB mRNA as reverse primers.

Primer design of the MYB65-GSP and MYB65-nested-GSP primers was performed downsteam of the predicted cleavage site using the Primer3Plus tool (https://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi; 25 September 2021) [77] with the following parameters complying with the recommendations of the GeneRacer kit: amplicon size 70–150 pb, primer length 23–28 nucleotides, primer melting temperature (Tm) 60–74 °C. The sequences of the resulting primers are shown in the Table S5.

## 4.8. Cloning, Colony PCR and Sequencing

After PCR product separation by gel electrophoresis (1.5% agarose gel in 1× Tris acetate EDTA buffer, migration 100 V/1h30) and visualization under ethidium bromide/UV light, amplification products with a size compatible with the occurrence of miRNA-directed cleavage were selected. They were purified with the DNA Clean & Concentrator kit (ZymoResearch, Irvine, CA, USA), then their blunt ends were changed into cohesive ends with dATP and Diamond Taq DNA Polymerase (Eurogentec, Seraing, Belgium), thanks to its terminal transferase activity allowing TA cloning. After another purification with the DNA Clean & Concentrator kit (ZymoResearch), cloning was performed using the TOPO-TA Cloning Kit for Sequencing (Invitrogen Life Technologies) following the manufacturer's instructions modified as follows: bacterial transformations were performed using 2.5 µL of ligation reaction with 50 µL of One Shot chemically competent cells. Transformed bacteria (50, 100 or 150 µL) were plated on a Luria-Bertani (LB) medium supplemented with 1.5% agar and with a final concentration of 100 µg/mL ampicillin.

Several produced colonies from each target gene were analyzed by colony PCR with M13 F/R primers followed by electrophoresis on a 1.5% agarose gel. Those of the PCR products displaying a size compatible with the predicted miRNA-directed cleavage were purified using the QIAquick PCR purification kit (Qiagen, Hilden, Germany). Sanger sequencing of the purified amplification products was performed by Eurofins Genomics (Köln, Germany).

# 5. Conclusions

Our results revealed that at least 17 miRNAs change their pattern of expression during dormancy in apple buds. Since most plant miRNAs control the activity of multigene TF families, the activity of these 17 miRNAs is expected to have a significant effect on bud development. Thus, the knowledge produced here can pave the way for using miRNAs as a biotechnological tool for improving apple tree traits related to dormancy and bud-break. However, this will undoubtedly require a step of miRNA validation and functional characterization that could be done, for example, using gene editing technologies.
**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/plants10122665/s1. Figure S1: A high resolution image of the GO enrichment map shown in Figure 3B. Table S1: List of small RNAs identified in the small RNA-seq experiment. Table S2: List of detected DE-miRNAs using EdgeR and Deseq2 programs. Table S3: List of potential miRNA target genes identified using psRNAtarget. Table S4: Results of the GO enrichment analysis done on the predicted DE-miRNA target genes. Table S5: List of primers used in this study.

**Author Contributions:** Conceptualization, F.A. and E.C.; methodology, J.G., E.D., J.E., K.L., B.H., G.S. and I.F.; investigation P.G., R.C.T., M.M.d.C.C., F.A., J.G., P.G. and G.S.; writing F.A., E.C. and J.G.; review and editing, J.L.; supervision, F.A. and E.C.; project administration, F.A.; funding acquisition, F.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by BAP INRAE and La région Occitanie, grant number 24001136. The APC was funded by INRAE SE.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The small RNA sequencing datasets supporting the results of this article are available in NCBI's Gene Expression Omnibus (GEO) repository through the GEO Series accession number GSE189658.

**Acknowledgments:** We thank Guillaume Perez (AGAPi, Montpellier, France) for technical support in the field. We are grateful to Estelle Jaligot for her precious advices about RLM 5-RACE experiments (AGAPi, Montpellier, France).

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Maurya, J.P.; Bhalerao, R.P. Photoperiod- and temperature-mediated control of growth cessation and dormancy in trees: A molecular perspective. *Ann. Bot.* **2017**, *120*, 351–360. [CrossRef]

2. Lang, G.A. Endo-, para-, and ecodormancy: Physiological terminology and classification for dormancy research. *HortScience* **1987**, *22*, 371–377.

3. da Silveira Falavigna, V.; Guitton, B.; Costes, E.; Andrés, F. I want to (bud) break free: The potential role of DAM and SVP-like genes in regulating dormancy cycle in temperate fruit trees. *Front. Plant Sci.* **2018**, *9*, 1990. [CrossRef]

4. Yamamoto, T.; Terakami, S. Genomics of pear and other rosaceae fruit trees. *Breed. Sci.* **2016**, *66*, 148–159. [CrossRef] [PubMed]

5. da Silveira Falavigna, V.; Severing, E.; Lai, X.; Estevan, J.; Farrera, I.; Hugouvieux, V.; Revers, L.F.; Zubieta, C.; Coupland, G.; Costes, E.; et al. Unraveling the role of MADS transcription factor complexes in apple tree dormancy. *New Phytol.* **2021**, *232*, 2071–2088. [CrossRef]

6. Bielenberg, D.G.; Wang, Y.; Fan, S.; Reighard, G.L.; Scorza, R.; Abbott, A.G. A deletion affecting several gene candidates is present in the evergrowing peach mutant. *J. Hered.* **2004**, *95*, 436–444. [CrossRef]

7. Vimont, N.; Fouché, M.; Campoy, J.A.; Tong, M.; Arkoun, M.; Yvin, J.-C.; Wigge, P.A.; Dirlewanger, E.; Cortijo, S.; Wenden, B. From bud formation to flowering: Transcriptomic state defines the cherry developmental phases of sweet cherry bud dormancy. *BMC Genom.* **2019**, *20*, 974. [CrossRef]

8. Bielenberg, D.G.; Rauh, B.; Fan, S.; Gasic, K.; Abbott, A.G.; Reighard, G.L.; Okie, W.R.; Wells, C.E. Genotyping by sequencing for SNP-based linkage map construction and QTL analysis of chilling requirement and bloom date in peach [*Prunus persica* (L.) Batsch]. *PLoS ONE* **2015**, *10*, e0139406. [CrossRef]

9. Abbott, A.G.; Zhebentyayeva, T.; Barakat, A.; Liu, Z. The genetic control of bud-break in trees. In *Advances in Botanical Research*; Plomion, C., Adam-Blondon, A.-F., Eds.; Academic Press: Cambridge, MA, USA, 2015; Volume 74, pp. 201–228, ISBN 978-0-12-398548-4.

10. Fadón, E.; Fernandez, E.; Behn, H.; Luedeling, E. A conceptual framework for winter dormancy in deciduous trees. *Agronomy* **2020**, *10*, 241. [CrossRef]

11. Chmielewski, F.; Gotz, K.; Homann, T.; Huschek, G.; Rawel, H. Identification of endodormancy release for cherries (*Prunus avium* L.) by abscisic acid and sugars. *J. Hortic.* **2017**, *4*, 3. [CrossRef]

12. Tuan, P.A.; Bai, S.; Saito, T.; Ito, A.; Moriguchi, T. Dormancy-associated MADS-Box (DAM) and the abscisic acid pathway regulate pear endodormancy through a feedback mechanism. *Plant Cell Physiol.* **2017**, *58*, 1378–1390. [CrossRef]

13. Tylewicz, S.; Petterle, A.; Marttila, S.; Miskolczi, P.; Azeez, A.; Singh, R.K.; Immanen, J.; Mähler, N.; Hvidsten, T.R.; Eklund, D.M.; et al. Photoperiodic control of seasonal growth is mediated by ABA acting on cell-cell communication. *Science* **2018**, *360*, 212–215. [CrossRef]
14. Henderson, I.R.; Zhang, X.; Lu, C.; Johnson, L.; Meyers, B.C.; Green, P.J.; Jacobsen, S.E. Dissecting *Arabidopsis thaliana* DICER function in small RNA processing, gene silencing and DNA methylation patterning. *Nat. Genet.* **2006**, *38*, 721–725. [CrossRef]

15. Wilson, R.C.; Doudna, J.A. Molecular mechanisms of RNA interference. *Annu. Rev. Biophys.* **2013**, *42*, 217–239. [CrossRef]

16. Chalker, D.L.; Yao, M.-C. DNA elimination in ciliates: Transposon domestication and genome surveillance. *Annu. Rev. Genet.* **2011**, *45*, 227–246. [CrossRef]

17. Matzke, M.A.; Mosher, R.A. RNA-directed DNA methylation: An epigenetic pathway of increasing complexity. *Nat. Rev. Genet.* **2014**, *15*, 394–408. [CrossRef]

18. Borges, F.; Martienssen, R.A. The Expanding world of small RNAs in plants. *Nat. Rev. Mol. Cell Biol.* **2015**, *16*, 727–741. [CrossRef]

19. D'Ario, M.; Griffiths-Jones, S.; Kim, M. Small RNAs: Big impact on plant development. *Trends Plant Sci.* **2017**, *22*, 1056–1068. [CrossRef]

20. Hyun, Y.; Richter, R.; Coupland, G. Competence to flower: Age-controlled sensitivity to environmental cues. *Plant Physiol.* **2017**, *173*, 36–46. [CrossRef]

21. Wang, J.-W.; Czech, B.; Weigel, D. MiR156-regulated SPL transcription factors define an endogenous flowering pathway in *Arabidopsis Thaliana*. *Cell* **2009**, *138*, 738–749. [CrossRef]

22. Ó'Maoiléidigh, D.S.; van Driel, A.D.; Singh, A.; Sang, Q.; Le Bec, N.; Vincent, C.; de Olalla, E.B.G.; Vayssières, A.; Romera Branchat, M.; Severing, E.; et al. Systematic analyses of the MIR172 family members of *Arabidopsis* define their distinct roles in regulation of APETALA2 during floral transition. *PLoS Biol.* **2021**, *19*, e3001043. [CrossRef]

23. Huo, H.; Wei, S.; Bradford, K.J. Delay of germination 1 (DOG1) regulates both seed dormancy and flowering time through microRNA pathways. *Proc. Natl. Acad. Sci. USA* **2016**, *113*, E2199–E2206. [CrossRef]

24. Niu, Q.; Li, J.; Cai, D.; Qian, M.; Jia, H.; Bai, S.; Hussain, S.; Liu, G.; Teng, Y.; Zheng, X. Dormancy-associated MADS-box genes and microRNAs jointly control dormancy transition in pear (*Pyrus pyrifolia* White Pear Group) flower bud. *J. Exp. Bot.* **2016**, *67*, 239–257. [CrossRef] [PubMed]

25. Rothkegel, K.; Sánchez, E.; Montes, C.; Greve, M.; Tapia, S.; Bravo, S.; Prieto, H.; Almeida, A.M. DNA methylation and small interference RNAs participate in the regulation of MADS-box genes involved in dormancy in sweet cherry (*Prunus Avium* L.). *Tree Physiol.* **2017**, *37*, 1–13. [CrossRef]

26. Bai, S.; Saito, T.; Ito, A.; Tuan, P.A.; Xu, Y.; Teng, Y.; Moriguchi, T. Small RNA and PARE Sequencing in flower bud reveal the involvement of SRNAs in endodormancy release of Japanese Pear (*Pyrus pyrifolia* 'Kosui'). *BMC Genom.* **2016**, *17*, 230. [CrossRef] [PubMed]

27. Zhang, Y.; Wang, Y.; Gao, X.; Liu, C.; Gai, S. Identification and characterization of microRNAs in tree peony during chilling induced dormancy release by high-throughput sequencing. *Sci. Rep.* **2018**, *8*, 4537. [CrossRef]

28. Smita, S.; Robben, M.; Deuja, A.; Accerbi, M.; Green, P.J.; Subramanian, S.; Fennell, A. Integrative analysis of gene expression and miRNAs reveal biological pathways associated with bud paradormancy and endodormancy in grapevine. *Plants* **2021**, *10*, 669. [CrossRef] [PubMed]

29. Yu, J.; Bennett, D.; Dardick, C.; Zhebentyayeva, T.; Abbott, A.G.; Liu, Z.; Staton, M.E. Genome-wide changes of regulatory non-coding RNAs reveal pollen development initiated at ecodormancy in peach. *Front. Mol. Biosci.* **2021**, *8*, 612881. [CrossRef]

30. Barakat, A.; Sriram, A.; Park, J.; Zhebentyayeva, T.; Main, D.; Abbott, A. Genome wide identification of chilling responsive microRNAs in *Prunus persica*. *BMC Genom.* **2012**, *13*, 481. [CrossRef]

31. Brentzinger, T.; Oberlin, S.; Schott, G.; Handler, D.; Svozil, J.; Barragan-Borrero, V.; Humbert, A.; Duharcourt, S.; Brennecke, J.; Voinnet, O. A universal method for the rapid isolation of all known classes of functional silencing small RNAs. *Nucleic Acids Res.* **2020**, *48*, e79. [CrossRef]

32. Tabuenca, M. Chilling requirements of apricot, peach and pear varieties. *Anal. Estac. Exp. Aula Dei* **1964**, *7*, 113–132.

33. Legave, J.-M.; Guédon, Y.; Malagi, G.; El Yaacoubi, A.; Bonhomme, M. Differentiated responses of apple tree floral phenology to global warming in contrasting climatic regions. *Front. Plant Sci.* **2015**, *6*, 1054. [CrossRef] [PubMed]

34. Axtell, M.J.; Meyers, B.C. Revisiting criteria for plant microRNA annotation in the era of big data. *Plant Cell* **2018**, *30*, 272–284. [CrossRef] [PubMed]

35. Moser, M.; Asquini, E.; Miolli, G.V.; Weigl, K.; Hanke, M.-V.; Flachowsky, H.; Si-Ammour, A. The MADS-box gene MdDAM1 controls growth cessation and bud dormancy in apple. *Front. Plant Sci.* **2020**, *11*, 1003. [CrossRef]

36. Li, J.; Xu, Y.; Niu, Q.; He, L.; Teng, Y.; Bai, S. Abscisic acid (ABA) promotes the induction and maintenance of pear (*Pyrus pyrifolia* White Pear Group) flower bud endodormancy. *Int. J. Mol. Sci.* **2018**, *19*, 310. [CrossRef] [PubMed]

37. Faust, M.; Erez, A.; Rowland, L.J.; Wang, S.Y.; Norman, H.A. Bud dormancy in perennial fruit trees: Physiological basis for dormancy induction, maintenance, and release. *HortScience* **1997**, *32*, 623–629. [CrossRef]

38. Reyes, J.L.; Chua, N.-H. ABA Induction of MiR159 controls transcript levels of two MYB factors during *Arabidopsis* seed germination. *Plant J. Cell Mol. Biol.* **2007**, *49*, 592–606. [CrossRef]

39. Song, C.; Zhang, D.; Zheng, L.; Zhang, J.; Zhang, B.; Luo, W.; Li, Y.; Li, G.; Ma, J.; Han, M. MiRNA and degradome sequencing reveal miRNA and their target genes that may mediate shoot growth in spur type mutant "Yanfu 6". *Front. Plant Sci.* **2017**, *8*, 441. [CrossRef]

40. Lee, Y.; Do, V.G.; Kim, S.; Kweon, H.; McGhie, T.K. Cold stress triggers premature fruit abscission through ABA-dependent signal transduction in early developing apple. *PLoS ONE* **2021**, *16*, e0249975. [CrossRef]
41. Chmielewski, F.M.; Baldermann, S.; Götz, K.P.; Homann, T.; Gödeke, K.; Schumacher, F.; Huschek, G.; Rawel, H.M. Abscisic acid related metabolites in sweet cherry buds (*Prunus avium* L.). *J. Hortic.* **2018**, *5*, 221. [CrossRef]

42. Wang, D.; Gao, Z.; Du, P.; Xiao, W.; Tan, Q.; Chen, X.; Li, L.; Gao, D. Expression of ABA metabolism-related genes suggests similarities and differences between seed dormancy and bud dormancy of peach (*Prunus persica*). *Front. Plant Sci.* **2016**, *6*, 1248. [CrossRef] [PubMed]

43. Miotto, Y.E.; Tessele, C.; Czermainski, A.B.C.; Porto, D.D.; da Silveira Falavigna, V.; Sartor, T.; Cattani, A.M.; Delatorre, C.A.; de Alencar, S.A.; da Silva-Junior, O.B.; et al. Spring is coming: Genetic analyses of the bud break date locus reveal candidate genes from the cold perception pathway to dormancy release in apple (*Malus × domestica* Borkh.). *Front. Plant Sci.* **2019**, *10*, 33. [CrossRef] [PubMed]

44. Samad, A.F.A.; Sajad, M.; Nazaruddin, N.; Fauzi, I.A.; Murad, A.M.A.; Zainal, Z.; Ismail, I. microRNA and transcription factor: Key players in plant regulatory network. *Front. Plant Sci.* **2017**, *8*, 565. [CrossRef]

45. Subramanian, S. Little RNAs go a long way: Long-distance signaling by microRNAs. *Mol. Plant* **2019**, *12*, 18–20. [CrossRef]

46. Conrad, A.O.; Yu, J.; Staton, M.E.; Audergon, J.-M.; Roch, G.; Decroocq, V.; Knagge, K.; Chen, H.; Zhebentyayeva, T.; Liu, Z.; et al. Association of the phenylpropanoid pathway with dormancy and adaptive trait variation in apricot (*Prunus armeniaca*). *Tree Physiol.* **2019**, *39*, 1136–1148. [CrossRef] [PubMed]

47. Fennell, A.Y.; Schlauch, K.A.; Gouthu, S.; Deluc, L.G.; Khadka, V.; Sreekantan, L.; Grimplet, J.; Cramer, G.R.; Mathiason, K.L. Short day transcriptomic programming during induction of dormancy in grapevine. *Front. Plant Sci.* **2015**, *6*, 834. [CrossRef] [PubMed]

48. Deluc, L.G.; Decendit, A.; Papastamoulis, Y.; Mérillon, J.-M.; Cushman, J.C.; Cramer, G.R. Water deficit increases stilbene metabolism in cabernet sauvignon berries. *J. Agric. Food Chem.* **2011**, *59*, 289–297. [CrossRef]

49. Liu, G.; Li, W.; Zheng, P.; Xu, T.; Chen, L.; Liu, D.; Hussain, S.; Teng, Y. Transcriptomic analysis of “Suli” Pear (*Pyrus pyrifolia* White Pear Group) buds during the dormancy by RNA-seq. *BMC Genom.* **2012**, *13*, 700. [CrossRef]

50. Brown, D.E.; Rashotte, A.M.; Murphy, A.S.; Normanly, J.; Tague, B.W.; Peer, W.A.; Taiz, L.; Muday, G.K. Flavonoids act as negative regulators of auxin transport in vivo in *Arabidopsis*. *Plant Physiol.* **2001**, *126*, 524–535. [CrossRef]

51. Peer, W.A.; Brown, D.E.; Tague, B.W.; Muday, G.K.; Taiz, L.; Murphy, A.S. Flavonoid accumulation patterns of transparent testa mutants of *Arabidopsis*. *Plant Physiol.* **2001**, *126*, 536–548. [CrossRef]

52. Buer, C.S.; Muday, G.K. The transparent Testa4 mutation prevents flavonoid synthesis and alters auxin transport and the response of *Arabidopsis* roots to gravity and light. *Plant Cell* **2004**, *16*, 1191–1205. [CrossRef]

53. Powell, L.E. The hormonal control of bud and seed dormancy in woody plants. In *Plant Hormones and Their Role in Plant Growth and Development*; Springer: Dordrecht, The Netherlands, 1987; pp. 539–552, ISBN 978-90-247-3498-6.

54. Liu, J.; Sherif, S.M. Hormonal orchestration of bud dormancy cycle in deciduous woody perennials. *Front. Plant Sci.* **2019**, *10*, 1136. [CrossRef]

55. Sharma, D.; Tiwari, M.; Pandey, A.; Bhatia, C.; Sharma, A.; Trivedi, P.K. microRNA858 is a potential regulator of phenylpropanoid pathway and plant development. *Plant Physiol.* **2016**, *171*, 944–959. [CrossRef] [PubMed]

56. Nesi, N.; Jond, C.; Debeaujon, I.; Caboche, M.; Lepiniec, L. The *Arabidopsis* TT2 gene encodes an R2R3 MYB domain protein that acts as a key determinant for proanthocyanidin accumulation in developing seed. *Plant Cell* **2001**, *13*, 2099–2114. [CrossRef] [PubMed]

57. Zhou, M.; Zhang, K.; Sun, Z.; Yan, M.; Chen, C.; Zhang, X.; Tang, Y.; Wu, Y. LNK1 and LNK2 corepressors interact with the MYB3 transcription factor in phenylpropanoid biosynthesis. *Plant Physiol.* **2017**, *174*, 1348–1358. [CrossRef] [PubMed]

58. Mallory, A.C.; Dugas, D.V.; Bartel, D.P.; Bartel, B. microRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. *Curr. Biol. CB* **2004**, *14*, 1035–1046. [CrossRef] [PubMed]

59. Raman, S.; Greb, T.; Peaucelle, A.; Blein, T.; Laufs, P.; Theres, K. Interplay of MiR164, CUP-SHAPED COTYLEDON genes and lateral suppressor controls axillary meristem formation in *Arabidopsis thaliana*. *Plant J. Cell Mol. Biol.* **2008**, *55*, 65–76. [CrossRef]

60. Wang, M.; Le Moigne, M.-A.; Bertheloot, J.; Crespel, L.; Perez-Garcia, M.-D.; Ogé, L.; Demotes-Mainard, S.; Hamama, L.; Davière, J.-M.; Sakr, S. BRANCHED1: A key hub of shoot branching. *Front. Plant Sci.* **2019**, *10*, 76. [CrossRef]

61. Zhan, J.; Chu, Y.; Wang, Y.; Diao, Y.; Zhao, Y.; Liu, L.; Wei, X.; Meng, Y.; Li, F.; Ge, X. The MiR164-GhCUC2-GhBRC1 module regulates plant architecture through abscisic acid in cotton. *Plant Biotechnol. J.* **2021**, *19*, 1839–1851. [CrossRef]

62. Sakaguchi, J.; Watanabe, Y. MiR165/166 and the development of land plants. *Dev. Growth Differ.* **2012**, *54*, 93–99. [CrossRef]

63. Zheng, C.; Halaly, T.; Acheampong, A.K.; Takebayashi, Y.; Jikumaru, Y.; Kamiya, Y.; Or, E. Abscisic acid (ABA) regulates grape bud dormancy, and dormancy release stimuli may act through modification of ABA metabolism. *J. Exp. Bot.* **2015**, *66*, 1527–1542. [CrossRef]

64. Alonso-Peral, M.M.; Li, J.; Li, Y.; Allen, R.S.; Schnippenkoetter, W.; Ohms, S.; White, R.G.; Millar, A.A. The microRNA159-regulated GAMYB-like genes inhibit growth and promote programmed cell death in *Arabidopsis*. *Plant Physiol.* **2010**, *154*, 757–771. [CrossRef]

65. Guo, C.; Xu, Y.; Shi, M.; Lai, Y.; Wu, X.; Wang, H.; Zhu, Z.; Poethig, R.S.; Wu, G. Repression of MiR156 by MiR159 regulates the timing of the juvenile-to-adult transition in *Arabidopsis*. *Plant Cell* **2017**, *29*, 1293–1304. [CrossRef] [PubMed]

66. Gao, T.; Meng, X.; Zhang, W.; Jin, W. *MiR-Island: An Ultrafast and Memory-Efficient Tool for Plant MiRNA Annotation and Expression Analysis*; Zhejiang Sci-Tech University: Zhejiang, China, 2019.
67. Daccord, N.; Celton, J.-M.; Linsmith, G.; Becker, C.; Choisne, N.; Schijlen, E.; van de Geest, H.; Bianco, L.; Micheletti, D.; Velasco, R.; et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. *Nat. Genet.* **2017**, *49*, 1099–1106. [CrossRef] [PubMed]

68. Griffiths-Jones, S. MiRBase: microRNA sequences and annotation. *Curr. Protoc. Bioinf.* **2010**, *29*, Unit 12.9.1–Unit 12.9.10. [CrossRef]

69. Yi, X.; Zhang, Z.; Ling, Y.; Xu, W.; Su, Z. PNRD: A plant non-coding RNA database. *Nucleic Acids Res.* **2015**, *43*, D982–D989. [CrossRef] [PubMed]

70. Chen, Y.; Lun, A.T.L.; Smyth, G.K. Differential expression analysis of complex RNA-seq experiments using EdgeR. In *Statistical Analysis of Next Generation Sequencing Data*; Datta, S., Nettleton, D., Eds.; Frontiers in Probability and the Statistical Sciences; Springer International Publishing: Cham, Switzerland, 2014; pp. 51–74, ISBN 978-3-319-07212-8.

71. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef]

72. Lorenz, R.; Bernhart, S.H.; Höner zu Siederdissen, C.; Tafer, H.; Flamm, C.; Stadler, P.F.; Hofacker, I.L. ViennaRNA package 2.0. *Algorithms Mol. Biol.* **2011**, *6*, 26. [CrossRef]

73. Dai, X.; Zhuang, Z.; Zhao, P.X. PsRNA Target: A plant small RNA target analysis server (2017 release). *Nucleic Acids Res.* **2018**, *46*, W49–W54. [CrossRef] [PubMed]

74. Maere, S.; Heymans, K.; Kuiper, M. BiNGO: A cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. *Bioinformatics* **2005**, *21*, 3448–3449. [CrossRef] [PubMed]

75. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. *J. R. Stat. Soc. Ser. B Methodol.* **1995**, *57*, 289–300. [CrossRef]

76. Ruijter, J.M.; Ramakers, C.; Hoogaars, W.M.H.; Karlen, Y.; Bakker, O.; van den Hoff, M.J.B.; Moorman, A.F.M. Amplification efficiency: Linking baseline and bias in the analysis of quantitative PCR data. *Nucleic Acids Res.* **2009**, *37*, e45. [CrossRef]

77. Untergasser, A.; Nijveen, H.; Rao, X.; Bisseling, T.; Geurts, R.; Leunissen, J.A.M. Primer3Plus, an enhanced web interface to Primer3. *Nucleic Acids Res.* **2007**, *35*, W71–W74. [CrossRef]

78. Nordgård, O.; Kvaløy, J.T.; Farmen, R.K.; Heikkilä, R. Error propagation in relative real-time reverse transcription polymerase chain reaction quantification models: The balance between accuracy and precision. *Anal. Biochem.* **2006**, *356*, 182–193. [CrossRef] [PubMed]

79. Perini, P.; Pasquali, G.; Margis-Pinheiro, M.; de Oliviera, P.R.D.; Revers, L.F. Reference genes for transcriptional analysis of flowering and fruit ripening stages in apple (*Malus × domestica* Borkh.). *Mol. Breed.* **2014**, *34*, 829–842. [CrossRef]

80. Gébelin, V.; Argout, X.; Engchuan, W.; Pitollat, B.; Duan, C.; Montoro, P.; Leclercq, J. Identification of novel microRNAs in hevea brasiliensis and computational prediction of their targets. *BMC Plant Biol.* **2012**, *12*, 18. [CrossRef] [PubMed]