# The identification of the Rosa S -locus provides new insights into the breeding and wild origins of continuous-flowering roses

Koji Kawamura, Yoshihiro Ueda, Shogo Matsumoto, Takanori Horibe, Shungo Otagaki, Li Wang, Guoliang Wang, Laurence Hibrand-Saint Oyant, Fabrice Foucher, Marcus Linde, et al.

► To cite this version:

Koji Kawamura, Yoshihiro Ueda, Shogo Matsumoto, Takanori Horibe, Shungo Otagaki, et al. The identification of the Rosa S -locus provides new insights into the breeding and wild origins of continuous-flowering roses. Horticulture research, 2022, 9, pp.uhac155. 10.1093/hr/uhac155. hal-03810227

HAL Id: hal-03810227

https://univ-angers.hal.science/hal-03810227

Submitted on 14 Oct 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>
Article

# The identification of the Rosa S-locus provides new insights into the breeding and wild origins of continuous-flowering roses

Koji Kawamura<fnref n="1" />*, Yoshihiro Ueda<fnref n="2" />·<fnref n="3" />, Shogo Matsumoto<fnref n="4" />, Takanori Horibe<fnref n="4" />·<fnref n="5" />, Shungo Otagaki<fnref n="4" />, Li Wang<fnref n="6" />, Guoliang Wang<fnref n="7" />·<fnref n="8" />, Laurence Hibrand-Saint Oyant<fnref n="9" />, Fabrice Foucher<fnref n="9" />, Marcus Linde<fnref n="10" /> and Thomas Debener<fnref n="10" />

Abstract

This study aims to: (i) identify the Rosa S-locus controlling self-incompatibility (SI); (ii) test the genetic linkage of the S-locus with other loci controlling important ornamental traits, such as the continuous-flowering (CF) characteristic; (iii) identify the S-alleles (S<sub>c</sub>) of old Chinese CF cultivars (e.g. Old Blush, Slater's Crimson China) and examine the changes in the frequency of cultivars with S<sub>c</sub> through the history of breeding; (iv) identify wild species carrying the S<sub>c</sub>-alleles to infer wild origins of CF cultivars. We identified a new S-RNase (S<sub>C2</sub>) of Rosa chinensis in a contig from a genome database that has not been integrated into one of the seven chromosomes yet. Genetic mapping indicated that S<sub>C2</sub> is allelic to the previously-identified S-RNase (S<sub>C1</sub>) in chromosome 3. Pollination experiments with half-compatible pairs of roses confirmed that they are the pistil-determinant of SI. The segregation analysis of an F<sub>1</sub>-population indicated genetic linkage between the S-locus and the floral repressor gene KSN. The non-functional allele ksn is responsible for the CF characteristic. A total of five S-alleles (S<sub>C1-5</sub>) were identified from old CF cultivars. The frequency of cultivars with S<sub>c</sub> dramatically increased after the introgression of ksn from Chinese to European cultivars and remains high (80%) in modern cultivars, suggesting that S-genotyping is helpful for effective breeding. Wild individuals carrying S<sub>c</sub> were found in Rosa multiflora (S<sub>C1</sub>), Rosa chinensis var. spontanea (S<sub>C3</sub>), and Rosa gigantea (S<sub>C2</sub>, S<sub>C4</sub>), supporting the hypothesis of hybrid origins of CF cultivars and providing a new evidence for the involvement of Rosa multiflora.

Introduction

The rose is one of the globally most popular ornamental
plants, with a long history of breeding and cultivation. It
is important not only economically, but also culturally.
More than 30 000 rose cultivars have been developed [1]
mainly by cross breeding. Patterns of inheritance are
quite difficult to predict for most traits, as roses are pre-
dominantly outcrossing and highly heterozygous plants
[2]. As a consequence, the success of cross breeding in the
rose has largely depended on chance and the experience
of breeders, requiring enormous efforts to make new
cultivars with desirable traits [3]. Advances in scientific
knowledge on rose genetics have been much awaited.
Researchers have developed molecular markers to con-
struct genetic linkage maps of the rose, clarifying the

underlying genetic mechanisms controlling ornamen-
tal traits [4] and, genome sequencing finished in 2018
[5–7]. Key genetic factors controlling important orna-
mental traits, such as scent production [8, 9], continuous-
flowering (CF) [10], and double-flower (DF) [5, 11] have
been identified.

Although identification of the S-locus controlling self-
incompatibility (SI) in roses is essential for improving
breeding at the diploid level, it has not been fully elu-
cidated. The gametophytic SI is controlled by a single
S-locus with multiple alleles, and when one of the two
S-alleles of the pistil matches those of pollen, the pollen
is recognized as self and is rejected [12]. The pistil S
gene of Rosaceae encodes an extracellular ribonuclease
called S-RNase, which acts as a cytotoxin in self-pollen

<footnote n="1">¹Department of Environmental Engineering, Osaka Institute of Technology, Japan</footnote>

<footnote n="2">²Gifu International Academy of Horticulture, Japan</footnote>

<footnote n="3">³Gifu World Rose Garden, Japan</footnote>

<footnote n="4">⁴Graduate School of Bioagricultural Sciences, Nagoya University, Japan</footnote>

<footnote n="5">⁵College of Bioscience and Biotechnology, Chubu University, Japan</footnote>

<footnote n="6">⁶College of Life Sciences, Sichuan University, China</footnote>

<footnote n="7">⁷Jiangsu Provincial Department of Agriculture and Rural Affairs, China</footnote>

<footnote n="8">⁸Agricultural University of Nanjing, China.</footnote>

<footnote n="9">⁹Univ Angers, INRAE, Institut Agro, IRHS, SFR QUASAV, F-49000 Angers, France</footnote>

<footnote n="10">¹⁰Leibniz Universität, Hannover, Germany</footnote>

<footnote>*Corresponding author. E-mail: koji.kawamura@oit.ac.jp</footnote>
tubes. The pollen S gene is involved in detoxifying non-self-S-RNases and encodes an *F-box* gene [12]. The “col-laborative non-self-recognition” model was proposed for SI in the Solanaceous plant, *Petunia* [13], where multiple *F-box* genes are involved in pollen specificity, and each targets a subset of non-self-S-RNases for detoxification. The SI of Rosaceae may adopt this “non-self-recognition by a multiple factors” system [12], except for the self-recognition by a single *F-box* gene system in *Prunus* [14]. In the genome databases of *R. chinensis* “Old Blush”, two different regions have been proposed as the S-locus [5, 15], both of which are located on the same chromosome 3 within a few Mbp of each other. Vieira *et al.* [15] concluded that their candidate region was the true S-locus because its S-RNase (hereafter called as $S_{C1}$ S-RNase) has a stronger similarity to the *Prunus* S-RNase than the S-RNase36 identified by Hibrand-Saint Oyant *et al.* [5]. Chen *et al.* [16] identified an S-locus-like region in the genome of *Rosa rugosa*, whose the S-RNase is orthologous to the $S_{C1}$ S-RNase of Old Blush. Du *et al.* [17] identified the S-RNase controlling SI in *Fragaria*, which is also more similar to the $S_{C1}$ S-RNase rather than S-RNase36. These results indicate that the $S_{C1}$ S-RNase is the true S gene controlling SI in Old Blush. However, validation of the *Rosa* S-locus by pollination experiments has been conducted with only a few individuals [15]. Pollination tests with many individuals are required to confirm the result. Furthermore, Vieira *et al.* [15] identified the $S_{C1}$ S-RNase of Old Blush from the genome database of Raymond *et al.* [6] but failed to identify its allelic gene in the other genome database of Hibrand-Saint Oyant *et al.* [5]. The S-RNase36 identified by Hibrand-Saint Oyant *et al.* [5] has a single intron and weak homology to the $S_{C1}$ S-RNase, which has two introns. As the two genome databases of Old Blush represent different haplotypes of chromosome 3 [18], there should be unidentified S-RNase in the genome database of Hibrand-Saint Oyant *et al.* [5]. Thus, the first objective of this study is to (i) identify the other allele of S-RNases in Old Blush and to confirm the *Rosa* S-locus with strong evidence from a number of pollination experiments.

The second objective of this study is to (ii) test the genetic linkages of the S-locus with other loci controlling important ornamental traits. In the chromosome 3 where the *Rosa* S-locus is assumed to be located, there are other important loci controlling valuable traits for ornamental plants, such as CF [19], double flower [20], thornlessness [21], and resistance against black spot disease [22]. These previous studies found strong skewness in the segregation of these traits and hypothesized that their underlying genes are genetically linked with the S-locus. The genetic linkage of the S-locus and other loci may have important consequences for rose breeding. For example, the CF characteristic is controlled by a single recessive locus on the chromosome 3, and the underlying gene is the non-functional allele (ksn) of RoKSN [5, 10], a floral repressor gene in roses [23]. Roses homozygous at the non-functional allele ksn will be targets in the cross

breeding of CF roses. But, if the ksn alleles of parental roses are strongly linked with the same S-allele, breeding of CF roses with the ksn-homozygote will be hampered by SI, where the pistil rejects the pollen with the same S-allele linked to ksn.

The third objective of this study is to (iii) determine the S-allele ($S_C$) of old Chinese CF cultivars and to clarify the changes in the frequency of rose cultivars with $S_C$ during the history of rose breeding. The ksn allele conferring the CF characteristic in modern rose cultivars originated from old Chinese cultivars around 200 years ago [10, 24]. The original Chinese cultivars introducing the CF characteristic into modern roses are thought to be the four China roses, i.e. Slater's Crimson China, Parsons' Pink China (Old Blush), Hume's Blush Tea-scented China, and Parks' Yellow Tea-scented China [25]. After the introduction of ksn from China to Europe, the frequency of rose cultivars carrying the Sc should increase but may progressively decrease during the 200 years' history of rose breeding due to the hybridizations with European cultivars. However, if Sc is genetically linked with ksn, the strong artificial selection on ksn [24] might result in a ksn-associated increase in the frequency of Sc during the history of rose breeding.

The fourth objective of this study is to (iv) identify candidate genotypes for wild ancestors of old Chinese CF cultivars by using $S_C$ as markers. The introduction of Chinese CF cultivars is one of the most revolutionary events throughout the history of rose breeding [24, 25], and the Chinese cultivars have profound effects on the genetic bases of modern rose cultivars [26]. However, the wild ancestral origin of the old Chinese CF cultivars have not been completely elucidated. Molecular phylogenetic studies have indicated their complex hybrid origins [27–32]. The whole genome sequencing will play a critical role in the elucidation of the hybrid origin of old Chinese cultivars but is too complicated to analyze the massive data for many candidate species. There is also a large intraspecific genetic variation within the species [33, 34], and there might be also natural hybrids between the species. Due to the high genetic divergences of the plant S-locus [35], we expect that the S-allele specific marker is useful to pinpoint the candidate wild ancestors of old cultivars.

# Results

## Genome-wide identification of candidate S-RNase

We performed a genome-wide search for candidate S-RNase genes in the genome databases of Old Blush [5, 6], *R. multiflora* [7], and *R. rugosa* [16, 39] (Supplementary data Table D2) and identified 16 genes located on four chromosomes (Fig. 1a). Of these, seven genes were expressed in the pistil according to the pistil transcriptomes of Old Blush, *R. multiflora*, and *R. rugosa* (Fig. 1b). These candidate genes were named by chromosome number and letter, e.g. 3D and 6A (Fig. 1a). The specific genes of Old Blush, *R. multiflora*, and *R. rugosa* were named with Rc, Rm,
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Genome-wide identification and genetic analyses of candidate S-RNase genes in the rose. (a) Distribution of the candidate genes in R. chinensis "Old Blush" genome. Numbers and genomic positions are shown based on the genome database of Raymond et al. [6]. Pistil-expressed genes are shown in red. (b) Pistil and stamen expression levels of the candidate genes in roses. FPKM values were calculated from RNA-seq data of the pistil and stamen. Averages and standard errors of three individuals are shown for R. chinensis "Old Blush". From the transcriptome data of eight R. multiflora plants, genes orthologous to the reference genes of Old Blush were identified, and their averages and standard errors of the FPKM values are shown for R. multiflora. (c) Molecular phylogenetic tree of the pistil-expressed, candidate genes, including Prunus S-RNases as references. The tree was rooted with T2-RNase of Arabidopsis thaliana RNS2 (NM129536) [36]. Bold branches have FastTree support values [37, 38] above 0.95 based on 1000 resamples.</figcaption></figure>

and Rg, respectively (Fig. 1c). The putative amino acid sequences were obtained from the cDNA sequences, and a molecular phylogenetic tree of the pistil-expressed, S-RNase-like genes was constructed, including the Prunus S-RNase (Fig. 1c). The results show that the 3D and 0A genes are the primary candidates for S-RNase controlling SI in the rose due to their (i) high expression levels in the pistil (Fig. 1b), (ii) similarity to Prunus S-RNase (Fig. 1c), and (iii) high genetic divergences (Fig. 1c).

The 3D S-RNase gene in Old Blush named as "Rc3D_Sc1" (Fig. 1c) is the same gene previously identified by Vieira et al. [15] as the candidate S-RNase in roses. We

identified nine additional 3D S-RNase-like genes from the genome databases and pistil transcriptomes of R. multiflora and R. rugosa (Fig. 1c). One of them, named as "Rg3D_S18", is the same gene previously reported as S-RNase (Chr4.718) in the S-locus of R. rugosa [16]. The high genetic divergence (64.8%) in the 3D S-RNase gene agrees well with the assumption that this gene encodes the true S-RNase controlling SI in roses (Table 1). One of the 3D genes of R. multiflora, named as "Rm3D_Sc1", is identical to the Rc3D_Sc1 of Old Blush, indicating that the R. multiflora and Old Blush share the same S-alleles (Sc1).
<table><caption>**Table 1.** The identification of candidate genes for the Rosa S-RNase. According to the phylogenetic tree (Fig. 1c), four genes (3D, 0A, 3A, 6C) with high similarities to Prunus S-RNases are considered to be candidate S-RNase genes</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Data source<sup>b</sup></th><th colspan="4">Candidate S-RNase genes<sup>a</sup></th></tr><tr><th>3D</th><th>0A</th><th>3A</th><th>6C</th></tr></thead><tbody><tr><td rowspan="2">R. chinensis "Old Blush"</td><td>Genome database [5]</td><td>×</td><td>Rc0A_S<sub>C2</sub></td><td>Rc3A</td><td>Rc6C_a</td></tr><tr><td>Genome database [6]</td><td>Rc3D_S<sub>C1</sub></td><td>×</td><td>Rc3A</td><td>Rc6C_a</td></tr><tr><td rowspan="3">R. rugosa</td><td>Genome database [16]</td><td>Rg3D_S<sub>18</sub></td><td>×</td><td>×</td><td>Rg6C_GD1</td></tr><tr><td>Genome database [39]</td><td>Rg3D_S<sub>20</sub></td><td>×</td><td>×</td><td>×</td></tr><tr><td>PT of the plant Rg46</td><td>Rg3D_S<sub>18</sub></td><td>Rg0A_S<sub>19</sub></td><td>×</td><td>×</td></tr><tr><td rowspan="12">R. multiflora</td><td>Genome database [7]</td><td>Rm3D_S<sub>15</sub></td><td>×</td><td>Rm3A</td><td>Rm6C</td></tr><tr><td>PT of the plant Rm08</td><td>Rm3D_S<sub>16</sub></td><td>Rm0A_S<sub>6</sub></td><td>Rm3A_Rm08</td><td>×</td></tr><tr><td>PT of the plant Rm09</td><td>Rm3D_S<sub>7</sub></td><td>Rm0A_S<sub>8</sub></td><td>Rm3A_Rm09</td><td>Rm6C_Rm09</td></tr><tr><td>PT of the plant Rm27</td><td>Rm3D_S<sub>9</sub></td><td>×</td><td>Rm3A_Rm27</td><td>Rm6C_Rm27</td></tr><tr><td>PT of the plant Rm28</td><td>Rm3D_S<sub>13</sub></td><td>Rm0A_S<sub>10</sub></td><td>×</td><td>Rm6C_Rm28</td></tr><tr><td>PT of the plant Rm33</td><td>Rm3D_S<sub>9</sub></td><td>Rm0A_S<sub>12</sub></td><td>×</td><td>Rm6C_Rm33</td></tr><tr><td>PT of the plant Rm3<br>(tetraploid)</td><td>Rm3D_S<sub>C1</sub></td><td>×</td><td>×</td><td>×</td></tr><tr><td>PT of the plant Rm2</td><td>Rm3D_S<sub>15</sub></td><td>×</td><td>×</td><td>×</td></tr><tr><td>PT of the plant Rm1</td><td>Rm3D_S<sub>C1</sub></td><td>×</td><td>×</td><td>×</td></tr><tr><td>PT of the plant Rm1</td><td>Rm3D_S<sub>14</sub></td><td>×</td><td>×</td><td>×</td></tr><tr><td colspan="2">Average pairwise identity of protein sequences</td><td>64.8%</td><td>62.9%</td><td>97.8%</td><td>94.9%</td></tr><tr><td colspan="2">Expression levels in pistil<sup>c</sup></td><td>***</td><td>***</td><td>*</td><td>**</td></tr></tbody></table>

In the same clade of the Prunus S-RNase and the Rosa 3D S-RNase, there were three other genes (6C, 3A, 0A; Fig. 1c). The 6C and 3A genes are located in different genomic positions from the 3D S-RNase (Fig. 1a). Therefore, they are not allelic genes to 3D genes. The 6C and 3A genes have low or undetectable expression levels in pistils (Fig. 1b) and low genetic divergences between alleles (>94% identity; Table 1). One allele of the 6C gene of Old Blush [5] and that of R. rugosa [39] are pseudogenes due to frameshift mutations (Supplementary Data Table D2). The 3A gene is homozygous in Old Blush and cannot be identified in two R. rugosa genomes. These results indicate that the 6C and 3A genes are not functional S-RNase genes.

The 0A gene is supposed to be the allelic gene to 3D S-RNase in roses. The 0A gene of Old Blush, named as "Rc0A_S<sub>C2</sub>", was identified in the chromosome 0 (i.e. a contig not assigned to one of the seven chromosomes) of the genome database of Hibrand-Saint Oyant et al. [5]. We identified five 0A S-RNase-like genes from the genome databases and pistil transcriptomes of R. multiflora and R. rugosa (Fig. 1c). These 0A S-RNase-like genes are strongly expressed in pistils (Fig. 1b) and exhibit high genetic divergence (62.9%; Table 1), like the 3D S-RNase gene. According to the qPCR analyses of different tissues and developmental stages (Supplementary information 2), the 0A and 3D S-RNase genes show pistil-specific expressions. The analyses of pistil transcriptomes of seven diploid R. multiflora plants showed that four plants (Rm08, Rm09, Rm28, Rm33) expressed one 3D gene and one 0A

gene, while the other three plants (Rm27, Rm2, Rm1) expressed no 0A genes but two 3D genes (Table 1). This also suggests that the 0A and 3D genes are allelic.

To test the hypothesis that the 0A S-RNase is allelic to the 3D S-RNase, we genetically mapped these genes using a diploid mapping population. The female parent of the diploid mapping population (FW), "The Fairy" (TF), has two 0A genes, one of which is identical to the Rc0A_S<sub>C2</sub> (hereafter called as S<sub>C2</sub> for simplicity) of Old Blush and the other called S<sub>21</sub> (Supplementary information 4). The TF has an S<sub>C2</sub>/S<sub>21</sub> genotype of the 0A S-RNase. In the male parent, RW, we identified one 3D S-RNase-like gene and named it S<sub>1w</sub> (Supplementary information 4); it has a sequence similar to the Rc3D_S<sub>C1</sub> of Old Blush (hereafter called as S<sub>C1</sub> for simplicity). The RW has an S<sub>1w</sub>/S<sub>x</sub> genotype of the 3D S-RNase (S<sub>x</sub> is an anonymous allele). The 0A and 3D S-RNase genes are mapped to homologous positions of chromosome 3 of the female (TF3) and male (RW3) genetic linkage maps, respectively (Fig. 2a). Estimation of genomic position of the S gene from map cM of the integrated map (IT3) indicates that the S gene is located on 41.5 Mbp of the chromosome 3 (RC3) in the genome database of Hibrand-Saint Oyant et al. [5] (Fig. 2b) and is located on 5.6 Mbp of the chromosome 3 (Chr3) of the other genome database of Raymond et al. [6] (Fig. 2c). In the latter genome database, the S<sub>C1</sub> S-RNase is located at 5.5 Mbp of the Chr3 close to the estimated position (Fig. 2c).

The two haploid genome databases of R. chinensis "Old Blush" are inversely oriented (Fig. 3a). The estimated

<footnote><sup>a</sup>Genes identified in genome databases and pistil transcriptomes. × = unidentified.</footnote>

<footnote><sup>b</sup>PT = Pistil transcriptome (RNaseq).</footnote>

<footnote><sup>c</sup>The average FPKM values of R. multiflora (Fig. 1b): ***, >100; **, 50–100; *10–50.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Genetic mapping of candidate S-RNase genes identified from the genome database of R. chinensis "Old Blush". (a) Genetic linkage maps of a diploid population FW [19]: The chromosome 3 of the female parent map TF3, the male parent map RW3, and their integrated map IT3. (b)-(c) Estimations of genomic positions of the S gene from the map cM of IT3. By using regression lines made by surrounding markers (green), the genomic positions of the S gene (red) were estimated in the genome database of (b) Hibrand-Saint Oyant et al. [5] and (c) Raymond et al. [6].</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Synteny of the two genome databases of R. chinensis "Old Blush" chromosome 3 and the estimated position of S-RNase. Dot plots of orthologous genes of two genome databases: Chr3 [6] and RC3 [5] are made using SynMap of CoGe [40]. The CDS of the two genome databases were used to generate the dot plots. (a) Whole chromosome scale, and (b) the focus of the region marked by a blue rectangle. The S<sub>C1</sub> (3D) S-RNase is located on the 5.5 Mbp position of Chr3, and there is a 500 kbp gap in the synteny at 41.4–41.9 Mbp of RC3. The S<sub>C2</sub> (0A) S-RNase is estimated to be located on this gap region.</figcaption></figure>

position of the S gene in RC3 (41.5 Mbp) is in an orthologous position to the S<sub>C1</sub> S-RNase position in Chr3 (5.5 Mbp) and is located in a gap of synteny between the two genome databases (Fig. 3b).

## The identification of the S-locus F-box genes

There should be pollen-expressed F-box genes flanking the S-RNase. This prediction was confirmed by the genome database search and transcriptomes of the stamens of Old Blush (Table 2). From 5.3 to 5.8 Mbp region in Chr3, we identified 12 F-box genes and one S-RNase (S<sub>C1</sub>). Homologous F-box genes were identified in 41.3–41.9 Mbp region in RC3, and three other F-box genes were also identified in the contig RC0 containing the S<sub>C2</sub> S-RNase (Table 2). All F-box genes were expressed in the stamen but not in the pistil. According to the

nomenclature of Kubo et al. [13], these F-box genes were termed as SLF or FBX (Supplementary Table S7-2). We also confirmed that the similar number of SLFs were identified in 500–600 kbp regions flanking to the 3D S-RNase-like genes of R. multiflora (S<sub>15</sub> S-RNase, S<sub>16</sub> S-RNase) and R. rugosa (S<sub>18</sub> S-RNase, S<sub>20</sub> S-RNase) genome databases, and they were expressed in the stamens but not in the pistils (Table S7-1a,b).

Figure 4 shows the structures of the putative S-locus regions identified in the Rosa genome databases. In the S-locus of Old Blush chromosome RC3, there is no S-RNase but there is a poly-N region, and its adjacent 10kbp regions (Block-a, Block-b) show close sequence identities with the end regions of the contig RC0. Block-a is 11220 bp with a 98.8% identity between the RC3 and RC0. Block-b is 11382 bp with a 97.5% identity. The S<sub>C2</sub>
<table><caption>**Table 2.** List of candidate S-locus F-box genes identified by in silico blast search on R. chinensis "Old Blush" genome databases</caption><thead><tr><th rowspan="2">Type<sup>c</sup></th><th rowspan="2">Gene name</th><th rowspan="2">S-genotype</th><th colspan="3">Genomic position<sup>a</sup></th><th rowspan="2">Length</th><th rowspan="2">Annotation<sup>d</sup></th><th colspan="2">FPKM<sup>b</sup></th></tr><tr><th>Chromosome</th><th>Start</th><th>End</th><th>Stamen</th><th>Pistil</th></tr></thead><tbody><tr><td>FBX12</td><td>S<sub>C1</sub>_FBX12_pseudo</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 391 421</td><td>5 390 102</td><td>1320</td><td>NA</td><td>Truncated</td><td></td></tr><tr><td>SLF1</td><td>S<sub>C1</sub>_SLF1</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 434 392</td><td>5 433 145</td><td>1248</td><td>Chr3g0455861</td><td>22 (1.6)</td><td>1 (0.4)</td></tr><tr><td>SLF2</td><td>S<sub>C1</sub>_SLF2</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 475 946</td><td>5 477 154</td><td>1209</td><td>Chr3g0455891m</td><td>18 (2.9)</td><td>1 (0.3)</td></tr><tr><td><strong>S-RNase</strong></td><td><strong>S<sub>C1</sub>_SRNase</strong></td><td><strong>S<sub>C1</sub></strong></td><td><strong>Chr3</strong></td><td><strong>5 498 387</strong></td><td><strong>5 488 716</strong></td><td><strong>9672</strong></td><td><strong>Chr3g0455911m</strong></td><td><strong>0 (0)</strong></td><td><strong>100 (32)</strong></td></tr><tr><td>SLF3</td><td>S<sub>C1</sub>_SLF3</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 559 472</td><td>5 558 219</td><td>1254</td><td>Chr3g0455931m</td><td>6 (1.9)</td><td>1 (0.1)</td></tr><tr><td>SLF4</td><td>S<sub>C1</sub>_SLF4</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 618 754</td><td>5 619 992</td><td>1239</td><td>Chr3g0455991</td><td>18 (2.8)</td><td>1 (0.2)</td></tr><tr><td>FBX1</td><td>S<sub>C1</sub>_FBX1</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 722 023</td><td>5 720 695</td><td>1329</td><td>Chr3g0456091</td><td>8 (0.8)</td><td>0 (0.2)</td></tr><tr><td>SLF5</td><td>S<sub>C1</sub>_SLF5</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 735 463</td><td>5 734 183</td><td>1281</td><td>Chr3g0456131</td><td>13 (2.5)</td><td>1 (0.2)</td></tr><tr><td>SLF6</td><td>S<sub>C1</sub>_SLF6</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 751 338</td><td>5 752 606</td><td>1269</td><td>Chr3g0456171</td><td>11 (0.5)</td><td>1 (0.2)</td></tr><tr><td>SLF8</td><td>S<sub>C1</sub>_SLF8</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 803 902</td><td>5 805 158</td><td>1257</td><td>Chr3g0456241</td><td>8 (2.4)</td><td>1 (0.3)</td></tr><tr><td>SLF9</td><td>S<sub>C1</sub>_SLF9</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 809 247</td><td>5 808 021</td><td>1227</td><td>Chr3g0456251</td><td>20 (2.6)</td><td>1 (0.4)</td></tr><tr><td>SLF10</td><td>S<sub>C1</sub>_SLF10</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 857 662</td><td>5 856 391</td><td>1272</td><td>Chr3g0456291</td><td>7 (0.6)</td><td>0 (0.3)</td></tr><tr><td>SLF11</td><td>S<sub>C1</sub>_SLF11</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 880 743</td><td>5 879 523</td><td>1221</td><td>Chr3g0456311</td><td>5 (2)</td><td>1 (0.2)</td></tr><tr><td>FBX4</td><td>S<sub>C1</sub>_FBX4_pseudo</td><td>S<sub>C1</sub></td><td>Chr3</td><td>5 890 721</td><td>5 889 444</td><td>1278</td><td>NA</td><td>Truncated</td><td></td></tr><tr><td>FBX2</td><td>S<sub>C2</sub>_FBX2</td><td>S<sub>C2</sub></td><td>RC0</td><td>27 704 970</td><td>27 703 705</td><td>1266</td><td>NA</td><td>5 (0.6)</td><td>1 (0.2)</td></tr><tr><td><strong>S-RNase</strong></td><td><strong>S<sub>C2</sub>_SRNase</strong></td><td><strong>S<sub>C2</sub></strong></td><td><strong>RC0</strong></td><td><strong>27 779 508</strong></td><td><strong>27 755 726</strong></td><td><strong>23 783</strong></td><td><strong>NA</strong></td><td><strong>0 (0)</strong></td><td><strong>44 (44)</strong></td></tr><tr><td>FBX3</td><td>S<sub>C2</sub>_FBX3</td><td>S<sub>C2</sub></td><td>RC0</td><td>27 852 918</td><td>27 854 177</td><td>1260</td><td>RC0G0207600</td><td>17 (1)</td><td>1 (0.4)</td></tr><tr><td>SLF7</td><td>S<sub>C2</sub>_SLF7</td><td>S<sub>C2</sub></td><td>RC0</td><td>27 963 602</td><td>27 962 373</td><td>1230</td><td>RC0G0208100m</td><td>4 (0.7)</td><td>0 (0.1)</td></tr><tr><td>FBX4</td><td>S<sub>C2</sub>_FBX4</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 290 222</td><td>41 291 499</td><td>1278</td><td>RC3G0342300</td><td>16 (2.4)</td><td>1 (0.5)</td></tr><tr><td>SLF11</td><td>S<sub>C2</sub>_SLF11</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 310 056</td><td>41 311 276</td><td>1221</td><td>RC3G0342500</td><td>5 (0.9)</td><td>0 (0.1)</td></tr><tr><td>SLF10</td><td>S<sub>C2</sub>_SLF10</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 344 314</td><td>41 345 582</td><td>1269</td><td>RC3G0343000</td><td>5 (0.5)</td><td>0 (0.2)</td></tr><tr><td>SLF9</td><td>S<sub>C2</sub>_SLF9</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 348 095</td><td>41 349 321</td><td>1227</td><td>RC3G0343100</td><td>18 (3.1)</td><td>1 (0.5)</td></tr><tr><td>SLF8</td><td>S<sub>C2</sub>_SLF8</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 352 234</td><td>41 350 978</td><td>1257</td><td>RC3G0343200</td><td>9 (2.6)</td><td>1 (0.2)</td></tr><tr><td>SLF6</td><td>S<sub>C2</sub>_SLF6_pseudo</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 369 680</td><td>41 368 410</td><td>1271</td><td>RC3G0343500m</td><td>Truncated</td><td></td></tr><tr><td>SLF5</td><td>S<sub>C2</sub>_SLF5</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 380 179</td><td>41 381 417</td><td>1239</td><td>RC3G0343800m</td><td>16 (2.8)</td><td>1 (0.3)</td></tr><tr><td>SLF4</td><td>S<sub>C2</sub>_SLF4_1</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 646 962</td><td>41 648 194</td><td>1233</td><td>RC3G0344900m</td><td>4 (0.5)</td><td>0 (0)</td></tr><tr><td>S<sub>C2</sub>_SLF4_2</td><td>S<sub>C2</sub></td><td>S<sub>C2</sub></td><td>RC3</td><td>41 695 111</td><td>41 696 343</td><td>1233</td><td>RC3G0345200</td><td>4 (1.7)</td><td>0 (0.1)</td></tr><tr><td>SLF3</td><td>S<sub>C2</sub>_SLF3</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 876 741</td><td>41 877 979</td><td>1239</td><td>RC3G0346100</td><td>10 (2.8)</td><td>1 (0.2)</td></tr><tr><td>SLF2</td><td>S<sub>C2</sub>_SLF2</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 896 100</td><td>41 897 320</td><td>1221</td><td>RC3G0346200m</td><td>13 (2.2)</td><td>1 (0.2)</td></tr><tr><td>SLF1</td><td>S<sub>C2</sub>_SLF1</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 907 657</td><td>41 908 904</td><td>1248</td><td>RC3G0346300</td><td>15 (3.1)</td><td>1 (0.2)</td></tr><tr><td>FBX5</td><td>S<sub>C2</sub>_FBX5</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 932 562</td><td>41 931 243</td><td>1320</td><td>RC3G0346700</td><td>4 (1.6)</td><td>1 (0.2)</td></tr><tr><td>FBX12</td><td>S<sub>C2</sub>_FBX12</td><td>S<sub>C2</sub></td><td>RC3</td><td>41 965 102</td><td>41 966 437</td><td>1336</td><td>RC3G0347000m</td><td>Truncated</td><td></td></tr></tbody></table>

S-RNase in the contig RC0 appears to be integrated into the poly-N region of RC3. Furthermore, SLF7 is present in the contig RC0, and its orthologous genes were identified in the S-locus regions of R. rugosa and R. multiflora (See also Fig.S7-2 for the alignment). Further, the co-segregations of SLF5 in RC3 and S<sub>C2</sub> S-RNase in RC0 were tested using the FW mapping population (Supplementary information 8) and found to be perfectly co-segregated (n = 97).

## Validation of the S-RNase based S-genotyping by pollination experiments

In order to validate these S-RNase genes as the pistil-determinant of SI, pollination experiments on the pairs of diploid roses that share the same S-alleles were performed (Table 3). For example, plants with the S-genotype S<sub>C1</sub>/S<sub>x</sub> (where S<sub>x</sub> is an anonymous genotype that includes several different S-RNase genotypes) were selected and pollinated by pollen collected from the plants with S<sub>C1</sub>/S<sub>y</sub> (where S<sub>y</sub> is another anonymous genotype). After seed

maturation, the S-genotypes of the seeds were determined, and the genotypic frequency of pollen fertilizing the ovule were calculated. A total 111 seeds were analyzed, and results show that no S<sub>C1</sub> pollen fertilized the ovules (i.e. all seeds are derived from pollen with S<sub>y</sub> genotypes). These pollination tests with half-compatible pairs were performed for R. multiflora plants sharing S<sub>C1</sub> S-RNase-like genes, i.e. Rm3D genes, S<sub>7</sub>, S<sub>9</sub>, S<sub>11</sub>, or S<sub>13</sub> S-RNase (Fig. 1c), providing the same results with S<sub>C1</sub> by analyzing total 174 seeds (Table 3). Furthermore, the pollination experiment of TF (S<sub>C2</sub>/S<sub>21</sub>) with the pollen of Old Blush (S<sub>C1</sub>/S<sub>C2</sub>) supported the prediction that the S<sub>C2</sub> pollen of Old Blush was rejected by the pistil of TF by analyzing 60 seeds (Table 3). The cross pollinations with R. multiflora plants sharing S<sub>C2</sub> S-RNase-like genes, i.e. Rm0A genes, S<sub>6</sub>, S<sub>10</sub>, S<sub>12</sub> S-RNase (Fig. 1c), providing the consistently same results for total 124 seeds (Table 3). These data strongly support the hypothesis that the 3D and 0A S-RNase genes are the pistil-determinant of SI in roses.

<footnote><sup>a</sup>Genome data source is Raymond et al. [6] for S<sub>C1</sub>, Hibrand-Saint Oyant et al. [5] for S<sub>C2</sub>.</footnote>

<footnote><sup>b</sup>FPKM (Fragments per kilobase of exon per million reads mapped) was calculated from the RNA-seq data of three individuals of R. chinensis "Old Blush" (OB15, OB20, OB75), and averages with standard errors in parenthesis are shown.</footnote>

<footnote><sup>c</sup>Type was defined as either SLF (S-locus linked F-box) or FBX, with numbers indicating groups with homologous protein sequences. We categorized the type as SLF when at least four S-haplotypes out of six (S<sub>C1</sub>, S<sub>C2</sub>, S<sub>15</sub>, S<sub>16</sub>, S<sub>18</sub>, S<sub>20</sub>) had the F-box genes, whereas we called them FBX only when less than three S-haplotypes had the F-box genes (Table S7-2).</footnote>

<footnote><sup>d</sup>NA = No Annotation in original database. The suffix "m" indicates "modified" annotation of original one.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 4.** The structure of the Rosa S-locus. The putative S-locus genomic regions were extracted from genome databases, and the positions of the S-RNase and F-box (SLF, FBX) genes were plotted. Homologous F-box genes are connected by dotted lines. Data are available in Table 1 for R. chinensis "Old Blush" and Supplementary Table S7-1a,b for R. rugosa and R. multiflora.</figcaption></figure>

## Linkage between the S-locus and important ornamental traits

To consider the effect of SI on the breeding of roses, the degree of genetic linkage between the S-locus and the genes underlying important ornamental characteristics (CF and DF) was tested. The KSN gene (controlling CF) has a 13.5Mbp distance from the S-locus, and the AP2-like gene (controlling DF) has a 9.0 Mbp distance from the S-locus (Table 4). By using two diploid F₁ hybrid populations [19, 41], the recombinant frequencies between the S-locus and these genes were estimated. The recombinant frequency between S and KSN is 20% in RW, and those between S and AP2-like is 13% in RW and 40% in TF in the FW population. The recombinant frequency between S and AP2-like is 22% in 93/1–119 in the 94/1 population. Except for the high recombination between the S and AP2-like loci in TF, there are significant genetic linkages between the S and the KSN loci and the S and the AP2-like loci. This may act as an internal constraint on rose breeding, which is discussed later.

## Sc-alleles of old Chinese CF cultivars

To identify the Chinese S-alleles ($S_C$) introduced into Europe in the 18th century and examine their fate during

the past 200 years' history of rose breeding, we first analysed the S-genotypes of old Chinese CF cultivars (Table 5). The S-genotyping showed that the old Chinese cultivars frequently shared $S_{C1}$ and $S_{C2}$ alleles with Old Blush (Table 5). Hume's Blush Tea-scented China has the same S-genotype ($S_{C1}/S_{C2}$) as Old Blush. The other CF cultivars, such as Slater's Crimson China, R. chinensis "Mutabilis", and R. chinensis, have either $S_{C1}$ or $S_{C2}$, indicating that they have unidentified S-alleles. These unidentified $S_C$-alleles were amplified with RT-PCR using mRNA prepared from their pistils with degenerate primers for S-RNase in roses and sequenced (Fig. S10-1). Partial ($S_{C3}$, $S_{C4}$) and full ($S_{C5}$) S-RNase cDNA sequences were obtained and analysed via BLASTX searches against a local protein database of the Rosa S-RNase and S-RNase-like proteins, shown in Figure 1c. The BLAST search confirmed that the $S_{C3-5}$ sequences are closest to the 0A S-RNase proteins (See Fig. S10-2 for the alignment). As a consequence, S-genotypes were estimated as follows: Mutabilis = $S_{C1}/S_{C5}$, R. chinensis = $S_{C1}/S_{C4}$, Slater's Crimson China = $S_{C1}/S_{C3}$ or $S_{C1}/S_{C2}/S_{C3}$ (triploid type), and Sanguinea = $S_{C2}/S_{C3}$ (Table 5).

The genotyping of the KSN gene responsible for the CF characteristic identified five KSN genotypes (Table 5). Three KSN alleles, i.e. KSNW (wild allele), ksnCopia (the
<table><caption>**Table 3.** Validation of the S-RNase-based genotyping by pollination experiments with half-compatible pairs. See Supplementary information 5 for full details of the pollination experiments</caption><thead><tr><th rowspan="2">Recipient plant (e)</th><th rowspan="2">S-genotype<sup>b</sup></th><th colspan="2">Pollen donor (e)</th><th colspan="2">Number of seeds (n) by S-genotype</th><th rowspan="2">S-genotype</th><th rowspan="2">n</th><th rowspan="2">Number of Pollen Genotypes fertilizing the ovule</th></tr><tr><th>Plant ID<sup>a</sup></th><th>S-genotype<sup>b</sup></th><th>S-genotype</th><th>n</th></tr></thead><tbody><tr><td rowspan="2">M</td><td>S<sub>x</sub></td><td>S<sub>C1</sub></td><td>M</td><td>S<sub>C1</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>C1</sub></td><td>0</td><td>S<sub>C1</sub> / S<sub>C1</sub></td><td>44</td></tr><tr><td>S<sub>8</sub></td><td>S<sub>7</sub></td><td>Rm08</td><td>S<sub>7</sub></td><td>S<sub>6</sub></td><td>S<sub>8</sub> / S<sub>7</sub></td><td>0</td><td>S<sub>7</sub> / S<sub>6</sub></td><td>16</td></tr><tr><td rowspan="2">M</td><td>S<sub>x</sub></td><td>S<sub>9</sub></td><td>M</td><td>S<sub>9</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>9</sub></td><td>0</td><td>S<sub>9</sub> / S<sub>9</sub></td><td>24</td></tr><tr><td>S<sub>x</sub></td><td>S<sub>11</sub></td><td>M</td><td>S<sub>11</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>11</sub></td><td>0</td><td>S<sub>11</sub> / S<sub>11</sub></td><td>21</td></tr><tr><td rowspan="2">M</td><td>S<sub>x</sub></td><td>S<sub>13</sub></td><td>M</td><td>S<sub>13</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>13</sub></td><td>0</td><td>S<sub>13</sub> / S<sub>13</sub></td><td>27</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">TF</td><td>S<sub>21</sub></td><td>S<sub>C2</sub></td><td>OB</td><td>S<sub>C2</sub></td><td>S<sub>C1</sub></td><td>S<sub>21</sub> / S<sub>C2</sub></td><td>0</td><td>S<sub>C2</sub> / S<sub>C2</sub></td><td>25</td></tr><tr><td>S<sub>x</sub></td><td>S<sub>6</sub></td><td>M</td><td>S<sub>6</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>6</sub></td><td>0</td><td>S<sub>6</sub> / S<sub>y</sub></td><td>23</td></tr><tr><td rowspan="2">M</td><td>S<sub>x</sub></td><td>S<sub>10</sub></td><td>M</td><td>S<sub>10</sub></td><td>S<sub>y</sub></td><td>S<sub>x</sub> / S<sub>10</sub></td><td>0</td><td>S<sub>10</sub> / S<sub>10</sub></td><td>17</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Rm32</td><td>S<sub>13</sub></td><td>S<sub>12</sub></td><td>Rm33</td><td>S<sub>12</sub></td><td>S<sub>y</sub></td><td>S<sub>13</sub> / S<sub>12</sub></td><td>0</td><td>S<sub>12</sub> / S<sub>9</sub></td><td>15</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

copia-retrotransposon inserted, non-functional allele) [10], and ksn<sup>null</sup> (the deletion, non-functional allele) [5] were distinguished by PCR. New primers were designed to amplify the ksn<sup>null</sup> allele (Fig.S10-3). Old Blush, Hume's Blush Tea-scented China, and Mutabilis are heterozygous for the two non-functional alleles (ksn<sup>copia</sup>/ksn<sup>null</sup>), whereas Slater's Crimson China and R. chinensis are homozygous for one non-functional allele (ksn<sup>copia</sup>/ksn<sup>null</sup>), and Sanguinea is homozygous for the other non-functional allele (ksn<sup>null</sup>/ksn<sup>null</sup>). Other old Chinese cultivars with once-flowering (OF) behaviour all have the wild KSN allele (KS<sup>N</sup>). The ap2 allele (a transposon-inserted allele) of AP2-like [5, 11] exists in all old Chinese cultivars with a double flower phenotype (Table 5).

## Introgression of Chinese S-alleles into European roses

153 rose cultivars with a variety of breeding histories were genotyped to test for the presence of ksn and ap2 and their associated S<sub>C</sub> (See Supplementary data Table D6 for the list of cultivars). The frequencies of rose cultivars with these genes were calculated, along with their breeding periods (Fig. 5). Only 25% of European roses bred before 1850 had either ksn<sup>copia</sup> or ksn<sup>null</sup>. In particular, roses in the Gallica, Damask, Centifolia, and Alba groups have neither ksn<sup>copia</sup> nor ksn<sup>null</sup> and show no signs of introgression from Chinese roses (Table S10-1). In contrast, some old roses in the Moss group, such as Mousseline (bred in 1855), had ksn alleles, indicating the onset of introgression from Chinese CF roses in this period. These roses also had S<sub>C1</sub> or S<sub>C2</sub>, demonstrating the parallel introgression of S<sub>C</sub> into European roses.

The frequency of roses carrying ksn alleles dramatically increased to 86% of roses bred from 1850–1900, gradually increasing to 100% in roses bred from 1980–2020 (Fig. 5). On the other hand, the frequency of roses with S<sub>C</sub> reached a peak (95%) in roses bred from 1850–1900, slightly decreasing to 83% in roses bred from 1900–1940, but not changing much until recently (Fig. 5). Due to the small number of samples, these changes in the frequency of rose cultivars after the period II are not statistically significant. The frequency of rose cultivars with ksn in the period II (86%) is not significantly different from that in the period V (100%) (Fisher's exact test, p = 0.0507). The frequency of rose cultivars with S<sub>C</sub> in the period II (95%) is not significantly different from that in the period V (82%) (Fisher's exact test, p = 0.2316).

The S<sub>C1</sub> and S<sub>C2</sub> are major S<sub>C</sub> alleles, and the frequency of roses with the S<sub>C3</sub> increased during the 19th century (Fig. 5c). In contrast, roses with the S<sub>C4</sub> and S<sub>C5</sub> were found with low frequency (<10%) in modern cultivars.

The frequency of roses with ap2 showed the same trend of frequency changes with breeding periods as ksn (Fig. 5), suggesting that the ap2 allele also originated from

<footnote><sup>a</sup>M= Multiple plants with different S-alleles of S<sub>x</sub> and S<sub>y</sub> were used for experiments. TF = The Fairy, OB = Rosa chinensis "Old Blush". Rm = R. multiflora.</footnote>

<footnote><sup>b</sup>The common S-alleles shared by parents are underlined. S<sub>C1</sub> is S-RNase in Chr3 of the Old Blush genome database [6], and S<sub>C2</sub> is the new S-RNase identified in RC0 of the Old Blush genome database [5]. Other S-alleles are S-RNase genes identified in R. multiflora, S<sub>7</sub>, S<sub>9</sub>, S<sub>11</sub>, and S<sub>13</sub> are orthologous to S<sub>C1</sub> and S<sub>6</sub>, S<sub>10</sub>, and S<sub>12</sub> are orthologous to S<sub>C2</sub> (Fig. 1c). S<sub>x</sub> and S<sub>y</sub> indicate anonymous S-alleles (S<sub>x</sub> ≠ S<sub>y</sub> in the same row).</footnote>
<table><caption>**Table 4.** Linkage between the S-locus and important ornamental traits in the rose</caption><thead><tr><th></th><th></th><th>S-locus</th><th>AP2-like<br>(Double-flower)</th><th>KSN<br>(Continuous-flowering)</th></tr></thead><tbody><tr><td>Genomic position (bp)<sup>a</sup></td><td></td><td>5 488 716 - 5 498 387</td><td>14 492 546 - 14 506 765</td><td>18 979 892 - 18 989 895</td></tr><tr><td>Gene annotation</td><td></td><td>Chr3g0455911</td><td>Chr3g0468481 - Chr3g0468491</td><td>Chr3g0473011 Chr3g0473021</td></tr><tr><td>Physical distance (Mbp)</td><td></td><td>0</td><td>9.0</td><td>13.5</td></tr><tr><td>Recombination<br>frequency (%)<sup>b</sup></td><td>RW</td><td>-</td><td>13% (13/97)</td><td>20% (19 × 97)</td></tr><tr><td></td><td>TF</td><td>-</td><td>40% (39/97)</td><td>NA</td></tr><tr><td></td><td>93/1–119</td><td>-</td><td>22% (11/50)</td><td>NA</td></tr></tbody></table>

<table><caption>**Table 5.** Genotypes of genes controlling important ornamental characteristics in old Chinese rose cultivars. The presence (●) or absence (×) of each gene was determined by PCR. Roses are classified into five types based on the KSN genotypes: Type 1 = ksn<sup>copia</sup> / ksn<sup>null</sup>; Type 2 = ksn<sup>copia</sup> / ksn<sup>copia</sup>; Type 3 = ksn<sup>null</sup> / ksn<sup>null</sup>; Type 4 = ksn<sup>copia</sup> / KSN<sup>w</sup>; Type 5 = KSN<sup>w</sup> / KSN<sup>w</sup>. Roses in Types 1–3 are continuous-flowering (CF), and those in Types 4–5 are once-flowering (OF).</caption><thead><tr><th rowspan="2">Type</th><th rowspan="2">Name</th><th colspan="5">S-allele</th><th colspan="3">KSN</th><th>AP2-like</th><th colspan="2">Phenotype<sup>†</sup></th></tr><tr><th>S<sub>C1</sub></th><th>S<sub>C2</sub></th><th>S<sub>C3</sub></th><th>S<sub>C4</sub></th><th>S<sub>C5</sub></th><th>ksn<sup>copia</sup></th><th>ksn<sup>null</sup></th><th>KSN<sup>w</sup></th><th>ap2</th><th>Bloom</th><th>Flower</th></tr></thead><tbody><tr><td rowspan="3">1</td><td>Rosa chinensis "Old Blush"</td><td>●</td><td>●</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>×</td><td>●</td><td>CF</td><td>D</td></tr><tr><td>Hume's Blush Tea-scented China</td><td>●</td><td>●</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>×</td><td>●</td><td>CF</td><td>D</td></tr><tr><td>Rosa chinensis "Mutabilis"</td><td>●</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>●</td><td>×</td><td>×</td><td>CF</td><td>S</td></tr><tr><td rowspan="3">2</td><td>Rosa chinensis</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>●</td><td>CF</td><td>D</td></tr><tr><td>Slater's Crimson China (×2)</td><td>●</td><td>×</td><td>●</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>●</td><td>CF</td><td>D</td></tr><tr><td>Slater's Crimson China (×3)</td><td>●</td><td>●</td><td>×</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>●</td><td>CF</td><td>D</td></tr><tr><td>3</td><td>Rosa chinensis "Sanguinea"</td><td>×</td><td>●</td><td>●</td><td>×</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>CF</td><td>S</td></tr><tr><td rowspan="3">4</td><td>Rosa chinensis "Single white-eye"</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>×</td><td>OF</td><td>S</td></tr><tr><td>Rosa chinensis "Narrow-leaflet"</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>×</td><td>×</td><td>×</td><td>OF</td><td>D</td></tr><tr><td>Rosa chinensis "Major"</td><td>×</td><td>×</td><td>●</td><td>●</td><td>×</td><td>●</td><td>×</td><td>●</td><td>×</td><td>OF</td><td>D</td></tr><tr><td rowspan="5">5</td><td>Fortune's Double Yellow</td><td>●</td><td>×</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>OF</td><td>D</td></tr><tr><td>Rosa odorata var. erubescens</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>OF</td><td>D</td></tr><tr><td>Parks' Yellow Tea-scented China</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>OF</td><td>D</td></tr><tr><td>Rosa multiflora "Carnea"</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>OF</td><td>D</td></tr><tr><td>Rosa odorata "Double Light Yellow"</td><td>×</td><td>●</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>●</td><td>●</td><td>OF</td><td>D</td></tr></tbody></table>

Chinese roses. Old European roses with the DF phenotype, such as Rosa gallica officinalis, Quatre Saisons, and Chapeau de Napoleon, had no ap2 allele (Table S10-1), indicating that there is another genetic factor for DF phenotype.

## The wild origin of S<sub>c</sub>-alleles

By screening a total of 95 plants from 25 wild Rosa species (Table S11-1) with S<sub>c</sub>-specific PCRs, putative wild ancestors of S<sub>c1-4</sub> were identified (Table 6). Positive PCR amplifications with S<sub>c1</sub>-specific primers were observed for two species, R. multiflora and Rosa brunonii (Table S11-2). The sequencing of the PCR products shows that R. multiflora sequences are 100% identical to the S<sub>c1</sub>-sequence of Old Blush, while R. brunonii has a 98.4% identity with S<sub>c1</sub> (8 SNPs per 500 bp). In the RNA-seq analysis of pistil-expressed genes, it was already found that three R. multiflora plants (Rm1, Rm2, Rm3) have S<sub>c1</sub> S-RNase (= Rm3D_S<sub>c1</sub> in Fig. 1c), with 100% identical amino

acid sequences to Old Blush S<sub>c1</sub> (= Rc3D_S<sub>c1</sub>). The cDNA sequence of Rm3D_S<sub>c1</sub> is also 100% identical to Rc3D_S<sub>c1</sub> (Supplementary data Table D4).

No wild roses were found with positive PCR amplifications from S<sub>c2</sub>-specific primers (Supplementary data Table D6). However, it was found that the genome re-sequencing individual of R. gigantea (SRR6175515) has 100% identical sequences to S<sub>c2</sub> S-RNase of Old Blush (Fig. S11-2). For the S<sub>c3</sub> S-RNase isolated from Slater's Crimson China, PCR amplifications with S<sub>c3</sub>-specific primers were found for three species, R. chinensis var. spontanea, R. gigantea, and R. multiflora var. cathayensis (Table 6). Sequencing of the PCR products showed that only R. chinensis var. spontanea has a 100% identical sequence with S<sub>c3</sub> of Slater's Crimson China (Table S11-3). For the S<sub>c4</sub> S-RNase isolated from R. chinensis, PCR amplifications with S<sub>c4</sub>-specific primers were found for three species, R. gigantea, Rosa soulieana, and Rosa helenae (Table 6). Sequencing of the PCR products showed that only one

<footnote><sup>a</sup>Chr3 in the homozygous genome of R. chinensis "Old Blush" [6]</footnote>

<footnote><sup>b</sup>Recombination frequencies in S–AP2 and S–KSN were estimated in two diploid F₁-hybrid populations: the FW population consists of 97 individuals derived from the pollen donor RW and the seed parent TF [19]. 94/1 population consists of 50 individuals derived from the pollen donor 93/1–117 and the seed parent 93/1–119 [41]. For the two populations, genotypes of S, AP2, and KSN were determined and used to calculate the recombination frequencies. See Supplementary information 9 for detailed methods for genotyping.</footnote>

<footnote>†, CF = Continuous-flowering, OF = Once-flowering, D = double-flower, S = single-flower.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Frequency of rose cultivars with genes that originated from old Chinese cultivars. (a) Frequency of roses carrying CF genes (ksn), associated S genes ($S_C$), and double-flower gene (ap2). ksn indicates the frequency with which roses have $ksn^{copia}$ and/or $ksn^{null}$; $S_C$ indicates the frequency with which roses have $S_{C1}$, $S_{C2}$, $S_{C3}$, $S_{C4}$, and/or $S_{C5}$. The number of rose cultivars studied are shown above breeding periods. (b) Individual frequency of roses that have $ksn^{copia}$ or $ksn^{null}$, and those roses that have $S_{C1}$, $S_{C2}$, $S_{C3}$, $S_{C4}$, or $S_{C5}$. Breeding periods: I = before 1850, II = 1850–1900, III = 1900–1940, IV = 1940–1980, V = 1980–2020. The data used to calculate frequencies are available in Supplementary data Table D6.</figcaption></figure>

<table><caption>**Table 6.** A survey of Chinese wild *Rosa* species with $S_C$-specific PCRs to identify candidate ancestral species of old continuous-flowering roses</caption><thead><tr><th rowspan="2">Section</th><th rowspan="2">Species</th><th rowspan="2">No. of individuals studied<sup>a</sup></th><th colspan="5">PCR tests for the presence/absence of specific S-RNase alleles<sup>b</sup></th></tr><tr><th>$S_{C1}$</th><th>$S_{C2}$</th><th>$S_{C3}$</th><th>$S_{C4}$</th><th>$S_{C5}$</th></tr></thead><tbody><tr><td rowspan="3">Chinenses</td><td>Rosa chinensis var. spontanea</td><td>34</td><td></td><td></td><td>★<sup>5</sup></td><td></td><td></td></tr><tr><td>Rosa gigantea<sup>c</sup></td><td>8</td><td></td><td>★<sup>1</sup></td><td>Δ<sup>2</sup></td><td>★<sup>1</sup></td><td></td></tr><tr><td>Rosa lucidissima</td><td>2</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="10">Synstylae</td><td>Rosa multiflora</td><td>10</td><td></td><td>★<sup>1</sup></td><td></td><td></td><td></td></tr><tr><td>Rosa multiflora var. cathayensis</td><td>9</td><td></td><td></td><td>Δ<sup>2</sup></td><td></td><td></td></tr><tr><td>Rosa brunonii</td><td>4</td><td></td><td>Δ<sup>2</sup></td><td></td><td></td><td></td></tr><tr><td>Rosa helenae</td><td>2</td><td></td><td></td><td></td><td>Δ<sup>2</sup></td><td></td></tr><tr><td>Rosa soulieana</td><td>2</td><td></td><td></td><td></td><td>Δ<sup>1</sup></td><td></td></tr><tr><td>Rosa rubus</td><td>4</td><td></td><td></td><td></td><td></td><td>Δ<sup>1</sup></td></tr><tr><td>Rosa luciae</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rosa moschata</td><td>2</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rosa anemoneflora</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

plant (pink flower type) of *R. gigantea* has a 100% identical sequence to the $S_{C4}$ (Table S11-4). For the $S_{C5}$ S-RNase isolated from Mutabilis, a positive PCR amplification with $S_{C5}$-specific primers was found only in a wild individual of *R. rubus* (Table 6), while the sequencing showed that the PCR product from *Rosa rubus* is not identical to $S_{C5}$ of Mutabilis (4 SNPs per 221 bp; Fig. S11-3). The putative genetic connections inferred from the shares of S-alleles are summarized in Figure 6.

## Discussion

### Identification and validation of the Rosa S-locus

The *Rosa* S-locus was identified and confirmed by clarifying some unresolved issues in previous studies. The genome sequencing of *R. chinensis* “Old Blush” proposed a candidate region of the S-locus [5]. This region corresponds to the region of 3C gene (Fig. 1a). We concluded that the 3C gene is not the true S-RNase because of (i) the low or undetectable expressions of 3C gene in the

<footnote><sup>a</sup>Information of sample source is available in Supplementary data Table D6.</footnote>

<footnote><sup>b</sup>★ = no. of individuals with 100% identical sequences to CF cultivars, Δ = no. of individuals with positive PCR amplifications, while the sequences were not identical to those of CF cultivars. Empty cells = no PCR amplifications for all individuals studied.</footnote>

<footnote><sup>c</sup>Results of one individual is based on a genome data (SRR6175515).</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Hypothetical wild ancestral species of old Chinese cultivars inferred from the shares of S-alleles. Five S-alleles (S<sub>C1</sub>, S<sub>C2</sub>, S<sub>C3</sub>, S<sub>C4</sub>, and S<sub>C5</sub>) were identified from old Chinese CF cultivars and can be traced to four wild ancestors. The uppermost roses are CF cultivars, and the others are once-flowering cultivars and species. The shares of S-alleles are shown by dotted line connections. Cultivars are classified into five types based on their KSN genotypes (Table 5), indicated by different colored boxes.</figcaption></figure>

pistil (Fig. 1b) and (ii) the low level of sequence divergence between alleles and among individuals (Fig. 1c). Vieira et al. [15] analyzed another genome of Old Blush [6] and firstly reported that the 3D gene is the Rosa S-RNase. The S-RNase “Rchinensis1_3-Rchinensis2_27” in Vieira et al. [15] corresponds to the 3D gene (= S<sub>C1</sub> S-RNase) in this study. Vieira et al. [15] validated the S-RNase based on the results of no fruit set of a few individuals pollinated by other individuals with the same S-genotypes. This cannot exclude the possibility that inbreeding depression results in no fruit set. We provide strong evidence that the 3D

S-RNase is the pistil determinant of the SI in roses by genotyping more than 400 seeds produced by a number of half-compatible pairs of roses (Table 3).

This study also identified the 0A gene (= S<sub>C2</sub> S-RNase) in the contig RC0 of the Old Blush genome database (Fig. 1c) and suggested via the mapping approach (Fig. 2,3) that it is an unidentified allele of the S<sub>C1</sub> S-RNase of Old Blush. The sequence analysis of the Rosa S-locus suggested that the contig RC0 can be integrated into a poly-N region of the S-locus of RC3 (Fig. 4). The perfect co-segregation of SLF5 in RC3 and of S<sub>C2</sub> S-RNase in RC0 in a mapping
population (n=97; See Supplementary information 8) supports this hypothesis. Furthermore, we identified S<sub>C2</sub> S-RNase-like genes from the pistil transcriptomes of wild individuals of R. multiflora and R. rugosa (Rm0A and Rg0A genes; Table 1). The pollination experiments using half-compatible pairs of roses that share S<sub>C2</sub> or S<sub>C2</sub> S-RNase-like genes indicated that these genes are the pistil-determinants of SI in the rose (Table 3). To determine the precise location of S<sub>C2</sub> S-RNase in the S-locus of Old Blush, a BAC library will need to be constructed and sequenced, as shown by Liang et al. [42] in their analyses of the S-locus structure of Citrus.

Based on the identification of SLFs flanking to the S-RNase, we estimated that the S-locus of Old Blush spans approximately 500 kbp (Fig. 4). Chen et al. [16] reported that the S-locus of R. rugosa spanned 667kbp, including one S-RNase and 19 F-box genes (Table S7-1b). In order to confirm the physical size of the Rosa S-locus, further co-segregation analyses of SLFs and S-RNase are required. The evolutionary divergence analysis of SLFs in comparison with S-RNase (Fig.S7-3) suggest that the SI of Rosa is controlled by the non-self recognition system as previously reported by Vieira et al. [15]. In the non-self recognition system, polyploidization will break down the SI [12, 13]. In support to the prediction, we found that the colchicine-induced chromosome doubling of a diploid rose resulted in self-compatible tetraploid (Table S6-1).

## Insights into the breeding of CF roses

The estimations of the degree of genetic linkage between the S-locus and the *ksn* controlling CF and the *ap2* controlling DF indicate weak but significant genetic linkages between them. Since the CF is a recessive characteristic, the CF rose (*ksn*-homozygote) cannot be created if *ksn* remains to link with a same S-allele. The 20% of recombination frequency (Table 4) suggests that when we cross roses carrying the *ksn* linked with a same S-allele, only 20% seedlings are expected to be CF. Therefore, information of S-alleles linked with *ksn* is helpful to make the CF rose breeding effective. The SI constraint on breeding would be lower in DF than in CF. Since DF is a dominant characteristic, one *ap2*-allele is enough to make DF roses [5, 11]. Furthermore, we suggest that there had been already other genes or alleles for DF in European roses before the introgression of *ap2* from China (Table S10-1).

Further studies are necessary to confirm and assess the degrees of genetic linkages between the S-locus and the loci controlling important ornamental traits, such as CF, DF, thornlessness, and resistance against black spot disease. Our estimation of the degree of genetic linkage depends on a few hybrid populations (Table 4), which cannot represent diverse rose cultivars used for breeding. In addition, the low recombination rates associated with regions adjacent to breakpoints in inversion heterozygotes [43] might result in a stronger genetic linkage of the S-locus with *ksn*null than with *ksn*copia.

As we expected, S<sub>C</sub> alleles were introduced into European roses in the 18th century with the CF allele *ksn*. However, the frequency of rose cultivars having S<sub>C</sub> remains high (> 80%) in modern cultivars, although it tended to decrease from 95% of cultivars from 1850–1900 (Fig. 5a). Due to the fact that many rose cultivars possess the same Sc alleles, S-genotyping in advance of breeding is helpful for making the diploid rose breeding effective. Furthermore, due to the fact that the two Sc alleles S<sub>C4</sub> and S<sub>C5</sub> are not common (< 10%) in modern cultivars (Fig. 5c), original Chinese CF cultivars of S<sub>C4</sub> (R. chinensis) and S<sub>C5</sub> (Mutabilis) are still useful breeding materials to introduce these rare S-alleles into CF cultivars.

As the SI can break down with polyploidization (Table S6-1), breeding between polyploid cultivars may not necessarily consider SI constraints. However, substantial portion (22–35%) of modern cultivars are estimated to be diploid [44, 45], and most wild species and Chinese old cultivars are diploid. Therefore, diploid rose breeding is still an important part of rose breeding.

## Insights into the wild origin of CF roses

We identified the putative ancestors of the S<sub>C</sub>-alleles of old Chinese cultivars by screening wild roses with the S<sub>C</sub>-specific primers (Table 6). The results confirmed the hybrid origin of the old Chinese cultivars with R. chinensis var. *spontanea* and R. gigantea in the section *Chinenses* (*Indicae*) and the introgression from R. multiflora in the section *Synstylae* (Fig. 6). The genome sequencing of Old Blush reported a sign of introgression from the section *Synstylae* [6], but which species of the section is involved in the formation of the Old Blush genome is not clarified. Yang et al. [30] discussed that R. multiflora is a candidate species contributing to the introgression, but there are many other candidates in the section. We investigated nine candidate species in the section *Synstylae* in southwestern China (Table 6) and demonstrated that only R. multiflora has an identical S<sub>C1</sub>-allele with Old Blush, providing a new evidence for the genetic link between R. multiflora and old Chinese cultivars.

## Materials and methods

### Genome-wide identification of candidate S-RNase

By using S-RNase of other Rosaceae crops, including *Malus domestica* (AAA79841.1), *Prunus dulcis* (AAL35960.2), *Prunus avium* (BAA36389.1, CAC27788.1), and *Prunus persica* (BAF42768.1) as queries, the Old Blush genome databases [5, 6] were searched using TBLASTN to identify candidate S-RNase. All genomic regions with significant hits (E-value <10<sup>-10</sup>), including the regions without any annotations, were listed and manually annotated to infer the coding regions. The other genome databases of *Rosa multiflora* [7] and R. rugosa [16, 39] were then searched using the candidate S-RNase genes of Old Blush as queries to identify their orthologous genes. For the
phylogenetic reconstruction of S-RNase-like genes in roses, deduced amino acid sequences were aligned using MUSCLE and then converted back to DNA sequences. The maximum-likelihood phylogenetic tree was constructed from the nucleotide protein-coding sequence alignment by FastTree [37, 38] using the Jukes-Cantor model of nucleotide evolution.

## RNA-seq analysis

To confirm the expression of candidate S-RNase in the pistil and to isolate new S-RNase alleles, RNA sequencing (RNA-seq) was used. Flower buds one or a few days before anthesis were collected from three individuals of Old Blush, eight individuals of *R. multiflora*, and one individual of *R. rugosa*. Pistils and stamens were collected from the buds and immediately frozen with liquid nitrogen. Total RNAs were extracted using a commercial kit according to the protocol described in Dubois *et al.* [46], and RNA-seq data was obtained through poly-A purification and the 150 bp paired-end method. The RNA-seq reads were mapped to the CDSs of specific genes and whole genome data to calculate FPKM (Fragments per kilobase of exon per million reads mapped) values. To analyze the data of wild *R. multiflora* and *R. rugosa* plants without any specific reference genome databases, the RNA-seq reads were assembled first, followed by the construction of a local database of mRNA sequences, and the database was blasted using the Old Blush genes as queries to identify orthologous genes. The assembly of RNA-seq data was conducted by using the velvet and tadpole algorithm, with default parameters in Geneious Prime 2020.

## Genetic mapping of the new S-RNase in chromosome zero

Genetic mapping of a newly-identified S-RNase in the contig that is not assigned to seven chromosomes was performed by using an F₁ diploid mapping population (FW) [19, 47, 48]. Three new genetic markers linked to the candidate S-RNase were added to the previous map to estimate the genomic position of the S-gene (Supplementary information 4).

## SLF identification

ORFs longer than 1kbp were extracted from the 1Mbp genomic regions surrounding the S-RNase in the Old Blush, *R. multiflora*, and *R. rugosa* genomes, and the F-box genes were identified by a Blast search of the ORFs.

## Validation of the S-RNase based S-genotyping by pollination experiments

A total of 20 pairs of diploid roses that share one S-allele (i.e. half-compatible) were used for pollination experiments to validate the candidate S-RNase gene. Old Blush (SC₁/SC₂), The Fairy (SC₂/S₂₁), *R. chinensis* "Single white-eye" (SC₁/S₁₂), and 11 wild individuals of *R. multiflora* were selected based on their S-genotypes. Before anthesis, flower buds were bagged to prevent open pollination. Petals and anthers were removed at the balloon stage,

and outcross pollen grains were put on the exposed stigma. Pollinations were carried out from April to May in 2018 and 2020, and matured fruits were collected from September to October of the same years. The fruits were opened in the laboratory, the achenes (seeds) were collected, and the S-genotypes of the seeds were determined by PCR (Supplementary information 5).

## Linkage between the S-locus and important ornamental traits

Two F₁-diploid mapping populations (FW [19] and 94/1 [41]) were used to estimate the recombination frequencies between the S-locus and the genes controlling CF (KSN) and DF (AP2-like). The genotyping of the S-locus, KSN, and AP2-like was performed by PCR (Supplementary information 9), and recombination frequencies were calculated.

## Introgression of the Chinese S-alleles into modern roses

A total of 153 rose cultivars were selected from a wide range of breeding ages and classifications, and young leaves were collected from 20 rose gardens and nurseries in Japan from 2014–2020 (Supplementary data Table D6). DNA was extracted from the young leaves by using the Nucleospin Plant II kit (Macherey-Nagel) according to the manufacturer's protocol. A PCR was performed to test whether the roses have specific S-alleles, mutated-alleles *ksn*, wild-allele KSN<sup>W</sup>, or the mutated-allele *ap2* by using EmeraldAmp PCR Master Mix (TaKaRa) with thermal cycling: (1) 2 min of 95°C; (2) 30 s of 95°C; (3) 30 s of T<sub>m</sub>°C; (4) 20–50 s of 72°C; and (5) go back to step (2) 29 times. Information on primer sequence and annealing temperature (T<sub>m</sub>) for PCR is available from Supplementary data Table D1. To identify other Chinese S-alleles originally linked with *ksn*, we used three old Chinese CF cultivars, Slater's Crimson China, Mutabilis, and Rosa chinensis, extracted RNAs from their pistils, and performed RT-PCR with primer sets designed on conserved sites of S-RNase. Partial or full sequences of three new S-alleles were determined and named as S<sub>C3</sub>, S<sub>C4</sub>, and S<sub>C5</sub>, and specific primers were designed for each (Supplementary information 10).

## Wild roses carrying the same S-alleles as old cultivars in China

By screening wild rose species carrying the same S-alleles as the old Chinese CF cultivars, the wild ancestral species were inferred. A total of 25 Rosa species were surveyed (Table S11-1), with a focus on 13 species in southwestern China (Sichuan and Yunnan Provinces; Fig. S11-1), where a wild type of *R. chinensis*, named *R. chinensis* var. *spon-tanea*, is naturally distributed. Field collections of wild plants in southwestern China were performed during the flowering seasons of 2018 and 2019. Forty wild plants of seven species were sampled, and S-genotyping was performed by PCR. The nucleotide sequences of the PCR products were determined to confirm the results.
# Acknowledgements

This study was supported by the JSPS KAKENHI No. 24688004 and 17H04616. We would like to express our appreciation for Cheng Tingting and Jin Tang in Sichuan University, and Ryosuke Sakamoto in Osaka Institute of Technology, for their helpful support during the fieldwork in China. We are also appreciative of the support from The Rose Culture Institute, Japan, and Drs. Masaki Ochiai and Kunio Yamada at Gifu University, for providing rose materials.

# Contributions

K.K. designed the project and wrote the manuscript. K.K., Y.U., W.L., and G.W. conducted fieldwork in China. K.K., S.M., M.L., and T.D. performed pollination experiments. K.K., T.H., and S.O. analyzed candidate genes, L.H.O., F.F., M.L., and T.D. provided mapping populations. All authors are involved in final manuscript editing.

# Data availability

Sequence data are available in Supplementary data Tables. NGS read data of RNA-seq are available from DDBJ, BioProject PRJDB12320, with accession numbers from DRR321244 to DRR321267.

# Conflict of interests

The authors declare no conflicts of interest associated with this manuscript.

# Supplementary data

Supplementary data is available at *Horticulture Research* online.

# References

1. Young MA, Schorr P, Baer R. *Modern Roses* 12. Shreveport: American Rose Society; 2007.

2. Leus L, Van Laere K, De Riek J *et al*. Ornamental crops. In: Van Huylenbroeck J (ed.), *Handbook of Plant Breeding* Vol. 11. Cham: Springer, 2018, 719–67.

3. Zlesak DC. Rose. In: Anderson NO, ed. *Flower Breeding and Genetics*. Springer: Dordrecht, 2006,695–740.

4. Debener T, Linde M. Exploring complex ornamental genomes, the rose as a model plant. *Crit Rev Plant Sci.* 2009;**28**:267–80.

5. Hibrand Saint-Oyant L, Ruttink T, Hamama L *et al*. A high-quality genome sequence of *Rosa chinensis* to elucidate ornamental traits. *Nat Plants*. 2018;**4**:473–84.

6. Raymond O, Gouzy J, Just J *et al*. The *Rosa* genome provides new insights into the domestication of modern roses. *Nat Genet.* 2018;**50**:772–7.

7. Nakamura N, Hirakawa H, Sato S *et al*. Genome structure of *Rosa multiflora*, a wild ancestor of cultivated roses. *DNA Res.* 2018;**25**: 113–21.

8. Scalliet G, Piola F, Douady CJ *et al*. Scent evolution in Chinese roses. *PNAS*. 2008;**105**:5927–32.

9. Magnard JL, Roccia A, Caissard JC *et al*. Biosynthesis of monoter- pene scent compounds in roses. *Science*. 2015;**349**:81–3.

10. Iwata H, Gaston A, Remay A *et al*. The TFL1 homologue KSN is a regulator of continuous flowering in rose and strawberry. *Plant J.* 2012;**69**:116–25.

11. François L, Verdenaud M, Fu X *et al*. A miR172 target-deficient AP2-like gene correlates with the double flower phenotype in roses. *Sci Rep.* 2018;**8**:12912.

12. Sassa H. Molecular mechanism of the S-RNase-based gametophytic self-incompatibility in fruit trees of Rosaceae. *Breed Sci.* 2016;**66**:116–21.

13. Kubo K, Entani T, Takara A *et al*. Collaborative non-self recognition system in S-RNase-based self-incompatibility. *Science*. 2010;**330**:796–9.

14. Akagi T, Henry IM, Morimoto T *et al*. Insights into the Prunus-specific S-RNase-based self-incompatibility system from a genome-wide analysis of the evolutionary radiation of S locus-related F-box genes. *Plant Cell Physiol.* 2016;**57**:1281–94.

15. Vieira J, Pimenta J, Gomes A *et al*. The identification of the Rosa S-locus and implications on the evolution of the *Rosaceae* gametophytic self-incompatibility systems. *Sci Rep.* 2021;**11**:3710.

16. Chen F, Su L, Hu S *et al*. A chromosome-level genome assembly of rugged rose (*Rosa rugosa*) provides insights into its evolution, ecology, and floral characteristics. *Hortic Res.* 2021;**8**:141.

17. Du J, Ge C, Li T *et al*. Molecular characteristics of S-RNase alleles as the determinant of self-incompatibility in the style of *Fragaria viridis*. *Hortic Res.* 2021;**8**:185.

18. Smulders MJM, Arens P, Bourke PM *et al*. In the name of the rose, a roadmap for rose research in the genome era. *Hortic Res.* 2019;**6**:65.

19. Kawamura K, Hibrand-Saint Oyant L, Crespel L *et al*. Quantitative trait loci for flowering time and inflorescence architecture in rose. *Theor Appl Genet.* 2011;**122**:661–75.

20. Debener T, Bretzke M, Dreier M *et al*. Genetic and molecular analyses of key loci involved in self incompatibility and floral scent in roses. *Acta Hortic.* 2010;**870**:183–90.

21. Zhou NN, Tang KX, Jeauffre J *et al*. Genetic determinism of prickles in rose. *Theor Appl Genet.* 2020;**133**:3017–35.

22. Lopez Arias DC, Chastellier A, Thouroude T *et al*. Characterization of black spot resistance in diploid roses with QTL detection, meta-analysis and candidate-gene identification. *Theor Appl Genet.* 2020;**133**:3299–321.

23. Randoux M, Davière JM, Jeauffre J *et al*. RoKSN, a floral repres- sor, forms protein complexes with RoFD and RoFT to regulate vegetative and reproductive development in rose. *New Phytol.* 2014;**202**:161–73.

24. Soufflet-Freslon V, Araou E, Jeauffre J *et al*. Diversity and selec- tion of the continuous-flowering gene, RoKSN, in rose. *Hortic Res.* 2021;**8**:76.

25. Hurst CC. Notes on the origin and evolution of our garden roses. *J Roy Hort Soc.* 1941;**66**: 73-82, 242-250, 282-289.

26. Liorzou M, Pernet A, Li S *et al*. Nineteenth century French rose (*Rosa* sp.) germplasm shows a shift over time from a European to an Asian genetic background. *J Exp Bot.* 2016;**67**: 4711–25.

27. Meng J, Fougère-Danezan M, Zhang LB *et al*. Untangling the hybrid origin of the Chinese tea roses, evidence from DNA sequences of single-copy nuclear and chloroplast genes. *Plant Syst Evol.* 2011;**297**:157–70.

28. Zhu ZM, Gao XF, Fougère-Danezan M. Phylogeny of *Rosa* sections *Chinenses* and *Synstylae* (*Rosaceae*) based on chloroplast and nuclear markers. *Mol Phylogenet Evol.* 2015;**87**:50–64.

29. Tan J, Wang J, Luo L *et al*. Genetic relationships and evolution of old Chinese garden roses based on SSRs and chromosome diversity. *Sci Rep.* 2017;**7**:15437.
30. Yang C, Ma Y, Cheng B et al. Molecular evidence for hybrid origin and phenotypic variation of Rosa section Chinenses. *Genes*. 2020;11:996.

31. Debray K, Le Paslier MC, Bérard A et al. Unveiling the patterns of reticulated evolutionary processes with phylogenomics: hybridization and polyploidy in the genus Rosa. *Syst Biol*. 2022;71:547–69.

32. Cui W-H, Du XY, Zhong MC et al. Complex and reticulate origin of edible roses (Rosa, Rosaceae) in China. *Hortic Res*. 2022;9:uhab051.

33. Meng J, He SL, Li DZ et al. Nuclear genetic variation of Rosa odorata var. gigantea (Rosaceae), population structure and conservation implications. *Tree Genet Genomes*. 2016;12:65.

34. Jian HY, Zhao L, Zhang H et al. Phylogeography and population genetics of Rosa chinensis var. spontanea and R. lucidissima Complex, the important ancestor of modern roses. *Fornt. Plant Sci*. 2022;13:851396.

35. Schueler S, Tusch A, Scholz F. Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites. *Mol Ecol*. 2006;15:3231–43.

36. Aguiar B, Vieira J, Cunha AE et al. Convergent evolution at the gametophytic self-incompatibility system in *Malus* and *Prunus*. *PLoS One*. 2015;10:e0126138.

37. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. *Mol Biol Evol*. 2009;26:1641–50.

38. Prince MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. *PLoS One*. 2010;5:e9490.

39. Zang F, Ma Y, Tu X et al. A high-quality chromosome-level genome of wild Rosa rugosa. *DNA Res*. 2021;28:dsab017.

40. Lyons E, Freeling M. How to usefully compare homologous plant genes and chromosomes as DNA sequences. *Plant J*. 2008;53:661–73.

41. Debener T, Mattiesch L. Construction of a genetic linkage map for roses using RAPD and AFLP markers. *Theor Appl Genet*. 1999;99:891–9.

42. Liang M, Cao Z, Zhu A et al. Evolution of self-compatibility by a mutant *S*m-*RNase* in citrus. *Nat Plants*. 2020;6:131–42.

43. Navarro A, Betran E, Barbadilla A et al. Recombination and gene flux caused by gene conversion and crossing over in inversion heterokaryotypes. *Genetics*. 1997;146:695–709.

44. Ueckert J, Byrne DH, Crosby K et al. The utilization of the polyploid nature of roses. *Acta Hortic*. 2015;1064:73–8.

45. Zlesak DC. Pollen diameter and guard cell length as predictors of ploidy in diverse rose cultivars, species, and breeding lines. *Flori Ornam Biotech*. 2009;3:53–70.

46. Dubois A, Raymond O, Maene M et al. Tinkering with the C-function, a molecular frame for the selection of double flowers in cultivated roses. *PLoS One*. 2010;5:e9288.

47. Kawamura K, Hibrand-Saint Oyant L, Thouroude T et al. Inheritance of garden rose architecture and its association with flowering behaviour. *Tree Genet Genomes*. 2015;11:22.

48. Kawamura K, Hibrand-Saint Oyant L, Foucher F et al. Kernel methods for phenotyping complex plant architecture. *J Theor Biol*. 2014;342:83–92.