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

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.

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.

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.

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.

An efficient CRISPR/Cas9 platform for targeted genome editing in rose (Rosa hybrida)

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-related nuclease 9 (Cas9) system enables precise, simple editing of genes in many animals and plants. However, this system has not been applied to rose (Rosa hybrida) due to the genomic complexity and lack of an efficient transformation technology for this plant. Here, we established a platform for screening single-guide RNAs (sgRNAs) with high editing efficiency for CRISPR/Cas9-mediated gene editing in rose using suspension cells. We used the Arabidopsis thaliana U6-29 promoter, which showed high activity for driving sgRNA expression, to modify the CRISPR/Cas9 system. We used our highly efficient optimized CRISPR/Cas9 system to successfully edit RhEIN2, encoding an indispensable component of the ethylene signaling pathway, resulting in ethylene insensitivity in rose. Our optimized CRISPR/Cas9 system provides a powerful toolbox for functional genomics, molecular breeding, and synthetic biology in rose.

Genome-Wide Identification of the Rose SWEET Gene Family and Their Different Expression Profiles in Cold Response between Two Rose Species

Sugars Will Eventually be Exported Transporter (SWEET) gene family plays indispensable roles in plant physiological activities, development processes, and responses to biotic and abiotic stresses, but no information is known for roses. In this study, a total of 25 RcSWEET genes were identified in Rosa chinensis 'Old Blush' by genome-wide analysis and clustered into four subgroups based on their phylogenetic relationships. The genomic features, including gene structures, conserved motifs, and gene duplication among the chromosomes of RcSWEET genes, were characterized. Seventeen types of cis-acting elements among the RcSWEET genes were predicted to exhibit their potential regulatory roles during biotic and abiotic stress and hormone responses. Tissue-specific and cold-response expression profiles based on transcriptome data showed that SWEETs play widely varying roles in development and stress tolerance in two rose species. Moreover, the different expression patterns of cold-response SWEET genes were verified by qRT-PCR between the moderately cold-resistant species R. chinensis 'Old Blush' and the extremely cold-resistant species R. beggeriana. Especially, SWEET2a and SWEET10c exhibited species differences after cold treatment and were sharply upregulated in the leaves of R. beggeriana but not R. chinensis 'Old Blush', indicating that these two genes may be the crucial candidates that participate in cold tolerance in R. beggeriana. Our results provide the foundation for function analysis of the SWEET gene family in roses, and will contribute to the breeding of cold-tolerant varieties of roses.

Genome-Wide Identification and Expression Analysis of the TIR-NBS-LRR Gene Family and Its Response to Fungal Disease in Rose (Rosa chinensis)

Roses, which are one of the world’s most important ornamental plants, are often damaged by pathogens, resulting in serious economic losses. As a subclass of the disease resistance gene family of plant nucleotide-binding oligomerization domain (NOD)-like receptors, TIR-NBS-LRR (TNL) genes play a vital role in identifying pathogen effectors and activating defense responses. However, a systematic analysis of the TNL gene family is rarely reported in roses. Herein, 96 intact TNL genes were identified in Rosa chinensis. Their phylogenies, physicochemical characteristics, gene structures, conserved domains and motifs, promoter cis-elements, microRNA binding sites, and intra- and interspecific collinearity relationships were analyzed. An expression analysis using transcriptome data revealed that RcTNL genes were dominantly expressed in leaves. Some RcTNL genes responded to gibberellin, jasmonic acid, salicylic acid, Botrytis cinerea, Podosphaera pannosa, and Marssonina rosae (M. rosae); the RcTNL23 gene responded significantly to three hormones and three pathogens, and exhibited an upregulated expression. Furthermore, the black spot pathogen was identified as M. rosae. After inoculating rose leaves, an expression pattern analysis of the RcTNL genes suggested that they act during different periods of pathogen infection. The present study lays the foundations for an in-depth investigation of the TNL gene function and the mining of disease resistance genes in roses.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

High Expression of ALTERNATIVE OXIDASE 2 in Latent Axillary Buds Suggests Its Key Role in Quiescence Maintenance in Rosebush

Most vegetative axes remain quiescent as dormant axillary buds until metabolic and hormonal signals, driven by environmental changes, trigger bud outgrowth. While the resumption of growth activity is well documented, the establishment and maintenance of quiescence is comparatively poorly understood, despite its major importance in the adaptation of plants to the seasonal cycle or in the establishment of their shape. Here, using the rosebush Rosa hybrida 'Radrazz' as a plant model, we highlighted that the quiescent state was the consequence of an internal and active energy control of buds, under the influence of hormonal factors previously identified in the bud outgrowth process. We found that the quiescent state in the non-growing vegetative axis of dormant axillary buds displayed a low energy state along with a high expression of the ALTERNATIVE OXIDASE 2 (AOX2) and the accumulation of the corresponding protein. Conversely, AOX2 expression and protein amount strongly decreased during bud burst as energy status shifted to a high state, allowing growth. Since AOX2 can deviate electrons from the cytochrome pathway in the mitochondrial respiratory chain, it could drastically reduce the formation of ATP, which would result in a low energy status unfavorable for growth activities. We provide evidence that the presence/absence of AOX2 in quiescent/growing vegetative axes of buds was under hormonal control and thus may constitute the mechanistic basis of both quiescence and sink strength manifestation, two important aspects of budbreak.

BOX38, a DNA Marker for Selection of Essential Oil Yield of Rosa × rugosa

Rosa rugosa L. was a famous aromatic plant whose cultivars (Rosa × rugosa) have been widely used in the perfume industry in Asia. The perfume market looks for rose cultivars bearing higher essential oil, while the oil yields of most R. × rugosa have not been evaluated due to limiting conditions, such as insufficient cultivation areas. Here, we tested the yield and the aroma components of essential oil of 19 R. × rugosa. The results indicated that the yields of nerol, citronellol, and geraniol could represent an alternative index of the total yield of essential oil. Sequence syntenic analysis indicated that the Rosa genus specific cis-element Box38 was highly polymorphic. The Box38 region isolation of Rosa × rugosa by flanked primers proved that Box38 repeat number was significantly positively correlated with the essential oil yield of the corresponding cultivar. In the breeding of Rosa × rugosa, six-Box38-repeat could be a robust threshold for selection of high-essential-oil roses. Together, we found that Box38 was a DNA marker for essential oil yield and that it would be helpful in the early selection and breeding of essential oil roses.

Chromosome-Level Assembly of Flowering Cherry (Prunus campanulata) Provides Insight into Anthocyanin Accumulation

The flowering cherries (genus Prunus, subgenus Cerasus) are popular ornamental trees in China, Japan, Korea, and elsewhere. Prunus campanulata Maxim. is an important species of flowering cherry native to Southern China, which is also distributed in Taiwan, the Ryukyu Islands of Japan, and Vietnam. It produces bell-shaped flowers with colors ranging from bright pink to crimson during the Chinese Spring Festival from January to March each year. We selected the P. campanulata cultivar "Lianmeiren", with only 0.54% of heterozygosity, as the focus of this study, and generated a high-quality chromosome-scale genome assembly of P. campanulata by combining Pacific Biosciences (PacBio) single-molecule sequencing, 10× Genomics sequencing, and high-throughput chromosome conformation capture (Hi-C) technology. We first assembled a 300.48 Mb genome assembly with a contig N50 length of 2.02 Mb. In total, 28,319 protein-coding genes were predicted from the genome, 95.8% of which were functionally annotated. Phylogenetic analyses indicated that P. campanulata diverged from a common ancestor of cherry approximately 15.1 million years ago. Comparative genomic analyses showed that the expanded gene families were significantly involved in ribosome biogenesis, diterpenoid biosynthesis, flavonoid biosynthesis, and circadian rhythm. Furthermore, we identified 171 MYB genes from the P. campanulata genome. Based on the RNA-seq of five organs at three flowering stages, expression analyses revealed that the majority of the MYB genes exhibited tissue-specific expression patterns, and some genes were identified as being associated with anthocyanin accumulation. This reference sequence is an important resource for further studies of floral morphology and phenology, and comparative genomics of the subgenera Cerasus and Prunus.

Multiomics analysis reveals the mechanisms underlying the different floral colors and fragrances of Rosa hybrida cultivars

The color and fragrance of rose flowers affect their commercial value. However, several rose varieties with new floral colors developed by the bud mutation method lost their fragrance during the breeding process, raising the question: Is there a relationship between floral color and aroma traits? Rose cultivar 'Yellow Island' (YI) with intensely aroma and yellow petals, while its bud mutant 'Past Feeling' (PF) with light aroma and pink petals mixing some yellow, two cultivars were used to explore this question using multiomics approaches. We investigated the genomic polymorphisms between PF and YI by whole-genome resequencing. 71 differentially abundant metabolites and 155 related differentially expressed genes identified in petals between PF and YI. From this, we constructed a model of metabolic changes affecting floral color and fragrance integrating shikimate, terpenoid, carotenoid, and green leaf volatile metabolites and predicted the associated key genes and transcription factors. This study provides a reference for understanding the molecular mechanism of variation in rose floral color and aroma traits.

Pedigree-based analysis in multi-parental diploid rose populations reveals QTLs for cercospora leaf spot disease resistance

Cercospora leaf spot (CLS) (Cercospora rosicola) is a major fungal disease of roses (Rosa sp.) in the southeastern U.S. Developing CLS-resistant cultivars offers a potential solution to reduce pesticide use. Yet, no work has been performed on CLS resistance. This study aimed to identify QTLs and to characterize alleles for resistance to CLS. The study used pedigree-based QTL analysis to dissect the genetic basis of CLS resistance using two multi-parental diploid rose populations (TX2WOB and TX2WSE) evaluated across five years in two Texas locations. A total 38 QTLs were identified across both populations and distributed over all linkage groups. Three QTLs on LG3, LG4, and LG6 were consistently mapped over multiple environments. The LG3 QTL was mapped in a region between 18.9 and 27.8 Mbp on the Rosa chinensis genome assembly. This QTL explained 13 to 25% of phenotypic variance. The LG4 QTL detected in the TX2WOB population spanned a 35.2 to 39.7 Mbp region with phenotypic variance explained (PVE) up to 48%. The LG6 QTL detected in the TX2WSE population was localized to 17.9 to 33.6 Mbp interval with PVE up to 36%. Also, this study found multiple degrees of favorable allele effects (q-allele) associated with decreasing CLS at major loci. Ancestors 'OB', 'Violette', and PP-M4-4 were sources of resistance q-alleles. These results will aid breeders in parental selection to develop CLS-resistant rose cultivars. Ultimately, high throughput DNA tests that target major loci for CLS could be developed for routine use in a DNA-informed breeding program.

Creating novel ornamentals via new strategies in the era of genome editing

Ornamental breeding has traditionally focused on improving novelty, yield, quality, and resistance to biotic or abiotic stress. However, achieving these goals has often required laborious crossbreeding, while precise breeding techniques have been underutilized. Fortunately, recent advancements in plant genome sequencing and editing technology have opened up exciting new frontiers for revolutionizing ornamental breeding. In this review, we provide an overview of the current state of ornamental transgenic breeding and propose four promising breeding strategies that have already proven successful in crop breeding and could be adapted for ornamental breeding with the help of genome editing. These strategies include recombination manipulation, haploid inducer creation, clonal seed production, and reverse breeding. We also discuss in detail the research progress, application status, and feasibility of each of these tactics.

Genomic Approaches for the Study of Flower Development in Floriculture Crops

The advances in genomics and bioinformatics have made possible the study in non-model plants of phenotypes associated to flower development. Floriculture crops are an interesting source of traits associated to flower development such as the transition between zygomorphic and actinomorphic flowers or the production of flowers with double and triple corollas. In this chapter, we summarize the material and methods for the use of floriculture crops to study flower development using genomic tools, from the sequencing and assembly of a reference genome to QTL and RNA-Seq analysis to search candidate genes associated to specific traits.

Developing best practices for genotyping-by-sequencing analysis in the construction of linkage maps

BACKGROUND

Genotyping-by-sequencing (GBS) provides affordable methods for genotyping hundreds of individuals using millions of markers. However, this challenges bioinformatic procedures that must overcome possible artifacts such as the bias generated by polymerase chain reaction duplicates and sequencing errors. Genotyping errors lead to data that deviate from what is expected from regular meiosis. This, in turn, leads to difficulties in grouping and ordering markers, resulting in inflated and incorrect linkage maps. Therefore, genotyping errors can be easily detected by linkage map quality evaluations.

RESULTS

We developed and used the Reads2Map workflow to build linkage maps with simulated and empirical GBS data of diploid outcrossing populations. The workflows run GATK, Stacks, TASSEL, and Freebayes for single-nucleotide polymorphism calling and updog, polyRAD, and SuperMASSA for genotype calling, as well as OneMap and GUSMap to build linkage maps. Using simulated data, we observed which genotype call software fails in identifying common errors in GBS sequencing data and proposed specific filters to better handle them. We tested whether it is possible to overcome errors in a linkage map using genotype probabilities from each software or global error rates to estimate genetic distances with an updated version of OneMap. We also evaluated the impact of segregation distortion, contaminant samples, and haplotype-based multiallelic markers in the final linkage maps. Through our evaluations, we observed that some of the approaches produce different results depending on the dataset (dataset dependent) and others produce consistent advantageous results among them (dataset independent).

CONCLUSIONS

We set as default in the Reads2Map workflows the approaches that showed to be dataset independent for GBS datasets according to our results. This reduces the number of required tests to identify optimal pipelines and parameters for other empirical datasets. Using Reads2Map, users can select the pipeline and parameters that best fit their data context. The Reads2MapApp shiny app provides a graphical representation of the results to facilitate their interpretation.

Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome

High-quality genome sequences help to elucidate the genetic basis of numerous biological processes and track species evolution. For flax (Linum usitatissimum L.)—a multifunctional crop, high-quality assemblies from Oxford Nanopore Technologies (ONT) data were unavailable, largely due to the difficulty of isolating pure high-molecular-weight DNA. This article proposes a scheme for gaining a contiguous L. usitatissimum assembly using Nanopore data. We developed a protocol for flax nuclei isolation with subsequent DNA extraction, which allows obtaining about 5 μg of pure high-molecular-weight DNA from 0.5 g of leaves. Such an amount of material can be collected even from a single plant and yields more than 30 Gb of ONT data in two MinION runs. We performed a comparative analysis of different genome assemblers and polishers on the gained data and obtained the final 447.1-Mb assembly of L. usitatissimum line 3896 genome using the Canu—Racon (two iterations)—Medaka combination. The genome comprised 1695 contigs and had an N50 of 6.2 Mb and a completeness of 93.8% of BUSCOs from eudicots_odb10. Our study highlights the impact of the chosen genome construction strategy on the resulting assembly parameters and its eligibility for future genomic studies.

Genome-wide identification of apple PPI genes and a functional analysis of the response of MxPPI1 to Fe deficiency stress

Iron (Fe) deficiency affects plant growth and development. The proton pump interactor (PPI) in plants responds to multiple abiotic stresses, although it has not been well characterized under Fe deficiency stress. In this study, we systematically identified and analyzed the PPI gene family in apple. Three PPI candidate genes were found, and they contained 318-1349 amino acids and 3-7 introns. Under Fe deficiency stress, we analyzed the expression of all the PPI genes in roots of apple rootstock Malus xiaojinensis. Expression of the gene MD11G1247800, designated PPI1, is obviously induced by Fe deficiency treatment in M. xiaojinensis. We first cloned MxPPI1 from M. xiaojinensis and determined its subcellular localization, which indicated that it is localized in the cell membrane and nucleus in tobacco. We found that the level of expression of the MxPPI1 protein increased significantly under Fe deficiency stress in apple calli. Moreover, overexpressing MxPPI1 in apple calli enhanced the activities of ferric chelate reductase and H⁺-ATPase, H⁺ secretion, MxHA2 gene expression and total Fe content when compared with the wild type calli. We further found that MxPPI1 interacted with MxHA2 using bimolecular fluorescence complementation and luciferase complementation assays. Overall, we demonstrated that MxPPI1 interacts with MxHA2 to enhance the activity of H⁺-ATPase to regulate Fe absorption in M. xiaojinensis.

Characterization and expression analysis of WRKY genes during leaf and corolla senescence of Petunia hybrida plants

Several families of transcription factors (TFs) control the progression of senescence. Many key TFs belonging to the WRKY family have been described to play crucial roles in the regulation of leaf senescence, mainly in Arabidopsis thaliana. However, little is known about senescence-associated WRKY members in floricultural species. Delay of senescence in leaves and petals of Petunia hybrida, a worldwide ornamental crop are highly appreciated traits. In this work, starting from 28 differentially expressed WRKY genes of A. thaliana during the progression of leaf senescence, we identified the orthologous in P. hybrida and explored the expression profiles of 20 PhWRKY genes during the progression of natural (age-related) leaf and corolla senescence as well as in the corollas of flowers undergoing pollination-induced senescence. Simultaneous visualization showed consistent and similar expression profiles of PhWRKYs during natural leaf and corolla senescence, although weak expression changes were observed during pollination-induced senescence. Comparable expression trends between PhWRKYs and the corresponding genes of A. thaliana were observed during leaf senescence, although more divergence was found in petals of pollinated petunia flowers. Integration of expression data with phylogenetics, conserved motif and cis-regulatory element analyses were used to establish a list of candidates that could regulate more than one senescence process. Our results suggest that several members of the WRKY family of TFs are tightly linked to the regulation of senescence in P. hybrida.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01243-y.

RrTTG1 promotes fruit prickle development through an MBW complex in Rosa roxburghii

Fruit prickles are widely distributed on the pericarp and exhibit polymorphic traits at different developmental stages. Although they are multicellular appendages that are well-known for helping plants defend against biotic and abiotic stresses, their origination and molecular mechanism are still less known. Here, we studied the origination and molecular mechanism of fruit prickles in Rosa roxburghii. Using morphological and histological observations, we found that the fruit prickle primordium of R. roxburghii originated from the ground meristem that underwent cell division to form flagelliform prickles, continued to enlarge, and finally lignified to form mature fruit prickles. We amplified a homolog of candidate gene TRANSPARENT TESTA GLABRA1 (TTG1) from R. roxburghii, named RrTTG1. RrTTG1 harbored four conserved WD-repeat domains and was exclusively nuclear-localized. Using qRT-PCR and in situ hybridization, we found that RrTTG1 was constitutively expressed and highly expressed during the initiation and cell expansion phases of fruit prickles. Ectopic expression analysis in Arabidopsis proved that RrTTG1 substantially enhanced the number of trichome and pigmentation production and inhibited root hair formation. Besides, RrTTG1 complemented the phenotypes of the ttg1 mutant in Arabidopsis, thus indicating that RrTTG1 played pleiotropic roles akin to AtTTG1. We demonstrated that the RrTTG1 only interacted with RrEGL3, a homolog of ENHANCER OF GLABRA3 (EGL3), via yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays. Briefly, RrTTG1 might positively regulate the initiation of fruit prickle primordium and cell enlargement by forming the RrTTG1-RrEGL3-RrGL1 complex in R. roxburghii. Therefore, our results help characterize the RrTTG1 in R. roxburghii and also elucidate the establishment of the prickles regulatory system in the Rosaceae plants.