Transcriptomic and metabolomic changes triggered by Macrosiphum rosivorum in rose (Rosa longicuspis).

BACKGROUND: Rose is one of the most popular flowers in the wold. Its field growth and quality are negatively affected by aphids. However, the defence mechanisms used by rose plants against aphids are unclear. Therefore, to understand the defence mechanism of rose under aphid stress, transcriptome and metabolome techniques were used to investigate the regulation mechanism in R. longicuspis infected with M. rosivorum. RESULT: In our study, after inoculation with M. rosivorum, M. rosivorum quickly colonized R. longicuspis. A total of 34,202 genes and 758 metabolites were detected in all samples. Under M. rosivorum stress, R. longicuspis responded by MAPK cascades, plant hormone signal transduction pathway activation, RlMYBs and RlERFs transcription factors expression and ROS production. Interestingly, the 'brassinosteroid biosynthesis' pathway was significantly enriched in A3 d-vs.-A5 d. Further analysis showed that M. rosivorum induced the biosynthesis of secondary metabolites such as terpenoids, tannins and phenolic acids, among others. Importantly, the 'glutathione metabolic' and 'glucosinolate biosynthesis' pathways were significantly enriched, which involved in the rose against aphids. CONCLUSION: Our study provides candidate genes and metabolites for Rosa defence against aphids. This study provides a theoretical basis for further exploring the molecular regulation mechanism of rose aphid resistance and aphid resistance breeding in the future.

RcTGA1 and glucosinolate biosynthesis pathway involvement in the defence of rose against the necrotrophic fungus Botrytis cinerea.

BACKGROUND: Rose is an important economic crop in horticulture. However, its field growth and postharvest quality are negatively affected by grey mould disease caused by Botrytis c. However, it is unclear how rose plants defend themselves against this fungal pathogen. Here, we used transcriptomic, metabolomic and VIGS analyses to explore the mechanism of resistance to Botrytis c. RESULT: In this study, a protein activity analysis revealed a significant increase in defence enzyme activities in infected plants. RNA-Seq of plants infected for 0 h, 36 h, 60 h and 72 h produced a total of 54 GB of clean reads. Among these reads, 3990, 5995 and 8683 differentially expressed genes (DEGs) were found in CK vs. T36, CK vs. T60 and CK vs. T72, respectively. Gene annotation and cluster analysis of the DEGs revealed a variety of defence responses to Botrytis c. infection, including resistance (R) proteins, MAPK cascade reactions, plant hormone signal transduction pathways, plant-pathogen interaction pathways, Ca2+ and disease resistance-related genes. qPCR verification showed the reliability of the transcriptome data. The PTRV2-RcTGA1-infected plant material showed improved susceptibility of rose to Botrytis c. A total of 635 metabolites were detected in all samples, which could be divided into 29 groups. Metabonomic data showed that a total of 59, 78 and 74 DEMs were obtained for T36, T60 and T72 (T36: Botrytis c. inoculated rose flowers at 36 h; T60: Botrytis c. inoculated rose flowers at 60 h; T72: Botrytis c. inoculated rose flowers at 72 h) compared to CK, respectively. A variety of secondary metabolites are related to biological disease resistance, including tannins, amino acids and derivatives, and alkaloids, among others; they were significantly increased and enriched in phenylpropanoid biosynthesis, glucosinolates and other disease resistance pathways. This study provides a theoretical basis for breeding new cultivars that are resistant to Botrytis c. CONCLUSION: Fifty-four GB of clean reads were generated through RNA-Seq. R proteins, ROS signalling, Ca2+ signalling, MAPK signalling, and SA signalling were activated in the Old Blush response to Botrytis c. RcTGA1 positively regulates rose resistance to Botrytis c. A total of 635 metabolites were detected in all samples. DEMs were enriched in phenylpropanoid biosynthesis, glucosinolates and other disease resistance pathways.

The Complete Chloroplast Genome of a Key Ancestor of Modern Roses, Rosa chinensis var. spontanea, and a Comparison with Congeneric Species.

Rosa chinensis var. spontanea, an endemic and endangered plant of China, is one of the key ancestors of modern roses and a source for famous traditional Chinese medicines against female diseases, such as irregular menses and dysmenorrhea. In this study, the complete chloroplast (cp) genome of R. chinensis var. spontanea was sequenced, analyzed, and compared to congeneric species. The cp genome of R. chinensis var. spontanea is a typical quadripartite circular molecule of 156,590 bp in length, including one large single copy (LSC) region of 85,910 bp and one small single copy (SSC) region of 18,762 bp, separated by two inverted repeat (IR) regions of 25,959 bp. The GC content of the whole genome is 37.2%, while that of LSC, SSC, and IR is 42.8%, 35.2% and 31.2%, respectively. The genome encodes 129 genes, including 84 protein-coding genes (PCGs), 37 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. Seventeen genes in the IR regions were found to be duplicated. Thirty-three forward and five inverted repeats were detected in the cp genome of R. chinensis var. spontanea. The genome is rich in SSRs. In total, 85 SSRs were detected. A genome comparison revealed that IR contraction might be the reason for the relatively smaller cp genome size of R. chinensis var. spontanea compared to other congeneric species. Sequence analysis revealed that the LSC and SSC regions were more divergent than the IR regions within the genus Rosa and that a higher divergence occurred in non-coding regions than in coding regions. A phylogenetic analysis showed that the sampled species of the genus Rosa formed a monophyletic clade and that R. chinensis var. spontanea shared a more recent ancestor with R. lichiangensis of the section Synstylae than with R. odorata var. gigantea of the section Chinenses. This information will be useful for the conservation genetics of R. chinensis var. spontanea and for the phylogenetic study of the genus Rosa, and it might also facilitate the genetics and breeding of modern roses.

Graft-accelerated virus-induced gene silencing facilitates functional genomics in rose flowers.

Rose has emerged as a model ornamental plant for studies of flower development, senescence, and morphology, as well as the metabolism of floral fragrances and colors. Virus-induced gene silencing (VIGS) has long been used in functional genomics studies of rose by vacuum infiltration of cuttings or seedlings with an Agrobacterium suspension carrying TRV-derived vectors. However, VIGS in rose flowers remains a challenge because of its low efficiency and long time to establish silencing. Here we present a novel and rapid VIGS method that can be used to analyze gene function in rose, called 'graft-accelerated VIGS', where axillary sprouts are cut from the rose plant and vacuum infiltrated with Agrobacterium. The inoculated scions are then grafted back onto the plants to flower and silencing phenotypes can be observed within 5 weeks, post-infiltration. Using this new method, we successfully silenced expression of the RhDFR1, RhAG, and RhNUDX1 in rose flowers, and affected their color, petal number, as well as fragrance, respectively. This grafting method will facilitate high-throughput functional analysis of genes in rose flowers. Importantly, it may also be applied to other woody species that are not currently amenable to VIGS by conventional leaf or plantlet/seedling infiltration methods.

Identification and Validation of Reference Genes for qRT-PCR Analysis of Petal-Color-Related Genes in Rosa praelucens.

The flower's color is regarded as one of the most outstanding features of the rose. Rosa praelucens Byhouwer, an endemic and critically endangered decaploid wild rose species, is abundant in phenotypic diversity, especially in flower color variation, from white to different degrees of pink. The mechanism underlying this variation, e.g., the level of petal-color-related genes, is worth probing. Seven candidate reference genes for qRT-PCR analysis, including tubulin α chain (TUBA), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone H2B (Histone2A), eukaryotic translation elongation factor 1-α (EEF1A), 60S ribosomal protein (RPL37), eukaryotic translation initiation factor 1-α (EIF1A), and aquaporins (AQP), were detected from the transcriptome datasets of full blooming flowers of white-petaled and pink-petaled individuals, and their expression stabilities were evaluated through qRT-PCR analysis. According to stability rankings analysis, EEF1A showed the highest stability and could be chosen as the most suitable reference gene. Moreover, the reliability of EEF1A was demonstrated via qRT-PCR analysis of six petal-color-related target genes, the expression patterns of which, through EEF1A normalization, were found to be consistent with the findings of transcriptome analysis. The result provides an optimal reference gene for exploring the expression level of petal-color-related genes in R. praelucens, which will accelerate the dissection of petal-color-variation mechanisms in R. praelucens.

Simulation of climate change effect on the global distribution of Rosa multiflora.

Rosa multiflora, originated from East Asia, is one of the original ancestors of modern roses. It is also an important genetic resource and rootstock for rose cultivation. Due to its high resistance and vigorous growth, R. multiflora has become an invasive species in some introduction sites, such as North America. To explore the correlation between the suitable habitat of R. multiflora and climate change, we predicted its potential geographic distribution with an optimized MaxEnt model based on 1246 distribution records and nine bioclimatic variables. The results showed that the mean temperature of the coldest quarter, minimum temperature of the coldest month, precipitation of the warmest quarter, and isothermality were significant bioclimatic variables affecting the potential geographic distribution of R. multiflora. Under current climate conditions, R. multiflora naturally distributed in the plains and hilly areas to the east and south of the Loess Plateau. The distribution pattern in the mid-holocene was similar to its current distribution, but the highly suitable distribution area was in the south of North China Plain, the Sichuan Basin, and parts of the Middle-Lower Yangtze Plain. During the last interglacial, the suitable areas generally contrac-ted southward, while the highly suitable areas significantly expanded and mainly located in the Sichuan Basin, the Middle-Lower Yangtze Plains, the Yunnan-Guizhou Plateau, and the Southeast Hills. Beyond its natural distribution in East Asia, R. multiflora had been introduced and spread to most parts of Europe and the central and eastern United States. The distribution area of R. multiflora would expand under three warming scenarios of different shared socioeconomic pathways (SSP1-2.6, SSP2-4.5, and SSP5-8.5) during 2041-2060 and 2081-2100. Its average distribution center (centroid) would shift towards higher latitude, indicating that the distribution of R. multiflora was closely related to climate change and that global warming might lead to an expansion of its distribution area. These results would improve our understanding of the ecological adaptability of R. multiflora, facilitate the predicting of its future distribution, and provide a theoretical basis for monitoring and early warning measures following its introduction.

Loss of Rose Fragrance under Chilling Stress Is Associated with Changes in DNA Methylation and Volatile Biosynthesis.

Rose plants are widely cultivated as cut flowers worldwide and have economic value as sources of natural fragrance and flavoring. Rosa 'Crimson Glory', whose petals have a pleasant fragrance, is one of the most important cultivars of edible rose plants. Flower storage at low-temperature is widely applied in production to maintain quality; however, chilling results in a decrease in aromatic volatiles. To determine the molecular basis underlying the changes in aromatic volatile emissions, we investigated the changes in volatile compounds, DNA methylation patterns, and patterns of the transcriptome in response to chilling temperature. The results demonstrated that chilling roses substantially reduced aromatic volatile emissions. We found that these reductions were correlated with the changes in the methylation status of the promoters and genic regions of the genes involved in volatile biosynthesis. These changes mainly occurred for CHH (H = A, T, or C) which accounted for 51% of the total methylation. Furthermore, transcript levels of scent-related gene Germacrene D synthase (RhGDS), Nudix hydrolase 1 (RhNUDX1), and Phenylacetaldehyde reductase (RhPAR) of roses were strikingly depressed after 24 h at low-temperature and remained low-level after 24 h of recovery at 20 °C. Overall, our findings indicated that epigenetic regulation plays an important role in the chilling tolerance of roses and lays a foundation for practical significance in the production of edible roses.

Phylogeography and Population Genetics of Rosa chinensis var. spontanea and R. lucidissima Complex, the Important Ancestor of Modern Roses.

Rosa chinensis var. spontanea and R. lucidissima complex are the morphologically very similar key ancestors of modern roses with high importance in rose research and breeding. Although widely distributed in subtropical central and southwestern China, these two taxa are highly endangered. We sampled a total of 221 specimens and 330 DNA samples from 25 populations across the two taxa's whole range. Leaf morphological traits were compared. Two chloroplast DNA intergenic spacers (trnG-trnS, petL-psbE) and ITS were used for population genetics and phylogenetic study to delimit the boundary between the two taxa, assess the genetic variation, uncover the possible evolutionary mechanism responsible for the differentiation within the complex, and make the conservation recommendations. The complex exhibited high levels of genetic variation (h TcpDNA = 0.768, h TITS = 0.726) and high population differentiation even over small geographic distance. We suggest R. chinensis var. spontanea and R. lucidissma be treated as independent taxa, and the northern populations around and within the Sichuan Basin being R. chinensis var. spontanea, having broader leaflets and paler full-blooming flowers, while those in the middle and southern Yunnan-Guizhou Plateau and the adjacent regions being R. lucidissma, having narrower leaflets and darker full-blooming flowers. Transitional areas between the southeastern Sichuan Basin and northeastern Guizhou are the contact or the hybridization zone of the two taxa. Ancestral haplotypes of the complex (R. lucidissma) evolved at about 1.21-0.86 Mya in southeastern Yunnan-Guizhou Plateau and its adjacent regions and survived there during the Quaternary Oscillation. Ancestral haplotypes of R. chinensis var. spontanea deviated from R. lucidissma at about 0.022-0.031 Mya at the transitional areas (Daloushan and Wulingshan Mountains) between the northeastern edge of Yunnan-Guizhou Plaeteau and the southeastern border of Sichuan Basin, where they survived the LGM. The evolution of the complex included spatial isolation and inter-species hybridization. The complex's endangered status might be the result of over-exploitation for its ornamental and medical value, or due to reforestation of some originally open habitats. We provide specific recommendations for the two taxa's in situ and ex situ conservation.

Isolation and identification of a putative scent-related gene RhMYB1 from rose.

Rose fragrances play an important role in attracting pollinators and commercial value. However, some genes involved in rose floral scent metabolism are less well understood. Here, wild-type scented rose (WR) and its spontaneous non-scented mutant rose (MR) were analyzed. SPME-GC/MS analysis showed that relative content of 1-ethenyl-4-methoxy-benzene represented was significantly different between WR and MR. We have isolated an EST encoding a MYB family of transcription factor from SSH libraries of the two roses in the previous studies, and designated RhMYB1. In the study, the full-length cDNA of RhMYB1 was identified and characterized by rapid amplification of cDNA ends (RACE). The RhMYB1 full-length cDNA was 1,125 bp containing an 882 bp open reading frame, which encodes a precursor protein of 294 amino acids. Sequence alignments revealed that RhMYB1 shared high similarity with other plants R-type MYB, and RhMYB1 contained a DNA binding domain. Northern blot analysis revealed that RhMYB1 was expressed specifically in flower petal, moreover, the expression level of RhMYB1 in WR increased along with scent emission, and decreased when the scent emission decreased. It is suggested that RhMYB1 might be a putative identification of gene involved in the biosynthesis of rose scent.

[Development of new SSR markers from EST of SSH cDNA libraries on rose fragrance].

The new SSR markers of rose related fragrance were developed based on the SSH cDNA libraries of rose floral scent mutant. In this study, 10 EST-SSRs (2.6%) from 391 ESTs in the libraries were identified. Six EST-SSRs primers were designed to sequence flanking SSRs. The primer pairs designed were screened on the wild-type Jinyindao, which has flowers full of pleasant scent, and the mutant-type Wangriqinghuai without perceivable floral scent. Five primer pairs were amplified effectively in Jinyindao and Wangriqinghuai, and 3 were polymorphic between Jinyindao and Wangriqinghuai. Eighteen rose cultivars including fragrant roses and nonfragrant roses were identified by the five prime pairs. These results proved that EST-SSR markers are effective markers to identify the polymorphism of the rose.

The development of a high-density genetic map significantly improves the quality of reference genome assemblies for rose.

Roses are important woody plants featuring a set of important traits that cannot be investigated in traditional model plants. Here, we used the restriction-site associated DNA sequencing (RAD-seq) technology to develop a high-density linkage map of the backcross progeny (BC1F1) between Rosa chinensis 'Old Blush' (OB) and R. wichuraiana 'Basyes' Thornless' (BT). We obtained 643.63 million pair-end reads and identified 139,834 polymorphic tags that were distributed uniformly in the rose genome. 2,213 reliable markers were assigned to seven linkage groups (LGs). The length of the genetic map was 1,027.425 cM in total with a mean distance of 0.96 cM per marker locus. This new linkage map allowed anchoring an extra of 1.21/23.14 Mb (12.18/44.52%) of the unassembled OB scaffolds to the seven reference pseudo-chromosomes, thus significantly improved the quality of assembly of OB reference genome. We demonstrate that, while this new linkage map shares high collinearity level with strawberry genome, it also features two chromosomal rearrangements, indicating its usefulness as a resource for understanding the evolutionary scenario among Rosaceae genomes. Together with the newly released genome sequences for OB, this linkage map will facilitate the identification of genetic components underpinning key agricultural and biological traits, hence should greatly advance the studies and breeding efforts of rose.

Yalongjiang River Has Had an Important Role in the Dispersal and Divergence of Rosa soulieana in the Hengduan Mountains of China.

The Hengduan Mountains are the core of the Sino-Himalayan Floristic Region. Rosa soulieana Crép. is an important wild rose species that is widely distributed in the Hengduan Mountains. To provide better future utilization of this wild rose, and also to add some possible proof of the effect of geomorphological and ecological characteristics of the Hengduan Mountains on the current spatial distribution and genetic diversity of local species, the genetic diversity and genetic structure of 556 individuals from 37 populations of R. soulieana were studied using fluorescent amplified fragment length polymorphisms (AFLPs). R. soulieana showed a moderately high level of genetic diversity and a high level of genetic differentiation at the species level. The total percentage of polymorphic loci, total heterozygosity (Ht), Shannon index (I), and heterozygosity value within populations (Hs) were 97.8%, 0.253, 0.339, and 0.139, respectively. More than half of the total genetic variation (54.0%) occurred within populations, and the overall gene differentiation coefficient (Gst) was 0.451. The genetic differentiation among populations was positively and significantly correlated with geographic distance. The neighbor-joining cluster and the Bayesian analysis divided all the populations and individuals into 3 groups, and did not support the morphology based intraspecific varieties. The results confirmed that the ancient R. soulieana of the third group survived in northwestern Yunnan and Yalongjiang valley and then moved upnorth along the valley. The spatial distribution of the other two groups was the result of allopatric divergence due to long period of adaptation to the different climatic conditions of its distribution at either side of the Yalongjiang River.

The Rosa chinensis cv. Viridiflora Phyllody Phenotype Is Associated with Misexpression of Flower Organ Identity Genes.

Phyllody is a flower abnormality in which leaf-like structures replace flower organs in all whorls. Here, we investigated the origin and the molecular mechanism of phyllody phenotype in Rosa chinensis cv. Viridiflora, an ancient naturally occurring Chinese mutant cultivar. Reciprocal grafting experiments and microscopy analyses, demonstrated that the phyllody phenotype in Viridiflora is not associated with phytoplasmas infection. Transcriptome comparisons by the mean of RNA-Seq identified 672 up-regulated and 666 down-regulated genes in Viridiflora compared to its closely related genotype R. chinensis cv. Old Blush. A fraction of these genes are putative homologs of genes known to be involved in flower initiation and development. We show that in flower whorl 2 of Viridiflora, a down-regulation of the floral organ identity genes RcPISTILLATA (RcPI), RcAPETALA3 (RcAP3) and RcSEPALLATA3 (RcSEP3), together with an up-regulation of the putative homolog of the gene SUPPRESSOR of OVEREXPRESSION of CONSTANS1 (RcSOC1) are likely at the origin of the loss of petal identity and leaf-like structures formation. In whorl 3 of Viridiflora, ectopic expression of RcAPETALA2 (RcAP2) along with the down regulation of RcPI, RcAP3, and RcSEP3 is associated with loss of stamens identity and leaf-like structures formation. In whorl 4, the ectopic expression of RcAP2 associated with a down-regulation of RcSEP3 and of the C-class gene RcAGAMOUS correlate with loss of pistil identity. The latter also suggested the antagonist effect between the A and C class genes in the rose. Together, these data suggest that modified expression of the ABCE flower organ identity genes is associated with the phyllody phenotype in the rose Viridiflora and that these genes are important for normal flower organs development.

Transcriptome and gene expression analysis during flower blooming in Rosa chinensis 'Pallida'.

Rosa chinensis 'Pallida' (Rosa L.) is one of the most important ancient rose cultivars originating from China. It contributed the 'tea scent' trait to modern roses. However, little information is available on the gene regulatory networks involved in scent biosynthesis and metabolism in Rosa. In this study, the transcriptome of R. chinensis 'Pallida' petals at different developmental stages, from flower buds to senescent flowers, was investigated using Illumina sequencing technology. De novo assembly generated 89,614 clusters with an average length of 428bp. Based on sequence similarity search with known proteins, 62.9% of total clusters were annotated. Out of these annotated transcripts, 25,705 and 37,159 sequences were assigned to gene ontology and clusters of orthologous groups, respectively. The dataset provides information on transcripts putatively associated with known scent metabolic pathways. Digital gene expression (DGE) was obtained using RNA samples from flower bud, open flower and senescent flower stages. Comparative DGE and quantitative real time PCR permitted the identification of five transcripts encoding proteins putatively associated with scent biosynthesis in roses. The study provides a foundation for scent-related gene discovery in roses.

Isolation, Molecular Characterization, and Mapping of Four Rose MLO Orthologs.

Powdery mildew is a major disease of economic importance in cut and pot roses. As an alternative to conventional resistance breeding strategies utilizing single-dominant genes or QTLs, mildew resistance locus o (MLO)-based resistance might offer some advantages. In dicots such as Arabidopsis, pea, and tomato, loss-of-function mutations in MLO genes confer high levels of broad-spectrum resistance. Here, we report the isolation and characterization of four MLO homologs from a large rose EST collection isolated from leaves. These genes are phylogenetically closely related to other dicot MLO genes that are involved in plant powdery mildew interactions. Therefore, they are candidates for MLO genes involved in rose powdery mildew interactions. Two of the four isolated genes contain all of the sequence signatures considered to be diagnostic for MLO genes. We mapped all four genes to three linkage groups and conducted the first analysis of alternative alleles. This information is discussed in regards to a reverse genetics approach aimed at the selection of rose plants that are homozygous for loss-of-function in one or more MLO genes.