# Haplotype-resolved genome assembly of the diploid Rosa chinensis provides insight into the mechanisms underlying key ornamental traits

Xiaoni Zhang, Quanshu Wu, Lan Lan, Dan Peng, Huilin Guan, Kaiqing Luo, Manzhu Bao, Mohammed Bendahmane, Xiaopeng Fu, Zhiqiang Wu

► To cite this version:

Xiaoni Zhang, Quanshu Wu, Lan Lan, Dan Peng, Huilin Guan, et al.. Haplotype-resolved genome assembly of the diploid Rosa chinensis provides insight into the mechanisms underlying key ornamental traits. Molecular Horticulture, 2024, 4 (1), pp.14. 10.1186/s43897-024-00088-1. hal-04794862

HAL Id: hal-04794862

https://hal.inrae.fr/hal-04794862v1

Submitted on 21 Nov 2024

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>
<figure><img src="image_2.png" /></figure>

# Haplotype-resolved genome assembly of the diploid *Rosa chinensis* provides insight into the mechanisms underlying key ornamental traits

Xiaoni Zhang<sup>1,2,3†</sup>, Quanshu Wu<sup>1†</sup>, Lan Lan<sup>2,3,4</sup>, Dan Peng<sup>2,3</sup>, Huilin Guan<sup>1</sup>, Kaiqing Luo<sup>2,3</sup>, Manzhu Bao<sup>1</sup>, Mohammed Bendahmane<sup>1,5*</sup>, Xiaopeng Fu<sup>1*</sup> and Zhiqiang Wu<sup>2,3*</sup>

## Abstract

Roses are consistently ranked at the forefront in cut flower production. Increasing demands of market and changing climate conditions have resulted in the need to further improve the diversity and quality of traits. However, frequent hybridization leads to highly heterozygous nature, including the allelic variants. Therefore, the absence of comprehensive genomic information leads to them making it challenging to molecular breeding. Here, two haplotype-resolved chromosome genomes for *Rosa chinensis* 'Chilong Hanzhu' (2n = 14) which is high heterozygous diploid old Chinese rose are generated. An amount of genetic variation (1,605,616 SNPs, 209,575 indels) is identified. 13,971 allelic genes show differential expression patterns between two haplotypes. Importantly, these differences hold valuable insights into regulatory mechanisms of traits. *RcMYB114b* can influence cyanidin-3-glucoside accumulation and the allelic variation in its promoter leads to differences in promoter activity, which as a factor control petal color. Moreover, gene family expansion may contribute to the abundance of terpenes in floral scents. Additionally, *RcANT1*, *RcDA1*, *RcAG1* and *RcSVP1* genes are involved in regulation of petal number and size under heat stress treatment. This study provides a foundation for molecular breeding to improve important characteristics of roses.

**Keywords** *Rosa chinensis*, Haplotype-resolved genome, Petal colour, Flower development

## Core

Two haplotype-resolved chromosome genomes for highly heterozygous *Rosa chinensis* were first reported. *RcMYB114b* can function as an activator of anthocyanin biosynthesis in petals. The regulatory mechanisms of flower development under high temperature were investigated.

## Gene and Accession Numbers

Sequence data and rose genome with annotation from this article can be found in the database of the National Center for Biotechnology Information (NCBI) under the accession numbers PRJNA932466.

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

<footnote>†Xiaoni Zhang and Quanshu Wu contributed equally to this work.</footnote>

<footnote>*Correspondence:
Mohammed Bendahmane
mohammed.bendahmane@ens-lyon.fr
Xiaopeng Fu
fuxiaopeng@mail.hzau.edu.cn
Zhiqiang Wu
wuzhiqiang@caas.cn</footnote>

<footnote n="1">¹ National Key Laboratory for Germplasm Innovation & Utilization of Horticultural Crops, College of Horticulture and Forestry Sciences, Huazhong Agricultural University, Wuhan 430070, China</footnote>

<footnote n="2">² Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China</footnote>

<footnote n="3">³ Kunpeng Institute of Modern Agriculture at Foshan, Foshan 528200, China</footnote>

<footnote n="4">⁴ College of Science, Health, Engineering and Education, Murdoch University, 90 South Street, Murdoch, WA 6150, Australia</footnote>

<footnote n="5">⁵ Laboratoire Reproduction Et Development Des Plantes, INRA-CNRS-Lyon1-ENS, Ecole Normale Supérieure de Lyon, 520074 Lyon, France</footnote>
## Introduction

Roses are the most popular cut flowers in the world due to their beautiful colours and rich fragrances. The cultivation of roses can be dated back 5,000 years among civilizations in Asia (Wang 2007). Notably, *Rosa chinensis* cultivars like ‘Old Blush’ (OB) and ‘Crimson China’ were highly significant old Chinese rose in the history of rose breeding that were introduced to Europe in the eighteenth century, and a large number of modern rose varieties were cultivated through frequent crosses of these cultivars (Huylenbroeck, 2018). Then, rose production entered a period of vigorous development, during which many of the most sought-after traits were cultivated (Young et al. 2007). Roses become the dominant species in cut flower production globally, owing to its appealing petal colours, intricate floral patterns, rich fragrances and high vase life.

To meet huge market demands and expedite breeding progress, many efforts have been invested in unraveling the molecular mechanisms underlying key traits. For petal colour, the genes involved in the anthocyanin biosynthesis pathway, *dihydroflavonol-4-reductase (DFR)* and *flavonol synthase (FLS)* (Suzuki et al. 2014), *MYB*, *basic helix-loop-helix (bHLH)* and *WD40*, were reported (Koes et al. 2005, Bendahmane et al. 2013, Raymond et al. 2018). For floral morphology, the differential expression of *AGAMOUS* (*RhAG*) in roses contributed to double-petal or simple-petal flower formation (Dubois et al. 2010), similar to findings in other species (Jung et al. 2014, Wollmann et al. 2010). Moreover, a microRNA172 (miR172) target-deficient *APETALA2-like (AP2L)* gene was correlated with the double-flower phenotype in rose (François et al. 2018), peach (Gattolin et al. 2018), and carnation (Wang et al. 2020, Zhang et al. 2022). For scent, different rose varieties exhibit different principal aroma components (Scalliet et al. 2008, Magnard et al. 2015, Zhou et al. 2020a). Additionally, with the ongoing change in global temperatures, and especially the frequent occurrence of high temperature (Seneviratne et al., 2021), it has a significant impact on flower yields and morphology in roses. The mechanism behind of these phenotypes are poorly deeply understood.

An increasing number of studies have provided preliminary evidence that allelic variations may play important roles during evolution and trait selection. In rose, the presence or absence of a retrotransposon insertion in the *Koushin* (*RcKSN*) allele determined the recurrent/once blooming character (Iwata et al. 2012). The duplication and specialization of *Nudix hydrolase 1* (*NUDX1*) in Rosaceae led to geraniol production in petals (Conart et al. 2022). In *Citrus*, one allele of *AbRuby2Full* acted as an anthocyanin activator, while other *CgRuby2Short* as repressor (Huang et al. 2018). Most of the ornamental plants are highly heterozygous, especially the rose. However comprehensive

genomic information at a haplotype-resolved level is still lacking for the published rose genome (Chen et al. 2021, Raymond et al. 2018, Zhong et al. 2021).

Here, we present a high-quality, haplotype-resolved genome of *R. chinensis* ‘Chilong Hanzhu’ (‘CH’), a variety of the ‘Crimson China’ cultivar, which is one of the two significant old Chinese rose cultivars (in addition to ‘OB’). We identified a lot of variations and observed that the allelic genes exhibited differential expression between two haplotypes. We performed the regulatory mechanisms involved in important ornamental traits, and identified key related regulatory genes. Our findings provide fundamental resources and insights for accelerating the genetic improvement of roses.

## Results

### Assembly and annotation of the haplotype genome of R. chinensis ‘Chilong Hanzhu’

We generated ~25 Gb of PacBio HiFi reads and ~34 Gb raw MGI short reads from *R. chinensis* 'CH' (2n=14) (Fig. 1a and Table S1). The genome size of 'CH' was estimated to be ~538 Mb, with a heterozygosity of 2.57% by survey (Fig. S1 and Table S2). We obtained the initial assembly containing two divergent haplotypes with Hifiasm software after removing organelle sequences, resulting in genomes of 518 Mb (hA) and 541 Mb (hB), with contig N50 values of 10.89 Mb and 18.05 Mb, respectively. In addition, we observed that the heterozygous regions were mostly collapsed into single-copy homozygous content of the two genomes by KAT software suggested the correctness of genome phasing (Fig. S2). Most (97.9% and 98.4%) of the genome contigs were respectively anchored to 7 pseudochromosomes by Hi-C data (Table S3), which yielded well BUSCO evaluation (Table 1 and Table S4). To confirm the accuracy of phasing, we used the following strategy: 14 chromosomes were placed in a Hi-C heatmap with strong interaction patterns (Fig. S3); Subsequent remapping efforts demonstrated that high mapping rates for HiFi reads (99.85% and 99.84%, respectively) and MGI reads (98.58% and 98.39%, respectively) (Table S6); GC content and second generational data coverage depth shown the genomes without external pollution (Fig. S4); The quality value were 47.084 (hapA) and 47.415 (hapB); LAIs were 19.93 (hA) and 21.24 (hB); In addition, good collinearity was observed between the two haplotype genomes and between each haplotype genome separately and the previous ‘OB’ genome (Fig. 1b and Fig. S5); The evaluated switch error rate is 5.30% in our assembled genome. Collectively, these results strongly indicated that we obtained two high-quality haplotyped genomes and hapB as the primary haplotype with higher quality.

We predicted 31,744 genes in hA and 32,340 in hB, and the BUSCO evaluation results showed completeness of
<figure><img src="image_4.png" /><figcaption>**Fig. 1** The features of *Rosa chinensis* 'CH' genome. **a** The plant of *R. chinensis* 'CH'. **b** Distribution of 'CH' genomic features. (I) Circular representation of the 14 chromosomes of two haplotype genomes, length in Mb, (II) gene density, (III) GC content in 500 kb windows, and (IV) each linking line in the center of the circle connects a pair of homologous genes between the two haplotype genomes of 'CH'</figcaption></figure>

<table><caption>**Table 1** Statistics for *R. chinensis* 'CH' Haplotype genomes</caption><thead><tr><th></th><th>HapA</th><th>HapB</th></tr></thead><tbody><tr><td>Contig N50</td><td>10.89 Mb</td><td>18.05 Mb</td></tr><tr><td>scaffolds N50</td><td>69.83 Mb</td><td>74.66 Mb</td></tr><tr><td>Longest contig</td><td>92,402,791 bp</td><td>92,644,500 bp</td></tr><tr><td>Chromosome number</td><td>7</td><td>7</td></tr><tr><td>Number of contigs</td><td>342</td><td>264</td></tr><tr><td>Number of scaffolds</td><td>193</td><td>176</td></tr><tr><td>Genome-size</td><td>518,698,571 bp</td><td>541,343,551 bp</td></tr><tr><td>Rate of anchoring (%)</td><td>97.9</td><td>98.4</td></tr><tr><td>GC content (%)</td><td>38.95</td><td>38.84</td></tr><tr><td>BUSCO_genome (%)</td><td>98.4</td><td>98.4</td></tr><tr><td>TE (%)</td><td>59.6</td><td>61.83</td></tr><tr><td>Number of predicted genes</td><td>31,744</td><td>32,340</td></tr><tr><td>BUSCO_protein (%)</td><td>98.6</td><td>99.0</td></tr></tbody></table>

98.6% and 99.0%, respectively (Table 1 and Table S6, 7, 8), proving the high quality of the gene structure prediction results, which were sufficient for further analyses. We found that the predicted genes were mainly distributed at the ends of the chromosomes (Fig. 1b). Among the TEs, 59.60% (hA) and 61.83% (hB) were predicted to be transposable elements, among which Copia LTR-RTs accounted for 14.27% and 13.46% and Gypsy accounted for 12.79% and 12.94%, respectively (Table S9). In addition, 710 (hA) and 677 (hB)

transfer RNA, 1,746 and 2,083 ribosomal RNA (rRNA), 667 and 687 small nuclear RNA, and 162 and 153 microRNA (miRNA) genes were predicted using Infernal.

## Genetic variation of the two haplotype genomes in highly heterozygous 'CH'

To evaluate the divergence between the two haplotypes of 'CH', we identified genetic polymorphisms between the 7 homologous chromosome pairs. The Syteny and Rearrangement Identifier (SyRI) tool revealed 201 syntenic blocks, 1,605,616 single nucleotide polymorphisms (SNPs), 209,575 indels (1–50 bp insertion and deletion), 59 inversions, 1,169 translocations, 5 tandem repeats, 1,425 duplications and 2,489 large insertion and large deletion (>50 bp) between the two haplotypes (Table S10–11 and Fig. S6a). What's more, the biggest inversion of chromosome 2 was supported by HiFi reads mapping (Fig. S6b). Among these variations, the lengths of inversions, translocations and duplications were mainly distributed in the around 1,000 bp (Fig. 2a–c). In addition, we identified ~83 Mb of unique regions in hB, ~69 Mb of unique regions in hA and highly diverse regions in the two 'CH' haplotypes (Fig. 2d–f, Table S11). Based on syteny and sequence similarity, 16,889 pairs (33,778 genes, 52.70% of all annotated genes) which contain at least one SNP variation were considered allelic genes, alongside 6,207 "single allele" with identical coding sequences
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Genetic variation of the haplotype genome in high heterozygous *R. chinensis* 'CH'. **a–c** The length distribution of duplications (**a**), inversions (**b**) and translocations (**c**) between the two haplotypes genomes. **d–e** The genomic structural variation characteristics between the two haplotype genomes. **hA** and **hB** mean the haplotype A and haplotype B of 'CH' genome, separately. **f** Collinearity and structural variation distribution between two haplotypes of 'CH'. **g** Differential expression ratios of alleles in four tissues of roots, stems, leaves, and flowers</figcaption></figure>

between the two haplotypes. To understand the expres-
sion landscape of allelic genes, we analyzed the tran-
scriptome data of four tissues (roots, stems, leaves, and
flowers) using the Kallisto pipeline (Fig. S7). A total of
13,971 of allelic genes were differentially expressed in dif-
ferent tissues. Among these expressed allelic genes, 1,342
(~10%) pairs unequal expression (P < 0.01) between two
alleles in all tissues exhibited, as did 3,949 (~28%) pairs
in at least one tissue (Fig. 2g and Fig. S8). These differen-
tially expression of alleles may be related to the structural
variations.

Phylogenetic evolution

To understand the phylogenetic evolution of 'CH', we
searched for conserved orthologous gene families in
the genomes of *Nymphaea colorata*, *Oryza sativa*,
*Arabidopsis thaliana*, *Fragaria vesca*, *Malus domestica*,

Prunus communis, *P. persica*, *Rubus occidentalis*, *R. chinensis* 'OB', and *R. rugosa*. Among these species, *N. colorata*, *O. sativa* and *A. thaliana* were selected as outgroups. A total of 681 single-copy orthologous groups were identified with OrthoFinder to construct a species tree of 'CH' (hapB/hB) and representative Rosaceae species. 'CH' hB and 'OB' constituted sister groups belonging to *R. chinensis*, and *R. rugosa* was the closest member of the species. The divergence time of 'CH' (hB) and 'OB' was estimated to be ~1.1 MYA. In addition, that of *R. chinensis* and *R. rugosa* was estimated to be ~2.7 MYA (Fig. 3a). After replace genome hB with hA and using the same pipeline, the similar phylogenetic relationship was obtained (Fig. S9). What's more, we identified (498 and 320 respectively) expanded gene families and (1,350 and 604) contracted gene families in 'CH' hB and hA (Fig. 3b). Compared to either 'CH' hA
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Evolution of the *R. chinensis* 'CH' and influence of family gene expansion on the main components of floral aroma. **a** Phylogenetic tree for 'CH' hB and ten other eudicot species contained three outgroups. Gene family expansions are indicated in green, and contractions in black; The estimated divergence time (million years ago, MYA) is indicated at each node; numbers in brackets are the 95% confidence intervals (each center is defined as mean value). The red dot represents a calibration point. **b–c** Gene Ontology (GO) enrichment analysis of expanded gene families in 'CH' (b) and 'OB' (c). **d** Accumulation histogram of four groups of volatile organic compounds of *R. chinensis* 'OB' and 'CH'</figcaption></figure>

or hB, there were more expanded gene families in 'OB', especially terpene-related gene families (Fig. 3b–c and Fig. S10). In order to confirm whether there was difference of the volatile organic compounds present between the two cultivars, the GC-MS experiment was performed. We found that there were more terpenoids and derivatives in 'OB' than in 'CH' (Fig. 3d). By analyzing the transcriptome data of 'OB' at different flower developmental stages, we found that eight related genes were expressed (FPKM > 1) with increasing pattern of expression during the flower development (Fig. S11). Overall, the eight expanded family genes above in 'OB' may contribute to more terpenoids and derivatives in 'OB' than in 'CH'.

## Identification of pigments and differentially expressed genes in red and pink petals

To explore the molecular mechanism underlying the colour difference between red and pink petals, the metabolites of these petals were identified. Cyanidin-3,5-diglucoside and cyanidin-3-glucoside were found to be the main anthocyanin components of the red and pink petals (Fig. 4a). Cyanidin-3,5-diglucoside was the predominant component in the pink petals of 'OB', while cyanidin-3-glucoside was the predominant component in those of 'CH' (Fig. 4a).

By analyzing the transcriptome data of petals (non-coloured petals vs. colouring petals, colouring petals vs. coloured petals, and red petal regions vs. the white
<figure><img src="image_7.png" /><figcaption>c</figcaption></figure>

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

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

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

<figure><img src="image_11.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>g</figcaption></figure>

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

**Fig. 4** Potential anthocyanin regulatory pathway of red or pink petals in *R. chinensis*. a The metabolites were identified by LC-MS in 'CH' and 'OB' petals. Peak 1–3 represent cyanidin-3,5-diglucoside, pelargonidin-3,5-O-diglucoside and cyanidin-3-glucoside, respectively. b Proposed biosynthetic pathway of anthocyanin synthesis in rose. c The different 35S:RcMYB114b transgenic lines of rose turned red. d Expression of RcMYB114b (left) and anthocyanin biosynthesis pathway genes (right) in different 35S:RcMYB114b transgenic lines of rose callus. The bars indicate the SD. Stars mean significant differences (*** represents *P* < 0.001, ** represents *P* < 0.01, * represents *P* < 0.05) analyzed by one-way ANOVA. e. LC-MS analysis of 35S:RcMYB114b transgenic rose callus and control. Peak 1–3 represent cyanidin-3,5-diglucoside, pelargonidin-3,5-O-diglucoside and cyanidin-3-glucoside, respectively. f. The transient expression of RcMYB114b in apple fruit. The black triangle represents the site of the injection. g. Expression of anthocyanin biosynthesis pathway genes in RcMYB114b transgenic lines of apple. The bars indicate the SD. Stars mean significant differences (*** represents *P* < 0.001, ** represents *P* < 0.01, * represents *P* < 0.05) analyzed by one-way ANOVA
heart petal regions) in 'CH' and 'OB', we found that the expression of anthocyanin biosynthesis pathway genes, including *phenylalanine ammonia-lyase* (*PAL*), *chal-cone synthase* (*CHS*) and flavanone 3-hydroxylase (*F3H*), was higher in red petals than in pink petals and in coloured petal regions than in the white heart region of petals (Fig. 4b). In addition, there were some related differentially expressed transcription factors (TFs) identified between different petal colouring stages, such as *MYB*, *bHLH*, ethylene response factor (*ERF*) (Fig. S12). Among these genes, the expression of *RcMYB114b*, encoding MYB transcription factors (Fig. 4b and Fig. S13), showed an upwards trend during the process of petal colouring development in 'CH', while, *RcMYB308a* and *RcMYB308b* mainly in 'CH' coloured stage. And their expression trend was confirmed by real-time quantitative PCR (qRT-PCR) in 'CH' and 'OB' (Fig. S14). Notably, based on the variation between identified alleles between two haplotypes and two materials (Table S10-S11 and Fig. S15), we found that there were many variants between alleles of *RcMYB114b* promoter, suggesting that *RcMYB114b* may contribute more to the formation of petal colour.

### RcMYB114b contributes to the formation of petal colour

To investigate the function of *RcMYB114b*, the coding sequences (CDSs) of *RcMYB114b* amplified from 'OB' and 'CH' were compared, and their amino acid sequences were identical. Overexpression of *RcMYB114b* in 'OB' calli resulted in a colour change, with the transgenic calli turning red due to cyanidin-3-glucoside accumulation, while the control calli did not change colour (Fig. 4c-e). Transcriptome sequencing of the transgenic lines of rose calli revealed 2,832 differentially expressed genes (DEGs) compared to the control calli and these DEGs were associated with various KEGG enrichment terms (Fig. S16). Among the DEGs, key genes involved in the flavonoid biosynthesis pathway, such as *CHS*, *F3H*, *flavonoid 3'-hydroxylase* (*F3'H*), *DFR*, and *leucoanthocyanidin reductase* (*LAR*), and transcription factors *RcMYB114b* and *RcbHLH95*, showed significant changes in expression level (Fig. S16c, d). Moreover, the qRT-PCR results validated the upregulation of gene expression for *RcCHS*, *chalcone isomerase* (*RcCHI*), *RcF3H*, *RcF3'H*, *RcDFR*, *anthocyanidin synthase* (*RcANS*), *RcLAR* and *anthocyanidin 3-O-glucosyltransferase* (*Rc3GT*) in transgenic calli, consistent with the RNA-seq results (Fig. 4d). Furthermore, the transient transformation of *RcMYB114b* into *M. pumila* (Rosaceae) demonstrated similar changes: the injection site turned red, and related gene (*MdCHS*, *MdCHI*, *MdF3'H*, *MdDFR*, *MdANS* and *lavonoid 3-O-glucosyltransferase* (*MdUFGT1*)) expression was upregulated, while the empty vector injection site showed

no change in colour (Fig. 4f-g). Collectively, these results indicate that *RcMYB114b* functions as an activator of anthocyanin biosynthesis, influencing cyanidin-3-glucoside accumulation in petals by activating the expression of genes involved in anthocyanin synthesis pathway and regulating the formation of red colouration in petals.

To clarify why *RcMYB114b* expression was different between 'OB' and 'CH', the promoter sequences of *RcMYB114b* in 'OB' and 'CH' were identified based on the haplotype genomes and amplified, and we found that there was only one promoter genotype of *RcMYB114b* in 'CH', named *RcMYB114bprolong*, while in 'OB' there were two genotypes. One was identical to that in 'CH' and the other had multiple SNPs, 6 indels and 2 large deletions, and was named *RcMYB114bproshort* (Fig. S17). Transient promoter activity assay was performed and the results showed that *RcMYB114bprolong* was more active than *RcMYB114bproshort* (Fig. S18), which was the reason for the different expression levels of *RcMYB114b* in 'CH' and 'OB'.

### Multi-TFs affect the response of rose petal number and size to high temperature

During the cultivation process, we found that the petal number of 'CH' plants was obviously affected by environmental factors. When 'CH' plants were treated at high (35 / 30 °C) and normal temperatures (25 / 20 °C), it was very interesting to see that the double-flower phenotype (more than 25 petals) of 'CH' plants observed at normal temperatures shifted to a single-flower phenotype (~5 petals), with smaller petals, at high temperatures (Fig. 5a-c). In 'OB', however, there was no significant reduction in petals (Fig. S19). While, stamen and carpel number were affected both in 'CH' and 'OB' under high temperature treatment.

In addition, we observed a significant decrease in the number of abaxial epidermis of petals under high temperature treatment (25.33 ± 1.67 to 20.33 ± 2.23), accompanied by a significant increase in petal cell area (107.23 ± 19.44 µm² to 157.82 ± 28.99 µm²) at FBGP, potentially contributing to the reduction in petal size (Fig. 5c). Then, we collected samples from the four stages for RNA sequencing. Through weighted correlation network analysis (WGCNA) (Fig. S20, 21), we identified four modules (Blue, Darkgrey, Lavenderblush3 and Darkseagreen) potentially associated with petal and stamen primordium development, which were further analyzed using Gene Ontology (GO) analysis (Fig. S22). Among the genes of the modules, the expression of *RcANT1* (homologous to *AIN-TEGUMENTA* in *Arabidopsis*), reported to regulate cell proliferation (Krizek 1999, Delgado-Benarroch
<figure><img src="image_14.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>d</figcaption></figure>

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

<figure><img src="image_19.png" /><figcaption>f</figcaption></figure>

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

**Fig. 5** Changes in the number and size of petals in *R. chinensis* 'CH' under high temperature treatment. **a** Flower morphology at different stages under high temperature (downer) and normal temperature (upper) treatment. **b** Statistics of the number and size of floral organs under different temperature treatments. Stars (***) mean significant differences (*P* < 0.001) analyzed by two-way ANOVA. **c** The cell density middle of petals under high temperature (right) and normal temperature (left) treatment. **d–e** The expression of *RcANT1*, *RcDA1*, *RcAG1*, *RcAG2* and *RcSVP1*. Stars mean significant differences (*** represents *P* < 0.001, ** represents *P* < 0.01, * represents *P* < 0.05) analyzed by T-test. **f** The coverage of *RcCAP2LWT* exon 10th with miR172 binding site under high temperature treatment. Stars mean significant differences (* represents *P* < 0.05) analyzed by T-test. **g** The dual-LUC assays. Red represents higher signal intensity and blue represents lower signal intensity. The bars indicate the SD. Stars (**) mean significant differences (*P* < 0.001) analyzed by one-way ANOVA.
et al. 2009) was downregulated after high-temperature treatment in FBGP, while the expression of the *RcDA1* gene (homologous to *DA* in *Arabidopsis*, which restricted the duration of the cell proliferation phase of floral organ growth (Li et al. 2008) was upregulated (Fig. 5d). The changes in *RcANT1* and *RcDA1* expression may contribute to the cell proliferation of petals in response to high temperature in rose.

Furthermore, we conducted a comparative analysis of DEGs between the two key cell proliferation groups (25 °C ISPS vs HBPS and 35 °C ISPS vs HBPS). There were 842 common DEGs in the two groups containing 13 MADS box genes related to flower development and 1 class A gene *RcAP2L* (Fig. S23 and Table S12). Among them, the expression of *RcAG1* increased with the formation of flower organ primordia and was higher in ISPS and HBPS under high-temperature treatment than under normal temperatures. *RcAG2* also had a similar expression pattern, but its expression level was lower than *RcAG1* (Fig. 5e). The expression of *RcSVP1* gene was significantly down-regulated. Dual-luciferase assays indicated that *RcSVP1* involved in inhibiting the expression of *RcAG1* and affect petal numbers under high temperature (Fig. 5g). However, it indicated that *RcSVP1* could not directly bind to the promoter of *RcAG1* (Fig. S24). The two *RcAP2L* alleles of ‘CH’ were obtained, one contained about 10 kb TE insertion, while the other did not. Both have miR172 binding sites in gene (Fig. S25). Based on transcriptome data, five splicing variants of *RcAP2L* in CH were obtained (Fig. S26). Notably, the level of reads coverage of miR172 binding site on exon 10th under high temperature environment was significantly higher than that under normal temperature at HBPS stage (Fig. 5f), which may result in stronger inhibition *RcAP2L* translation in high-temperature environments by miR172.

## Discussion

### Haplotype genomic information provided a new insight to trait diversity of R. chinensis

To maintain consistent diversity/quality of ornamental traits, rose is usually propagated vegetatively to maintain a state of high genomic heterozygosity. At present, there are some excellent germplasm resources for rose cultivars. However, if we seek to further improve these germplasms through hybridization, the original excellent traits may be lost, which severely limits the improvement of rose. Previous studies have found that molecular breeding using genomic information can improve breeding efficiency (Li and Zhang 2013, Liu et al. 2022). Especially, the haplotype genome is obtained to fully preserve the integrity of genetic information and greatly accelerate the breeding process such as in potato (Zhou et al. 2020b) and Bletilla

*striata* (Jiang et al. 2022). In this study, two haplotype-resolved genomes of *R. chinensis* 'CH' were first assembled. The genetic variation and differential expression of alleles between the two haplotypes were evaluated. All of these results provide guidance for exploring the regulatory mechanisms of important traits and assessing haplotype diversity to accelerate the improvement of rose in the future. To date, the accuracy phasing and accuracy assembly of highly heterozygous species have still posed significant challenges with current sequencing technologies. In the future, a 'CH' genetic map will be useful to validate the phasing. With advancements in sequencing technology, and the update of assembly software such as hifiasm, these advancements pave the way possibilities in assembling the genomes of highly heterozygous flowers or polyploid species, enabling the acquisition of more comprehensive and precise typing genome information.

Allelic variation confers the imbalanced expression of alleles, and contribute to trait diversity. For example, in apple, the presence of the TE insertion in the allele of *MYB* was positively associated with ASE and petal color (Tian et al. 2022). In Lettuce, the alleles of *Red Lettuce Leaves 1 (RLL1)* gene, a *bHLH* transcription factor, with 5 bp base deletion or not, regulated the leaf colour red or green (Su et al. 2020). In our study, we propose that *RcMYB114b* is an anthocyanin activator that is similar to those found in many plants, such as citrus, sweet cherry, and snapdragon (Schwinn et al. 2006, Jin et al. 2016, Huang et al. 2018). But it's interesting that in rose *RcMYB114b* was found to control the petal colour change from pink to red due to variations in promoter sequences. It was further proven that the allelic sequences may have different functions, so more comprehensive sets of clean haplotype genome information are needed to elucidate the molecular mechanism of important ornamental traits.

### Complex and diverse mechanisms for regulating flower color and pattern of roses

Anthocyanins are important metabolites that regulate the formation of petal color. Previous studies have reported that S6 subgroup of MYB involved in anthocyanin accumulation, contained 5 genes (*RcMYB113a*, *RcMYB113b*, *RcMYB113c*, *RcMYB114a* and *RcMYB114b*) (Li et al., 2021). Three of them (*MYB113a*, *RcMYB113c* called *RcMYB114* in previous reports, and *RcMYB114a*), have been reported (Lin-Wang et al., 2010; Li et al. 2022, Yan et al. 2023). In our study, another new gene *RcMYB114b* was identified. These reported showed that although these five genes belong to the same subfamily, they not only regulate different metabolite accumulation, but also have different expression sites/periods. For example, *RcMYB114a* (Yan et al. 2023) and *RcMYB114b* promoted cyanidin-3-glucoside accumulation. While *RcMYB113c* promoted
cyanidin-3-O-sophoroside accumulation (Li et al. 2022). In contrast, another homolog of *RcMYB113a*, *FaMYB10*, has been reported in strawberry, and overexpression of *FaMYB10* promoted the accumulation of cyanidin-glucoside, pelargonidin-glucoside and pelargonidin (Lin-Wang et al., 2010). In addition, *RcMYB114a* effected the accumulation of anthocyanin in response to light/dark, while *RcMYB114b* cannot. Interestingly, their expression pattern also was significantly different. The expression of *RcMYB114a* was up-regulated at the flower in early opening stage and highest at the fully open flower stage, whereas the expression of *RcMYB114b* was highest at the flower in early opening stage and down-regulated at the fully open flower stage. Different from above, *RcMYB113c* was largely unexpressed (Yan et al. 2023). In addition, different expression patterns of homologous genes were also found in strawberry, which may affect fruit color (Lin-Wang et al., 2010). All these founds indicate that the S6 subgroup genes may have been newly functionalized to regulate the accumulation of different substances, and subfunctionalized to regulate the same substances accumulation, but at different sites and periods. The discovery of a new gene *RcMYB114b* in our study provides new insights into the functional differentiation of the S6 subgroup genes in rose.

High temperatures have an important impact on flower yields and morphology of rose. In our study, under high-temperature treatment, we found that not only the number of petals but also petal size changed in the proliferation stage of petal cells. Previous studies not only showed *DA1* (Li et al. 2008, Gao et al. 2021), *ANT1* (Krizek 1999, Randall et al. 2015), and Class -A (*AP2* and *AP1*), -B (*AP3* and *PISTILLATA*), -C (AG) genes in volved in petal development, but also showed gene epigenetic regulation play important roles in petal number changes (Law and Jacobsen 2010, Jones 2012). For example, reduced expression of *AG* promoted the formation of double petal flowers, and increased methylation of the *RhAG* promoter at low temperature (15 / 5 °C) has been linked to a higher number of petals in roses (Ma et al. 2015). In addition, the transcription may also regulate the petal number (François et al. 2018). Consider the miR172 binding site was the key factor effecting petal number rather than one isoform's expression. In our results, we detected the level of reads coverage of miR172 binding site on exon 10th was higher under high temperature environment. As previous study found that the more expression of miR172 binding site, it could more decrease the *RcAP2L* protein to effect petal number (Chen 2004). And interestingly, under high temperature treatment, ‘CH’ and ‘OB’ showed different phenotypes which containing different genotypes. In a nutshell, all

these indicated that flower development under high-temperature treatment involves a very complex regulatory network, and more studies are required to clarify this in the future.

## Materials and methods

### Plant materials

*R. chinensis* 'CH' and 'OB' are widely recognized old rose cultivars. Varieties selected from these two cultivars for this study were growing in the same experimental field under natural conditions at the Flower Garden of Huazhong agricultural University. The leaves of 'CH' were subjected to DNA extraction for genome sequencing, and extracts of its young stems, flower buds, leaves and young roots were used for genome annotation.

To explore the formation of petal colour, samples of petals in three stages (noncoloured petals, colouring petals, coloured petals) and two petal parts (the red region and the white heart region of coloured petals) were collected for RNA-seq. A total of five samples were collected from 'CH' for analysis of the anthocyanin regulatory pathway. The different petal stages selected were consistent with the published petal sampling stages of 'OB' (Han et al. 2017). The petals of six flowers at each developmental stage were collected as one biological replicate, with a total of three biological replicates for RNA-seq, anthocyanin extraction and liquid chromatography-mass spectrometry (LC-MS) analysis.

To analyze the effect of high-temperature treatment on rose flower development, sixty individual 'CH' plants were pruned and then randomly placed into two incubators with two different parameter settings: (1) 16 / 8 h day/night at 35 °C / 30 °C and (2) 16 / 8 h day/night at 25 °C / 20 °C. After 50 d, the four different stages of flower buds: the vegetative meristem stage (DBO); the initiation stages of petals/petal-like structures and stamens/stamen-like structures (ISPS); the stage in which the hypanthium starts to sink below the perianth and stamens (HBPS); and the stage of flower buds with young noncoloured petals (FBGP), were observed individually and were collected under a microscope (SOPTOP SZN, Yuyao, China) in three biological replicates for RNA-seq.

The numbers of sepals, petals, stamens and carpels were counted in at least thirty flowers under heat treatment and control conditions. Statistical analysis was performed using two-way ANOVA along with the least significant difference test using GraphPad Prism version 9.0, and *P* < 0.05 was considered significant. Transgenic callus of the same line was snap frozen in liquid nitrogen and stored at -80 °C in the refrigerator until transcriptionome sequencing, quantitative detection, pigment extraction and LC-MS.
## Microscopic examination and cell counting

Cell photography and counting were performed as described by previous reports (Dewitte et al. 2007) with some modifications. In short, petal samples at stage of flower buds with young noncoloured petals (FBGP) in the above two incubators were taken as slices at 50% of the petal length from the petal top. The slices were embedded in 3% agarose gel (Biosharp, Beijing, China) and 60 µm sections were cut on a semi-automatic Leica VT1200 vibrating blade microtome (Weztlar, Germany). The abaxial epidermis cells of the slices were photographed using a Novel digital microscope (NLCD500, Nanjing, China). Fifteen flowers were used in each treatment. Numbers of cells were counted in a visual field of 50×50 µm². The area of cell was measured in ImageJ software. Statistical analysis was performed using T-test by GraphPad Prism version 9.0, and *P*<0.05 was considered significant.

## Genome sequencing, assembly and anchoring

MGI data sequencing. The healthy young leaves of *R. chinensis* 'CH' were quickly frozen in liquid nitrogen and stored in a -80 °C freezer for MGI and HiFi sequencing. A library with an insert size of 150 bp was produced using MGI to survey the genome.

HiFi data sequencing. Genomic DNA was prepared by purification with a QIAGEN® Genomic kit (QIA-GEN) according to the manufacturer's instructions. More than 10 µg of genomic DNA was used to construct ~ 20-kb-target-size Sequel SMRTbell libraries in circular consensus sequencing (CCS) mode for PacBio HiFi sequencing.

Hi-C data sequencing. A Hi-C library was constructed from young leaves of 'CH', and sequenced via the Illumina NovaSeq platform. Low-quality sequences and adaptor sequences were filtered out, and the unique mapped paired-end reads were extracted by mapping the clean paired-end reads to the draft assembled sequence.

To estimate the size of the genome, the 21-mer frequency distribution was analyzed using Jellyfish (v2.3.0) (https://github.com/jamesturk/jellyfish). Hifiasm (v0.16.1-r375) software (Cheng et al. 2021) was used to assemble the genome of two haplotypes of *R. chinensis* 'CH' using HiFi reads by diploid phased assembly using the parameter "hifiasm -o asm -h1 HIC_1.fq.gz -h2 HIC_2.fq.gz genome.fastq.gz". The resulting contigs were aligned to rose plastid genomes using BLASTn (v2.9.0), and the contigs in which more than 95% identity with plastid genome were filtered. Then, the clean data were used for two rounds of error correction by Pilon (v1.23). Hi-C reads were mapped to the assembled scaffolds using Juicer (v1.6)

(Durand et al. 2016) with the default parameters, and 3D-DNA (Dudchenko et al. 2017) was used to cluster and order them into 14 chromosomes with parameters "-r 0". After manual adjustment by Juicebox v1.11.08, we obtained the two haplotype genomes (hA and hB), which were fully resolved at the chromosomal level. We used the following strategies to assess the quality of the genomes: (1) The BUSCO (v4.0.6) database (eudicotyledons_odb10 database) was used to assess the continuity of the assembly with parameters "-m geno". (2) The heterozygous regions were detected by KAT software and the MGI and HiFi reads were mapped back to the genome. (3) The mapping rate as well as genome coverage of sequencing reads were assessed using samtools v1.7. (4) The Quality Value score was used to evaluate this genome by Merqury (Rhie et al. 2020). (5) Whole-genome with fourteen chromosomes Hi-C contact heatmaps were plotted using ALLHIC. (6) The chromosomes of 'CH' were aligned to the OB reference genome (Raymond et al. 2018) using MUMmer v3.1 with default parameter ('OB' genome as reference genome and 'CH' genome as query genome). (7) The switch errors were evaluated by the calc_switchErr (https://github.com/tangerzhang/calc_switchErr).

## Annotation of repetitive sequences and protein-coding genes

TEs in 'CH' were identified with the an extensive de novo TE annotator pipeline (Ou et al. 2019) with parameter "-anno 1 -force 1 -debug 1 -sensitive 1", which combines the widely used RepeatModeler package with Repbase and seven other structure-based methods for de novo LTR identification.

Protein-coding genes for each haplotype genome were predicted using EVidenceModeller (Haas et al. 2008) by integrating three strategies: ab initio gene predictions; RNA-seq read mapping; and protein homology. Ab initio gene predictions were generated by AUGUSTUS v3.3.3 (Stanke et al. 2006), SNAP v2013-02-16 (Korf 2004), and GlimmerHMM v3.0.4 (Majoros et al. 2004). For homology-based prediction, the protein sequences from eight species, including *Arabidopsis thaliana*, *Oryza sativa*, *Vitis vinifera* and the Rosaceae members *R. chinensis* 'OB', *R. chinensis* 'OB' _DH, *P. persica*. v2, *Pyrus betulifolia*, and *F. vesca* v4 (https://www.rosaceae.org/tools/jbrowse), were integrated by Exonerate v2.2.0 (http://www.ebi.ac.uk/~guy/exonerate/) and AUGUSTUS. For transcriptome-based prediction, RNA-seq data were mapped using HiSAT2 v2.1.0 (Kim et al. 2015) and then assembled with PASA v 2.4.1 (Haas et al. 2003). The final annotation of protein completeness was evaluated with BUSCO.
All predicted proteins were aligned against the UniProt SwissProt and NR databases. The functions of the best-matched proteins were assigned to the predicted proteins. Functional annotation was performed using InterProScan v5.21 (Jones et al. 2014). GO terms were assigned according to the InterPro classification.

## The overall view of the comparisons and diversity analysis between the two haplotypes

To identify the variation between the two genomes, we used SyRI coordinates to align them pairwise and identify variations such as candidate inversions, intra- and inter-chromosomal translocations (Goel et al. 2019). To identify the homologous regions between two haplotypes, the syntenic blocks were constructed based on well-aligned genes (CDSs) using MCScanX package. We screened for allelic genes based on the following criteria: (1) paired regions must be located on homologous chromosomes within syntenic blocks; (2) The gene has another best homologous gene on another haplotype and each of the gene pairs should be a reciprocal blast best hit; (3) The coding sequence alignment needed to exhibit at least one SNP. Genes that satisfied these specific criteria were classified as allelic genes. If genes had identical coding sequences within syntenic blocks between the two haplotypes, they were labeled as “single allele”.

## Allelic genes expression analysis in different tissues

We performed allelic gene expression analysis using the allele-specific mapping of Kallisto which was used in the comparison of homologous expression in polyploid wheat (Ramírez-González et al. 2018) and potato (Zhou et al. 2020b), and we used the same software to obtain the expression levels, in transcripts per million (TPM), of genes on both haplotypes based on RNA-seq data from different tissues (roots, stems, leaves and flowers). These genes with average TPM (value > 1) for any of tissues were selected as expressed genes. We then tested for expression differences between allelic genes using the T-test with a significance threshold of *P* < 0.01.

## Phylogenetic evolutionary analysis

The species phylogenetic tree of *N. colorata*, *O. sativa*, *A. thaliana*, *F. vesca*, *M. domestica*, *P. communis*, *P. persica*, *R. occidentalis*, *R. chinensis* 'OB' and *R. rugosa* was constructed based on single-copy homologous genes with OrthoFinder (Emms and Kelly 2015). The hA and hB genomes of *R. chinensis* 'CH' was separately selected for the construction of the species phylogenetic tree. The time calibrations from the TimeTree database (http://www.timetree.org/) were used to estimate divergence times. The

single-copy gene sequences were aligned with MAFFT (Katoh and Standley 2013). They were subsequently used to construct a maximum likelihood tree, and the 95% confidence intervals of divergence times were obtained with MCMCtree in the PAML package (Yang 2007).

## Transcriptome analysis

Total RNA was isolated using the EASYspin Plant RNA extraction kit (Aidlab, Beijing, China) following the manufacturer's protocol. All libraries (three biological replicates for every sample) were sequenced on the BGI platform. Transcriptome data for nine samples of different petal stages of 'OB' were downloaded from NCBI (PRJNA351281) (Han et al. 2017). After quality control and the removal of low-quality and adapter-containing reads with Trimmomatic-0.36, the clean data were aligned to the hB genome of 'CH' and 'OB' using HiSAT2 v2.1.0 (Kim et al. 2015). For accurate quantification of homologous gene expression, only the unique mapping reads were kept for further analysis. Transcript levels were normalized using the fragments per kilobase per million mapped reads (FPKM) method and coverages of exon were calculated with StringTie v1.3.5 (Pertea et al. 2016).

Differential expression analysis of two conditions/groups was performed using the DESeq2 (Love et al. 2014). The resulting *P*-values were adjusted using the Benjamini and Hochberg's approach for controlling the false discovery rate. Genes with an adjusted *P*-value < 0.05 (log2FoldChange > 2 | log2FoldChange < -2) found by DESeq2 were assigned as differentially expressed. The Venn diagram and kyoto encyclopedia of genes and genomes (KEGG) enrichment diagram were performed using R software.

## Volatile organic compound collection and analysis

Headspace solid-phase microextraction–gas chromatography–mass spectrometry (HS-SPME–GC–MS) technology was used to extract and identify volatile organic compounds. Scent collection was performed in this study according to previous study (Zhao et al. 2016). Freshly picked rose flowers samples of approximately 0.7 g were weighted and then added to a collection bottle. Five independent replicates were performed at the “coloured” stage of *R. chinensis* 'OB' and 'CH' petals. After water bath and extraction, the fibre was inserted into a GC injection port and then subjected to GC–MS analysis. For chromatographic separation, He (99.999%) was used as the carrier gas. The starting temperature of the column was 40 °C. After a 4 min hold, the temperature was increased to 160 °C at 2 °C/min. Subsequently, the temperature was increased to 250 °C at 15 °C/min with a final holding time of 5 min. The temperature of the transfer line to the mass spectrometer was set at 220 °C.
The mass range was 50 to 650 m/z, and full-scanning mode was employed. Xcalibur software (Thermo, Waltham, US) was used for the analysis of mass spectra. M/Z values were matched with those of the National Institute of Standards and Technology (NIST) library. The substances were classified according to the published classification (Zhou et al. 2020a), and the peak area of each type of substance was quantified to compare the proportions of the contents of each substance type.

## Genetic transformation in rose calli

Calli of the 'OB' rose cultivar were used for genetic transformation, and the CDS of *RcMYB114b* was amplified and cloned into the pRI 101-AN vector using corresponding primer sequences (Table S13). The constructs of 35S:*RcMYB114b* and vector-control were transformed into *Agrobacterium tumefaciens* GV3101 and then transformed into rose calli as previously reported (Vergne et al. 2010) but with minor changes. Briefly, embryos were bulk wounded with sterile sand. The cells were co-cultivated for 1 h at 100 rpm at room temperature. 2 d later, the embryos were washed several times. Approximately 10 d after transformation, embryos were transferred to embryogenic callus maintenance medium supplemented with kanamycin (60 mg·L⁻¹) and Timentin (150 mg·l⁻¹) and subcultured every 3–4 weeks.

## Transient expression assay in apple fruit

A transient gene expression assay was employed for the overexpression of *RcMYB114b* in apple fruit. The overexpression vector was obtained by cloning the sequence of *RcMYB114b* from the cDNA of 'CH' using primer sequences (Table S13), and fused into the pGreenII 62-SK vector. The empty SK and RcMYB114b-SK constructs were then introduced to *A. tumefaciens* strain GV3101 cells. The injections were carried out using a 1 ml needleless syringe as described previously (Li et al. 2012). The injected apples were kept in the dark at 23 °C overnight and were subsequently treated with constant light (100 µmol·m⁻²·s⁻¹) for colour development.

## qRT-PCR analysis

Total RNA from the 'CH' and 'OB' flower samples, the apple and rose calli of control and transgenic lines was extracted for qRT-PCR. TRUEscript RT MasterMix (Aidlab, Beijing, China) was employed for reverse transcription using 1 µg of RNA. The generated cDNA was a template for qRT-PCR assays with gene-specific rose (Table S13) and apple (Han et al. 2010, An et al. 2018, Kang et al. 2022) primers using Taq Pro Universal SYBR aPCR Master Mix (Vazume, Nanjing, China). The *RcACTIN* (Dubois et al.

2012) and *MdACTIN* (Li et al. 2012) genes were used as internal quantitative controls to normalize samples. The relative expression values were calculated using the comparative CT (2⁻ΔΔCT) method. Statistical analysis was performed via one-way ANOVA along with the least significant difference test using GraphPad Prism version 9.0, and *P*<0.05 was considered significant.

## Anthocyanin extraction and LC-MS analysis

The anthocyanin metabolites of rose 'CH' and 'OB' petals and transgenic calli were measured. The samples were powdered in liquid nitrogen. Then, a 0.2 g sample was weighed, placed in an ultrasonicator bath for 30 min with 2 ml of acidic methanol solution (70:27:2:1; v/v, CH₃OH: H₂O: HCOOH: CF₃COOH), and maintained in the dark at 4 °C overnight) (Gao and Mazza 1994, Liu et al. 2016). The extract was centrifuged (12,000 rpm, 10 min, 20 °C), and the supernatant was filtered through 0.22 µm syringe filters (Coolwind, Guangzhou, China) into a 2 ml amber screw autosampler vial (Jiedao, Jiangsu, China). LC-MS analysis was carried out on an Accurate-Mass Q-TOF LC/MS 6520 (Agilent, California, USA). A 150 mm×4.6 mm, i.d., 4.6 µm, HC-C18(2) column (Agilent, California, USA) was used. Referring to a reported method (Liu et al. 2016), 0.5% aqueous formic acid (A) and acetonitrile (B) were used as mobile phases. The gradients, column temperature, injection volume, and flow rate were programmed as described in a previous report. The ionization of anthocyanin was achieved with an ESI source in positive mode, and the parameters were set followed a previous report (Wan et al. 2018). The anthocyanin chromatograms were extracted at 540 nm.

## Dual-luciferase assays

The promoter fragments of *RcAG1* (1,445 bp upstream of start codon) were inserted into the pGreenII 0800-LUC vector to construct *RcAG1pro:LUC*. The open reading frame (ORF) of *RcSVP1* was inserted into the pGreenII 62-SK vector to construct 35Spro:*RcSVP1* (Table S13). The dual-luciferase assay followed the protocol described previously in *Nicotiana benthamiana* (Hellens et al. 2005, Sparkes et al. 2006). Firefly luciferase and Renilla luciferase were assayed 3 d after injection using a Dual-Luciferase Reporter Assay Kit (Promega, Madison, USA). Three independent experiments with a minimum of six replicates were performed to verify the results. Photographs were taken by using an IVIS Spectrum in vivo imaging system (PerkinElmer, Waltham, USA).

Remaining method descriptions are included to Supplementary Notes.
## Abbreviations

<table><tr><td>AP2L</td><td>APETALA2-like</td></tr><tr><td>ASE</td><td>Allele-specific expression</td></tr><tr><td>BUSCO</td><td>Benchmarking universal single-copy ortholog</td></tr><tr><td>CCS</td><td>Circular consensus sequencing</td></tr><tr><td>CDS</td><td>Coding sequences</td></tr><tr><td>CHS</td><td>Chalcone synthase</td></tr><tr><td>DEGs</td><td>Differentially expressed genes</td></tr><tr><td>DFR</td><td>Dihydroflavonol 4-reductase</td></tr><tr><td>ERF</td><td>Ethylene response factor</td></tr><tr><td>F3H</td><td>Flavanone 3-hydroxylase</td></tr><tr><td>F3'H</td><td>Flavonoid 3'-hydroxylase</td></tr><tr><td>FLS</td><td>Flavonol synthase</td></tr><tr><td>GO</td><td>Gene ontology</td></tr><tr><td>LAI</td><td>The long terminal repeat (LTR)-assembly index</td></tr><tr><td>LAR</td><td>Leucoanthocyanidin reductase</td></tr><tr><td>LC-MS</td><td>Liquid chromatography-mass spectrometry</td></tr><tr><td>LTR-RTs</td><td>Long terminal repeat retrotransposons</td></tr><tr><td>miR172</td><td>MicroRNA172</td></tr><tr><td>ORF</td><td>Open reading frame</td></tr><tr><td>PAL</td><td>Phenylalanine ammonia-lyase</td></tr><tr><td>qRT-PCR</td><td>Real-time quantitative PCR</td></tr><tr><td>RLL1</td><td>Red Lettuce Leaves 1</td></tr><tr><td>SNP</td><td>Single nucleotide polymorphism</td></tr><tr><td>SyRI</td><td>Synteny and rearrangement identifier</td></tr><tr><td>TEs</td><td>Transposable elements</td></tr><tr><td>TFs</td><td>Transcription factors</td></tr><tr><td>TPM</td><td>Transcript per million</td></tr><tr><td>GUS</td><td>β-Glucuronidase</td></tr><tr><td>WGCNA</td><td>Weighted correlation network analysis</td></tr></table>

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s43897-024-00088-1.

**Additional file 1: Fig. S1.** K-mer analysis for estimating the genome size of R. chinensis 'CH' indicating high heterozygosity. **Fig. S2.** The plot depicts the K-mer spectra of the haplotype assembly. **Fig. S3.** Hi-C interaction heatmaps of the two haplotypes (hap A and hap B) in R. chinensis 'CH'. **Fig. S4.** Sequence depth and GC content for the two haplotypes in R. chinensis 'CH'. **Fig. S5.** Dot plots with respect to the 'OB' genome reference for the two haplotypes assemblies of 'CH' respectively. **Fig. S6.** Collinearity of genomes from 'CH' and 'OB' by SyRI. **Fig. S7.** Expression of alleles in different tissues of R. chinensis 'CH'. **Fig. S8.** Statistics of ASE gene in different tissues. **Fig. S9.** Phylogenetic tree for 'CH' hA and ten other eudicot species contained three outgroups. **Fig. S10.** Gene Ontology (GO) enrichment analysis of expanded gene families in 'CH' hA (a) and 'OB' (b). **Fig. S11.** Expression of expanded family genes which Gene Ontology associated with terpene synthesis. **Fig. S12.** The differentially expressed genes between noncoloured petals and colouring petals of 'CH' hB. **Fig. S13.** Phylogenetic tree analysis of MYB genes from R. chinensis 'CH', 'OB' and Arabidopsis. **Fig. S14.** The differential genes' expressions form transcriptome was validated by qRT-PCR. **Fig. S15.** The structural variation of alleles was observed between two haplotypes and two materials by SyRI. **Fig. S16.** The differential expression gene between transgenic roses callus and control roses callus. **Fig. S17.** Promoter sequence of RcMYB114b in R. chinensis 'OB' and 'CH'. **Fig. S18.** GUS staining and GUS enzyme detection. **Fig. S19.** Morphological analysis of 'OB' flowers in high temperature treatment and control group. **Fig. S20.** Hierarchical cluster dendrogram showing co-expressed modules identified by weighted gene co-expression network analysis for the rose RNA-seq data. **Fig. S21.** Module-stage association analysis. **Fig. S22.** Gene Ontology (GO) enrichment analysis of four modules genes related to flower organ primordium formation. **Fig. S23.** The DEGs in two comparisons (25°C ISPS vs_HBPS and 35°C ISPS vs_HBPS). **Fig. S24.** Yeast one-hybrid (Y1H) assay and electrophoretic mobility shift assay (EMSA). **Fig. S25.** Genotype information of RcAP2L. **Fig. S26.** RcAP2L variable spliceosome sequence alignment analysis.

**Additional file 2: Table S1.** Summary of genome sequencing and read out for R. chinensis 'CH'. **Table S2.** Survey statistic results of R. chinensis 'CH'. **Table S3.** Summary of Hi-C reads mapping to the R. chinensis. **Table S4.** Result of the evaluation of the R. chinensis genome by BUSCO. **Table S5.** Genome sequence data comparison statistics. **Table S6.** Result of the evaluation of the R. chinensis coding genes by BUSCO. **Table S7.** Gene numbers and features of R. chinensis 'CH'. **Table S8.** The statistical results of gene function annotation of R. chinensis 'CH'. **Table S9.** Summary of repeat contents in R. chinensis 'CH' genome. **Table S10.** Structural variation between two haplotypes of R. chinensis 'CH'. **Table S11.** SNP variation between two haplotypes of R. chinensis 'CH'. **Table S12.** Genes related to flower development in DEGs **Table S13.** List of primers used.

## Acknowledgements

We thank members of the Wu laboratory and Fu laboratory. We thank Prof. Kang Chunying (Huazhong Agricultural University) for providing the pRI 101-AN vector.

## Authors' contributions

The ideas for the paper were conceived by X.F., Zh.W. and X.Z.. X.Z. and D.P. performed the assembly and annotation of the genome. X.Z. and L.L. conducted phylogenetic analysis. Q.W. and H.G. performed construction of plastid, transgenic experiments and the qRT-PCR. Q.W. participated in metabolite profiling. K.L. planted the plant materials. X.Z. and Q.W. drafted the manuscript with the help of X.F., Z.W., M.B and Mo.B. All authors read and approved the final manuscript.

## Funding

This work was supported by funding from the National key research and development program of China (2019YFD1000403 & 2021YFD1200205-03). The Scientific Research Foundation for Principle Investigator, Kunpeng Institute of Modern Agriculture at Foshan (KIMA-QD2022004), the Funding of Major Scientific Research Tasks, Kunpeng Institute of Modern Agriculture at Foshan (KIMA-ZDKY2022004), this work was also supported in part by the Chinese Academy of Agricultural Sciences Elite Youth Program to Zhiqiang Wu.

## Availability of data and materials

All data supporting the results of this study are included in the manuscript and its additional files. Genome assembly and annotations were deposited in the NCBI BioProject under accessions PRJNA932466.

## Declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors have no conflicts of interest to declare.

Received: 27 December 2023 Accepted: 19 February 2024

Published online: 16 April 2024

## References

* An JP, An XH, Yao JF, Wang XN, You CX, Wang XF, et al. BTB protein MdBT2 inhibits anthocyanin and proanthocyanidin biosynthesis by triggering MdMYB9 degradation in apple. Tree Physiol. 2018;38:1578-87.

* Bendahmane M, Dubois A, Raymond O, Bris ML. Genetics and genomics of flower initiation and development in roses. J Exp Bot. 2013;64:847-57.

* Chen X. A MicroRNA as a Translational Repressor of APETALA2 in Arabidopsis Flower Development. Science. 2004;303:2022-5.

* Chen F, Su LY, Hu SY, Xue JY, Liu H, Liu GH, 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.
Cheng HY, Concepcion GT, Feng XW, Zhang HW, Li H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021;18:170–5.

Conart C, Saclier N, Foucher F, Goubert C, Rius-Bony A, Paramita SN, et al. Duplication and Specialization of *NUDX1* in Rosaceae Led to Geraniol Production in Rose Petals. Mol Biol Evol. 2022;39:msac002.

Delgado-Benarroch L, Causier B, Weiss J, Egea-Cortines M. FORMOSA controls cell division and expansion during floral development in *Antirrhinum majus*. Planta. 2009;229:1219–29.

Dewitte W, Scofield S, Alcasabas AA, Maughan SC, Menges M, Braun N, et al. *Arabidopsis* CYCD3 D-type cyclins link cell proliferation and endocycles and are rate-limiting for cytokinin responses. Proc Natl Acad Sci USA. 2007;104:14537–42.

Dubois A, Raymond O, Maene M, Baudino S, Langlade NB, Boltz V, et al. Tinkering with the C-Function: A Molecular Frame for the Selection of Double Flowers in Cultivated Roses. PLoS ONE. 2010;5:e9288.

Dubois A, Carrere S, Raymond O, Pouvreau B, Cottret L, Roccia A, et al. Transcriptome database resource and gene expression atlas for the rose. BMC Genom. 2012;13:638.

Dudchenko O, Batra SS, Omer AD, Nyquist SK, Hoeger M, Durand NC, et al. De novo assembly of the *Aedes aegypti* genome using Hi-C yields chromosome-length scaffolds. Science. 2017;356:92–5.

Durand NC, Shamim MS, Machol I, Rao SS, Huntley MH, Lander ES, et al. Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments. Cell Syst. 2016;3:95–8.

Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16:157.

François L, Verdenaud M, Fu XP, Ruleman D, Dubois A, Vandenbussche M, et al. A miR172 target-deficient AP2-like gene correlates with the double flower phenotype in roses. Sci Rep. 2018;8:e12912.

Gao L, Mazzza G. Quantitation and Distribution of Simple and Acylated Anthocyanins and Other Phenolics in Blueberries. J Food Sci. 1994;59:1057–9.

Gao Q, Zhang N, Wang WQ, Shen SY, Bai C, Song XJ. The ubiquitin-interacting motif-type ubiquitin receptor HDR3 interacts with and stabilizes the histone acetyltransferase GW6a to control the grain size in rice. Plant Cell. 2021;33:3331–47.

Gattolin S, Cirilli M, Pacheco I, Ciacciulli A, Linge CD, Mauroux JB, et al. Deletion of the miR172 target site in a TOE-type gene is a strong candidate variant for dominant double-flower trait in Rosaceae. Plant J. 2018;96:358–71.

Goel M, Sun H, Jiao WB, Schneeberger K. SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol. 2019;20:277.

Haas BJ, Salzberg SL, Zhu W, Pertea M, Allen JE, Orvis J, et al. Automated eukaryotic gene structure annotation using eEvidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol. 2008;9:R7.

Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK, Hannick LI, et al. Improving the *Arabidopsis* genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003;31:5654–66.

Han Y, Vimolmangkang S, Soria-Guerra RE, Rosales-Mendoza S, Zheng D, Lygin AV, et al. Ectopic Expression of Apple *F3'H* Genes Contributes to Anthocyanin Accumulation in the *Arabidopsis* *tt7* Mutant Grown Under Nitrogen Stress. Plant Physiol. 2010;153:806–20.

Han Y, Wan H, Cheng T, Wang J, Yang W, Pan H, et al. Comparative RNA-seq analysis of transcriptome dynamics during petal development in *Rosa chinensis*. Sci Rep. 2017;7:43382.

Hellens RP, Allan AC, Friel EN, Bolitho K, Grafton K, Templeton MD, et al. Transient expression vectors for functional genomics, quantification of promoter activity and RNA silencing in plants. Plant Methods. 2005;1:13.

Huang D, Wang X, Tang Z, Yuan Y, Xu Y, He J, et al. Subfunctionalization of the *Ruby2-Ruby1* gene cluster during the domestication of citrus. Nat Plants. 2018;4:930–41.

Huylenbroeck J. Ornamental Crops Preface. Handb Plant Breed. 2018;11:V-Vi. https://doi.org/10.1007/978-3-319-90698-0.

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

Jiang L, Lin MF, Wang H, Song H, Zhang L, Huang QY, et al. Haplotype-resolved genome assembly of *Bletilla striata* (Thunb.) Reichb.f. to elucidate medicinal value. Plant J. 2022;111:1340–53.

Jin W, Wang H, Li M, Wang J, Yang Y, Zhang X, et al. The R2R3 MYB transcription factor *PavMYB10.1* involves in anthocyanin biosynthesis and determines fruit skin colour in sweet cherry (*Prunus avium* L.). Plant Biotechnol J. 2016;14:2120–33.

Jones P, Binns D, Chang HY, Fraser M, Li WZ, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.

Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13:484–92.

Jung JH, Lee S, Yun J, Lee M, Park CM. The miR172 target TOE3 represses *AGAMOUS* expression during *Arabidopsis* floral patterning. Plant Sci. 2014;215:29–38.

Kang H, Zhang TT, Li YY, Lin-Wang K, Espley RV, Du YP, et al. The apple BTB protein MdBT2 positively regulates MdCOP1 abundance to repress anthocyanin biosynthesis. Plant Physiol. 2022;190:305–18.

Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.

Koes R, Verweij W, Quattrochio F. Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005;10:236–42.

Korf I. Gene finding in novel genomes. BMC Bioinform. 2004;5:59.

Krizek BA. Ectopic expression of *AINTEGUMENTA* in *Arabidopsis* plants results in increased growth of floral organs. Dev Genet. 1999;25:224–36.

Law JA, Jacobsen SE. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11:204–20.

Li M, Zhang H, Yang Y, Wang H, Xue Z, Fan Y, et al. Rosa1, a Transposable Element-Like Insertion, Produces Red Petal Coloration in Rose Through Altering RcMYB114 Transcription. Front Plant Sci. 2022;13:857684.

Li M, Zhang X, Lin S, Wu Q, Bao M, Fu X. Genome-wide identification and analysis of the R2R3-MYB transcription factor in Rosa chinensis. Mol Plant Breed. 2021;19:7018–29.

Lin-Wang K, Bolitho K, Grafton K, Kortstee A, Karunairetnam S, McGhie TK, et al. An R2R3 MYB transcription factor associated with regulation of the anthocyanin biosynthetic pathway in Rosaceae. BMC Plant Biol. 2010;10:50.

Li Y, Zheng L, Corke F, Smith C, Bevan MW. Control of final seed and organ size by the DA1 gene family in *Arabidopsis thaliana*. Genes Dev. 2008;22:1331–6.

Li YY, Mao K, Zhao C, Zhao XY, Zhang HL, Shu HR, et al. MdCOP1 ubiquitin E3 ligases interact with MdMYB1 to regulate light-induced anthocyanin biosynthesis and red fruit coloration in apple. Plant Physiol. 2012;160:1011–22.

Li ZK, Zhang F. Rice breeding in the post-genomics era: from concept to practice. Curr Opin Plant Biol. 2013;16:261–9.

Liu H, Sun M, Du D, Pan H, Cheng T, Wang J, et al. Whole-transcriptome analysis of differentially expressed genes in the ray florets and disc florets of *Chrysanthemum morifolium*. BMC Genom. 2016;17:398.

Liu WW, Zhang YL, He H, He GM, Deng XW. From hybrid genomes to heterotic trait output: Challenges and opportunities. Curr Opin Plant Biol. 2022;66:102193.

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

Ma N, Chen W, Fan T, Tian Y, Zhang S, Zeng D, et al. Low temperature-induced DNA hypermethylation attenuates expression of *RhAG*, an *AGAMOUS* homolog, and increases petal number in rose (*Rosa hybrida*). BMC Plant Biol. 2015;15:237.

Magnard JL, Roccia A, Caissard JC, Vergne P, Sun P, Hecquet R, et al. Biosynthesis of monoterpane scent compounds in roses. Science. 2015;349:81–3.

Majoros WH, Pertea M, Salzberg SL. igrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004;20:2878–9.

Ou SJ, Su WJ, Liao Y, Chougule K, Agda JRA, Hellinga AJ, et al. Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. Genome Biol. 2019;20:275.

Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT. StringTie and Ballgown Nat Protoc. 2016;11:1650–67.
Ramírez-González RH, Borrill P, Lang D, Harrington SA, Brinton J, Venturini L, et al. The transcriptional landscape of polyploid wheat. Science. 2018;361:earar6089.

Randall RS, Sornay E, Dewitte W, Murray JA. AINTEGUMENTA and the D-type cyclin CYCD3;1 independently contribute to petal size control in *Arabidopsis*: evidence for organ size compensation being an emergent rather than a determined property. J Exp Bot. 2015;66:3991–4000.

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

Rhie A, Walenz BP, Koren S, Phillippy AM. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol. 2020;21:245.

Scalliet G, Piola F, Douady CJ, Rety S, Raymond O, Baudino S, et al. Scent evolution in Chinese roses. Proc Natl Acad Sci USA. 2008;105:5927–32.

Schwinn K, Venail J, Shang Y, Mackay S, Alm V, Butelli E, et al. A small family of *MYB*-regulatory genes controls floral pigmentation intensity and patterning in the genus *Antirrhinum*. Plant Cell. 2006;18:831–51.

Seneviratne SI, X Zhang, M Adnan, W Badi, C Dereczynski, A Di Luca, et al. Weather and Climate Extreme Events in a Changing Climate. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. https://www.ipcc.ch/report/ar6/wg1/chapter/chapter-11/. Accessed 29 June 2023.

Sparkes IA, Runions J, Kearns A, Hawes C. Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants. Nat Protoc. 2006;1:2019–25.

Stanke M, Schöffmann O, Morgenstern B, Waack S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinform. 2006;7:62.

Su W, Tao R, Liu W, Yu C, Yue Z, He S, et al. Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication. Plant Biotechnol J. 2020;18:479–90.

Suzuki K, Tsuda S, Fukui Y, Fukuchi-Mizutani M, Yonekura-Sakakibara K, Tanaka Y, et al. Molecular characterization of rose flavonoid biosynthesis genes and their application in petunia. Biotechnol Biotec Eq. 2014;14:56–62.

Tian Y, Thrimawithana A, Ding T, Guo J, Gleave A, Chagné D, et al. Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.). Plant Biotechnol J. 2022;20:1285–97.

Vergne P, Maene M, Gabant G, Chauvet A, Debener T, Bendahmane M. Somatic embryogenesis and transformation of the diploid *Rosa chinensis* cv Old Blush. Plant Cell Tiss Organ Cult. 2010;100:73–81.

Wan H, Yu C, Han Y, Guo X, Ahmad S, Tang A, et al. Flavonols and Carotenoids in Yellow Petals of Rose Cultivar (*Rosa* 'Sun City'): A Possible Rich Source of Bioactive Compounds. J Agric Food Chem. 2018;66:4171–81.

Wang GL. A study on the history of Chinese Roses from ancient works and images. Acta Hortic. 2007;751:347–56.

Wang QJ, Zhang XN, Lin SN, Yang SZ, Yan XL, Bendahmane M, et al. Mapping a double flower phenotype-associated gene *DcAP2L* in *Dianthus chinensis*. J Exp Bot. 2020;71:1915–27.

Wollmann H, Mica E, Todesco M, Long JA, Weigel D. On reconciling the interactions between *APETALAZ*, *miR172* and *AGAMOUS* with the ABC model of flower development. Development. 2010;137:3633–42.

Yan Y, Zhao J, Lin S, Li M, Liu J, Raymond O, et al. Light-mediated anthocyanin biosynthesis in rose petals involves a balanced regulatory module comprising transcription factors RhHY5, RhMYB114a, and RhMYB3b. J Exp Bot. 2023;74:5783–804.

Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.

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

Zhang X, Lin S, Peng D, Wu Q, Liao X, Xiang K, et al. Integrated multi-omic data and analyses reveal the pathways underlying key ornamental traits in carnation flowers. Plant Biotechnol J. 2022;20:1182–96.

Zhao CY, Xue J, Cai XD, Guo J, Li B, Wu S. Assessment of the key aroma compounds in rose-based products. J Food Drug Anal. 2016;24:471–6.

Zhong MC, Jiang XD, Yang GQ, Cui WH, Suo ZQ, Wang WJ, et al. Rose without prickle: genomic insights linked to moisture adaptation. Natl Sci Rev. 2021;8:nwab092.

Zhou LJ, Yu C, Cheng BX, Wan HH, Luo L, Pan HT, et al. Volatile compound analysis and aroma evaluation of tea-scented roses in China. Ind Crop Prod. 2020a;155: 112735.

Zhou Q, Tang D, Huang W, Yang ZM, Zhang Y, Hamilton JP, et al. Haplotype-resolved genome analyses of a heterozygous diploid potato. Nat Genet. 2020b;52:1018–23.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.