Genomic approach to study floral development genes in Rosa sp

Analysis on morphological characteristics and identification of candidate genes during the flowering development of alfalfa

Flower development is a crucial and complex process in the reproductive stage of plants, which involves the interaction of multiple endogenous signals and environmental factors. However, regulatory mechanism of flower development was unknown in alfalfa (Medicago sativa). In this study, the three stages of flower development of 'M. sativa cv. Gannong No. 5' (G5) and its early flowering and multi flowering mutant (MG5) were comparatively analyzed by transcriptomics. The results showed that compared with late bud stage (S1), 14287 and 8351 differentially expressed genes (DEGs) were identified at early flower stage (S2) in G5 and MG5, and 19941 and 19469 DEGs were identified at late flower stage (S3). Compared with S2, 9574 and 10870 DEGs were identified at S3 in G5 and MG5, respectively. Venn analysis revealed that 547 DEGs were identified among the three comparison groups. KEGG pathway enrichment analysis showed that these genes were involved in the development of alfalfa flowers through redox pathways and plant hormone signaling pathways. Key candidate genes including SnRK2, BSK, GID1, DELLA and CRE1, for regulating the development from buds to mature flowers in alfalfa were screened. In addition, differential expression of transcription factors such as MYB, AP2, bHLH, C2C2, MADS-box, NAC, bZIP, B3 and AUX/IAA also played an important role in this process. The results laid a theoretical foundation for studying the molecular mechanisms of the development process from buds to mature flowers in alfalfa.

Transcriptional Changes in Damask Rose Suspension Cell Culture Revealed by RNA Sequencing

Damask roses (Rosa x damascena) are widely used in cosmetics and pharmaceutics. Here, we established an in vitro suspension cell culture for calli derived from damask rose petals. We analyzed rose suspension cell transcriptomes obtained at two different time points by RNA sequencing to reveal transcriptional changes during rose suspension cell culture. Of the 580 coding RNAs (1.3%) highly expressed in the suspension rose cells, 68 encoded cell wall-associated proteins. However, most RNAs encoded by the chloroplasts and mitochondria are not expressed. Many highly expressed coding RNAs are involved in translation, catalyzing peptide synthesis in ribosomes. Moreover, the amide metabolic process producing naturally occurring alkaloids was the most abundant metabolic process during the propagation of rose suspension cells. During rose cell propagation, coding RNAs involved in the stress response were upregulated at an early stage, while coding RNAs associated with detoxification and transmembrane transport were upregulated at the late stage. We used transcriptome analyses to reveal important biological processes and molecular mechanisms during rose suspension cell culture. Most non-coding (nc) RNAs were not expressed in the rose suspension cells, but a few ncRNAs with unknown functions were highly expressed. The expression of ncRNAs and their target coding RNAs was highly correlated. Taken together, we revealed significant biological processes and molecular mechanisms occurring during rose suspension cell culture using transcriptome analyses.

Genetic and Biochemical Aspects of Floral Scents in Roses

Floral scents possess high ornamental and economic values to rose production in the floricultural industry. In the past two decades, molecular bases of floral scent production have been studied in the rose as well as their genetic inheritance. Some significant achievements have been acquired, such as the comprehensive rose genome and the finding of a novel geraniol synthase in plants. In this review, we summarize the composition of floral scents in modern roses, focusing on the recent advances in the molecular mechanisms of floral scent production and emission, as well as the latest developments in molecular breeding and metabolic engineering of rose scents. It could provide useful information for both studying and improving the floral scent production in the rose.

Comprehensive Genome-Wide Analysis of Histone Acetylation Genes in Roses and Expression Analyses in Response to Heat Stress

Roses have high economic values as garden plants and for cut-flower and cosmetics industries. The growth and development of rose plants is affected by exposure to high temperature. Histone acetylation plays an important role in plant development and responses to various stresses. It is a dynamic and reversible process mediated by histone deacetylases (HDAC) and histone acetyltransferases (HAT). However, information on HDAC and HAT genes of roses is scarce. Here, 23 HDAC genes and 10 HAT genes were identified in the Rosa chinensis ‘Old Blush’ genome. Their gene structures, conserved motifs, physicochemical properties, phylogeny, and synteny were assessed. Analyses of the expression of HDAC and HAT genes using available RNAseq data showed that these genes exhibit different expression patterns in different organs of the three analyzed rose cultivars. After heat stress, while the expression of most HDAC genes tend to be down-regulated, that of HAT genes was up-regulated when rose plants were grown at high-temperature conditions. These data suggest that rose likely respond to high-temperature exposure via modification in histone acetylation, and, thus, paves the way to more studies in order to elucidate in roses the molecular mechanisms underlying rose plants development and flowering.

Transcriptomic and physiological analysis reveals interplay between salicylic acid and drought stress in citrus tree floral initiation

Salicylic acid (SA) and drought stress promote more flowering in sweet orange. The physiological response and molecular mechanism underlying stress-induced floral initiation were discovered by transcriptome profiling. Numerous flowering-regulated genes were identified, and ectopically expressed CsLIP2A promotes early flowering in Arabidopsis. Floral initiation is a critical developmental mechanism associated with external factors, and citrus flowering is mainly regulated by drought stress. However, little is known about the intricate regulatory network involved in stress-induced flowering in citrus. To understand the molecular mechanism of floral initiation in citrus, flower induction was performed on potted Citrus sinensis trees under the combined treatment of salicylic acid (SA) and drought (DR). Physiological analysis revealed that SA treatment significantly normalized the drastic effect of drought stress by increasing antioxidant enzyme activities (SOD, POD, and CAT), relative leaf water content, total chlorophyll, and proline contents and promoting more flowering than drought treatment. Analysis of transcriptome changes in leaves from different treatments showed that 1135, 2728 and 957 differentially expressed genes (DEGs) were revealed in response to DR, SD (SA + DR), and SA (SA + well water) treatments in comparison with the well watered plants, respectively. A total of 2415, 2318 and 1933 DEGs were expressed in DR, SD, and SA in comparison with water recovery, respectively. Some key flowering genes were more highly expressed in SA-treated drought plants than in DR-treated plants. GO enrichment revealed that SA treatment enhances the regulation and growth of meristem activity under drought conditions, but no such a pathway was found to be highly enriched in the control. Furthermore, we focused on various hormones, sugars, starch metabolism, and biosynthesis-related genes. The KEGG analysis demonstrated that DEGs enriched in starch sucrose metabolism and hormonal signal transduction pathways probably account for stress-induced floral initiation in citrus. In addition, a citrus LIPOYLTRANSFERSAE 2A homologous (LIP2A) gene was upregulated by SD treatment. Ectopic expression of CsLIP2A exhibited early flowering in transgenic Arabidopsis. Taken together, this study provides new insight that contributes to citrus tree floral initiation under the SA-drought scenario as well as an excellent reference for stress-induced floral initiation in woody trees.

Detection of Reproducible Major Effect QTL for Petal Traits in Garden Roses

The detection of QTL by association genetics depends on the genetic architecture of the trait under study, the size and structure of the investigated population and the availability of phenotypic and marker data of sufficient quality and quantity. In roses, we previously demonstrated that major QTL could already be detected in small association panels. In this study, we analyzed petal number, petal size and fragrance in a small panel of 95 mostly tetraploid garden rose genotypes. After genotyping the panel with the 68 K Axiom WagRhSNP chip we detected major QTL for all three traits. Each trait was significantly influenced by several genomic regions. Some of the QTL span genomic regions that comprise several candidate genes. Selected markers from some of these regions were converted into KASP markers and were validated in independent populations of up to 282 garden rose genotypes. These markers demonstrate the robustness of the detected effects independent of the set of genotypes analyzed. Furthermore, the markers can serve as tools for marker-assisted breeding in garden roses. Over an extended timeframe, they may be used as a starting point for the isolation of the genes underlying the QTL.

Characterization of floral morphoanatomy and identification of marker genes preferentially expressed during specific stages of cotton flower development

Characterization of anther and ovule developmental programs and expression analyses of stage-specific floral marker genes in Gossypium hirsutum allowed to build a comprehensive portrait of cotton flower development before fiber initiation. Gossypium hirsutum is the most important cotton species that is cultivated worldwide. Although cotton reproductive development is important for fiber production, since fiber is formed on the epidermis of mature ovules, cotton floral development remains poorly understood. Therefore, this work aims to characterize the cotton floral morphoanatomy by performing a detailed description of anther and ovule developmental programs and identifying stage-specific floral marker genes in G. hirsutum. Using light microscopy and scanning electron microscopy, we analyzed anther and ovule development during 11 stages of flower development. To better characterize the ovule development in cotton, we performed histochemical analyses to evaluate the accumulation of phenolic compounds, pectin, and sugar in ovule tissues. After identification of major hallmarks of floral development, three key stages were established in G. hirsutum floral development: in stage 1 (early-EF), sepal, petal, and stamen primordia were observed; in stage 2 (intermediate-IF), primordial ovules and anthers are present, and the differentiating archesporial cells were observed, marking the beginning of microsporogenesis; and in stage 6 (late-LF), flower buds presented initial anther tapetum degeneration and microspore were released from the tetrad, and nucellus and both inner and outer integuments are developing. We used transcriptome data of cotton EF, IF and LF stages to identify floral marker genes and evaluated their expression by real-time quantitative PCR (qPCR). Twelve marker genes were preferentially expressed in a stage-specific manner, including the putative homologs for AtLEAFY, AtAPETALA 3, AtAGAMOUS-LIKE 19 and AtMALE STERILITY 1, which are crucial for several aspects of reproductive development, such as flower organogenesis and anther and petal development. We also evaluated the expression profile of B-class MADS-box genes in G. hirsutum floral transcriptome (EF, IF, and LF). In addition, we performed a comparative analysis of developmental programs between Arabidopsis thaliana and G. hirsutum that considered major morphoanatomical and molecular processes of flower, anther, and ovule development. Our findings provide the first detailed analysis of cotton flower development.

Genome-wide expression of low temperature response genes in Rosa hybrida L

Plants respond to low temperature stress during cold acclimation, a complex process involving changes in physiological and biochemical modifications. The rose serves as a good model to investigate low temperature responses in perennial ornamentals. In this study, a heterologous apple microarray is used to investigate genome-wide expression profiles in Rosa hybrida subjected to low temperature dark treatment. Transcriptome profiles are determined in floral buds at 0h, 2h, and 12h of low temperature treatment (4 °C). It is observed that a total of 134 transcripts are up-regulated and 169 transcripts are down-regulated in response to low temperature. Interestingly, a total of eight up-regulated genes, including those coding for two cytochrome P450 proteins, two ankyrin repeat family proteins, two metal ion binding proteins, and two zinc finger protein-related transcription factors, along with a single down-regulated gene, coding for a dynamin-like protein, are detected. Transcript profiles of 12 genes known to be involved in cold stress response are also validated using qRT-PCR. Furthermore, expression patterns of the AP2/ERF gene family of transcription factors are investigated in both floral buds and leaves. Overall, AP2/ERFs genes are more rapidly induced in leaves than in floral buds. Moreover, differential expression of several AP2/ERF genes are detected earlier in vegetative rather than in reproductive tissues. These findings highlight important roles of various low temperature response genes in mediating cold acclimation, thereby allowing roses to adapt to low temperatures, but without adversely affecting flower bud development and subsequent flowering, while vegetative tissues undergo early adaptation to low temperatures.

Biosynthesis of 2-Phenylethanol in Rose Petals Is Linked to the Expression of One Allele of RhPAAS

Floral scent is one of the most important characters in horticultural plants. Roses (Rosa spp.) have been cultivated for their scent since antiquity. However, probably by selecting for cultivars with long vase life, breeders have lost the fragrant character in many modern roses, especially the ones bred for the cut flower market. The genetic inheritance of scent characters has remained elusive so far. In-depth knowledge of this quantitative trait is thus very much needed to breed more fragrant commercial cultivars. Furthermore, rose hybrids harbor a composite genomic structure, which complexifies quantitative trait studies. To understand rose scent inheritance, we characterized a segregating population from two diploid cultivars, Rosa × hybrida cv H190 and Rosa wichurana, which have contrasting scent profiles. Several quantitative trait loci for the major volatile compounds in this progeny were identified. One among these loci contributing to the production of 2-phenylethanol, responsible for the characteristic odor of rose, was found to be colocalized with a candidate gene belonging to the 2-phenylethanol biosynthesis pathway: the PHENYLACETALDEHYDE SYNTHASE gene RhPAAS An in-depth allele-specific expression analysis in the progeny demonstrated that only one allele was highly expressed and was responsible for the production of 2-phenylethanol. Unexpectedly, its expression was found to start early during flower development, before the production of the volatile 2-phenylethanol, leading to the accumulation of glycosylated compounds in petals.

Comparative transcriptomics identifies patterns of selection in roses

BACKGROUND

Roses are important plants for human beings with pivotal economical and biological traits like continuous flowering, flower architecture, color and scent. Due to frequent hybridization and high genome heterozygosity, classification of roses and their relatives remains a big challenge.

RESULTS

Here, to identify potential markers for phylogenetic reconstruction and to reveal the patterns of natural selection in roses, we generated sets of high quality and comprehensive reference transcriptomes for Rosa chinensis 'Old Blush' (OB) and R. wichuriana 'Basye's Thornless' (BT), two species exhibiting contrasted traits of high economical importance. The assembled reference transcriptomes showed transcripts N50 above 2000 bp. Two roses shared about 10,073 transcripts (N50 = 2282 bp), in which a set of 5959 transcripts was conserved within genera of Rosa. Further comparison with species in Rosaceae identified 4447 transcripts being common (Rosaceae-common) in Rosa, Malus, Prunus, Rubus, and Fragaria, while a pool of 164 transcripts being specific for roses (Rosa-specific). Among the Rosaceae-common transcripts, 409 transcripts showed a signature of positive selection and a clustered expression in different tissues. Interestingly, nine of these rapidly evolving genes were related to DNA damage repair and responses to environmental stimulus, a potential associated with genome confliction post hybridization. Coincident with this fast evolution pattern in rose genes, 24 F-box and four TMV resistant proteins were significantly enriched in the Rosa-specific genes.

CONCLUSIONS

We expect that these Rosaceae-common and Rosa-specific transcripts should facilitate the phylogenetic analysis of Rosaceae plants as well as investigations of Rosa-specific biology. The data reported here could provide fundamental genomic tools and knowledge critical for understanding the biology and domestication of roses and for roses breeding.

A miR172 target-deficient AP2-like gene correlates with the double flower phenotype in roses

One of the well-known floral abnormalities in flowering plants is the double-flower phenotype, which corresponds to flowers that develop extra petals, sometimes even containing entire flowers within flowers. Because of their highly priced ornamental value, spontaneous double-flower variants have been found and selected for in a wide range of ornamental species. Previously, double flower formation in roses was associated with a restriction of AGAMOUS expression domain toward the centre of the meristem, leading to extra petals. Here, we characterized the genomic region containing the mutation associated with the switch from simple to double flowers in the rose. An APETALA2-like gene (RcAP2L), a member of the Target Of EAT-type (TOE-type) subfamily, lies within this interval. In the double flower rose, two alleles of RcAP2L are present, one of which harbours a transposable element inserted into intron 8. This insertion leads to the creation of a miR172 resistant RcAP2L variant. Analyses of the presence of this variant in a set of simple and double flower roses demonstrate a correlation between the presence of this allele and the double flower phenotype. These data suggest a role of this miR172 resistant RcAP2L variant in regulating RcAGAMOUS expression and double flower formation in Rosa sp.

Developmental transcriptome analysis of floral transition in Rosa odorata var. gigantea

Expression analyses revealed that floral transition of Rosa odorata var. gigantea is mainly regulated by VRN1, COLs, DELLA and KSN, with contributions by the effects of phytohormone and starch metabolism. Seasonal plants utilize changing environmental and developmental cues to control the transition from vegetative growth to flowering at the correct time of year. This study investigated global gene expression profiles at different developmental stages of Rosa odorata var. gigantea by RNA-sequencing, combined with phenotypic characterization and physiological changes. Gene ontology enrichment analysis of the differentially expressed genes (DEGs) between four different developmental stages (vegetative meristem, pre-floral meristem, floral meristem and secondary axillary buds) indicated that DNA methylation and the light reaction played a large role in inducing the rose floral transition. The expression of SUF and FLC, which are known to play a role in delaying flowering until vernalization, was down-regulated from the vegetative to the pre-floral meristem stage. In contrast, the expression of VRN1, which promotes flowering by repressing FLC expression, increased. The expression of DELLA proteins, which function as central nodes in hormone signaling pathways, and probably involve interactions between GA, auxin, and ABA to promote the floral transition, was well correlated with the expression of floral integrators, such as AGL24, COL4. We also identified DEGs associated with starch metabolism correlated with SOC1, AGL15, SPL3, AGL24, respectively. Taken together, our results suggest that vernalization and photoperiod are prominent cues to induce the rose floral transition, and that DELLA proteins also act as key regulators. The results summarized in the study on the floral transition of the seasonal rose lay a foundation for further functional demonstration, and have profound economic and ornamental values.

The Rosa genome provides new insights into the domestication of modern roses

Roses hold high cultural and economic importance as ornamentals and for the perfume industry. We report the rose whole genome sequencing and assembly and resequencing of major genotypes that contributed to rose domestication. We generated a homozygous genotype from a heterozygous diploid modern roses progenitor, Rosa chinensis ‘Old Blush’. Using Single Molecule Real-Time sequencing and a meta-assembly approach we obtained one of the most complete plant genomes to date. Diversity analyses highlighted the mosaic origin of ‘La France’, one of the first hybrids combining the growth vigor of European species and recurrent blooming of Chinese species. Genomic segments of Chinese ancestry revealed new candidate genes for recurrent blooming. Reconstructing regulatory and secondary metabolism pathways allowed us to propose a model of interconnected regulation of scent and flower color. This genome provides a foundation for understanding the mechanisms governing rose traits and will accelerate improvement in roses, Rosaceae and ornamentals.

Petal senescence: a hormone view

Flowers are highly complex organs that have evolved to enhance the reproductive success of angiosperms. As a key component of flowers, petals play a vital role in attracting pollinators and ensuring successful pollination. Having fulfilled this function, petals senesce through a process that involves many physiological and biochemical changes that also occur during leaf senescence. However, petal senescence is distinct, due to the abundance of secondary metabolites in petals and the fact that petal senescence is irreversible. Various phytohormones are involved in regulating petal senescence, and are thought to act both synergistically and antagonistically. In this regard, there appears to be developmental point during which such regulatory signals are sensed and senescence is initiated. Here, we review current understanding of petal senescence, and discuss associated regulatory mechanisms involving hormone interactions and epigenetic regulation.

MIKCC-type MADS-box genes in Rosa chinensis: the remarkable expansion of ABCDE model genes and their roles in floral organogenesis

MIKC^C-type MADS-box (MIKC^C) genes encode transcription factors that have crucial roles in controlling floral organogenesis and flowering time in plants. Although this gene family has been well characterized in many plant species, its evolutionary and comprehensive functional analysis in rose is lacking. In this study, 58 non-redundant MIKC^C uni-transcripts were extensively identified from rose transcriptomes. Phylogenetic analysis placed these genes into 12 clades with their Arabidopsis and strawberry counterparts, and revealed that ABCDE model (including AP1/FUL, AP3/PI, AG, and SEP clades), and SOC1 and AGL6 clade genes have remarkably expanded in Rosa chinensis, whereas genes from the FLC and AGL17 clades were undetectable. Sequence alignments suggest that the AP3/PI clade may contribute to more specific functions in rose due to a high variation of amino acid residues within its MADS-box domains. A comparative analysis of gene expression in specific floral organ differentiation stages and floral organs between R. chinensis cv. Old Blush and the closely related mutant genotype R. chinensis cv. Viridiflora (floral organs mutated into leaf-like structures) further revealed the roles of ABCDE model genes during floral organogenesis in rose. Analysis of co-expression networks provided an overview of the regulatory mechanisms of rose MIKC^C genes and shed light on both the prominent roles of AP3/PI clade genes in floral organogenesis and the roles of RcAGL19, RcAGL24, and RcSOC1 in regulating floral transition in rose. Our analyses provide an overall insight of MIKC^C genes in rose and their potential roles in floral organogenesis.

Comparative transcriptome analysis of the floral transition in Rosa chinensis 'Old Blush' and R. odorata var. gigantea

The floral transition is a crucial developmental event, but little is known about the underlying regulatory networks in seasonally and continuously flowering roses. In this study, we compared the genetic basis of flowering in two rose species, Rosa chinensis ‘Old Blush’, which flowers continuously, and R. odorata var. gigantea, which blooms in early spring. Gene ontology (GO) terms related to methylation, light reaction, and starch metabolism were enriched in R. odorata var. gigantea and terms associated with sugar metabolism were enriched in R. chinensis ‘Old Blush’ during the floral transition. A MapMan analysis revealed that genes involved in hormone signaling mediate the floral transition in both taxa. Furthermore, differentially expressed genes (DEGs) involved in vernalization, photoperiod, gibberellin (GA), and starch metabolism pathways converged on integrators, e.g., LFY, AGL24, SOC1, CAL, and COLs, to regulate the floral transition in R. odorata var. gigantea, while DEGs related to photoperiod, sugar metabolism, and GA pathways, including COL16, LFY, AGL11, 6PGDH, GASA4, and BAM, modulated the floral transition in R. chinensis ‘Old Blush.’ Our analysis of the genes underlying the floral transition in roses with different patterns of flowering provides a basis for further functional studies.

Transcriptome profiling of litchi leaves in response to low temperature reveals candidate regulatory genes and key metabolic events during floral induction

Background

Litchi (Litchi chinensis Sonn.) is an economically important evergreen fruit tree widely cultivated in subtropical areas. Low temperature is absolutely required for floral induction of litchi, but its molecular mechanism is not fully understood. Leaves of litchi played a key role during floral induction and could be the site of low temperature perception. Therefore, leaves were treated under different temperature (15 °C/25 °C), and high-throughput RNA sequencing (RNA-Seq) performed with leaf samples for the de novo assembly and digital gene expression (DGE) profiling analyses to investigate low temperature-induced gene expression changes.

Results

83,107 RNA-Seq unigenes were de novo assembled with a mean length of 1221 bp and approximately 61% of these unigenes (50,345) were annotated against public protein databases. Differentially-expressed genes (DEGs) under low temperature treatment in comparison with the control group were the main focus of our study. Hierarchical clustering analysis arranged 2755 DEGs into eight groups with three significant expression clusters (p-value ≤ 0.05) during floral induction. With the increasing contents of sugars and starch, the expression of genes involved in metabolism of sugars increased dramatically after low temperature induction. One FT gene (Unigene0025396, LcFT1) which produces a protein called ‘florigen’ was also detected among DEGs of litchi. LcFT1 exhibited an apparent specific tissue and its expression was highly increased after low temperature induction, GUS staining results also showed GUS activity driven by LcFT1 gene promoter can be induced by low temperature, which indicated LcFT1 probably played a pivotal role in the floral induction of litchi under low temperature.

Conclusions

Our study provides a global survey of transcriptomes to better understand the molecular mechanisms underlying changes of leaves in response to low temperature induction in litchi. The analyses of transcriptome profiles and physiological indicators will help us study the complicated metabolism of floral induction in the subtropic evergreen plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3747-x) contains supplementary material, which is available to authorized users.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

BACKGROUND

Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae).

RESULTS

The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species.

CONCLUSIONS

The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Transcriptome of the floral transition in Rosa chinensis 'Old Blush'

BACKGROUND

The floral transition plays a vital role in the life of ornamental plants. Despite progress in model plants, the molecular mechanisms of flowering regulation remain unknown in perennial plants. Rosa chinensis 'Old Blush' is a unique plant that can flower continuously year-round. In this study, gene expression profiles associated with the flowering transition were comprehensively analyzed during floral transition in the rose.

RESULTS

According to the transcriptomic profiles, 85,663 unigenes and 1,637 differentially expressed genes (DEGs) were identified, among which 32 unigenes were involved in the circadian clock, sugar metabolism, hormone, and autonomous pathways. A hypothetical model for the regulation of floral transition was proposed in which the candidate genes function synergistically the floral transition process. Hormone contents and biosynthesis and metabolism genes fluctuated during the rose floral transition process. Gibberellins (GAs) inhibited rose floral transition, the content of GAs gradually decreased and GA2ox and SCL13 were upregulated from vegetative (VM) meristem to floral meristem (FM). Auxin plays an affirmative part in mediating floral transition, auxin content and auxin-related gene expression levels were gradually upregulated during the floral transition of the rose. However, ABA content and ABA signal genes were gradually downregulated, suggesting that ABA passively regulates the rose floral transition by participating in sugar signaling. Furthermore, sugar content and sugar metabolism genes increased during floral transition in the rose, which may be a further florigenic signal that activates floral transition. Additionally, FRI, FY, DRM1, ELIP, COP1, CO, and COL16 are involved in the circadian clock and autonomous pathway, respectively, and they play a positively activating role in regulating floral transition. Overall, physiological changes associated with genes involved in the circadian clock or autonomous pathway collectively regulated the rose floral transition.

CONCLUSIONS

Our results summarize a valuable collective of gene expression profiles characterizing the rose floral transition. The DEGs are candidates for functional analyses of genes affecting the floral transition in the rose, which is a precious resource that reveals the molecular mechanism of mediating floral transition in other perennial plants.

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.

PLANT VOLATILES. Biosynthesis of monoterpene scent compounds in roses

The scent of roses (Rosa x hybrida) is composed of hundreds of volatile molecules. Monoterpenes represent up to 70% percent of the scent content in some cultivars, such as the Papa Meilland rose. Monoterpene biosynthesis in plants relies on plastid-localized terpene synthases. Combining transcriptomic and genetic approaches, we show that the Nudix hydrolase RhNUDX1, localized in the cytoplasm, is part of a pathway for the biosynthesis of free monoterpene alcohols that contribute to fragrance in roses. The RhNUDX1 protein shows geranyl diphosphate diphosphohydrolase activity in vitro and supports geraniol biosynthesis in planta.

Transcriptomic analysis of floral initiation in litchi (Litchi chinensis Sonn.) based on de novo RNA sequencing

Comparative transcriptome analysis of litchi ( Litchi chinensis Sonn.) buds at two developmental stages revealed multiple processes involving various phytohormones regulating floral initiation, and expression of numerous flowering-related genes. Floral initiation is a critical and complicated plant developmental process involving interactions of numerous endogenous and environmental factors, but little is known about the complex network regulating floral initiation in litchi (Litchi chinensis Sonn.). Illumina second-generation sequencing is an efficient method for obtaining massive transcriptional information resulting from phase changes in plant development. In this study, comparative transcriptomic analysis was performed with resting and emerging panicle stage buds, to gain further understanding of the molecular mechanisms involved in floral initiation in litchi. Abundance analysis identified 5,928 unigenes exhibiting at least twofold differences in expression between the two bud stages. Of these, 4,622 unigenes were up-regulated and 1,306 were down-regulated in panicle-emerging buds compared with resting buds. KEGG pathway enrichment analysis revealed that unigenes exhibiting differential expression were involved in the metabolism and signal transduction of various phytohormones. The expression levels of unigenes annotated as auxin, cytokinin, jasmonic acid, and salicylic acid biosynthesis were up-regulated, whereas those unigenes annotated as abscisic acid biosynthesis were down-regulated during floral initiation. In addition, 188 unigenes exhibiting sequence similarities to known flowering-related genes from other plants were differentially expressed during floral initiation. Thirteen genes were selected for confirmation of expression levels using quantitative-PCR. Our results provide abundant sequence resources for studying mechanisms underlying floral initiation in litchi and establish a platform for further studies of litchi and other evergreen fruit trees.

Anchoring linkage groups of the Rosa genetic map to physical chromosomes with tyramide-FISH and EST-SNP markers

In order to anchor Rosa linkage groups to physical chromosomes, a combination of the Tyramide-FISH technology and the modern molecular marker system based on High Resolution Melting (HRM) is an efficient approach. Although, Tyramide-FISH is a very promising technique for the visualization of short DNA probes, it is very challenging for plant species with small chromosomes such as Rosa. In this study, we successfully applied the Tyramide-FISH technique for Rosa and compared different detection systems. An indirect detection system exploiting biotinylated tyramides was shown to be the most suitable technique for reliable signal detection. Three gene fragments with a size of 1100 pb–1700 bp (Phenylalanine Ammonia Lyase, Pyrroline-5-Carboxylate Synthase and Orcinol O-Methyl Transferase) have been physically mapped on chromosomes 7, 4 and 1, respectively, of Rosa wichurana. The signal frequency was between 25% and 40%. HRM markers of these 3 gene fragments were used to include the gene fragments on the existing genetic linkage map of Rosa wichurana. As a result, three linkage groups could be anchored to their physical chromosomes. The information was used to check for synteny between the Rosa chromosomes and Fragaria.

Transcriptomic analysis of flower development in wintersweet (Chimonanthus praecox)

Wintersweet (Chimonanthus praecox) is familiar as a garden plant and woody ornamental flower. On account of its unique flowering time and strong fragrance, it has a high ornamental and economic value. Despite a long history of human cultivation, our understanding of wintersweet genetics and molecular biology remains scant, reflecting a lack of basic genomic and transcriptomic data. In this study, we assembled three cDNA libraries, from three successive stages in flower development, designated as the flower bud with displayed petal, open flower and senescing flower stages. Using the Illumina RNA-Seq method, we obtained 21,412,928, 26,950,404, 24,912,954 qualified Illumina reads, respectively, for the three successive stages. The pooled reads from all three libraries were then assembled into 106,995 transcripts, 51,793 of which were annotated in the NCBI non-redundant protein database. Of these annotated sequences, 32,649 and 21,893 transcripts were assigned to gene ontology categories and clusters of orthologous groups, respectively. We could map 15,587 transcripts onto 312 pathways using the Kyoto Encyclopedia of Genes and Genomes pathway database. Based on these transcriptomic data, we obtained a large number of candidate genes that were differentially expressed at the open flower and senescing flower stages. An analysis of differentially expressed genes involved in plant hormone signal transduction pathways indicated that although flower opening and senescence may be independent of the ethylene signaling pathway in wintersweet, salicylic acid may be involved in the regulation of flower senescence. We also succeeded in isolating key genes of floral scent biosynthesis and proposed a biosynthetic pathway for monoterpenes and sesquiterpenes in wintersweet flowers, based on the annotated sequences. This comprehensive transcriptomic analysis presents fundamental information on the genes and pathways which are involved in flower development in wintersweet. And our data provided a useful database for further research of wintersweet and other Calycanthaceae family plants.

An NAC transcription factor controls ethylene-regulated cell expansion in flower petals

Cell expansion is crucial for plant growth. It is well known that the phytohormone ethylene functions in plant development as a key modulator of cell expansion. However, the role of ethylene in the regulation of this process remains unclear. In this study, 2,189 ethylene-responsive transcripts were identified in rose (Rosa hybrida) petals using transcriptome sequencing and microarray analysis. Among these transcripts, an NAC (for no apical meristem [NAM], Arabidopsis transcription activation factor [ATAF], and cup-shaped cotyledon [CUC])-domain transcription factor gene, RhNAC100, was rapidly and dramatically induced by ethylene in the petals. Interestingly, accumulation of the RhNAC100 transcript was modulated by ethylene via microRNA164-dependent posttranscriptional regulation. Overexpression of RhNAC100 in Arabidopsis (Arabidopsis thaliana) substantially reduced the petal size by repressing petal cell expansion. By contrast, silencing of RhNAC100 in rose petals using virus-induced gene silencing significantly increased petal size and promoted cell expansion in the petal abaxial subepidermis (P < 0.05). Expression analysis showed that 22 out of the 29 cell expansion-related genes tested exhibited changes in expression in RhNAC100-silenced rose petals. Moreover, of those genes, one cellulose synthase and two aquaporin genes (Rosa hybrida Cellulose Synthase2 and R. hybrida Plasma Membrane Intrinsic Protein1;1/2;1) were identified as targets of RhNAC100. Our results suggest that ethylene regulates cell expansion by fine-tuning the microRNA164/RhNAC100 module and also provide new insights into the function of NAC transcription factors.

Validation of reference transcripts in strawberry (Fragaria spp.)

Contemporary methods to assay gene expression depend on a stable set of reference transcripts for accurate quantitation. A lack of well-tested reference genes slows progress in characterizing gene expression in high-value specialty crops. In this study, a set of strawberry (Fragaria spp.) constitutively expressed reference genes has been identified by merging digital gene expression data with expression profiling. Constitutive reference candidates were validated using quantitative PCR and hybridization. Several transcripts have been identified that show improved stability across tissues relative to traditional reference transcripts. Results are similar between commercial octoploid strawberry and the diploid model. Our findings also show that while some never-before-used references are appropriate for most applications, even the most stable reference transcripts require careful assessment across the diverse tissues and fruit developmental states before being adopted as controls.

RcRR1, a Rosa canina type-A response regulator gene, is involved in cytokinin-modulated rhizoid organogenesis

In vitro, a new protocol of plant regeneration in rose was achieved via protocorm-like bodies (PLBs) induced from the root-like organs named rhizoids that developed from leaf explants. The development of rhizoids is a critical stage for efficient regeneration, which is triggered by exogenous auxin. However, the role of cytokinin in the control of organogenesis in rose is as yet uncharacterized. The aim of this study was to elucidate the molecular mechanism of cytokinin-modulated rhizoid formation in Rosa canina. Here, we found that cytokinin is a key regulator in the formation of rhizoids. Treatment with cytokinin reduced callus activity and significantly inhibited rhizoid formation in Rosa canina. We further isolated the full-length cDNA of a type-A response regulator gene of cytokinin signaling, RcRR1, from which the deduced amino acid sequence contained the conserved DDK motif. Gene expression analysis revealed that RcRR1 was differentially expressed during rhizoid formation and its expression level was rapidly up-regulated by cytokinin. In addition, the functionality of RcRR1 was tested in Arabidopsis. RcRR1 was found to be localized to the nucleus in GFP-RcRR1 transgenic plants and overexpression of RcRR1 resulted in increased primary root length and lateral root density. More importantly, RcRR1 overexpression transgenic plants also showed reduced sensitivity to cytokinin during root growth; auxin distribution and the expression of auxin efflux carriers PIN genes were altered in RcRR1 overexpression plants. Taken together, these results demonstrate that RcRR1 is a functional type-A response regulator which is involved in cytokinin-regulated rhizoid formation in Rosa canina.

Integrative analysis of miRNA and mRNA profiles in response to ethylene in rose petals during flower opening

MicroRNAs play an important role in plant development and plant responses to various biotic and abiotic stimuli. As one of the most important ornamental crops, rose (Rosa hybrida) possesses several specific morphological and physiological features, including recurrent flowering, highly divergent flower shapes, colors and volatiles. Ethylene plays an important role in regulating petal cell expansion during rose flower opening. Here, we report the population and expression profiles of miRNAs in rose petals during flower opening and in response to ethylene based on high throughput sequencing. We identified a total of 33 conserved miRNAs, as well as 47 putative novel miRNAs were identified from rose petals. The conserved and novel targets to those miRNAs were predicted using the rose floral transcriptome database. Expression profiling revealed that expression of 28 known (84.8% of known miRNAs) and 39 novel (83.0% of novel miRNAs) miRNAs was substantially changed in rose petals during the earlier opening period. We also found that 28 known and 22 novel miRNAs showed expression changes in response to ethylene treatment. Furthermore, we performed integrative analysis of expression profiles of miRNAs and their targets. We found that ethylene-caused expression changes of five miRNAs (miR156, miR164, miR166, miR5139 and rhy-miRC1) were inversely correlated to those of their seven target genes. These results indicate that these miRNA/target modules might be regulated by ethylene and were involved in ethylene-regulated petal growth.

Genetics and genomics of flower initiation and development in roses

Roses hold high symbolic value and great cultural importance in different societies throughout human history. They are widely used as garden ornamental plants, as cut flowers, and for the production of essential oils for the perfume and cosmetic industries. Domestication of roses has a long and complex history, and the rose species have been hybridized across vast geographic areas such as Europe, Asia, and the Middle East. The domestication processes selected several flower characters affecting floral quality, such as recurrent flowering, double flowers, petal colours, and fragrance. The molecular and genetic events that determine some of these flower characters cannot be studied using model species such as Arabidopsis thaliana, or at least only in a limited manner. In this review, we comment on the recent development of genetic, genomic, and transcriptomic tools for roses, and then focus on recent advances that have helped unravel the molecular mechanisms underlying several rose floral traits.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars

Background

Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants—making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars.

Results

We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: ‘Vital’, ‘Maroussia’, and ‘Sympathy’ and Rosa rugosa Thunb. , respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes.

Conclusions

In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

Transcriptome database resource and gene expression atlas for the rose

BACKGROUND

For centuries roses have been selected based on a number of traits. Little information exists on the genetic and molecular basis that contributes to these traits, mainly because information on expressed genes for this economically important ornamental plant is scarce.

RESULTS

Here, we used a combination of Illumina and 454 sequencing technologies to generate information on Rosa sp. transcripts using RNA from various tissues and in response to biotic and abiotic stresses. A total of 80714 transcript clusters were identified and 76611 peptides have been predicted among which 20997 have been clustered into 13900 protein families. BLASTp hits in closely related Rosaceae species revealed that about half of the predicted peptides in the strawberry and peach genomes have orthologs in Rosa dataset. Digital expression was obtained using RNA samples from organs at different development stages and under different stress conditions. qPCR validated the digital expression data for a selection of 23 genes with high or low expression levels. Comparative gene expression analyses between the different tissues and organs allowed the identification of clusters that are highly enriched in given tissues or under particular conditions, demonstrating the usefulness of the digital gene expression analysis. A web interface ROSAseq was created that allows data interrogation by BLAST, subsequent analysis of DNA clusters and access to thorough transcript annotation including best BLAST matches on Fragaria vesca, Prunus persica and Arabidopsis. The rose peptides dataset was used to create the ROSAcyc resource pathway database that allows access to the putative genes and enzymatic pathways.

CONCLUSIONS

The study provides useful information on Rosa expressed genes, with thorough annotation and an overview of expression patterns for transcripts with good accuracy.