# Genomic approach to study floral development genes in Rosa sp.

Annick Dubois, Arnaud A. Remay, Olivier Raymond, Sandrine Balzergue, Aurélie Chauvet, Marion Maene, Yann Pécrix, Shu-Hua Yang, Julien Jeauffre, Tatiana Thouroude, et al.

► To cite this version:

Annick Dubois, Arnaud A. Remay, Olivier Raymond, Sandrine Balzergue, Aurélie Chauvet, et al. Ge-
nomic approach to study floral development genes in Rosa sp.. PLoS ONE, 2011, 6 (12), 10.1371/jour-
nal.pone.0028455 . hal-00753396

HAL Id: hal-00753396

https://inria.hal.science/hal-00753396v1

Submitted on 8 Oct 2018

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 CC BY 4.0 - Attribution - International License
# Genomic Approach to Study Floral Development Genes in *Rosa* sp.

Annick Dubois¹*, Arnaud Remay⁴*, Olivier Raymond¹*, Sandrine Balzergue²*, Aurélie Chauvet¹, Marion Maene¹, Yann Pécrix³, Shu-Hua Yang¹, Julien Jeauffre⁴, Tatiana Thouroude⁴, Véronique Boltz¹, Marie-Laure Martin-Magniette², Stéphane Janczarski¹, Fabrice Legeai⁵, Jean-Pierre Renou²,⁴, Philippe Vergne¹, Manuel Le Bris³, Fabrice Foucher⁴, Mohammed Bendahmane¹*

¹ Laboratoire Reproduction et Développement des Plantes, Institut Nationale de la Recherche Agronomique, Centre National de la Recherche Scientifique, Ecole Normale Supérieure, Lyon, France, ² Unité de Recherche en Génomique Végétale, Institut Nationale de la Recherche Agronomique, Centre National de la Recherche Scientifique, Evry, France, ³ Institut Méditerranéen d'Ecologie et de Paléoécologie, Centre National de la Recherche Scientifique, Université Paul Cézanne-Aix-Marseille III, Marseille, France, ⁴ UMR Génétique et Horticulture, Institut Nationale de la Recherche Agronomique, Agrocampus Ouest, Université d'Angers, Beaucouzé, France, ⁵ UMR Bio3P IRISA Equipe Symbiose Campus de Beaulieu, Institut Nationale de la Recherche Agronomique, Rennes, France

## Abstract

Cultivated for centuries, the varieties of rose have been selected based on a number of flower traits. Understanding the genetic and molecular basis that contributes to these traits will impact on future improvements for this economically important ornamental plant. In this study, we used scanning electron microscopy and sections of meristems and flowers to establish a precise morphological calendar from early rose flower development stages to senescing flowers. Global gene expression was investigated from floral meristem initiation up to flower senescence in three rose genotypes exhibiting contrasted floral traits including continuous versus once flowering and simple versus double flower architecture, using a newly developed Affymetrix microarray (Rosa1_Affyarray) tool containing sequences representing 4765 unigenes expressed during flower development. Data analyses permitted the identification of genes associated with floral transition, floral organs initiation up to flower senescence. Quantitative real time PCR analyses validated the mRNA accumulation changes observed in microarray hybridizations for a selection of 24 genes expressed at either high or low levels. Our data describe the early flower development stages in *Rosa* sp, the production of a rose microarray and demonstrate its usefulness and reliability to study gene expression during extensive development phases, from the vegetative meristem to the senescent flower.

**Citation:** Dubois A, Remay A, Raymond O, Balzergue S, Chauvet A, et al. (2011) Genomic Approach to Study Floral Development Genes in *Rosa* sp.. PLoS ONE 6(12): e28455. doi:10.1371/journal.pone.0028455

**Editor:** Miguel A. Blazquez, Instituto de Biología Molecular y Celular de Plantas, Spain

**Received:** November 2, 2011; **Accepted:** November 8, 2011; **Published:** December 14, 2011

**Copyright:** © 2011 Dubois et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was funded by the "Biologie Végétale" and the "Génétique et Amélioration des Plantes" Departments of the French Institut National de la Recherche Agronomique, and by the Région Rhônes-Alpes. Dr. Maene, Dr. Pécrix and Dr. Remay were supported by funds from the Région Rhône Alpes (Dr. Maene), The Region PACA (Dr. Pécrix) and by a joint grant from Région Pays de la Loire and the French 'Institut National de la Recherche Agronomique" (Dr. Remay). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

Roses are widely used as garden ornamental plants and cut flowers. A few flowering traits of roses are essential for the plants commercial value. Examples of these traits are plant architecture, continuous flowering, flower development, function and senescence, scent biosynthesis, reproduction and resistance to biotic and abiotic stresses. However, little is known about the molecular mechanisms that control these traits. This dearth of information limits the scope of rational selection to improve the ornamental plants. During the past decade, using model species such as *Arabidopsis thaliana*, tobacco, *Brachypodium distachyon*, rice or maize, researchers significantly enhanced our understanding of the various aspects of plant development and resistance to biotic and abiotic stresses, and of the molecular and genetic pathways associated with these aspects. However, these model species are not suitable for the studies of other flowering traits such as

recurrent blooming, scent production and double flower character. Rose represents an interesting ornamental model species to address some of these aspects.

Cultivated roses have a very ancient history. The two major areas of rose domestication were China and the peri-mediterranean area encompassing part of Europe and Middle East, where *Rosa chinensis* Jacq. and *R. gallica* L. (respectively) were bred and contributed predominantly to the subsequent selection process. Artificial crossing between Asian and European roses gave birth to "modern rose cultivars". Although testimonies and historical records have documented major crosses that led to modern roses, the genetic basis on which the modern rose cultivars are established is still poorly understood [1]. It has been reported that about 8 to 20 species out of about 200 wild species have contributed to the origin of present cultivars [2,3,4].

In *Rosa* sp., EST sequencing has identified novel genes whose expression is associated with several rose traits [5,6] such as the

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

<footnote>* E-mail: mohammed.bendahmane@ens-lyon.fr</footnote>

<footnote>¶ These authors contributed equally to this work.</footnote>
scent associated genes O-methyltransferases and alcohol acetyl-
transferase and floral associated genes [6,7,8,9,10,11,12,13]. EST
sequences were also used to generate a rose DNA microarray
comprising 350 selected ESTs [14]. Using this microarray,
researchers discovered several novel floral initiation genes and
flower scent-related candidate genes (i.e. germacrene D-synthase
encoding genes) [15]. However, this array contains only a limited
number of sequences that represent genes expressed at late petal
development stages.

With publicly available rose gene sequences, we generated a microarray and studied the gene expression throughout floral development, from the initial floral transition to floral senescence. We created an annotated flower EST database corresponding to 4834 genes and used the sequences to develop an Affymetrix microarray. With this microarray, we compared the transcriptome at different floral development stages. We found a good correlation between the microarray data and real time quantitative RT-PCR (qPCR) data for selected genes whose expression coincides with early, mid and late flower development stages. This dataset can help identify new rose genes associated with floral initiation, flower development and senescence.

**Results and Discussion**

Staging the floral transition and flower development in
*Rosa* sp

Understanding the genetic basis of flower formation in
ornamental plants such as roses is particularly important for
future cultivar improvement. We first analyzed the visible
morphological modifications during the floral process, from the
vegetative meristem to the senescent flower using three rose
cultivars, *Rosa wichurana*, *R. chinensis* cv. Old Blush and *R. x hybrida*
cv. Félicité et Perpétue. *Rosa wichurana* and *R. chinensis* cv. Old
Blush, two diploid roses, are among the few roses genotypes that
were used in the numerous crossings and hybridizations to create
the modern roses [2,16]. For example *R. chinensis* cv. Old Blush
contributed major traits, like recurrent flowering and components
of the characteristic 'tea scent' of modern roses [5,9,17], and *R.
wichurana* is a non recurrent flowering rose that contributed the
climbing trait for some garden roses [17]. The third rose, *R. x
hybrida* cv. Félicité et Perpétue (FP) is a cultivated hybrid. These
three cultivars were chosen because they have very different
flowering habits. For example *R. chinensis* cv. Old Blush was chosen
to study floral organogenesis, maturation and senescence, as it
flowers all year long in our greenhouse at ENS, Lyon. However,
continuing flowering limits our ability to sample enough vegetative
meristems for transcriptome analyses. Therefore, to collect
sufficient number of meristems, we also chose non recurrent
flowering roses, *R. wichurana* and *R. x hybrida* cv. Félicité et Perpétue
in greenhouse and field conditions at INRA, Angers.

Rose flowers are composed of four organ types arranged in
whorls, from the outer to the inner sepals, petals, stamens and
carpels. Flower development stages have been determined for
model plants such as *A. thaliana* [18]. However, these development
stages cannot be directly applied to the rose flower development.
In contrast to *A. thaliana* flowers that are composed of four
concentric whorls, rose flowers are composed of one whorl of 5
sepals and multiple whorls of petals, of stamens and of carpels.
Furthermore, the floral architecture of modern roses differs from
that of wild-type roses. For instance, modern rose varieties exhibit
double flower character of high number of petals and modified
numbers of stamens and carpels, whereas wild-type roses have 5
petals. Scanning electron microscopy (SEM) was used to image
floral initiation in *Rosa* sp (Figure 1). Based on these imaging data,

we divided the floral initiation process into three stages. After bud
outgrowth, the vegetative meristem is dome-shaped and narrow
with leaf primordia on its flanks (Stage VM1 for vegetative
meristem; Figure 1A, a, d). This structure is typical of a vegetative
meristem as previously described [19]. Rapidly, when the new
stems have acquired three fully expanded leaves, the meristem
enlarges, emerges and leaf primordia are now invisible (Stage
VM2, Figure 1A, b, e). We defined this VM2 stage as "pre-floral
stage". Then, the meristem becomes floral characterized by a flat,
large and doming structure (Stage FM for floral meristem;
Figure 1A, c, f). These morphological changes were similar in
the non-recurrent flowering roses, *R. wichurana* and *R. x hybrida* cv.
Félicité et Perpétue. Similar enlargement and doming of the
meristem were observed during the floral initiation in other related
Rosaceae [20].

Sections of floral meristem and young flower buds (Figure 1A,
g-k) were used to define the floral organogenesis steps in *R.
chinensis* cv. Old Blush. Five morphologically distinct developmen-
tal stages were easily distinguished under a dissecting microscope.
At flower development stage 1, the floral bud is surrounded by
bracts, the floral meristem is flat and five sepal primordia are
visible. Floral organs subsequently form following a radial gradient
so that the most external organs are the more differentiated. At
stage 2, petal primordia are apparent on the flank of the
hynanthium. At development stage 3 stamens primordia appear
on the flank of the hynanthium while petal primordia continue
developing. At stage 4, carpel primordia are the last organs that
appear in the center of the hynanthium, while the other organs
continue developing. At stage 5, all floral organs are apparent, and
the hynanthium starts to sink below the perianth and stamens.
During the onward development stages the hynanthium continues
to form and the flower becomes clearly visible (Figure 1 B). The
four types of floral organs continue developing and flowers start
opening (VP stage for visible petals) (Figure 1 B). Then the flower
fully opens (OF stage for open flower), and finally senesces (SF
stage for senescing flowers).

Rose EST database creation and Rosa1_Affymetrix
custom array design

We collected the available rose genes sequences (ESTs and
mRNA) and built a comprehensive database. Using sequence
clustering, we generated a dataset comprising 4765 unique
sequences (clusters and singletons) and deposited them in http://
urgi.versailles.inra.fr/GnpSeq.

For most of the clusters, one representative EST was chosen
based the following criteria. Its sequence is larger than 600
nucleotides and preferably corresponding to the 5' end gene
sequence. Because the rose is highly heterozygous, such strategy
should prevent using chimerical sequences that might have been
obtained during the clustering process. However, 343 clusters did
not meet the criteria above. For these 343 clusters, two or more
ESTs representing the unique sequence were used. In total, 5175
unique rose EST sequences representing 4765 unique sequences
were used for the Rosa1_Affymetrix array design and a total of
6,289 probe sets including Affymetrix control probesets were
designed. The arrays were manufactured by Affymetrix (http://
www.affymetrix.com).

Array sequences annotation

We used the Blastx algorithm against the nr database to identify
the best protein hits for the 5175 unique rose sequences, and
analyzed these results using Blast2go software [21]. 3959
sequences (76.5%) produced a significant match with one or more
entry in the database. Among the 3959 sequences, 222 (5.6%)

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Rose flower development stages. A. (a) to (f): Morphology of the floral transition in one-time flowering roses (R. wichurana) Schematic representation of the different stages observed during the floral transition in spring is shown in the upper panel from a vegetative meristem (VM) to a floral meristem (FM). a to c: Light microscopy of cross section of meristems. d to f: Environmental scanning electron microscopy images. Black bar: 10 µm. (g) to (k): Rose flower organogenesis stages. Cross sections of floral meristem and young flower buds. Images representing initiation of sepals (stages 1, 2, g), petals (stage 2, i), stamens (stages 3, h) and carpels (stage 4, j). k: hypanthium starts introverting below the floral organs (stages 5). Black bar: 50 µm (g,h,i); 200 µm (j,k). B. Visible rose flower stages. Pictures of rose flowers at flower bud with visible petals (stage VP), open flower stage (OF) and senescing Flower stage (SF).
doi:10.1371/journal.pone.0028455.g001</figcaption></figure>

could not be mapped with GO terms and 3737 had at least 1 GO term. For 1439 sequences, full automatic annotations were obtained. Analysis of GO biological process mapping showed that out of these 1439 sequences, 700 (48.6% of mapped sequences) were annotated as involved in primary metabolism processes and only 43 were annotated as putative secondary metabolism genes. 120 sequences (8.33% of mapped sequences) were mapped with the GO:0010468 annotation corresponding to regulation of gene expression. GO molecular function analysis showed that 38 sequences (2.6% of mapped sequences) had putative transcription factor activity (GO:00037000). The complete list of these sequences represented in the array, giving the first Blastx hit, the Blast2go computed annotation and gene ontology, is shown in Table S1. About 23.5% of the rose sequences produced no significant Blast hit in the gene databases. It is likely that the sequences of these genes have diverged far enough to render the annotation difficult. These highly divergent genes may have evolved functions that are be specific to the Rosa genus or Rosaceae family and are therefore of particular interest.

## Gene expression associated with rose floral initiation

We analyzed the transcriptomes of *R. wichurana* (*Rw*) and *R. x hybrida* *cv.* Félicité et Perpétue (FP) during floral initiation. Specifically, we compared vegetative (VM1) to pre-floral (VM2) stages and pre-floral to floral (FM) stages (Figure 2A). Such comparisons can uncover on genes potentially involved in the control of floral initiation. The rationale is that the genes up-regulated between vegetative and pre-floral buds are expected to be putative floral activators. Conversely, genes repressed between vegetative and pre-floral stages are expected to be putative floral inhibitors.

824 genes in *R. wichurana* and 652 genes in *R. x hybrida* *cv.* Félicité et Perpétue had a dynamic expression pattern between vegetative meristem (VM1) and pre-floral meristem (VM2) (Tables S2 and S3). Between VM1, VM2 and floral meristem (FM) stages, 302 (*Rw*) and 104 (FP) of these genes continued to be differentially expressed. During the VM1 to VM2 transition, 336 (*Rw*) and 301 (FP) genes were up-regulated between vegetative and floral stages, hence they represent candidates associated with floral initiation. 488 (*Rw*) and 351 (FP) genes were down-regulated and they are thus potential floral initiation repressors (Tables S2 and S3). To increase the confidence in the discovery of genes associated with floral induction, the overlapping genes from both datasets (*Rw* and FP) were selected. 258 differentially expressed genes during the VM1 to VM2 transition were common between FP and *Rw* samples. Among these genes, 222 out of 258 (86%) presented the following expression pattern. 131 genes are down-regulated between VM1 and VM2 stages and are thus putative floral repressors (top list in Table 1 and complete list in Table S4A). 91 gene are up-regulated between VM1 and VM2 stages and are thus putative floral activators (top list in Table 1 and complete list in Table S4B). Altogether, these genes are interesting candidates for studying floral initiation in *Rosa* sp.

Among the putative rose floral activators, the expression of the putative rose homologues of *SOC1* (*RhSOC1*) and *APETALA1* (*RhAPI*) were induced during the floral initiation both in *R. wichurana* and in *R. x hybrida* *cv.* Félicité et Perpétue (Tables S2 and S3; Figure 3), in agreement with previously reported data [13]. Therefore, like in *Arabidopsis* [22,23], in *Rosa* sp the expression of *RhSOC1* and *RhAPI* suggests that these genes may have similar function as floral integrator and floral meristem identity regulator, respectively. Among the genes that were differentially expressed in

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Description of the comparisons performed using micrarrays. **A.** To identify genes associated with floral initiation in *Rosa* using *R. wichurana* (*Rw*), *R. x hybrida* *cv.* Félicité et Perpétue (*FP*); Comparisons were done in the 2 genotypes; *VM1*: vegetative meristem stage; *VM2*: pre-floral meristems; *MF*: floral meristem. **B.** Schematic representation showing the rose flower development stages from flower organogenesis (stage 1) to onset of senescing flowers (stage SF). Arrows indicates the different transcriptome comparisons. *VP*: flower bud with visible petals; *OF*: open flower; *SF*: Senescing flower.
doi:10.1371/journal.pone.0028455.g002</figcaption></figure>

both roses during floral initiation, six (BI978989, BI978732, BI978794, EC589388, BQ104046, EC586448) showed similarities to genes involved in auxin transport or auxin signalling. Two auxin-repressed homologues (BI978989 and BI978794) were down-regulated and two auxin-induced homologues (BI978732 and BI978794) were up-regulated during the floral initiation process in *Rw* and *FP*, suggesting dynamic auxin signalling in the rose apex during the floral initiation and the organogenesis of the inflorescence meristems. Auxin and ethylene often interact synergistically [24–25]. We found genes involved in ethylene signalling were down-regulated during floral initiation in *Rw* and *FP*. These genes (EC586386 and AY919867) showed similarities with *EIN* and *EIL* genes [26]. *EIN* and *EIL* transcription factors are positive regulators of the ethylene signalling [27]. In *Arabidopsis*, ethylene delayed flowering as *acs* mutant flowered later [28]. In addition, during the floral initiation in *Rw*, two genes showing similarity with ethylene synthesis gene, ACC oxydase (AF441282) and ACC synthase (BQ105189) are down-regulated. Therefore, during the floral initiation, decrease in ethylene production may lead to diminution of EIN/EIL transcription factor and reduction of the ethylene signalling. These expression data suggest that ethylene and auxin may be involved in floral initiation process in rose although further experiments will be necessary to validate these hypotheses.

## Gene expression associated with rose floral development

We harvested six pools of samples corresponding to different flower development stages in *R. chinensis* *cv.* Old Blush (Figure 1) and compared the transcriptome in successive stages (Figure 2B). We found three distinct groups with common genes (T-test). These groups corresponded to early, mid and late floral development (Figure 2B). A total of 135, 401 and 456 sequences appeared

significantly and differentially regulated at least once during early, mid and late flower development stages, respectively.

To validate and evaluate the accuracy of the microarray data, we performed quantitative real-time PCR (qPCR). Twenty four genes were selected from the microarray transcriptomics comparisons based on previous bibliographic reports and/or deregulation levels, then, using qPCR, we further characterized the expression profiles (Figure 3; Figure S1). The correlation between the microarray results and those obtained by qPCR was assessed by calculating the Pearson's product moment correlation coefficient [47,48] (Table S5). Pearson's correlation coefficient was calculated between each pair of fold change as estimated by microarray and qPCR experiments. The statistical significance of each Pearson's correlation coefficient was assessed using the cor.test routine in R. A global correlation coefficient of 0.858 calculated by the average of every gene was observed. These results indicate that our microarrays are able to detect consistently both low and high fold-changes with high accuracy in different experimental conditions (Table S5).

## Transcriptome analyses during early flower development

135 genes were differentially expressed at during early floral organogenesis. Among these genes, 46 were found differentially expressed between stages 1+2 and 3+4 and 105 genes were differentially expressed between stages 3+4 and 5 (Table 2 and Table S6). An ACC synthase (AY803737) putative homologue was among the highly up-regulated genes between stages 1+2 and 3+4. In *Arabidopsis*, there are nine ACC synthases, many of which are expressed in the flower [29,30]. The floral organ identity MADS-box encoding genes [31,32,33], such as an *APETALA3* homologue (*RhTM6/MASAKOB3 AB055966*, Figure 3), the *AGAMOUS* ortholog (*RhAG*, *AB025645*, Figure 3), or the rose *PISTILLATA*

<figure><img src="image_7.png" /></figure>
<table><caption>**Table 1.** Top list of putative floral repressors and activators shared between *R. x. wichurana* and "Félicité et Perpétue".</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th colspan="2">R. x. wichurana</th><th colspan="2">"Félicité et Perpétue"</th></tr><tr><th>Log(ratio) VM2/WM1</th><th>Log(ratio) FM/WM2</th><th>Log(ratio) VM2/WM1</th><th>Log(ratio) FM/WM2</th></tr></thead><tbody><tr><td colspan="6"><strong>repressors</strong></td></tr><tr><td>EC58630</td><td>0,00</td><td>-2,77</td><td>-0,95</td><td>-1,03</td><td>-0,71</td></tr><tr><td>BQ104485</td><td>(Q5NE18) Formate dehydrogenase</td><td>-2,75</td><td>1,05</td><td>-1,14</td><td>-</td></tr><tr><td>EC589917</td><td>(Q3T923) Fra a 1 allergen (Fra a 1-A allergen)</td><td>-2,08</td><td>-0,34</td><td>-4,93</td><td>0,20</td></tr><tr><td>CF349421</td><td>(Q7XHM6) Hypothetical protein OSJNBb0095H08.9</td><td>-1,98</td><td>1,13</td><td>-0,75</td><td>-0,29</td></tr><tr><td>CF349812</td><td>(Q8H7G2) Hypothetical protein (Q8H7G2_ARATH)</td><td>-1,97</td><td>1,22</td><td>-1,00</td><td>0,32</td></tr><tr><td>EC587235</td><td>0,00</td><td>-1,81</td><td>0,76</td><td>-0,84</td><td>0,46</td></tr><tr><td>EC587845</td><td>(Q8L5Z1) Hypothetical protein At1g33810 (Q8L5Z1_ARATH)</td><td>-1,79</td><td>0,03</td><td>-1,65</td><td>0,41</td></tr><tr><td>CF349438</td><td>(O81644) Villin-2 (VIL12_ARATH)</td><td>-1,78</td><td>2,21</td><td>-0,87</td><td>0,07</td></tr><tr><td>CF349322</td><td>(Q1RST0) Peptidase S1 and S6</td><td>-1,76</td><td>2,40</td><td>-2,10</td><td>0,01</td></tr><tr><td>CF349916</td><td>(Q564G6) Galactomannan galactosyltransferase</td><td>-1,72</td><td>0,60</td><td>-1,01</td><td>-0,24</td></tr><tr><td>EC587239</td><td>0,00</td><td>-1,71</td><td>1,57</td><td>-1,68</td><td>0,36</td></tr><tr><td>BQ105944</td><td>ATP-dependent Clp protease ATP-binding subunit</td><td>-1,70</td><td>1,69</td><td>-1,08</td><td>-0,13</td></tr><tr><td>BQ105308</td><td>(Q533V0) Phospholipase D alpha (EC 3.1.4.4)</td><td>-1,69</td><td>1,85</td><td>-0,84</td><td>0,31</td></tr><tr><td>CF349664</td><td>(Q1S2R3) GIGANTEA protein (Q1S2R3_MEDTR)</td><td>-1,68</td><td>-0,13</td><td>-1,27</td><td>0,16</td></tr><tr><td>BI977439</td><td>(Q8L553) SCARECROW transcriptional regulator-like</td><td>-1,67</td><td>1,13</td><td>-1,29</td><td>-0,19</td></tr><tr><td>EC586479</td><td>0,00</td><td>-1,66</td><td>1,40</td><td>-1,24</td><td>-0,24</td></tr><tr><td>BQ104603</td><td>0,00</td><td>-1,65</td><td>2,31</td><td>-1,19</td><td>1,00</td></tr><tr><td>EC588955</td><td>0,00</td><td>-1,61</td><td>0,19</td><td>-1,01</td><td>0,67</td></tr><tr><td>BQ106662</td><td>(O04136) Homeobox protein knotted-1-like 3 (KNAP3))</td><td>-1,58</td><td>1,09</td><td>-1,02</td><td>-0,17</td></tr><tr><td>CF349422</td><td>(Q9SWH0) Plasma membrane proton ATPase</td><td>-1,53</td><td>1,54</td><td>-1,03</td><td>0,13</td></tr><tr><td>BQ106489</td><td>0,00</td><td>-1,52</td><td>0,98</td><td>-2,22</td><td>-</td></tr><tr><td>EC586088</td><td>(Q41695) Pectinacetylesterase precursor</td><td>-1,48</td><td>1,58</td><td>-1,32</td><td>0,86</td></tr><tr><td>EC588764</td><td>(Q2HTG1) GTP-binding signal recognition particle SRP54</td><td>-1,48</td><td>1,09</td><td>-1,26</td><td>-0,37</td></tr><tr><td>BI977401</td><td>(Q2AAC8) Cysteine proteinase</td><td>-1,46</td><td>0,60</td><td>-1,50</td><td>0,13</td></tr><tr><td>BQ104821</td><td>(Q2R3E0) Alpha-mannosidase</td><td>-1,46</td><td>1,43</td><td>-0,74</td><td>-0,11</td></tr><tr><td>BQ103923</td><td>(Q84V96) Aldehyde dehydrogenase 1 precursor</td><td>-1,46</td><td>0,74</td><td>-1,64</td><td>0,58</td></tr><tr><td>BQ104041</td><td>0,00</td><td>-1,45</td><td>1,98</td><td>-1,47</td><td>0,43</td></tr><tr><td>EC587517</td><td>(Q71BZ1) Type-B response regulator (Q71BZ1_CATRO)</td><td>-1,42</td><td>1,58</td><td>-0,70</td><td>-</td></tr><tr><td>EC588090</td><td>0,00</td><td>-1,42</td><td>-</td><td>-2,18</td><td>0,48</td></tr><tr><td>BQ105490</td><td>0,00</td><td>-1,37</td><td>1,83</td><td>-5,92</td><td>-0,84</td></tr><tr><td>BQ103990</td><td>(Q8RWI9) Hypothetical protein At3g21090 (Q8RWI9_ARATH)</td><td>-1,36</td><td>1,32</td><td>-1,01</td><td>-0,26</td></tr><tr><td>EC586608</td><td>0,00</td><td>-1,30</td><td>0,91</td><td>-0,99</td><td>0,38</td></tr><tr><td>BI978794</td><td>(Q05349) Auxin-repressed 12.5 kDa protein (12KD_FRAAN)</td><td>-1,30</td><td>0,49</td><td>-2,62</td><td>-0,80</td></tr><tr><td>CF349291</td><td>(Q9SGU9) Similar to O-succinylhomoserine sulfhydrylase</td><td>-1,26</td><td>0,89</td><td>-0,96</td><td>0,01</td></tr><tr><td>EC586448</td><td>(Q94BT2) Auxin-induced in root cultures protein 12-1,26 precursor</td><td>-1,26</td><td>-0,57</td><td>-0,77</td><td>-0,72</td></tr><tr><td colspan="6"><strong>activators</strong></td></tr><tr><td>BI978967</td><td>(Q6Z2K3) Putative Avr9/Cf-9 rapidly elicited</td><td>1,06</td><td>0,27</td><td>0,75</td><td>-0,04</td></tr><tr><td>BI977621</td><td>(Q8L5J6) Expansin 3 (Q8L5J6_MALDO)</td><td>1,08</td><td>-0,50</td><td>1,01</td><td>0,98</td></tr><tr><td>BQ104361</td><td>(Q650W6) Putative nucleic acid-binding prot.</td><td>1,08</td><td>-</td><td>1,40</td><td>-0,29</td></tr><tr><td>EC588171</td><td>(Q1SZF1) Allergen V5/Tpx-1 related</td><td>1,11</td><td>-0,75</td><td>1,29</td><td>0,79</td></tr></tbody></table>

<figure><img src="image_8.png" /></figure>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th colspan="2">R. x. wichurana</th><th colspan="2">"Félicité et Perpétue"</th></tr><tr><th>Log(ratio) VM2/WM1</th><th>Log(ratio) FM/WM2</th><th>Log(ratio) VM2/WM1</th><th>Log(ratio) FM/WM2</th></tr></thead><tbody><tr><td>EC586116</td><td>0,00</td><td>1,12</td><td>-1,17</td><td>1,14</td><td>0,34</td></tr><tr><td>EC589388</td><td>(Q1SHH7) Auxin responsive SAUR protein</td><td>1,14</td><td>-0,02</td><td>1,52</td><td>0,47</td></tr><tr><td>BI978946</td><td>(Q93Z01) AT5g58730</td><td>1,20</td><td>-1,03</td><td>1,55</td><td>-0,17</td></tr><tr><td>RoAGL20</td><td>(Q7Y137) MADS-box protein PTM5</td><td>1,22</td><td>-</td><td>0,97</td><td>0,14</td></tr><tr><td>BQ105514</td><td>0,00</td><td>1,23</td><td>-0,74</td><td>0,72</td><td>-0,02</td></tr><tr><td>BI977348</td><td>(Q94AQ7) Hypothetical protein At5g11280</td><td>1,23</td><td>0,50</td><td>0,73</td><td>-0,47</td></tr><tr><td>BQ103904</td><td>(Q41696) Cysteine protease precursor</td><td>1,25</td><td>-2,35</td><td>2,55</td><td>0,20</td></tr><tr><td>EC589855</td><td>0,00</td><td>1,27</td><td>-0,03</td><td>0,71</td><td>-0,44</td></tr><tr><td>BI978115</td><td>(Q84W81) Hypothetical protein At5g49800</td><td>1,27</td><td>0,34</td><td>1,85</td><td>-0,57</td></tr><tr><td>EC586690</td><td>Q2QXK7) F-box domain, putative</td><td>1,31</td><td>-0,91</td><td>1,36</td><td>-0,02</td></tr><tr><td>EC588294</td><td>(Q1S0D0) Glyoxalase/bleomycin resistance protein</td><td>1,32</td><td>-0,46</td><td>0,77</td><td>1,19</td></tr><tr><td>EC588783</td><td>(Q9LUC1) Putative protein At3g14740</td><td>1,34</td><td>0,27</td><td>1,19</td><td>-0,31</td></tr><tr><td>RoAP1a</td><td>(Q283Q1) APETALA1 protein</td><td>1,38</td><td>-0,40</td><td>2,06</td><td>0,88</td></tr><tr><td>EC587486</td><td>0,00</td><td>1,42</td><td>-0,62</td><td>1,38</td><td>-0,35</td></tr><tr><td>BI978732</td><td>(P32293) Auxin-induced protein 22A</td><td>1,47</td><td>-0,36</td><td>1,27</td><td>1,28</td></tr><tr><td>BQ104100</td><td>0,00</td><td>1,55</td><td>-0,94</td><td>1,49</td><td>1,03</td></tr><tr><td>BQ105108</td><td>(O65744) GDP dissociation inhibitor</td><td>1,63</td><td>-3,09</td><td>2,13</td><td>0,17</td></tr></tbody></table>

Log(ratio) of intensities are represented, italicized numbers represent ratios for which the p-value of the Bonferroni test was higher than 0.05. -: no value could be calculated.

doi:10.1371/journal.pone.0028455.t001

ortholog (RhPI/MASAKO BP, AB038462), were among the genes whose expression was up-regulated between stages 1+2 and 3+4 or between stages 3+4 and 5. Interestingly, genes that are predicted to have functions in cell wall remodelling, such as putative extracellular lipases (BQ106293, EC586717, EC588243, BI978064, BI977386, BQ105800), xyloglucan endotransglucosylase/hydrolase 2 (XTH2, DQ320658) [34], expansins (BI977621, EC589557), putative pectin esterase (BQ105504) and pectate lyase (BQ103887, BQ105987) were up-regulated between stages 3+4 and 5. This result supports the idea that very active cell wall remodeling coincides with the beginning of organ elongation that occurs mainly at stage 5. A putative gibberellin 2-oxidase (BQ105545) was up-regulated early during flower development. In *Arabidopsis*, a similar up-regulation of genes implicated in gibberellins metabolism and signaling have been described at early floral development [35,36]. In agreement with previously published data, our microarray analysis suggests that gibberellins are important during early floral development of rose plants [13,37]. Among the genes that showed strong down-regulation between stages 1+2 and 3+4, we found the putative orthologues of *PERIANTHIA* (PAN), *API* and *SOC1* (AGL20). In *Arabidopsis*, *PAN*, *API* and *SOC1* are expressed in the floral meristem, but their expression is down-regulated in the subsequent steps during floral organs differentiation [36,38,39,40], hence in agreement with the observed down-regulation of the rose homologues between flower development stages 1+2 and 3+4.

## Early to late floral development transition

Sequences corresponding to 401 genes were detected as differentially regulated between stages 5 and VP. Among these genes, 233 were down-regulated and 168 were up-regulated (see Table 3 for a selection of genes and Table S7 for full list). Genes

that exhibit strong similarities to genes involved in carotene, flavonoid and anthocyanin biosynthesis are up-regulated between stages 5 and VP. Among these genes, putative phytoene synthase (BI979026), zeta carotene desaturase (CF349648), lycopene beta-cyclase (BQ105122) are likely to be involved in carotenoid biosynthesis. The expression of UDP-glucose anthocyanidin-o-glucosyltransferase (AB201048/RhGT1), previously involved in anthocyanin synthesis [41], was strongly up-regulated. A similar strong up-regulation was observed for genes encoding putative phenylalanine ammonia-lyase (BQ105227), chalcone synthase (EC587811), flavonol synthase (AB038247) and anthocyanidin synthase (BI977949) (Figure 3). Altogether, these genes are likely good candidates involved in anthocyanins biosynthesis in rose petals.

Interestingly, genes predicted to encode five putative cyclins (EC586028, EC586517, EC587578, EC588351, and EC588489) and a putative cyclin dependent kinase (EC589228) are strongly down-regulated during floral organ morphogenesis. This down-regulation may reflect the transition from mitotic growth to post-mitotic growth where floral organs grow through cell expansion. Recently, Vanneste *et al.* showed that the transcriptional down-regulation of A2 type cyclins is a direct link between developmental programming and cell-cycle exit in *Arabidopsis thaliana* [42].

Fifteen genes encoding putative transcription factors were up-regulated, while nine were down-regulated. Among the up-regulated transcription factors, we found the putative orthologue of *SHP* (AB025643) [32] and a putative NAC domain protein (BI978992, Figure 3). BI978992 is homologous to *Arabidopsis NAC2*, a gene expressed in ovule integuments. The differential expression of *NAC2* between stages 5 and visible petals (VP) suggests its putatively conserved function with the *Arabidopsis NAC2*. Three putative MYB transcription factors were also up-

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Real time quantitative RT-PCR (qPCR) analysis of six selected differentially expressed genes during rose floral organogenesis, floral opening and senescence in *R. chinensis* cv. Old Blush. qPCR data (black histograms) are compared to the microarray hybridization data (white histograms). Microarray data is presented regardless of Bonferroni test success. Each pair of histograms represent successive comparisons between floral development stages 1+2, 3+4, 5, visible petals (VP), open flower (OF) and senescing flowers (SF). doi:10.1371/journal.pone.0028455.g003</figcaption></figure>

regulated (CF349636, BQ104100 and BI978095, Figure 3). These rose MYBs may be involved in organ elongation, as they share about 67% protein sequence similarity with AtMYB21, known to be involved in gibberellins/jasmonate-mediated control of stamen filament elongation [43].

## Late floral development

456 genes were differentially regulated at least once during the late phases of floral development, i.e. from visible petal (VP) stage to senescent flower (SF) stage. Most of these genes showed similar expression pattern when we compared stages VP to OF (open

flower) or stages VP to SF (See Table 4 for top list, and Table S8 for full data). This result indicates that the transcriptome becomes less dynamic at senescence stages and thus not so many differences are detected when comparing samples OF and SF to the VP sample. Gene ontology analysis showed that among the up-regulated genes, the three GO terms chlorophyll catabolic process, heterocycle catabolic process and cellular nitrogen compound catabolic process were significantly overrepresented as compared to the whole annotated set; the four GO terms nucleus, macromolecule biosynthetic process, intracellular non-membrane-bounded organelle and ribonucleoprotein complex were

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 2.** List of selected floral organogenesis associated genes in *R. chinensis* cv Old Blush.</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th colspan="2">R. chinensis cv Old Blush</th></tr><tr><th>Stages 3+4 vs 1+2</th><th>Stages 5 vs 3+4</th></tr></thead><tbody><tr><td>AY803737</td><td>Rosa hybrid cultivar 1-aminocyclopropane-1-carboxylase synthase 2 (ACS2)</td><td>2,99</td><td>-1,21</td></tr><tr><td>AB055966</td><td>Rosa rugosa MASAKO B3 mRNA for MADS-box protein,</td><td>2,67</td><td>1,04</td></tr><tr><td>AB025645</td><td>Rosa rugosa MASAKO C2 mRNA for MADS-box protein,</td><td>2,94</td><td>1,37</td></tr><tr><td>CF349463</td><td>(Q1S9M3) Lipase, active site (Q1S9M3_MEDTR)</td><td>2,68</td><td>-</td></tr><tr><td>BI978064</td><td>(Q9M8Y5) Putative GDSL-motif lipase/acylhydrolase (Q9M8Y5_ARATH)</td><td>2,10</td><td>1,17</td></tr><tr><td>BI977386</td><td>(Q9M8Y5) Putative GDSL-motif lipase/acylhydrolase (Q9M8Y5_ARATH)</td><td>1,99</td><td>1,08</td></tr><tr><td>EC586717</td><td>(Q1S3U7) Lipolytic enzyme, GDSL-L (Q1S3U7_MEDTR)</td><td>1,69</td><td>1,16</td></tr><tr><td>BQ105800</td><td>(Q1SAY6) Lipolytic enzyme, GDSL (Q1SAY6_MEDTR)</td><td>2,76</td><td>0,91</td></tr><tr><td>DQ320658</td><td>Rosa x borboniana xyloglucan endotransglucosylase/hydrolase 2 (Xth2)</td><td>2,55</td><td>0,74</td></tr><tr><td>BI977621</td><td>(Q8L5J6) Expansin 3 (Q8L5J6_MALDO)</td><td>-0,89</td><td>1,26</td></tr><tr><td>EC589557</td><td>(Q9SBT1) Expansin (Q9SBT1_FRAAN)</td><td>0,48</td><td>1,20</td></tr><tr><td>BQ105987</td><td>(Q94FT6) Pectate lyase B (Fragment) (Q94FT6_FRAAN)</td><td>0,75</td><td>1,27</td></tr><tr><td>BQ103887</td><td>(Q52PJ2) Ripening-related pectate lyase (Q52PJ2_MANIN)</td><td>1,21</td><td>1,11</td></tr><tr><td>BQ105504</td><td>(Q7X9B1) Pectinesterase (EC 3.1.1.11) (Q7X9B1_FRAAN)</td><td>1,44</td><td>1,35</td></tr><tr><td>BQ105545</td><td>(Q4W8C3) Gibberellin 2-oxidase (Q4W8C3_PHAAN)</td><td>0,42</td><td>-1,35</td></tr><tr><td>RoPAN</td><td>(Q9SX27) Putative bZIP transcription factor, PERIANTHIA (Q9SX27_ARATH)</td><td>-1,94</td><td>-2,08</td></tr><tr><td>RoAGL20</td><td>(Q7Y137) POPTM (Q7Y137) MADS-box protein PTM5</td><td>-2,77</td><td>-0,03</td></tr><tr><td>RoAP1b</td><td>(Q2XUP6) MADS-box protein</td><td>-0,98</td><td>-3,15</td></tr></tbody></table>

Log(ratio) of intensities are represented, italicized numbers represent ratios for which the p-value of the Bonferroni test was higher than 0.05.

doi:10.1371/journal.pone.0028455.t002

underrepresented. We could identify two genes encoding stay-green protein homologues (BI978267 and BQ106457) that are strongly up-regulated upon petal elongation and remain highly expressed throughout the final petal senescing process. Stay-green proteins have a major role in chlorophyll and photosynthetic pigments degradation and have been repeatedly described to be associated with the processes of fruit ripening and organ senescence [44]. Surprisingly, no gene related to ethylene biosynthesis or signaling was detected as differentially expressed during late floral development. However the *RbXTH1* and *RbEXPA1* genes, both induced during ethylene-triggered and field abscission [34,45], were strongly up-regulated between VP and OF stages and remained as such in senescing flowers. Among the down-regulated genes, the two GO terms protein metabolic process and plasma membrane were underrepresented as compared to the whole set (whole microarray GO terms) and the eight GO terms acyltransferase activity, acyl-carrier-protein biosynthetic process, acyl carrier activity, cellular carbohydrate metabolic process, polysaccharide metabolic process, fatty acid biosynthetic process, lipase activity and defense response to fungus were overrepresented (Table 5). The enrichment in the latter set may represent a slowdown of general metabolic pathways at the onset of flower senescence. Similar results were reported in *A. thaliana* during organs senescence where a down-regulation of the photosynthetic machinery accompanied by a reduction in expression of many cell wall biosynthetic genes reflecting a cessation of growth during senescence [46].

## Conclusions

We established a calendar of the floral initiation and development for the rose and developed a rose microarray that harbors sequence from genes expressed during the floral transition and whole floral development process in *Rosa* sp, from initiation up

to senescing flowers. This microarray and the floral development calendar were successfully used to identify genes whose expression correlated with different flower development stages. These multiple datasets represent an extensive study of rose floral development. This resource can be helpful to select candidate genes potentially involved in different horticultural traits, such as flowering, floral architecture, scent production and emission, senescence and abscission. We used the microarray developed herein to identify genes whose expression is associated with some of these rose important traits, such as flower initiation, development and senescence. Rosa1_Affyarray harbors sequences from ESTs found in petals of different rose genotypes [5,14] (http://urgi.versailles.inra.fr/GnpSeq) and thus may be helpful to identify genes associated with other rose traits such as scent biosynthesis and/or emission genes. The rose is among the species that exhibit the highest scent complexity [47–48] [12] and some scent biosynthesis pathways are unique to the rose or not yet identified in other model species including other members of the Rosaceae genus [11,49]. QTLs have been identified to be associated to several important traits of the rose [50]. However, the heterozygous genome of the rose complicates the breeding programs to select for several traits simultaneously. The identification of genes whose expression correlates with important ornamental traits can facilitate and accelerate candidate gene identification for rose breeding by marker assisted selection or genomic selection. For example, this dataset can provide researchers with a useful resource on the expression of candidate genes within a given mapping interval. Furthermore, the rapidly progressing high throughput sequencing technologies should allow the generation of precise genetic maps for the rose that could be combined to refined transcriptomics approaches to identify the genes responsible for important horticultural traits in the rose, and allow subsequent marker-assisted selection.

<figure><img src="image_12.png" /></figure>
<table><caption>**Table 3.** List of selected genes associated with early to late flower development in *R. chinensis* cv Old Blush.</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th>R. chinensis cv Old Blush</th></tr><tr><th>5 vs PA</th></tr></thead><tbody><tr><td>BI978095</td><td>(P93474) Myb26</td><td>8,00</td></tr><tr><td>BI978992</td><td>(Q50J79) NAM-like protein</td><td>5,28</td></tr><tr><td>AB038247</td><td>Rosa hybrid cultivar 'Kardinal' FLS mRNA for flavonol synthase</td><td>4,67</td></tr><tr><td>BQ105122</td><td>(Q9SEA0) Lycopene beta-cyclase</td><td>4,30</td></tr><tr><td>EC587811</td><td>(Q84UT9) Chalcone synthase</td><td>3,27</td></tr><tr><td>BQ104100</td><td>MYB domain class transcription factor</td><td>3,01</td></tr><tr><td>AB025643</td><td>Rosa rugosa MASAKO D1 mRNA for MADS-box protein.</td><td>3,00</td></tr><tr><td>CF349648</td><td>(Q5W5X6) Zeta-carotene desaturase ZDS2</td><td>2,97</td></tr><tr><td>BI979026</td><td>(Q2VEY1) Putative phytoene synthase</td><td>2,89</td></tr><tr><td>BI977949</td><td>(Q5UL09) Anthocyanidin synthase</td><td>2,18</td></tr><tr><td>AB201048</td><td>RhGT1 UDP-glucose: anthocyanidin 5,3-O-glucosyltransferase,</td><td>2,18</td></tr><tr><td>CF349636</td><td>(Q9ATD1) GHMYB9</td><td>2,12</td></tr><tr><td>BQ105227</td><td>(Q9M567) Phenylalanine ammonia-lyase 2</td><td>2,11</td></tr><tr><td>EC586028</td><td>(Q9SNV1) Cyclin D3a (Fragment)</td><td>-2,12</td></tr><tr><td>AB201051</td><td>RhGT4 mRNA UDP-glucose: flavonol 3-O-glucosyltransferase</td><td>-2,18</td></tr><tr><td>EC587392</td><td>(Q85342) Putative anthocyanidine rhamnosyl-transferase</td><td>-2,45</td></tr><tr><td>EC587578</td><td>(Q6T2Z6) Cyclin d3</td><td>-3,83</td></tr><tr><td>EC586734</td><td>(Q08733) Aquaporin PIP1.3</td><td>-4,52</td></tr><tr><td>RhCyc2</td><td>(Q9SBQ4) CYCB1-1 protein</td><td>-4,65</td></tr><tr><td>EC588351</td><td>(Q9SBQ4) CYCB1-1 protein</td><td>-4,72</td></tr><tr><td>EC58848</td><td>(P93557) Mitotic cyclin</td><td>-4,77</td></tr><tr><td>EC589228</td><td>(Q94EX2) At1g76540/cyclin dependent kinase</td><td>-5,05</td></tr><tr><td>EC586517</td><td>(Q4JF78) Cyclin-dependent kinase B</td><td>-5,26</td></tr></tbody></table>

Log(ratio) of intensities are represented, for all ratios the p-value of the Bonferroni test was lower than 0.05.
doi:10.1371/journal.pone.0028455.t003

# Materials and Methods

## Plant material

*R. wichurana* was obtained from 'Jardin de Bagatelle' (Paris, France) and *R. x hybrida* cv. Félicité et Perpétue from the Loubert Nursery (Rosier sur Loire, France). Plants were grown outdoors on their own roots as previously described [13]. In spring, at different time points (see results), terminal parts of the growing shoot were harvested and partly dissected (removal of young leaves). *R. chinensis* cv. Old Blush was propagated by cuttings from the Lyon Botanical Garden. Plants were grown in the greenhouse with 16 h/8 h day/night and 25°C/14°C day/night temperature. No specific permits were required for the described filed studies, no specific permissions were required for these locations, the location is not privately owned or protected, and the field studies did not involve endangered or protected species.

## Light microscopy and SEM imaging of meristems and early flower development

Samples were dissected under a binocular stereomicroscope and then fixed in 4% glutaraldehyde (v/v) in 0.1 M phosphate buffer (pH 7.2) for 2 h at 4°C under vacuum. Samples were dehydrated in a graded ethanol series and embedded in Technovit 7100 [51]. Sections of 1.5 to 2.0 µm (Leica RM 2165 microtome) were stained with toluidine blue and examined under an Olympus BH2-RFC microscope coupled to a 3CCD Sony camera.

For scanning electron microscopy, terminal part of the shoot was carefully dissected. After a fixation in 4% glutaraldehyde (v/v), followed by post-fixation with osmium tetroxide, the sample was dehydrated in a graded alcohol series and in acetone. Dehydration was completed by critical point drying. Sample were then coated with gold (MED 020 BALTEC) and observed with a JEOL JSM-63017 scanning electron microscope.

## RNA samples preparation

Two independent biological replicates were produced for each samples at different stages. For each biological repetition and each point, RNA samples were obtained by pooling vegetative or floral tissue from at least five different plants. For *R. chinensis* cv. Old Blush samples, meristems or flowers were dissected and collected individually on plants at developmental growth stages, cultivated in greenhouse conditions as previously described [55]. For *R. wichurana* and *R. x hybrida* cv. Félicité et Perpétue, RNA was extracted from non-dissected buds, including either the vegetative meristem and its surrounding leaves or the pre-floral/floral meristem and its surrounding leaves and bracts. Total RNA was extracted using RNeasy Plant Mini Kit (Qiagen) according to the supplier's instructions.

## AFFYMETRIX Array hybridization

RNA samples were checked for their integrity on The Agilent 2100 bioanalyzer according to the Agilent Technologies (Wald-broon, Germany).

<figure><img src="image_13.png" /></figure>
<table><caption>**Table 4.** List of selected floral maturation and senescence associated genes in *R. chinensis* cv Old Blush.</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th colspan="2">R. chinensis cv Old Blush</th></tr><tr><th>PA vs FE</th><th>PA vs FS</th></tr></thead><tbody><tr><td>BI977502</td><td>Brassinosteroid-regulated protein BRU1 precursor</td><td>8,82</td><td>8,79</td></tr><tr><td>BI978598</td><td>Early light-induced protein</td><td>8,38</td><td>7,99</td></tr><tr><td>EC587309</td><td>0,00</td><td>7,70</td><td>5,82</td></tr><tr><td>BI977376</td><td>Putative zinc finger protein At1g68190</td><td>7,35</td><td>7,13</td></tr><tr><td>BI978143</td><td>0,00</td><td>6,61</td><td>7,00</td></tr><tr><td>EC587486</td><td>0,00</td><td>6,31</td><td>4,81</td></tr><tr><td>BI978750</td><td>0,00</td><td>6,30</td><td>5,55</td></tr><tr><td>BI978596</td><td>Hypothetical protein</td><td>6,29</td><td>6,32</td></tr><tr><td>BQ104828</td><td>0,00</td><td>6,24</td><td>5,51</td></tr><tr><td>BQ105724</td><td>0,00</td><td>6,09</td><td>5,65</td></tr><tr><td>EC586975</td><td>Glycosyltransferase NTGT5a</td><td>6,07</td><td>6,28</td></tr><tr><td>BQ106572</td><td>0,00</td><td>6,02</td><td>5,87</td></tr><tr><td>BI977926</td><td>18.5 kDa class I heat shock protein</td><td>5,83</td><td>-</td></tr><tr><td>EC588495</td><td>0,00</td><td>5,60</td><td>5,54</td></tr><tr><td>BI978508</td><td>Hypothetical protein</td><td>5,51</td><td>-</td></tr><tr><td>BQ104475</td><td>0,00</td><td>5,50</td><td>6,25</td></tr><tr><td>BQ103870</td><td>0,00</td><td>5,26</td><td>5,12</td></tr><tr><td>BI977873</td><td>Hypothetical protein At5g63130</td><td>5,25</td><td>5,18</td></tr><tr><td>BQ103973</td><td>Tryptophan synthase alpha chain</td><td>5,24</td><td>4,61</td></tr><tr><td>BI977634</td><td>Aux/IAA protein</td><td>5,21</td><td>5,20</td></tr><tr><td>BQ105490</td><td>0,00</td><td>5,07</td><td>5,85</td></tr><tr><td>BQ106477</td><td>Protein WUSCHEL-like</td><td>5,03</td><td>5,54</td></tr><tr><td>BQ106330</td><td>0,00</td><td>5,02</td><td>5,22</td></tr><tr><td>BQ106091</td><td>Hypothetical protein At2g42570</td><td>4,91</td><td>4,52</td></tr><tr><td>BI977302</td><td>Bzip transcription factor</td><td>4,82</td><td>4,83</td></tr><tr><td>BI978926</td><td>AT5g11580</td><td>4,81</td><td>5,46</td></tr><tr><td>CF349316</td><td>Putative NADH dehydrogenase</td><td>4,81</td><td>5,01</td></tr><tr><td>EC589818</td><td>Putative calmodulin-related protein</td><td>4,80</td><td>-</td></tr><tr><td>BI978132</td><td>0,00</td><td>4,74</td><td>4,55</td></tr><tr><td>BQ105726</td><td>Expansin-like protein</td><td>4,73</td><td>3,90</td></tr><tr><td>BQ104701</td><td>0,00</td><td>4,68</td><td>4,81</td></tr><tr><td>EC586683</td><td>0,00</td><td>4,68</td><td>5,52</td></tr><tr><td>DQ320657</td><td>Rosa×borboniana expansin protein (ExpA1) mRNA</td><td>4,17</td><td>4,38</td></tr><tr><td>EC589229</td><td>Probable xyloglucan endotransglucosylase/hydrolase protein 8 precursor</td><td>3,20</td><td>-</td></tr><tr><td>BI978267</td><td>Senescence-inducible chloroplast stay-green protein 1</td><td>2,71</td><td>2,42</td></tr><tr><td>BQ106457</td><td>Senescence-inducible chloroplast stay-green protein 2</td><td>-</td><td>2,34</td></tr><tr><td>BQ104919</td><td>Pectate lyase</td><td>-4,04</td><td>-6,18</td></tr><tr><td>EC588897</td><td>Laccase</td><td>-4,04</td><td>-4,86</td></tr><tr><td>EC588483</td><td>Alpha-D-xylosidase precursor</td><td>-4,06</td><td>-3,72</td></tr><tr><td>EC587152</td><td>0,00</td><td>-4,08</td><td>-4,46</td></tr><tr><td>EC589569</td><td>Serine carboxypeptidase, putative</td><td>-4,10</td><td>-3,48</td></tr><tr><td>EC586717</td><td>Lipolytic enzyme, GDSL</td><td>-4,13</td><td>-2,59</td></tr><tr><td>EC587284</td><td>RNA-binding region RNP-1</td><td>-4,15</td><td>-3,42</td></tr><tr><td>BI977461</td><td>0,00</td><td>-4,21</td><td>-4,40</td></tr><tr><td>EC586015</td><td>Cold-regulated LTCOR12</td><td>-4,22</td><td>-4,44</td></tr><tr><td>BI978116</td><td>Putative cell-wall P4 protein</td><td>-4,24</td><td>-</td></tr><tr><td>BI977348</td><td>Hypothetical protein At5g11280</td><td>-4,27</td><td>-4,99</td></tr><tr><td>BQ106043</td><td>Isoamylase isoform 3</td><td>-4,27</td><td>-4,01</td></tr><tr><td>EC589137</td><td>Putative alpha-glucosidase</td><td>-4,33</td><td>-4,07</td></tr></tbody></table>

<figure><img src="image_14.png" /></figure>
<table><caption>**Table 4. Cont.**</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">annotation</th><th colspan="2">R. chinensis cv Old Blush</th></tr><tr><th>PA vs FE</th><th>PA vs FS</th></tr></thead><tbody><tr><td>EC587785</td><td>0,00</td><td>-4,37</td><td>-3,91</td></tr><tr><td>EC586984</td><td>Putative beta-expansin</td><td>-4,38</td><td>-4,99</td></tr><tr><td>CF349422</td><td>Plasma membrane proton ATPase</td><td>-4,41</td><td>-3,64</td></tr><tr><td>BI977751</td><td>0,00</td><td>-4,49</td><td>-4,43</td></tr><tr><td>EC589098</td><td>Phosphoethanolamine N-methyltransferase 1</td><td>-4,49</td><td>-5,79</td></tr><tr><td>BQ106293</td><td>GDSL-motif lipase/hydrolase-like protein</td><td>-4,51</td><td>-6,32</td></tr><tr><td>CF349724</td><td>Glucosyltransferase-like protein</td><td>-4,74</td><td>-4,10</td></tr><tr><td>BQ106328</td><td>Pathogenesis-related transcriptional factor and ERF</td><td>-5,10</td><td>-</td></tr><tr><td>EC586884</td><td>Proline-rich protein APG-like</td><td>-5,11</td><td>-3,72</td></tr><tr><td>CF349712</td><td>Senescence-inducible gene protein</td><td>-5,12</td><td>-5,47</td></tr><tr><td>BI978135</td><td>Plant lipid transfer protein/Par allergen</td><td>-5,16</td><td>-7,04</td></tr><tr><td>CF349692</td><td>Putative alcohol oxidase</td><td>-5,17</td><td>-6,66</td></tr><tr><td>EC588080</td><td>Hypothetical protein (At2g35760/T20F21.5)</td><td>-5,32</td><td>-5,20</td></tr><tr><td>BI977262</td><td>Putative lipase</td><td>-5,48</td><td>-6,04</td></tr><tr><td>CF349791</td><td>Globulin-like protein (Fragment)</td><td>-5,67</td><td>-5,53</td></tr><tr><td>AB121046</td><td>phloroglucinol O-methyltransferase, complete cds</td><td>-6,13</td><td>-5,06</td></tr><tr><td>BI978206</td><td>0,00</td><td>-8,51</td><td>-</td></tr><tr><td>BI978064</td><td>Putative GDSL-motif lipase/acylhydrolase</td><td>-10,47</td><td>-10,51</td></tr><tr><td>BI977386</td><td>Putative GDSL-motif lipase/acylhydrolase</td><td>-11,60</td><td>-11,42</td></tr></tbody></table>

Log(ratio) of intensities are represented, for all ratios the p-value of the Bonferroni test was lower than 0.05.
doi:10.1371/journal.pone.0028455.t004

Two µg of total RNA were used to synthesize biotin-labeled cRNAs with the One-cycle cDNA synthesis kit (Affymetrix, Santa Clara, CA). Superscript II reverse transcriptase and T7-oligo (dT) primers were used to synthesize the single strand of cDNA at 42°C during 1 hour followed by the synthesis of the double stranded cDNA by using DNA ligase, DNA polymerase I and RNaseH during 2 hours at 16°C. Clean up of the double-stranded cDNA was performed with Sample Cleanup Module (Affymetrix) followed by *in vitro* transcription (IVT) in presence of biotin-labeled UTP using GeneChip® IVT labelling Kit (Affymetrix). Quantity of the labelled-cRNA with RiboGreen® RNA Quantification Reagent (Turner Biosystems, Sunnyvale, CA) was determined after cleanup by the Sample Cleanup Module (Affymetrix). Fragmentation of 10 µg of labelled-cRNA was carried out for 35 minutes at 94°C, followed by hybridization during 16 hours at 45°C to Affymetrix GeneChip® Rosa1 Genome Array representing approximately 4869 genes. After hybridization, the arrays were washed with 2 different buffers (stringent: 6× SSPE, 0.01% Tween-20 and non-stringent: 100 mM MES, 0.1 M [Na⁺], 0.01% Tween-20) and stained with a complex solution including Streptavidin R-Phycoerythrin conjugate (Invitrogen/molecular probes, Carlsbad, CA) and anti Streptavidin biotinylated antibody (Vectors laboratories, Burlingame, CA). The washing and staining steps were performed in a GeneChip® Fluidics Station 450 (Affymetrix). The Affymetrix GeneChip® Rosa1 Genome Arrays were finally scanned with the GeneChip® Scanner 3000 7G piloted by the GeneChip® Operating Software (GCOS).

## Statistical Analysis of Microarray Data

The data were normalized with the gcrma algorithm [52], available in the Bioconductor package [53]. To determine differen-

tially expressed genes, we performed a usual two group t-test that assumes equal variance between groups. The variance of the gene expression per group is a homoscedastic variance, where genes displaying extremes of variance (too small or too large) were excluded. The raw P values were adjusted by the Bonferroni method, which controls the Family Wise Error Rate (FWER) [54]. A gene is declared differentially expressed if the Bonferroni P-Value is less than 0.05.

## Data Deposition

All this steps were performed on Affymetrix plateform at INRA-URGV, Evry. The raw. CEL files were imported in R software for data analysis. All raw and normalized data are available through the CATdb database (AFFY_PetalDvt_Lyon_Rose, [55]) and from the Gene Expression Omnibus (GEO) repository at the National Center for Biotechnology Information (NCBI) [56], accession number GSE18357.

## Validation of genes expression using quantitative real-time PCR

Only genes that were involved in floral development were analyzed for microarray data validation. One microgram total RNA (treated with DNase) was used in a reverse transcription assay with RevertAid M-MuLV Reverse Transcriptase (Fermentas, Burlington, Ontario). Target cDNAs were quantified by qPCR using FastStart universal SYBR green master (Roche, Basel, Switzerland) on a Step-OnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA USA). Expression levels were normalized with *RhαTubuline*, *RhGAPDH* and *RhEF1α* reference genes. These genes were validated as reference genes using the GeNorm application [57]. Three independent biological replicates (pools of dissected flowers from at least 5 different plants) were used for each experiment and two qPCR technical replicates were

<figure><img src="image_15.png" /></figure>
<table><caption>**Table 5.** Gossip analysis of GO terms enrichment in late flower development dataset (genes that are differentially expressed at least once during floral maturation and senescence).</caption><thead><tr><th></th><th>GO Term</th><th>Name</th><th>FDR</th><th>FWER</th><th>single test p-Value</th><th># in test group</th><th># in reference group</th><th># non annotated test</th><th># non annotated reference group</th><th>Over/Under</th></tr></thead><tbody><tr><td rowspan="6">Late floral development upregulated genes</td><td>GO:0005634</td><td>nucleus</td><td>0.0</td><td>0.0</td><td>0.012</td><td>0</td><td>107</td><td>56</td><td>1290</td><td>under</td></tr><tr><td>GO:0009059</td><td>macromolecule biosynthetic process</td><td>0.0</td><td>0.0</td><td>0.028</td><td>0</td><td>88</td><td>56</td><td>1309</td><td>under</td></tr><tr><td>GO:0043232</td><td>intracellular non-membrane-bounded organelle</td><td>0.0</td><td>0.0</td><td>0.030</td><td>0</td><td>86</td><td>56</td><td>1311</td><td>under</td></tr><tr><td>GO:0030529</td><td>ribonucleoprotein complex</td><td>0.0</td><td>0.0</td><td>0.030</td><td>0</td><td>84</td><td>56</td><td>1313</td><td>under</td></tr><tr><td>GO:0015996</td><td>chlorophyll catabolic process</td><td>0.008</td><td>0.008</td><td>5.43E-5</td><td>3</td><td>0</td><td>53</td><td>1397</td><td>over</td></tr><tr><td>GO:0046700</td><td>heterocycle catabolic process</td><td>0.028</td><td>0.062</td><td>5.13E-4</td><td>3</td><td>2</td><td>53</td><td>1395</td><td>over</td></tr><tr><td rowspan="6">Late floral development downregulated genes</td><td>GO:0019538</td><td>protein metabolic process</td><td>0.0</td><td>0.0</td><td>0.009</td><td>5</td><td>225</td><td>73</td><td>1150</td><td>under</td></tr><tr><td>GO:0005886</td><td>plasma membrane</td><td>0.012</td><td>0.006</td><td>0.013</td><td>4</td><td>192</td><td>74</td><td>1183</td><td>under</td></tr><tr><td>GO:0008415</td><td>acyltransferase activity</td><td>0.016</td><td>0.025</td><td>1.62E-4</td><td>7</td><td>17</td><td>71</td><td>1358</td><td>over</td></tr><tr><td>GO:0042967</td><td>acyl-carrier-protein biosynthetic process</td><td>0.016</td><td>0.028</td><td>2.15E-4</td><td>7</td><td>18</td><td>71</td><td>1357</td><td>over</td></tr><tr><td>GO:0000036</td><td>acyl carrier activity</td><td>0.027</td><td>0.079</td><td>5.73E-4</td><td>3</td><td>1</td><td>75</td><td>1374</td><td>over</td></tr><tr><td>GO:0044262</td><td>cellular carbohydrate metabolic process</td><td>0.027</td><td>0.084</td><td>6.73E-4</td><td>15</td><td>99</td><td>63</td><td>1276</td><td>over</td></tr><tr><td></td><td>GO:0005976</td><td>polysaccharide metabolic process</td><td>0.027</td><td>0.086</td><td>7.08E-4</td><td>10</td><td>48</td><td>68</td><td>1327</td><td>over</td></tr></tbody></table>

The reference group that was used corresponds to the full annotated sequences (sequences with GO terms) of the microarray. doi:10.1371/journal.pone.0028455.t005

performed for each biological replicate. Primer sequences are available in Table S9. The correlation between the microarray results, and those obtained by qPCR was assessed by calculating the Pearson's product moment correlation coefficient [58,59].

## Supporting Information

**Figure S1** Real time quantitative RT-PCR (qPCR) analysis of 18 selected differentially expressed genes during rose floral organogenesis and senescence in *R. chinensis* cv Old Blush.
(TIFF)

**Table S1** Full array sequences annotation and ontology.
(XLSX)

**Table S2** Genes differentially expressed during floral initiation in *R. wichurana*.
(XLSX)

**Table S3** Genes differentially expressed during floral initiation in *R. x hybrida* cv. Félicité et Perpétue.
(XLSX)

**Table S4** List of genes repressed (A) or activated (B) during flower initiation.
(XLSX)

**Table S5** Microarray and qRT-PCR results of 25 selected genes with their replicate-level Pearson correlation.
(DOCX)

**Table S6** Genes differentially expressed during early floral organogenesis in *R. chinensis* cv. Old Blush.
(XLSX)

**Table S7** Genes differentially expressed during floral organ elongation in *R. chinensis* cv. Old Blush.
(XLSX)

**Table S8** Genes differentially expressed during flower opening and senescence in *R. chinensis* cv. Old Blush.
(XLSX)

**Table S9** Primers used in this study.
(DOC)

## Acknowledgments

We thank Judit Szecsi and Sylvie Baudino for critical reading of the manuscript. We thank Alexis Lacroix, Isabelle Desbouchages, Priscilla Angelot and N. Dousset and J. Chameau taking care of the plants, M. Thellier and Michel Chevalier for the histological analysis, S/Georgeault and R. Filmont for the SEM studies.

<figure><img src="image_16.png" /></figure>
## Author Contributions

Conceived and designed the experiments: MB AD OR. Performed the experiments: AD AR OR SB AC MM YP SHY JJ TT VB MLMM SJ JPR

## References

1. Martin M, Piola F, Chessel D, Jay M, Heizmann P (2001) The domestication process of the Modern Rose: genetic structure and allelic composition of the rose complex. Theoretical-and-Applied-Genetics 102: 398-404.

2. De Vries DP, Dubois L (1996) Rose breeding: past, present, prospects. Acta Horticulturae 424: 241-248.

3. Gudin S (2001) Rose breeding technologies. Acta Horticulturae 547: 23-26.

4. Reynders-Aloisi S, Boliereau P (1996) Characterisation of genetic diversity in genus Rosa by Randomly Amplified Polymorphic DNA. Acta Horticulturae 424: 253-259.

5. Channeliere S, Riviere S, Scalliet G, Szecsi J, Jullien F, et al. (2002) Analysis of gene expression in rose petals using expressed sequence tags. FEBS Lett 515: 35-38.

6. Foucher F, Chevalier M, Corre C, Soufflet-Freslon V, Legeci F, et al. (2008) New resources for studying the rose flowering process. Genome 51: 827-837.

7. Guterman I, Masci T, Chen X, Negre F, Pichersky E, et al. (2006) Generation of phenylpropanoid pathway-derived volatiles in transgenic plants: rose alcohol acetyltransferase produces phenylethyl acetate and benzyl acetate in petunia flowers. Plant Mol Biol 60: 555-563.

8. Lavid N, Wang J, Shalit M, Guterman I, Bar E, et al. (2002) O-methyltransferases involved in the biosynthesis of volatile phenolic derivatives in rose petals. Plant Physiol 129: 1899-1907.

9. Scalliet G, Journot N, Jullien F, Baudino S, Magnard JL, et al. (2002) Biosynthesis of the major scent components 3,5-dimethoxytoluene and 1,3,5-trimethoxybenzene by novel rose O-methyltransferases. Febs Letters 523: PII S0014-5793(0002)02956-02953.

10. Scalliet G, Lionnet C, Le Becch M, Dutron L, Magnard JL, et al. (2006) Role of petal-specific orcinol O-methyltransferases in the evolution of rose scent. Plant Physiology 140: 18-28.

11. Scalliet G, Piola F, Douady CJ, Rety S, Raymond O, et al. (2008) Scent evolution in Chinese roses. Proceedings Of The National Academy Of Sciences Of The United States Of America 105: 5927-5932.

12. Shalit M, Guterman I, Volpin H, Bar E, Tamari T, et al. (2003) Volatile ester formation in roses. Identification of an acetyl-coenzyme A. Geraniol/Citronellol acetyltransferase in developing rose petals. Plant Physiol 131: 1868-1876.

13. Remay A, Lalanne D, Thouroude T, Le Couviour F, Hibrand-Saint Oyant L, et al. (2009) A survey of flowering genes reveals the role of gibberellins in floral control in rose. Theor Appl Genet 119: 767-781.

14. Guterman I, Shalit M, Menda N, Piestun D, Dafny-Yelin M, et al. (2002) Rose scent: Genomics approach to discovering novel floral fragrance-related genes. Plant Cell 14: 2325-2338.

15. Guterman I, Masci T, Chen XL, Negre F, Pichersky E, et al. (2006) Generation of phenylpropanoid pathway-derived volatiles in transgenic plants: Rose alcohol acetyltransferase produces phenylethyl acetate and benzyl acetate in petunia flowers. Plant Molecular Biology 60: 555-563.

16. Gudin S (2000) Rose: Genetics and breeding. In: Plant Breeding Reviews 17: 159-189.

17. Krussmann G (1981) The complete book of roses. Portland: Timber Press. xii, 436 p.

18. Smyth DR, Bowman JL, Meyerowitz EM (1990) Early flower development in Arabidopsis. Plant Cell 2: 755-767.

19. Chimonidou D (2003) Morphology and Anatomy: 'Flower Development and Abscission zone'. ENCYCLOPEDIA OF ROSE SCIENCE. Amsterdam: Elsevier Academic Press. pp 504-512.

20. Foster T, Johnston R, Seleznyova A (2003) A morphological and quantitative characterization of early floral development in apple (Malus × domestica Borkh.). Ann Bot 92: 199-206.

21. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420-3435.

22. Amasino R (2010) Seasonal and developmental timing of flowering. Plant J 61: 1001-1013.

23. Fornara F, de Montaigu A, Coupland G (2010) Snapshot: Control of flowering in Arabidopsis. Cell 141: 550, 550 e551-552.

24. Bennett MJ, Swarup R, Perry P, Hagenbeek D, Van Der Straeten D, et al. (2007) Ethylene upregulates auxin biosynthesis in Arabidopsis seedlings to enhance inhibition of root cell elongation. Plant Cell 19: 2186-2196.

25. Kang BG, Newcomb W, Burg SP (1971) Mechanism of Auxin-Induced Ethylene Production. Plant Physiology 47: 504-.

26. Ma N, Tan H, Liu X, Xue J, Li Y, et al. (2006) Transcriptional regulation of ethylene receptor and CTR genes involved in ethylene-induced flower opening in cut rose (Rosa hybrida) cv. Samantha. J Exp Bot 57: 2763-2773.

27. Helariutta Y, Bishopp A, Mahonen AP (2006) Signs of change: hormone receptors that regulate plant development. Development 133: 1857-1869.

28. Theologis A, Tsuchisaka A, Yu GX, Jin HL, Alonso JM, et al. (2009) A Combinatorial Interplay Among the 1-Aminocyclopropane-1-Carboxylate

PV MLB FF. Analyzed the data: AD OR SB MLMM PV MLB FF MB. Contributed reagents/materials/analysis tools: FL. Wrote the paper: AD MB.

Isoforms Regulates Ethylene Biosynthesis in Arabidopsis thaliana. Genetics 183: 979-1003.

29. Tsuchisaka A, Theologis A (2004) Heterodimeric interactions among the 1-amino-cyclopropane-1-carboxylate synthase polypeptides encoded by the Arabidopsis gene family. Proc Natl Acad Sci U S A 101: 2275-2280.

30. Yamagami T, Tsuchisaka A, Yamada K, Haddon WF, Harden LA, et al. (2003) Biochemical diversity among the 1-amino-cyclopropane-1-carboxylate synthase isozymes encoded by the Arabidopsis gene family. J Biol Chem 278: 49102-49112.

31. Dubois A, Raymond O, Maene M, Baudino S, Langlade NB, et al. (2010) Tinkering with the C-function: a molecular frame for the selection of double flowers in cultivated roses. PLoS One 5: e9288.

32. Kitahara K, Hibino Y, Aida R, Matsumoto S (2004) Ectopic expression of the rose AGAMOUS-like MADS-box genes 'MASAKO C1 and D1' causes similar homeotic transformation of sepal and petal in Arabidopsis and sepal in Torenia. Plant Science 166: 1245-1252.

33. Kitahara K, Hirai S, Fukui H, Matsumoto S (2001) Rose MADS-box genes 'MASAKO BP and B3' homologous to class B floral identity genes. Plant Science 161: 549-557.

34. Singh AP, Tripathi SK, Nath P, Sane AP (2011) Petal abscission in rose is associated with the differential expression of two ethylene-responsive xyloglucan endotransglucosylase/hydrolase genes, RbXTH1 and RbXTH2. J Exp Bot.

35. Kaufmann K, Wellmer F, Muino JM, Ferrier T, Wuest SE, et al. (2010) Orchestration of floral initiation by APETALA1. Science 328: 85-89.

36. Wellmer F, Alves-Ferreira M, Dubois A, Riechmann JL, Meyerowitz EM (2006) Genome-wide analysis of gene expression during early Arabidopsis flower development. PLoS Genet 2: e117.

37. Roberts AV, Blake PS, Lewis R, Taylor JM, Dunstan DI (1999) The Effect of Gibberellins on Flowering in Roses. J Plant Growth Regul 18: 113-119.

38. Borner R, Kampmann G, Chandler J, Gleissner R, Wisman E, et al. (2000) A MADS domain gene involved in the transition to flowering in Arabidopsis. Plant J 24: 591-599.

39. Das P, Ito T, Wellmer F, Vernoux T, Dedieu A, et al. (2009) Floral stem cell termination involves the direct regulation of AGAMOUS by PERIANTHIA. Development 136: 1605-1611.

40. Maier AT, Stehling-Sun S, Wollmann H, Demar M, Hong RL, et al. (2009) Dual roles of the bZIP transcription factor PERIANTHIA in the control of floral architecture and homeotic gene expression. Development 136: 1613-1620.

41. Ogata J, Kanno Y, Itoh Y, Tsugawa H, Suzuki M (2005) Plant biochemistry: anthocyanin biosynthesis in roses. Nature 435: 757-758.

42. Vanneste S, Coppens F, Lee E, Donner TJ, Xie Z, et al. (2011) Developmental regulation of CYCA2s contributes to tissue-specific proliferation in Arabidopsis. EMBO J 30: 3430-3441.

43. Cheng H, Song S, Xiao L, Soo HM, Cheng Z, et al. (2009) Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLoS Genet 5: e1000440.

44. Hortensteiner S (2009) Stay-green regulates chlorophyll and chlorophyll-binding protein degradation during senescence. Trends Plant Sci 14: 155-162.

45. Nath P, Sane AP, Tripathi SK (2007) Petal abscission in rose (Rosa bourboniana var Gruss an Teplitz) is associated with the enhanced expression of an alpha expansin gene, RbEXPA1. Plant Science 172: 481-487.

46. Stead A, Rogers HJ, Roberts JA, Wagstaff C (2007) Programmed cell death during floral development and senescence. Comparative Biochemistry and Physiology a-Molecular & Integrative Physiology 146: S199-S200.

47. Kovats ES (1987) Composition of essential oils: Part 7. Bulgarian oil of rose (Rosa Damascena mill.). Journal of Chromatography A 406: 185-222.

48. Nakamura S (1987) Scent and component analysis of the Hybrid Tea Rose. Perfumer & flavorist 112: 43-45.

49. Kaminaga Y, Schneppe J, Peel G, Kish CM, Ben-Nissan G, et al. (2006) Plant phenylacetaldehyde synthase is a bifunctional homotetrameric enzyme that catalyzes phenylalanine decarboxylation and oxidation. J Biol Chem 281: 23357-23366.

50. Spiller M, Berger RG, Debener T (2010) Genetic dissection of scent metabolic profiles in diploid rose populations. Theor Appl Genet 120: 1461-1471.

51. Baayen RP, Kroes GMLW, Lange W (1998) Histology of root rot of flax seedlings (Linum usitatissimum) infected by Fusarium oxysporum f.sp. lini. European Journal of Plant Pathology 104: 725-736.

52. Irizarry RA, Ooi SL, Wu Z, Boeke JD (2003) Use of mixture models in a microarray-based screening procedure for detecting differentially represented yeast mutants. Stat Appl Genet Mol Biol 2: Article1.

53. Gentleman R, Carey V (2002) Bioconductor. RNews 2: 1116.

54. Ge Y, Dudoit S, Speed TP (2003) Resampling-based multiple testing for microarray data analysis. TEST 12: 1-44.

55. Gagnot S, Tamby JP, Martin-Magniette ML, Bitton F, Taconnat L, et al. (2008) CATdb: a public access to Arabidopsis transcriptome data from the URGV-CATMA platform. Nucleic Acids Research 36: D986-D990.

<figure><img src="image_17.png" /></figure>
56. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, et al. (2007) NCBI GEO: mining tens of millions of expression profiles–database and tools update. Nucleic Acids Res 35: D760–765.

57. Speleman F, Vandesompele J, De Preter K, Pattyn F, Poppe B, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3.

58. Coppack SW (1990) Limitations of the Pearson product-moment correlation. Clin Sci (Lond) 79: 287.

59. Liu Y, Meng Q, Chen R, Wang J, Jiang S, et al. (2004) A new method to evaluate the similarity of chromatographic fingerprints: weighted pearson product-moment correlation coefficient. J Chromatogr Sci 42: 545–550.

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