A high-quality genome sequence of Rosa chinensis to elucidate ornamental traits

A chromosome-scale and haplotype-resolved genome assembly of tetraploid blackberry (Rubus L. subgenus Rubus Watson)

Blackberries (Rubus spp.) are globally consumed and well known for their rich anthocyanin and antioxidant content and distinct flavors. Improving blackberries has been challenging due to genetic complexity of traits and limited genomic resources. The blackberry genome has been particularly challenging to assemble due to its polyploid nature. Here, we present the first chromosome-scale and haplotype-phased assembly for the primocane-fruiting, thornless tetraploid blackberry selection BL1 (Rubus L. subgenus Rubus Watson). The genome assembly was generated using Oxford Nanopore Technology and Hi-C scaffolding, resulting in a 919 Mb genome distributed across 27 pseudochromosomes, with an N50 of 35.73 Mb. This assembly covers >92% of the genome length and contains over 98% of complete BUSCOs. Approximately, 58% of the assembly consists of repetitive sequences, with long terminal repeats being the most abundant class. A total of 87,968 protein-coding genes were predicted, of which, 82% were functionally annotated. Genome mining and RNA-Seq analyses identified possible candidate genes and transcription factors related to thornlessness and the key structural genes and transcription factors for anthocyanin biosynthesis. Activator genes including PAP1 and TTG1 and repressor genes such as ANL2 and MYBPA1 play an important role in the fine tuning of anthocyanin production during blackberry development. Resequencing of seven tetraploid blackberry cultivars/selections with different horticultural characteristics revealed candidate genes that could impact fruiting habit and disease resistance/susceptibility. This tetraploid reference genome should provide a valuable resource for accelerating genetic analysis of blackberries and facilitating the development of new improved cultivars with enhanced horticultural and nutritional traits.

Chromosomal level genome assembly of medicinal plant Rosa laevigata

Rosa laevigata Michx., an endemic perennial herbaceous plant in China, possesses significant medicinal value in traditional herbal practices. However, the absence of a reference genome has hindered its development and utilization. In this study, we present a chromosome-level de novo genome assembly of two haplotypes (Hap1 and Hap2) of R. laevigata by integration of Hifi long reads, BGI short reads, and Hi-C reads. The assembled Hap1 genome spans 493 Mb, Hap2 genome spans 479 Mb, and both of them were assigned to 7 chromosomes each. The mapping rate of BGI short reads to the two haplotypes genome is approximately 99.26% and 99.23%, and BUSCO assessment reveals that 98.6% and 98.7% of the genes are complete. Furthermore, we predicted 43,480 and 41,251 protein-coding genes in two haplotype genomes, respectively. The chromosome-level genome of R. laevigata enhances the genetic resources available for Rosa species and lays the groundwork for subsequent medicinal development.

Divergent Mechanisms of Internode Elongation in Response to Far-Red in Two Rose Genotypes

The quality of potted ornamental plants depends on their architecture, which should be compact and branched. Among the techniques for controlling this architecture, LED lighting, by manipulating light quality, offers an effective means of regulating elongation and branching. In rose, the addition of far-red (FR) light stimulated branching but induced excessive stem elongation, i.e., internode elongation. However, some varieties remained insensitive to this effect, demonstrating phenotypic stability. This study investigated the underlying mechanisms of internode elongation in response to FR in two rose cultivars, 'The Fairy' (TF) and Knock Out® Radrazz (KO), selected for their respective architectural plasticity and stability to FR. In TF, exposure to FR induced elongation of internodes, driven by cell division, with an increase in gibberellin A4 (GA4) level and a reduction in defense hormones (salicylic acid and jasmonic acid; JA). In contrast, in KO, exposure to FR did not induce internode elongation but caused cell elongation. This effect was accompanied by a reduction in cell number, modulated by hormonal changes (particularly GA4 and JA) and the inhibition of Block of cell proliferation 1, thereby limiting cell division. A deeper understanding of the mechanisms underlying architectural stability might lead to developing strategies to produce compact, branched plants, regardless of environmental conditions.

Phenotypic and genomic signatures across wild Rosa species open new horizons for modern rose breeding

The cultivation and domestication of roses reflects cultural exchanges and shifts in aesthetics that have resulted in today's most popular ornamental plant group. However, the narrow genetic foundation of cultivated roses limits their further improvement. Wild Rosa species harbour vast genetic diversity, yet their utilization is impeded by taxonomic confusion. Here we generated a phased and gap-free reference genome of Rosa persica for phylogenetic and population genomic analyses of a large collection of Rosa samples. The robust nuclear and plastid phylogenies support most of the morphology-based traditional taxonomy of Rosa. Population genomic analyses disclosed potential genetic exchanges among sections, indicating the northwest and southwest of China as two independent centres of diversity for Rosa. Analyses of domestication traits provide insights into selection processes related to flower colour, fragrance, double flower and resistance. This study provides a comprehensive understanding of rose domestication and lays a solid foundation for future re-domestication and innovative breeding efforts using wild resources.

A chromosomal-scale reference genome for Rosa hugonis

Rosa hugonis is widely distributed in the Hengduan Mountains, Qinling Mountains, and northern China. It is an important candidate species for ecological restoration, given its good adaptability. Here, we present the first high-quality chromosome-level assembly of R. hugonis based on HiFi reads and Hi-C data. The sequencing data were then assembled onto seven pseudochromosomes of R. hugonis. The genome sizes of R. hugonis is 337.92 Mb, with contig N50 length of 26.84 Mb. We annotated 36,218 protein-coding genes in R. hugonis. In summary, the high-quality genome sequences of R. hugonis provide a genetic roadmap for the study of its genetics and species relationships. This will facilitate future genomic comparative studies across more species within Rosa.

A high-quality chromosome-scale genome assembly of the Cherokee rose (Rosa laevigata)

Rosa laevigata is an excellent rose germplasm, highly resistant to aphid, and immune to both rose black spot and powdery mildew disease. It is also a well-known edible plant with a long history of medicinal use in China, having the effects of improving kidney function, inhibiting arteriosclerosis, and reducing inflammation. In this study, we assembled a high-quality chromosome-scale genome for R. laevigata by combining Illumina, PacBio, and Hi-C data, which has a length of approximately 494.2 Mb with a scaffold N50 of 68.6 Mb. A total of 493.2 Mb (99.8%) of the draft genome sequences were anchored on seven pseudochromosomes and two gapless pseudochromosomes were included in the final genome assembly. A total of 37,117 protein-coding genes were predicted, 34,047 of which were functionally annotated. Repeat annotation revealed 659,558 (285.6 Mb) repeat elements, accounting for 57.8% of the genome. The chromosome-scale genome provides valuable information to facilitate comparative genomic analysis of rose family and will accelerate genome-guided breeding and germplasm improvement of both R. laevigata itself and modern roses.

The telomere-to-telomere genome of flowering cherry (Prunus campanulata) reveals genomic evolution of the subgenus Cerasus

BACKGROUND

Prunus campanulata, a species of ornamental cherry, holds significant genetic and horticultural value. Despite the availability of various cherry genomes, a fully resolved telomere-to-telomere (T2T) assembly for this species has been lacking. Recent advancements in long-read sequencing technologies have made it possible to generate gap-free genome assemblies, providing comprehensive insights into genomic structures that were previously inaccessible.

FINDINGS

We present the first T2T genome assembly for P. campanulata "Lianmeiren" (v2.0), achieved through the integration of PacBio HiFi, ultra-long Oxford Nanopore Technologies, Illumina, and Hi-C sequencing. The assembly resulted in a highly contiguous genome with a total size of 266.23 Mb and a contig N50 of 31.6 Mb. The genome exhibits remarkable completeness (98.9% BUSCO) and high accuracy (quality value of 48.75). Additionally, 13 telomeres and putative centromere regions were successfully identified across the 8 pseudochromosomes. Comparative analysis with the previous v1.0 assembly revealed 336,943 single nucleotide polymorphisms, 107,521 indels, and 1,413 structural variations, along with the annotation of 1,402 new genes.

CONCLUSIONS

This T2T genome assembly of P. campanulata "Lianmeiren" provides a critical reference for understanding the genetic architecture of the species. It enhances our ability to study structural variations, gene function, and evolutionary biology within the Prunus genus.

Comparative Analysis of Ca2+/Cation Antiporter Gene Family in Rosa roxburghii and Enhanced Calcium Stress Tolerance via Heterologous Expression of RrCAX1a in Tobacco

Rosa roxburghii, a calciphilic species native to the mountainous regions of Southwest China, is renowned for its high vitamin C and bioactive components, making it valuable for culinary and medicinal uses. This species exhibits remarkable tolerance to the high-calcium conditions typical of karst terrains. However, the underlying mechanisms of this calcium resilience remain unclear. The Ca2+/cation antiporter (CaCA) superfamily plays a vital role in the transport of Ca2+ and other cations and is crucial for plant tolerance to metal stress. However, the roles and evolutionary significance of the CaCA superfamily members in R. roxburghii remain poorly understood. This study identified 22 CaCA superfamily genes in R. roxburghii, categorized into four subfamilies. The gene structures of these RrCaCAs show considerable conservation across related species. Selection pressure analysis revealed that all RrCaCAs are subject to purifying selection. The promoter regions of these genes contain numerous hormone-responsive and stress-related elements. qRT-PCR analyses demonstrated that H+/cation exchanger (CAX) RrCAX1a and RrCAX3a were highly responsive to Ca2+ stress, cation/Ca2+ exchanger (CCX) RrCCX4 to Mg2+ stress, and RrCCX11a to Na+ stress. Subcellular localization indicated that RrCAX1a is localized to the plant cell membrane, and its stable transformation in tobacco confirmed its ability to confer enhanced resistance to heavy Ca2+ stresses, highlighting its crucial role in the high-calcium tolerance mechanisms of R. roxburghii. This research establishes a foundation for further molecular-level functional analyses of the adaptation mechanisms of R. roxburghii to high-calcium environments.

Sugar Transport and Signaling in Shoot Branching

The source-sink relationship is critical for proper plant growth and development, particularly for vegetative axillary buds, whose activity shapes the branching pattern and ultimately the plant architecture. Once formed from axillary meristems, axillary buds remain dormant or become active to grow into new branches. This transition is notably driven by the regulation of the bud sink strength, which is reflected in the ability to unload, metabolize and store photoassimilates. Plants have so far developed two main mechanisms for unloading sugars (sucrose) towards sink organs, a symplasmic pathway and an apoplasmic pathway, but so far limited investigations have been reported about the modes of sugar uptake during the transition from the dormant to the active outgrowth state of the bud. The available data indicate that the switch from dormant bud to active outgrowing state, requires sugar and is shortly preceded by an increase in bud metabolic activity and a remobilization of the stem starch reserves in favor of growing buds. This activation of the bud sink strength is accompanied by an up-regulation of the main markers of apoplasmic unloading, such as sugar transporters (sucrose transporters-SUTs; sugar will eventually be exported transporters-SWEETs), sucrose hydrolyzing enzymes (cell wall invertase-CWINV) and sugar metabolic pathways (glycolysis/tricarboxylic cycle-TCA; oxidative pentose phosphate pathway-OPPP). As these results are limited to a few species, they are not sufficient to provide a complete and accurate picture of the mode(s) of sugar unloading toward axillary buds and deserve to be complemented by additional studies in a wide variety of plants using systems integration, combining genetic, molecular and immunolocalization approaches. Altogether, we discuss here how sugar is a systemic regulator of shoot branching, acting both as an energy-rich molecule and a signaling entity in the establishment of the bud sink strength.

Unraveling the genetic basis of Rhizobium rhizogenes-mediated transformation and hairy root formation in rose using a genome-wide association study

KEY MESSAGE

Multiple QTLs reveal the polygenic nature of R. rhizogenes-mediated transformation and hairy root formation in roses, with five key regions explaining 12.0-26.9% of trait variability and transformation-related candidate genes identified. Understanding genetic mechanisms of plant transformation remains crucial for biotechnology. This is particularly relevant for roses and other woody ornamentals that exhibit recalcitrant behavior in transformation procedures. Rhizobium rhizogenes-mediated transformation leading to hairy root (HR) formation provides an excellent model system to study transformation processes and host-pathogen interactions. Therefore, this study aimed to identify quantitative trait loci (QTLs) associated with HR formation and explore their relationship with adventitious root (AR) formation in rose as a model for woody ornamentals. A diversity panel of 104 in vitro grown rose genotypes was transformed with R. rhizogenes strain ATCC 15834 carrying a green fluorescent protein reporter gene. Phenotypic data on callus and root formation were collected for laminae and petioles. A genome-wide association study using 23,419 single-nucleotide polymorphism markers revealed significant QTLs on chromosomes one and two for root formation traits. Five key genomic regions explained 12.0-26.9% of trait variability, with some peaks overlapping previously reported QTLs for AR formation. This genetic overlap was supported by weak to moderate correlations between HR and AR formation traits, particularly in petioles. Candidate gene identification through literature review and transcriptomic data analysis revealed ten candidate genes involved in bacterial response, hormone signaling, and stress responses. Our findings provide new insights into the genetic control of HR formation in roses and highlight potential targets for improving transformation efficiency in ornamental crops, thereby facilitating future research and breeding applications.

Development of a robust SNP marker set for genotyping diverse gene bank collections of polyploid roses

BACKGROUND

Due to genetic depletion in nature, gene banks play a critical role in the long-term conservation of plant genetic resources and the provision of a wide range of plant genetic diversity for research and breeding programs. Genetic information on accessions facilitates gene bank management and can help to conserve limited resources and to identify taxonomic misclassifications or mislabelling. Here, we developed SNP markers for genotyping 4,187 mostly polyploid rose accessions from large rose collections, including the German Genebank for Roses.

RESULTS

We filtered SNP marker information from the RhWag68k Axiom SNP array using call rates, uniformity of the four allelic dosage groups and chromosomal position to improve genotyping efficiency. After conversion to individual PACE® markers and further filtering, we selected markers with high discriminatory power. These markers were used to analyse 4,187 accessions with a mean call rate of 91.4%. By combining two evaluation methods, the mean call rate was increased to 95.2%. Additionally, the robustness against the genotypic groups used for calling was evaluated, resulting in a final set of 18 markers. Analyses of 94 pairs of assumed duplicate accessions included as controls revealed unexpected differences for eight pairs, which were confirmed using SSR markers. After removing the duplicates and filtering for accessions that were robustly called with all 18 markers, 141 out of the 1,957 accessions showed unexpected identical marker profiles with at least one other accession in our PACE® and SSR analysis. Given the attractiveness of NGS technologies, 13 SNPs from the marker set were also analysed using amplicon sequencing, with 76% agreement observed between PACE® and amplicon markers.

CONCLUSIONS

Although sampling error cannot be completely excluded, this is an indication that mislabelling occurs in rose collections and that molecular markers may be able to detect these cases. In future applications, our marker set could be used to develop a core reference set of representative accessions, and thus optimise the selection of gene bank accessions.

Haplotype-resolved genome assembly and resequencing provide insights into the origin and breeding of modern rose

Modern rose (Rosa hybrida) is a recently formed interspecific hybrid and has become one of the most important and widely cultivated ornamentals. Here we report the haplotype-resolved chromosome-scale genome assembly of the tetraploid R. hybrida 'Samantha' ('JACmantha') and a genome variation map of 233 Rosa accessions involving various wild species, and old and modern cultivars. Homologous chromosomes of 'Samantha' exhibit frequent homoeologous exchanges. Population genomic and genomic composition analyses reveal the contributions of wild Rosa species to modern roses and highlight that R. odorata and its derived cultivars are important contributors to modern roses, much like R. chinensis 'Old Blush'. Furthermore, selective sweeps during modern rose breeding associated with major agronomic traits, including continuous and recurrent flowering, double flower, flower senescence and disease resistance, are identified. This study provides insights into the genetic basis of modern rose origin and breeding history, and offers unprecedented genomic resources for rose improvement.

Genome-wide association study (GWAS) analyses of early anatomical changes in rose adventitious root formation

Adventitious root (AR) formation is a genetically complex trait with high genotypic variability. Therefore, only a limited range of cultivars are currently propagated by cuttings in rose. In this study, we analysed the anatomy of in vitro shoots, the early formation of root primordia (RP) and the formation of ARs in a diverse set of 106 rose genotypes. Correlation analysis indicated that the growth in shoot diameter and the vasculature dimensions after 1 week of rooting contributed to successful AR formation. Using phenotypic data for genome-wide association studies (GWAS) analyses, nine significantly associated single nucleotide polymorphisms (SNPs) and genomic regions contributing to various RP and AR formation traits were identified. The contribution of genomic regions to trait variation was notably greater for traits associated with earlier processes than for traits associated with later developmental stages. The combination of RP and AR data allowed the detection of regions by GWAS that contain factors that potentially limit RP emergence. Homologues of 47 genes known to be involved in AR formation from the literature could be assigned to the identified peaks. Further studies are needed to investigate the suitability of SNPs exhibiting strong effects as allele-specific PCR markers for use in breeding.

Multi-omics analyzes of Rosa gigantea illuminate tea scent biosynthesis and release mechanisms

Rose is an important ornamental crop cultivated globally for perfume production. However, our understanding of the mechanisms underlying scent production and molecular breeding for fragrance is hindered by the lack of a reference genome for tea roses. We present the first complete telomere-to-telomere (T2T) genome of Rosa gigantea, with high quality (QV > 60), including detailed characterization of the structural features of repetitive regions. The expansion of genes associated with phenylpropanoid biosynthesis may account for the unique tea scent. We uncover the release rhythm of aromatic volatile organic compounds and their gene regulatory networks through comparative genomics and time-ordered gene co-expression networks. Analyzes of eugenol homologs demonstrate how plants attract pollinators using specialized phenylpropanoids in specific tissues. This study highlights the conservation and utilization of genetic diversity from wild endangered species through multi-omics approaches, providing a scientific foundation for enhancing rose fragrance via de novo domestication.

Genetic dissection of stem and leaf rachis prickles in diploid rose using a pedigree-based QTL analysis

Introduction

Prickles are often deemed undesirable traits in many crops, including roses (Rosa sp.), and there is demand for rose cultivars with no or very few prickles. This study aims to identify new and/or validate reported quantitative trait loci (QTLs) associated with stem and leaf rachis prickle density, characterize the effects of functional haplotypes for major QTLs, and identify the sources of QTL-alleles associated with increased/decreased prickle density in roses.

Methods

QTL mapping using pedigree-based analysis (PBA), and haplotype analysis were conducted on two multi-parental diploid rose populations (TX2WOB and TX2WSE).

Results and discussion

Twelve QTLs were identified on linkage groups (LGs) 2, 3, 4, and 6. The major QTLs for the stem prickle density were located between 42.25 and 45.66 Mbp on chromosome 3 of the Rosa chinensis genome assembly, with individual QTLs explaining 18 to 49% of phenotypic variance (PVE). The remaining mapped QTLs were minor. As for the rachis prickle density, several QTLs were detected on LG3, 4, and 6 with PVE 8 to 17%. Also, this study identified that ancestors R. wichurana 'Basye's Thornless', 'Old Blush', and the pollen parent of M4-4 were common sources of favorable alleles (q) associated with decreased prickle density, whereas 'Little Chief' and 'Srche Europy' were the source of unfavorable alleles (Q) in the TX2WOB and TX2WSE populations, respectively. The outcomes of this work complement other studies to locate factors that affect prickle density. These results can also be utilized to develop high-throughput DNA tests and apply parental selection to develop prickle-free rose cultivars.

Genomic data provides insights into the evolutionary history and adaptive differentiation of two tetraploid strawberries

Over the decades, evolutionists and ecologists have shown intense interest in the role of polyploidization in plant evolution. Without clear knowledge of the diploid ancestor(s) of polyploids, we would not be able to answer fundamental ecological questions such as the evolution of niche differences between them or its underlying genetic basis. Here, we explored the evolutionary history of two Fragaria tetraploids, Fragaria corymbosa and Fragaria moupinensis. We de novo assembled five genomes including these two tetraploids and three diploid relatives. Based on multiple lines of evidence, we found no evidence of subgenomes in either of the two tetraploids, suggesting autopolyploid origins. We determined that Fragaria chinensis was the diploid ancestor of F. corymbosa while either an extinct species affinitive to F. chinensis or an unsampled population of F. chinensis could be the progenitor of F. moupinensis. Meanwhile, we found introgression signals between F. chinensis and Fragaria pentaphylla, leading to the genomic similarity between these two diploids. Compared to F. chinensis, gene families related to high ultraviolet (UV)-B and DNA repair were expanded, while those that responded towards abiotic and biotic stresses (such as salt stress, wounding, and various pathogens) were contracted in both tetraploids. Furthermore, the two tetraploids tended to down-regulate defense response genes but up-regulate UV-B response, DNA repairing, and cell division gene expression compared to F. chinensis. These findings may reflect adaptions toward high-altitude habitats. In summary, our work provides insights into the genome evolution of wild Fragaria tetraploids and opens up an avenue for future works to answer deeper evolutionary and ecological questions regarding the strawberry genus.

Genome-Wide Analyses of CCHC Family Genes and Their Expression Profiles under Drought Stress in Rose (Rosa chinensis)

CCHC-type zinc finger proteins (CCHC-ZFPs), ubiquitous across plant species, are integral to their growth, development, hormonal regulation, and stress adaptation. Roses (Rosa sp.), as one of the most significant and extensively cultivated ornamentals, account for more than 30% of the global cut-flower market. Despite its significance, the CCHC gene family in roses (Rosa sp.) remains unexplored. This investigation identified and categorized 41 CCHC gene members located on seven chromosomes of rose into 14 subfamilies through motif distribution and phylogenetic analyses involving ten additional plant species, including Ginkgo biloba, Ostreococcus lucimarinus, Arabidopsis thaliana, and others. This study revealed that dispersed duplication likely plays a crucial role in the diversification of the CCHC genes, with the Ka/Ks ratio suggesting a history of strong purifying selection. Promoter analysis highlighted a rich presence of cis-acting regulatory elements linked to both abiotic and biotic stress responses. Differential expression analysis under drought conditions grouped the 41 CCHC gene members into five distinct clusters, with those in group 4 exhibiting pronounced regulation in roots and leaves under severe drought. Furthermore, virus-induced gene silencing (VIGS) of the RcCCHC25 member from group 4 compromised drought resilience in rose foliage. This comprehensive analysis lays the groundwork for further investigations into the functional dynamics of the CCHC gene family in rose physiology and stress responses.

Genome assembly of Hibiscus sabdariffa L. provides insights into metabolisms of medicinal natural products

Hibiscus sabdariffa L. is a widely cultivated herbaceous plant with diverse applications in food, tea, fiber, and medicine. In this study, we present a high-quality genome assembly of H. sabdariffa using more than 33 Gb of high-fidelity (HiFi) long-read sequencing data, corresponding to ∼20× depth of the genome. We obtained 3 genome assemblies of H. sabdariffa: 1 primary and 2 partially haplotype-resolved genome assemblies. These genome assemblies exhibit N50 contig lengths of 26.25, 11.96, and 14.50 Mb, with genome coverage of 141.3, 86.0, and 88.6%, respectively. We also utilized 26 Gb of total RNA sequencing data to predict 154k, 79k, and 87k genes in the respective assemblies. The completeness of the primary genome assembly and its predicted genes was confirmed by the benchmarking universal single-copy ortholog analysis with a completeness rate of 99.3%. Based on our high-quality genomic resources, we constructed genetic networks for phenylpropanoid and flavonoid metabolism and identified candidate biosynthetic genes, which are responsible for producing key intermediates of roselle-specific medicinal natural products. Our comprehensive genomic and functional analysis opens avenues for further exploration and application of valuable natural products in H. sabdariffa.

RcOST1L phosphorylates RcPIF4 for proteasomal degradation to promote flowering in rose

Flowering is a vital agronomic trait that determines the economic value of most ornamental plants. The flowering time of rose (Rosa spp.) is photoperiod insensitive and is thought to be tightly controlled by light intensity, although the detailed molecular mechanism remains unclear. Here, we showed that rose plants flower later under low-light (LL) intensity than under high-light (HL) intensity, which is mainly related to the stability of PHYTOCHROME-INTERACTING FACTORs (RcPIFs) mediated by OPEN STOMATA 1-Like (RcOST1L) under different light intensity regimes. We determined that HL conditions trigger the rapid phosphorylation of RcPIFs before their degradation. A yeast two-hybrid screen identified the kinase RcOST1L as interacting with RcPIF4. Moreover, RcOST1L positively regulated rose flowering and directly phosphorylated RcPIF4 on serine 198 to promote its degradation under HL conditions. Additionally, phytochrome B (RcphyB) enhanced RcOST1L-mediated phosphorylation of RcPIF4 via interacting with the active phyB-binding motif. RcphyB was activated upon HL and recruited RcOST1L to facilitate its nuclear accumulation, in turn leading to decreased stability of RcPIF4 and flowering acceleration. Our findings illustrate how RcPIF abundance safeguards proper rose flowering under different light intensities, thus uncovering the essential role of RcOST1L in the RcphyB-RcPIF4 module in flowering.

PE (Prickly Eggplant) encoding a cytokinin-activating enzyme responsible for the formation of prickles in eggplant

Eggplant is one of the most important vegetables worldwide, with some varieties displaying prickles. These prickles, present on the leaves, stems, and fruit calyxes, posing challenges during cultivation, harvesting, and transportation, making them an undesirable agronomic trait. However, the genetic mechanisms underlying prickle morphogenesis in eggplant remain poorly understood, impeding genetic improvements. In this study, genetic analyses revealed that prickle morphogenesis is governed by a single dominant nuclear gene, termed PE (Prickly Eggplant). Subsequent bulk segregant RNA-sequencing (BSR-seq) and linkage analysis preliminarily mapped PE to chromosome 6. This locus was then fine mapped to a 9233 bp interval in a segregating population of 1109 plants, harboring only one candidate gene, SmLOG1, which encodes a LONELY GUY (LOG)-family cytokinin biosynthetic enzyme. Expression analyses via transcriptome and qRT-PCR demonstrate that SmLOG1 is predominantly expressed in immature prickles. CRISPR-Cas9 knockout experiments targeting SmLOG1 in prickly parental line 'PI 381159' abolished prickles across all tissues, confirming its critical role in prickle morphogenesis. Sequence analysis of SmLOG1 pinpointed variations solely within the non-coding region. We developed a cleaved amplified polymorphic sequences (CAPS) marker from a distinct SNP located at -735-bp within the SmLOG1 promoter, finding significant association with prickle variation in 190 eggplant germplasms. These findings enhance our understanding of the molecular mechanisms governing prickle development in eggplant and facilitate the use of marker-assisted selection (MAS) for breeding prickleless cultivars.

Chromosome-level genome assembly of a cliff plant Taihangia rupestris var. ciliata provides insights into its adaptation and demographic history

BACKGROUND

Cliffs are recognized as one of the most challenging environments for plants, characterized by harsh conditions such as drought, infertile soil, and steep terrain. However, they surprisingly host ancient and diverse plant communities and play a crucial role in protecting biodiversity. The Taihang Mountains, which act as a natural boundary in eastern China, support a rich variety of plant species, including many unique to cliff habitats. However, it is little known how cliff plants adapt to harsh habitats and the demographic history in this region.

RESULTS

To better understand the demographic history and adaptation of cliff plants in this area, we analyzed the chromosome-level genome of a representative cliff plant, T. rupestris var. ciliata, which has a genome size of 769.5 Mb, with a scaffold N50 of 104.92 Mb. The rapid expansion of transposable elements may have contributed to the increasing genome and its ability to adapt to unique and challenging cliff habitats. Comparative analysis of the genome evolution between Taihangia and non-cliff plants in Rosaceae revealed a significant expansion of gene families associated with oxidative phosphorylation, which is likely a response to the abiotic stresses faced by cliff plants. This expansion may explain the long-term adaptation of Taihangia to harsh cliff environments. The effective population size of the two varieties has continuously decreased due to climatic fluctuations during the Quaternary period. Furthermore, significant differences in gene expression between the two varieties may explain the varied leaf phenotypes and adaptations to harsh conditions in different natural distributions.

CONCLUSION

Our study highlights the extraordinary adaptation of T. rupestris var. ciliata, shedding light on the evolution of cliff plants worldwide.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

Molecular and genetic regulation of petal number variation

Floral forms with an increased number of petals, also known as double-flower phenotypes, have been selected and conserved in many domesticated plants, particularly in ornamentals, because of their great economic value. The molecular and genetic mechanisms that control this trait are therefore of great interest, not only for scientists, but also for breeders. In this review, we summarize current knowledge of the gene regulatory networks of flower initiation and development and known mutations that lead to variation of petal number in many species. In addition to the well-accepted miR172/AP2-like module, for which many questions remain unanswered, we also discuss other pathways in which mutations also lead to the formation of extra petals, such as those involved in meristem maintenance, hormone signalling, epigenetic regulation, and responses to environmental signals. We discuss how the concept of 'natural mutants' and recent advances in genomics and genome editing make it possible to explore the molecular mechanisms underlying double-flower formation, and how such knowledge could contribute to the future breeding and selection of this trait in more crops.

The transcription factor RhMYB17 regulates the homeotic transformation of floral organs in rose (Rosa hybrida) under cold stress

Low temperatures affect flower development in rose (Rosa hybrida), increasing petaloid stamen number and reducing normal stamen number. We identified the low-temperature-responsive R2R3-MYB transcription factor RhMYB17, which is homologous to Arabidopsis MYB17 by similarity of protein sequences. RhMYB17 was up-regulated at low temperatures, and RhMYB17 transcripts accumulated in floral buds. Transient silencing of RhMYB17 by virus-induced gene silencing decreased petaloid stamen number and increased normal stamen number. According to the ABCDE model of floral organ identity, class A genes APETALA 1 (AP1) and AP2 contribute to sepal and petal formation. Transcription factor binding analysis identified RhMYB17 binding sites in the promoters of rose APETALA 2 (RhAP2) and APETALA 2-LIKE (RhAP2L). Yeast one-hybrid assays, dual-luciferase reporter assays, and electrophoretic mobility shift assays confirmed that RhMYB17 directly binds to the promoters of RhAP2 and RhAP2L, thereby activating their expression. RNA sequencing further demonstrated that RhMYB17 plays a pivotal role in regulating the expression of class A genes, and indirectly influences the expression of the class C gene. This study reveals a novel mechanism for the homeotic transformation of floral organs in response to low temperatures.

Indole-3 acetic acid negatively regulates rose black spot disease resistance through antagonizing the salicylic acid signaling pathway via jasmonic acid

MAIN CONCLUSION

IAA cooperates with JA to inhibit SA and negatively regulates rose black spot disease resistance. Black spot disease caused by the fungus Marssonina rosae is the most prevalent and severe ailment in rose cultivation, leading to the appearance of black spots on leaves and eventual leaf fall, significantly impacting the utilization of roses in gardens. Salicylic acid (SA) and jasmonic acid (JA) are pivotal hormones that collaborate with indole-3 acetic acid (IAA) in regulating plant defense responses; however, the detailed mechanisms underlying the induction of black spot disease resistance by IAA, JA, and SA remain unclear. In this study, transcript analysis was conducted on resistant (R13-54) and susceptible (R12-26) lines following M. rosae infection. In addition, the impact of exogenous interference with IAA on SA- and JA-mediated disease resistance was examined. The continuous accumulation of JA, in synergy with IAA, inhibited activation of the SA signaling pathway in the early infection stage, thereby negatively regulating the induction of effective resistance to black spot disease. IAA administration alleviated the inhibition of SA on JA to negatively regulate the resistance of susceptible strains by further enhancing the synthesis and accumulation of JA. However, IAA did not contribute to the negative regulation of black spot resistance when high levels of JA were inhibited. Virus-induced gene silencing of RcTIFY10A, an inhibitor of the JA signaling pathway, further suggested that IAA upregulation led to a decrease in disease resistance, a phenomenon not observed when the JA signal was inhibited. Collectively, these findings indicate that the IAA-mediated negative regulation of black spot disease resistance relies on activation of the JA signaling pathway.

Advances in genomics and genome editing for improving strawberry (Fragaria ×ananassa)

The cultivated strawberry, Fragaria ×ananassa, is a recently domesticated fruit species of economic interest worldwide. As such, there is significant interest in continuous varietal improvement. Genomics-assisted improvement, including the use of DNA markers and genomic selection have facilitated significant improvements of numerous key traits during strawberry breeding. CRISPR/Cas-mediated genome editing allows targeted mutations and precision nucleotide substitutions in the target genome, revolutionizing functional genomics and crop improvement. Genome editing is beginning to gain traction in the more challenging polyploid crops, including allo-octoploid strawberry. The release of high-quality reference genomes and comprehensive subgenome-specific genotyping and gene expression profiling data in octoploid strawberry will lead to a surge in trait discovery and modification by using CRISPR/Cas. Genome editing has already been successfully applied for modification of several strawberry genes, including anthocyanin content, fruit firmness and tolerance to post-harvest disease. However, reports on many other important breeding characteristics associated with fruit quality and production are still lacking, indicating a need for streamlined genome editing approaches and tools in Fragaria ×ananassa. In this review, we present an overview of the latest advancements in knowledge and breeding efforts involving CRISPR/Cas genome editing for the enhancement of strawberry varieties. Furthermore, we explore potential applications of this technology for improving other Rosaceous plant species.

Cotyledonary somatic embryo is one kind of intermediate material similar to callus in the process of in vitro tissue culture from Rosa hybrida 'John F. Kennedy'

BACKGROUND

Rose is recognized as an important ornamental plant worldwide, and it is also one of the most widely used flowers in gardens. At present, the improvement of rose traits is still difficult and uncertain, and molecular breeding can provide new ideas for the improvement of modern rose varieties. Somatic embryos are quite good receptors for genetic transformation. However, little is known about the molecular mechanisms underlying during the regeneration process of rose somatic embryos. To elucidate the molecular regulation mechanism of somatic embryo plantlet regeneration, the relationship between the differences in traits of the two different regenerated materials and the significantly differentially expressed genes (DEGs) related to phytohormone pathways in the process of regeneration were be investigated.

RESULTS

These representative two regenerated samples from single-piece cotyledonary somatic embryo (SPC) culture of Rosa hybrida 'John F. Kennedy', were harvested for transcriptome analysis, with the SPC explants at the initial culture (Day 0) as the control. The differentially expressed genes (DEGs) in the materials from two different types for regeneration approach (SBF type: the regeneration approach type of single bud formed from SPC explants; MBF type: the regeneration approach type of multiple buds formed from SPC explants) were be screened by means of the transcriptome sequencing technology. In this study, a total of about 396.24 million clean reads were obtained, of which 78.95-82.92% were localized to the reference genome, compared with the initial material (CK sample), there were 5594 specific genes in the material of SBF type and 6142 specific genes in the MBF type. The DEGs from the SBF type material were mainly concentrated in the biological processes of GO terms such as phytohormones, substance transport, cell differentiation, and redox reaction. The KEGG enrichment analysis revealed these DEGs were more active in ubiquinone and other terpenoid-quinone biosynthesis, fatty acid elongation, steroid biosynthesis, and glycosphingolipid biosynthesis-globo and isoglobo series. In contrast, the DEGs induced by the MBF type material were mainly associated with the biological processes such as phytohormones, phosphorylation, photosynthesis and signal transduction. According to KEGG analysis, these DEGs of MBF type were significantly enriched in the porphyrin and chlorophyll metabolism, brassinosteroid biosynthesis, carotenoid biosynthesis, and peroxisome. Furthermore, the results from the phytohormone pathways analysis showed that the auxin-responsive factor SAUR and the cell wall modifying enzyme gene XTH were upregulated for expression but the protein phosphatase gene PP2C was downregulated for expression in SBF type; the higher expression of the ethylene receptor ETR, the ethylene transduction genes EBF1/2, the transcription factor EIN3, and the ethylene-responsive transcription factor ERF1/2 were induced by MBF type.

CONCLUSIONS

According to the GO and KEGG analysis, it indicated the DEGs between two different regenerated materials from somatic embryos were significantly different which might be causing morphological differences. That was somatic embryos from Rosa hybrida 'John F. Kennedy' could regenerate plantlet via both classic somatic embryogenesis (seed-like germination) and organogenesis, cotyledonary somatic embryos should be considered as one kind of intermediate materials similiar to callus, rather than the indicator materials for somatic embryogenesis.

Genome-Wide Analysis of the WRKY Transcription Factor Family in Roses and Their Putative Role in Defence Signalling in the Rose-Blackspot Interaction

WRKY transcription factors are important players in plant regulatory networks, where they control and integrate various physiological processes and responses to biotic and abiotic stresses. Here, we analysed six rose genomes of 5 different species (Rosa chinensis, R. multiflora, R. roxburghii, R. sterilis, and R. rugosa) and extracted a set of 68 putative WRKY genes, extending a previously published set of 58 WRKY sequences based on the R. chinensis genome. Analysis of the promoter regions revealed numerous motifs related to induction by abiotic and, in some cases, biotic stressors. Transcriptomic data from leaves of two rose genotypes inoculated with the hemibiotrophic rose black spot fungus Diplocarpon rosae revealed the upregulation of 18 and downregulation of 9 of these WRKY genes after contact with the fungus. Notably, the resistant genotype exhibited the regulation of 25 of these genes (16 upregulated and 9 downregulated), while the susceptible genotype exhibited the regulation of 20 genes (15 upregulated and 5 downregulated). A detailed RT-qPCR analysis of RcWRKY37, an orthologue of AtWRKY75 and FaWRKY1, revealed induction patterns similar to those of the pathogenesis-related (PR) genes induced in salicylic acid (SA)-dependent defence pathways in black spot inoculation experiments. However, the overexpression of RcWRKY37 in rose petals did not induce the expression of any of the PR genes upon contact with black spot. However, wounding significantly induced the expression of RcWRKY37, while heat, cold, or drought did not have a significant effect. This study provides the first evidence for the role of RcWRKY37 in rose signalling cascades and highlights the differences between RcWRKY37 and AtWRKY75. These results improve our understanding of the regulatory function of WRKY transcription factors in plant responses to stress factors. Additionally, they provide foundational data for further studies.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Evaluation of aphid resistance on different rose cultivars and transcriptome analysis in response to aphid infestation

BACKGROUND

The rose is one of the most important ornamental flowers in the world for its aesthetic beauty but can be attacked by many pests such as aphids. Aphid infestation causes tremendous damage on plant tissues leading to harmed petals and leaves. Rose cultivars express different levels of resistance to aphid infestation yet the information remains unclear. Not only that, studies about the transcriptional analysis on defending mechanisms against aphids in rose are limited so far.

RESULTS

In this study, the aphid resistance of 20 rose cultivars was evaluated, and they could be sorted into six levels based on the number ratio of aphids. And then, a transcriptome analysis was conducted after aphid infestation in one high resistance (R, Harmonie) and one highly susceptibility (S, Carefree Wonder) rose cultivar. In open environment the majority of rose cultivars had the highest aphid number at May 6th or May 15th in 2020 and the resistance to infestation could be classified into six levels. Differential expression analysis revealed that there were 1,626 upregulated and 767 downregulated genes in the R cultivar and 481 upregulated and 63 downregulated genes in the S cultivar after aphid infestation. Pathway enrichment analysis of the differentially expressed genes revealed that upregulated genes in R and S cultivars were both enriched in defense response, biosynthesis of secondary metabolites (phenylpropanoid, alkaloid, and flavonoid), carbohydrate metabolism (galactose, starch, and sucrose metabolism) and lipid processing (alpha-linolenic acid and linolenic acid metabolism) pathways. In the jasmonic acid metabolic pathway, linoleate 13S-lipoxygenase was specifically upregulated in the R cultivar, while genes encoding other crucial enzymes, allene oxide synthase, allene oxide cyclase, and 12-oxophytodienoate reductase were upregulated in both cultivars. Transcription factor analysis and transcription factor binding search showed that WRKY transcription factors play a pivotal role during aphid infestation in the R cultivar.

CONCLUSIONS

Our study indicated the potential roles of jasmonic acid metabolism and WRKY transcription factors during aphid resistance in rose, providing clues for future research.

Comparative transcriptome and metabolome analysis revealed diversity in the response of resistant and susceptible rose (Rosa hybrida) varieties to Marssonina rosae

Rose black spot disease caused by Marssonina rosae is among the most destructive diseases that affects the outdoor cultivation and production of roses; however, the molecular mechanisms underlying the defensive response of roses to M. rosae have not been clarified. To investigate the diversity of response to M. rosae in resistant and susceptible rose varieties, we performed transcriptome and metabolome analyses of resistant (KT) and susceptible (FG) rose varieties and identified differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in response to M. rosae at different time points. In response to M. rosae, DEGs and DAMs were mainly upregulated compared to the control and transcription factors were concentrated in the WRKY and AP2/ERF families. Gene Ontology analysis showed that the DEGs of FG were mainly enriched in biological processes, such as the abscisic acid-activated signaling pathway, cell wall, and defense response, whereas the DEGs of KT were mainly enriched in Golgi-mediated vesicle transport processes. Kyoto Encyclopedia of Genes and Genomes analysis showed that the DEGs of both varieties were concentrated in plant-pathogen interactions, plant hormone signal transduction, and mitogen-activated protein kinase signaling pathways, with the greatest number of DEGs associated with brassinosteroid (BR) in the plant hormone signal transduction pathway. The reliability of the transcriptome results was verified by qRT-PCR. DAMs of KT were significantly enriched in the butanoate metabolism pathway, whereas DAMs of FG were significantly enriched in BR biosynthesis, glucosinolate biosynthesis, and tryptophan metabolism. Moreover, the DAMs in these pathways were significantly positively correlated with the DEGs. Disease symptoms were aggravated when FG leaves were inoculated with M. rosae after 24-epibrassinolide treatment, indicating that the response of FG to M. rosae involves the BR signaling pathway. Our results provide new insights into the molecular mechanisms underlying rose response to M. rosae and lay a theoretical foundation for formulating rose black spot prevention and control strategies and cultivating resistant varieties.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

A 49-bp deletion of PmAP2L results in a double flower phenotype in Prunus mume

The double flower is an important trait with substantial ornamental value. While mutations in PETALOSA TOE-type or AG (AGAMOUS) genes play a crucial role in enhancing petal number in ornamental plants, the complete mechanism underlying the formation of double flowers remains to be fully elucidated. Through the application of bulked segregant analysis (BSA), we identified a novel gene, APETALA2-like (PmAP2L), characterized by a 49-bp deletion in double-flowered Prunus mume. β-Glucuronidase (GUS) staining and luciferase reporter assays confirmed that the 49-bp deletion in PmAP2L reduced its binding with Pmu-miRNA172a. Phylogenetic analysis and microsynteny analysis suggested that PmAP2L was not a PETALOSA TOE-type gene, and it might be a new gene controlling the formation of double flower in P. mume. Subsequently, overexpression of PmAP2L-D in tobacco led to a significant rise in the number of stamens and the conversion of stamens to petals. Furthermore, silencing of the homologue of RC5G0530900 in rose significantly reduced the number of petals. Using transient gene expression in P. mume flower buds, we determined the functional differences between PmAP2L-D and PmAP2-S in controlling flower development. Meanwhile, DNA-affinity purification sequencing (DAP-seq), yeast hybrid assays and luciferase reporter assays indicated that PmAP2L negatively regulated the floral organ identity genes by forming a repressor complex with PmTPL and PmHDA6/19. Overall, these findings indicate that the variation in PmAP2L is associated with differences in the regulation of genes responsible for floral organ identity, providing new insights into the double-flower trait and double-flower breeding in plants.

Integrated proteomic analysis reveals interactions between phosphorylation and ubiquitination in rose response to Botrytis infection

As two of the most abundant post-translational modifications, phosphorylation and ubiquitination play a significant role in modulating plant-pathogen interactions and increasing evidence indicates their crosstalk in plant immunity. Rose (Rosa sp.) is one of the most important ornamental plants and can be seriously infected by Botrytis cinerea. Here, integrated proteomics analysis was performed to detect global proteome, phosphorylation, and ubiquitination changes in rose upon B. cinerea infection and investigate the possible phosphorylation and ubiquitination crosstalk. A total of 6165 proteins, 11 774 phosphorylation and 10 582 ubiquitination sites, and 77 phosphorylation and 13 ubiquitination motifs were identified. Botrytis cinerea infection resulted in 169 up-regulated and 122 down-regulated proteins, 291 up-regulated and 404 down-regulated phosphorylation sites, and 250 up-regulated and 634 down-regulated ubiquitination sites. There were 12 up-regulated PR10 proteins and half of them also showed reduced ubiquitination. A lot of kinases probably involved in plant pattern-triggered immunity signaling were up-regulated phosphoproteins. Noticeably, numerous kinases and ubiquitination-related proteins also showed a significant change in ubiquitination and phosphorylation, respectively. A cross-comparison of phosphoproteome and ubiquitylome indicated that both of two post-translational modifications of 104 proteins were dynamically regulated, and many putative pattern-triggered immunity signaling components in the plant plasma membrane were co-regulated. Moreover, five selected proteins, including four PR10 proteins and a plasma membrane aquaporin, were proven to be involved in rose resistance to B. cinerea. Our study provides insights into the molecular mechanisms underlying rose resistance to B. cinerea and also increases the database of phosphorylation and ubiquitination sites in plants.

Gene-flow investigation between garden and wild roses planted in close distance

Rose is a major ornamental plant, and a lot of cultivars with attractive morphology, color and scent have been generated by classical breeding. Recent progress of genetic modification produces a novel cultivar with attractive features. In both cases, a major problem is the gene-flow from cultivated or genetically modified (GM) plants to wild species, causing reduction of natural population. To investigate whether gene-flow occurs in wild species, molecular analysis with DNA markers with higher efficient technique is useful. Here we investigated the gene-flow from cultivated roses (Rosa×hybrida) to wild rose species planted in close distance in the field. The overlapping flowering periods and visiting insects suggest that pollens were transported by insects between wild and cultivated roses. We examined the germination ratio of seeds from wild species, and extracted DNA and checked with KSN and APETALA2 (AP2) DNA markers to detect transposon insertions. Using two markers, we successfully detected the outcross between wild and cultivated roses. For higher efficiency, we established a bulking method, where DNA, leaves or embryos were pooled, enabling us to that check the outcross of many plants. Our results suggest that wild species and garden cultivars can cross in close distance, so that they should be planted in distance, and checked the outcross with multiple DNA markers.

Comparison of nocturnal and diurnal metabolomes of rose flowers and leaves

INTRODUCTION

Roses are one of the most essential ornamental flowers and are commonly used in perfumery, cosmetics, and food. They are rich in bioactive compounds, which are of interest for therapeutic effects.

OBJECTIVES

The objective of this study was to understand the kinds of changes that occur between the nocturnal and diurnal metabolism of rose and to suggest hypotheses.

METHODS

Reversed-phase ultrahigh-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry or triple quadrupole mass spectrometry (TQ MS/MS) was used for nontargeted metabolomics and hormonal profiling respectively. For metabolite annotation, accurate mass spectra were compared with those in databases.

RESULTS

The hormonal profile of flowers showed an increase in jasmonate at night, while that of leaves indicated an increase in the salicylic acid pathway. Nontargeted analyses of the flower revealed a switch in the plant's defense mechanisms from glycosylated metabolites during the day to acid metabolites at night. In leaves, a significant decrease in flavonoids was observed at night in favor of acid metabolism to maintain a level of protection. Moreover, it might be possible to place back some of the annotated molecules on the shikimate pathway.

CONCLUSION

The influence of day and night on the metabolome of rose flowers and leaves has been clearly demonstrated. The hormonal modulations occurring during the night and at day are consistent with the plant circadian cycle. A proposed management of the sesquiterpenoid and triterpenoid biosynthetic pathway may explain these changes in the flower. In leaves, the metabolic differences may reflect night-time regulation in favor of the salicylic acid pathway.

Robust markers associated with floral traits in roses are suitable for marker-assisted selection across gene pools

We investigated the potential of markers associated with floral traits for parental selection in a cut rose breeding program. We analysed six Kompetitive Allele Specific PCR (KASP) markers for three important floral traits, petal length, petal number and scent, derived from experiments in a garden rose population. The six markers were applied to genotype a collection of 384 parental genotypes used for commercial cut rose breeding. We phenotyped a selection of progeny derived from pairs of parents having either high or low dosages of (contrasting) marker alleles associated with these traits. Significant differences were found between the contrasting progeny groups for each of the traits, although parents with the optimal allele dosage combinations could not always be used for the crosses. This not only supports the robustness of these marker‒trait associations but also demonstrates their potential for commercial rose breeding. It also demonstrates the use of marker information generated in garden rose populations for cut rose breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01438-5.

A High-Quality Reference Genome Assembly of Prinsepia uniflora (Rosaceae)

This study introduces a meticulously constructed genome assembly at the chromosome level for the Rosaceae family species Prinsepia uniflora, a traditional Chinese medicinal herb. The final assembly encompasses 1272.71 megabases (Mb) distributed across 16 pseudochromosomes, boasting contig and super-scaffold N50 values of 2.77 and 79.32 Mb, respectively. Annotated within this genome is a substantial 875.99 Mb of repetitive sequences, with transposable elements occupying 777.28 Mb, constituting 61.07% of the entire genome. Our predictive efforts identified 49,261 protein-coding genes within the repeat-masked assembly, with 45,256 (91.87%) having functional annotations, 5127 (10.41%) demonstrating tandem duplication, and 2373 (4.82%) classified as transcription factor genes. Additionally, our investigation unveiled 3080 non-coding RNAs spanning 0.51 Mb of the genome sequences. According to our evolutionary study, P. uniflora underwent recent whole-genome duplication following its separation from Prunus salicina. The presented reference-level genome assembly and annotation for P. uniflora will significantly facilitate the in-depth exploration of genomic information pertaining to this species, offering substantial utility in comparative genomics and evolutionary analyses involving Rosaceae species.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

Designing of future ornamental crops: a biotechnological driven perspective

With a basis in human appreciation of beauty and aesthetic values, the new era of ornamental crops is based on implementing innovative technologies and transforming symbols into tangible assets. Recent advances in plant biotechnology have attracted considerable scientific and industrial interest, particularly in terms of modifying desired plant traits and developing future ornamental crops. By utilizing omics approaches, genomic data, genetic engineering, and gene editing tools, scientists have successively explored the underlying molecular mechanism and potential gene(s) behind trait regulation such as floral induction, plant architecture, stress resistance, plasticity, adaptation, and phytoremediation in ornamental crop species. These signs of progress lay a theoretical and practical foundation for designing and enhancing the efficiency of ornamental plants for a wide range of applications. In this review, we briefly summarized the existing literature and advances in biotechnological approaches for the improvement of vital traits in ornamental plants. The future ornamental plants, such as light-emitting plants, biotic/abiotic stress detectors, and pollution abatement, and the introduction of new ornamental varieties via domestication of wild species are also discussed.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

GWAS of adventitious root formation in roses identifies a putative phosphoinositide phosphatase (SAC9) for marker-assisted selection

Rose propagation by cuttings is limited by substantial genotypic differences in adventitious root formation. To identify possible genetic factors causing these differences and to develop a marker for marker-assisted selection for high rooting ability, we phenotyped 95 cut and 95 garden rose genotypes in a hydroponic rooting system over 6 weeks. Data on rooting percentage after 3 to 6 weeks, root number, and root fresh mass were highly variable among genotypes and used in association mappings performed on genotypic information from the WagRhSNP 68 K Axiom SNP array for roses. GWAS analyses revealed only one significantly associated SNP for rooting percentage after 3 weeks. Nevertheless, prominent genomic regions/peaks were observed and further analysed for rooting percentage after 6 weeks, root number and root fresh mass. Some of the SNPs in these peak regions were associated with large effects on adventitious root formation traits. Very prominent were ten SNPs, which were all located in a putative phosphoinositide phosphatase SAC9 on chromosome 2 and showed very high effects on rooting percentage after 6 weeks of more than 40% difference between nulliplex and quadruplex genotypes. SAC9 was reported to be involved in the regulation of endocytosis and in combination with other members of the SAC gene family to regulate the translocation of auxin-efflux PIN proteins via the dephosphorylation of phosphoinositides. For one SNP within SAC9, a KASP marker was successfully derived and used to select genotypes with a homozygous allele configuration. Phenotyping these homozygous genotypes for adventitious root formation verified the SNP allele dosage effect on rooting. Hence, the presented KASP derived from a SNP located in SAC9 can be used for marker-assisted selection in breeding programs for high rooting ability in the future.

Pedigree-based QTL analysis of flower size traits in two multi-parental diploid rose populations

Rose (Rosa spp.) is one of the most economically important ornamental species worldwide. Flower diameter, flower weight, and the number of petals and petaloids are key flower-size parameters and attractive targets for DNA-informed breeding. Pedigree-based analysis (PBA) using FlexQTL software was conducted using two sets of multi-parental diploid rose populations. Phenotypic data for flower diameter (Diam), flower weight (fresh (FWT)/dry (DWT)), number of petals (NP), and number of petaloids (PD) were collected over six environments (seasons) at two locations in Texas. The objectives of this study were to 1) identify new and/or validate previously reported QTL(s); 2) identify SNP haplotypes associated with QTL alleles (Q-/q-) of a trait and their sources; and 3) determine QTL genotypes for important rose breeding parents. Several new and previously reported QTLs for NP and Diam traits were identified. In addition, QTLs associated with flower weight and PD were identified for the first time. Two major QTLs with large effects were mapped for all traits. The first QTL was at the distal end of LG1 (60.44-60.95 Mbp) and was associated with Diam and DWT in the TX2WOB populations. The second QTL was consistently mapped in the middle region on LG3 (30.15-39.34 Mbp) and associated with NP, PD, and flower weight across two multi-parent populations (TX2WOB and TX2WSE). Haplotype results revealed a series of QTL alleles with differing effects at important loci for most traits. This work is distinct from previous studies by conducting co-factor analysis to account for the DOUBLE FLOWER locus while mapping QTL for NP. Sources of high-value (Q) alleles were identified, namely, 'Old Blush' and Rosa wichuraiana from J14-3 for Diam, while 'Violette' and PP-J14-3 were sources for other traits. In addition, the source of the low-value (q) alleles for Diam was 'Little Chief', and Rosa wichuraiana through J14-3 was the source for the remaining traits. Hence, our results can potentially inform parental/seedling selections as means to improve ornamental quality in roses and a step towards implementing DNA-informed techniques for use in rose breeding programs.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

RcSPL1-RcTAF15b regulates the flowering time of rose (Rosa chinensis)

Rose (Rosa chinensis), which is an economically valuable floral species worldwide, has three types, namely once-flowering (OF), occasional or re-blooming (OR), and recurrent or continuous flowering (CF). However, the mechanism underlying the effect of the age pathway on the duration of the CF or OF juvenile phase is largely unknown. In this study, we observed that the RcSPL1 transcript levels were substantially upregulated during the floral development period in CF and OF plants. Additionally, accumulation of RcSPL1 protein was controlled by rch-miR156. The ectopic expression of RcSPL1 in Arabidopsis thaliana accelerated the vegetative phase transition and flowering. Furthermore, the transient overexpression of RcSPL1 in rose plants accelerated flowering, whereas silencing of RcSPL1 had the opposite phenotype. Accordingly, the transcription levels of floral meristem identity genes (APETALA1, FRUITFULL, and LEAFY) were significantly affected by the changes in RcSPL1 expression. RcTAF15b protein, which is an autonomous pathway protein, was revealed to interact with RcSPL1. The silencing and overexpression of RcTAF15b in rose plants led to delayed and accelerated flowering, respectively. Collectively, the study findings imply that RcSPL1-RcTAF15b modulates the flowering time of rose plants.

A cytosolic bifunctional geranyl/farnesyl diphosphate synthase provides MVA-derived GPP for geraniol biosynthesis in rose flowers

Geraniol derived from essential oils of various plant species is widely used in the cosmetic and perfume industries. It is also an essential trait of the pleasant smell of rose flowers. In contrast to other monoterpenes which are produced in plastids via the methyl erythritol phosphate pathway, geraniol biosynthesis in roses relies on cytosolic NUDX1 hydrolase which dephosphorylates geranyl diphosphate (GPP). However, the metabolic origin of cytosolic GPP remains unknown. By feeding Rosa chinensis "Old Blush" flowers with pathway-specific precursors and inhibitors, combined with metabolic profiling and functional characterization of enzymes in vitro and in planta, we show that geraniol is synthesized through the cytosolic mevalonate (MVA) pathway by a bifunctional geranyl/farnesyl diphosphate synthase, RcG/FPPS1, producing both GPP and farnesyl diphosphate (FPP). The downregulation and overexpression of RcG/FPPS1 in rose petals affected not only geraniol and germacrene D emissions but also dihydro-β-ionol, the latter due to metabolic cross talk of RcG/FPPS1-dependent isoprenoid intermediates trafficking from the cytosol to plastids. Phylogenetic analysis together with functional characterization of G/FPPS orthologs revealed that the G/FPPS activity is conserved among Rosaceae species. Site-directed mutagenesis and molecular dynamic simulations enabled to identify two conserved amino acids that evolved from ancestral FPPSs and contribute to GPP/FPP product specificity. Overall, this study elucidates the origin of the cytosolic GPP for NUDX1-dependent geraniol production, provides insights into the emergence of the RcG/FPPS1 GPPS activity from the ancestral FPPSs, and shows that RcG/FPPS1 plays a key role in the biosynthesis of volatile terpenoid compounds in rose flowers.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Meta-Analysis of Rose Rosette Disease-Resistant Quantitative Trait Loci and a Search for Candidate Genes

Rose rosette disease (RRD), caused by the rose rosette emaravirus (RRV), is a major viral disease in roses (Rosa sp.) that threatens the rose industry. Recent studies have revealed quantitative trait loci (QTL) for reduced susceptibility to RRD in the linkage groups (LGs) 1, 5, 6, and 7 in tetraploid populations and the LGs 1, 3, 5, and 6 in diploid populations. In this study, we seek to better localize and understand the relationship between QTL identified in both diploid and tetraploid populations. We do so by remapping the populations found in these studies and performing a meta-analysis. This analysis reveals that the peaks and intervals for QTL using diploid and tetraploid populations co-localized on LG 1, suggesting that these are the same QTL. The same was seen on LG 3. Three meta-QTL were identified on LG 5, and two were discovered on LG 6. The meta-QTL on LG 1, MetaRRD1.1, had a confidence interval (CI) of 10.53 cM. On LG 3, MetaRRD3.1 had a CI of 5.94 cM. MetaRRD5.1 had a CI of 17.37 cM, MetaRRD5.2 had a CI of 4.33 cM, and MetaRRD5.3 had a CI of 21.95 cM. For LG 6, MetaRRD6.1 and MetaRRD6.2 had CIs of 9.81 and 8.81 cM, respectively. The analysis also led to the identification of potential disease resistance genes, with a primary interest in genes localized in meta-QTL intervals on LG 5 as this LG was found to explain the greatest proportion of phenotypic variance for RRD resistance. The results from this study may be used in the design of more robust marker-based selection tools to track and use a given QTL in a plant breeding context.