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

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.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

MicroRNAs Roles in Plants Secondary Metabolism

Plant growth and development is dependent on the regulation of classes of microRNAs (miRNAs) that have emerged as important gene regulators. These miRNAs can regulate plant gene expression to function. They play an important roles in biological homeostasis and environmental response controls. A wide range of plant biological and metabolic processes, including developmental timing, tissues specific development, and differentiation, depends on miRNAs. They perpetually regulate secondary metabolite functions in different plant family lines. Mapping of molecular phylogenies shows the distribution of secondary metabolism in the plant territory. More importantly, a lot of information related to miRNA regulatory processes in plants is revealed, but the role of miRNAs in secondary metabolism regulation and functions of the metabolites are still unclear. In this review, we pinnacle some potential miRNAs regulating the secondary metabolite biosynthesis activities in plants. This will provide an alternative knowledge for functional studies of secondary metabolism.

Novel approaches on identification of conserved miRNAs for broad-spectrum Potyvirus control measures

Potyviridae comprises more than 200 ssRNA viruses, many of which have a broad host range and geographical distributions. Potyvirids (members of Potyviridae) infect several economically important plants such as saffron, cardamom, cucumber, pepper, potato, tomato, yam, etc. Cumulatively, potyvirids cause a substantial economic loss. The major bottleneck in developing an efficient antiviral strategy is that viruses quickly evade host immunity owing to their higher mutation and recombination rates. Due to this reason, the emergence of newer and improved broad-spectrum approaches to combat viral infections is essential. The use of microRNA's (miRNA) to circumvent viral infection against animal viruses has been successfully employed. Fewer studies reported the development of efficient miRNA-based antivirus resistant strategies against plant viruses and none focused on multiple virus resistance. We focused on potyviruses since studies are limited and identification of conserved miRNAs among various host plants would be an initiative to design broad-spectrum antivirus strategies. In this study, we predicted evolutionarily conserved miRNAs by BLAST searching of reported miRNAs from 15 plants against the GSS and EST sequences of banana. A total of nine miRNAs were predicted and screened in nine diverse potyvirids' hosts (Banana, Tomato, Green gram, Jasmine, Chilli, Coriander, Onion, Rose and Colocasia) belonging to eight different orders (Zingiberales, Solanales, Fabales, Lamiales, Apiales, Asperagales, Rosales and Alismatales). Results suggested that miR168 and miR162 are conserved among all the selected plants. This comprehensive study laid the foundations to design broad-spectrum antivirus resistance using miRNAs. To conclude miR168 and miR162 are conserved among many plants and play a crucial role in evading virus infection which could be used as a potential candidate for developing antiviral strategies against potyvirid infections.

Genes and genome editing tools for breeding desirable phenotypes in ornamentals

We review the main genes underlying commercial traits in cut flower species and critically discuss the possibility to apply genome editing approaches to produce novel variation and phenotypes. Promoting flowering and flower longevity as well as creating novelty in flower structure, colour range and fragrances are major objectives of ornamental plant breeding. The novel genome editing techniques add new possibilities to study gene function and breed new varieties. The implementation of such techniques, however, relies on detailed information about structure and function of genomes and genes. Moreover, improved protocols for efficient delivery of editing reagents are required. Recent results of the application of genome editing techniques to elite ornamental crops are discussed in this review. Enabling technologies and genomic resources are reviewed in relation to the implementation of such approaches. Availability of the main gene sequences, underlying commercial traits and in vitro transformation protocols are provided for the world's best-selling cut flowers, namely rose, lily, chrysanthemum, lisianthus, tulip, gerbera, freesia, alstroemeria, carnation and hydrangea. Results obtained so far are described and their implications for the improvement of flowering, flower architecture, colour, scent and shelf-life are discussed.

Integrated sRNAome and RNA-Seq analysis reveals miRNA effects on betalain biosynthesis in pitaya

BACKGROUND

MicroRNAs (miRNAs) and their regulatory functions in anthocyanin, carotenoid, and chlorophyll accumulation have been extensively characterized in many plant species. However, the miRNA regulatory mechanism in betalain biosynthesis remains mostly unknown.

RESULTS

In this study, 126 conserved miRNAs and 41 novel miRNAs were first isolated from Hylocereus monacanthus, among which 95 conserved miRNAs belonged to 53 miRNA families. Thirty-four candidate miRNAs related to betalain biosynthesis were differentially expressed. The expression patterns of those differential expressed miRNAs were analyzed in various pitaya tissues by RT-qPCR. A significantly negative correlation was detected between the expression levels of half those miRNAs and corresponding target genes. Target genes of miRNAs i.e. Hmo-miR157b-HmSPL6-like, Hmo-miR160a-Hpcyt P450-like3, Hmo-miR6020-HmCYP71A8-like, Hmo-novel-2-HmCYP83B1-like, Hmo-novel-15-HmTPST-like, Hmo-miR828a-HmTT2-like, Hmo-miR858-HmMYB12-like, Hmo-miR858-HmMYBC1-like and Hmo-miR858-HmMYB2-like were verified by 5'RACE and transient expression system in tobacco.

CONCLUSIONS

Hmo-miR157b, Hmo-miR160a, Hmo-miR6020 Hmo-novel-2, Hmo-novel-15, Hmo-miR828a and Hmo-miR858 play important roles in pitaya fruit coloration and betalain accumulation. Our findings provide new insights into the roles of miRNAs and their target genes of regulatory functions involved in betalain biosynthesis of pitaya.

Hybrid-Transcriptome Sequencing and Associated Metabolite Analysis Reveal Putative Genes Involved in Flower Color Difference in Rose Mutants

Gene mutation is a common phenomenon in nature that often leads to phenotype differences, such as the variations in flower color that frequently occur in roses. With the aim of revealing the genomic information and inner mechanisms, the differences in the levels of both transcription and secondary metabolism between a pair of natural rose mutants were investigated by using hybrid RNA-sequencing and metabolite analysis. Metabolite analysis showed that glycosylated derivatives of pelargonidin, e.g., pelargonidin 3,5 diglucoside and pelargonidin 3-glucoside, which were not detected in white flowers (Rosa 'Whilte Mrago Koster'), constituted the major pigments in pink flowers. Conversely, the flavonol contents of petal, such as kaempferol-3-glucoside, quercetin 3-glucoside, and rutin, were higher in white flowers. Hybrid RNA-sequencing obtained a total of 107,280 full-length transcripts in rose petal which were annotated in major databases. Differentially expressed gene (DEG) analysis showed that the expression of genes involved in the flavonoid biosynthesis pathway was significantly different, e.g., CHS, FLS, DFR, LDOX, which was verified by qRT-PCR during flowering. Additionally, two MYB transcription factors were found and named RmMYBAN2 and RmMYBPA1, and their expression patterns during flowering were also analyzed. These findings indicate that these genes may be involved in the flower color difference in the rose mutants, and competition between anthocyanin and flavonol biosynthesis is a primary cause of flower color variation, with its regulation reflected by transcriptional and secondary metabolite levels.

Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosarugosa Thunb

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles by regulating other genes. Rosa rugosa Thunb. is an important ornamental and edible plant, yet there are only a few studies on the miRNAs and their functions in R. rugosa.

RESULTS

We sequenced 10 samll RNA profiles from the roots, petals, pollens, stamens, and leaves and 4 RNA-seq profiles in leaves and petals to analysis miRNA, phasiRNAs and mRNAs in R. rugosa. In addition, we acquired a degradome sequencing profile from leaf of R. rugosa to identify miRNA and phasiRNA targets using the SeqTar algorithm. We have identified 321 conserved miRNA homologs including primary transcripts for 25 conserved miRNAs, and 22 novel miRNAs. We identified 592 putative targets of the conserved miRNAs or tasiRNAs that showed significant accumulations of degradome reads. We found differential expression patterns of conserved miRNAs in five different tissues of R. rugosa. We identified three hundred and thirty nine 21 nucleotide (nt) PHAS loci, and forty nine 24 nt PHAS loci, respectively. Our results suggest that miR482 triggers generations of phasiRNAs by targeting nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes in R. rugosa. Our results also suggest that the deregulated genes in leaves and petals are significantly enriched in GO terms and KEGG pathways related to metabolic processes and photosynthesis.

CONCLUSIONS

These results significantly enhanced our knowledge of the miRNAs and phasiRNAs, as well as their potential functions, in R. rugosa.

Tissue-specific transcriptional biomarkers in medicinal plants: Application of large-scale meta-analysis and computational systems biology

Biosynthesis of secondary metabolites in plant is a complex process, regulated by many genes and influenced by several factors. In recent years, the next-generation sequencing (NGS) technology and advanced statistical analysis such as meta-analysis and computational systems biology have provided novel opportunities to overcome biological complexity. Here, we performed a meta-analysis on publicly available transcriptome datasets of twelve economically significant medicinal plants to identify differentially expressed genes (DEGs) between shoot and root tissues and to find the key molecular features which may be effective in the biosynthesis of secondary metabolites. Meta-analysis identified a total of 880 genes with differential expression between two tissues. Functional enrichment and KEGG pathway analysis indicated that the functions of those DEGs are highly associated with the developmental process, starch metabolic process, response to stimulus, porphyrin and chlorophyll metabolism, biosynthesis of secondary metabolites and phenylalanine metabolism. In addition, systems biology analysis of the DEGs was applied to find protein-protein interaction network and discovery of significant modules. The detected modules were associated with hormone signal transduction, transcription repressor activity, response to light stimulus and epigenetic processes. Finally, analysis was extended to search for putative miRNAs that are associated with DEGs. A total of 31 miRNAs were detected which belonged to 16 conserved families. The present study provides a comprehensive view to better understand the tissue-specific expression of genes and mechanisms involved in secondary metabolites synthesis and may provide candidate genes for future researches to improve yield of secondary metabolites.

Genome engineering in ornamental plants: Current status and future prospects

Ornamental plants, like roses, carnations, and chrysanthemums, are economically important and are sold all over the world. In addition, numerous cut and garden flowers add colors to homes and gardens. Various strategies of plant breeding have been employed to improve traits of many ornamental plants. These approaches span from conventional techniques, such as crossbreeding and mutation breeding, to genetically modified plants. Recently, genome editing has become available as an efficient means for modifying traits in plant species. Genome editing technology is useful for genetic analysis and is poised to become a common breeding method for ornamental plants. In this review, we summarize the benefits and limitations of conventional breeding techniques and genome editing methods and discuss their future potential to accelerate the rate breeding programs in ornamental plants.

Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna

Development of flower color in plants is a complex process. Among others, it is an important trait for ornamental flowering plants. Canna is a flowering ornamental plant of family Cannaceae. To understand the molecular mechanism of flower color development in Canna, RNA sequencing from flower tissues of two contrasting flower color cultivars, Red President (RP) and Tropical Sunrise (TS) was performed. More than 27.0 million and 19.0 million clean reads were obtained from RP and TS, respectively. The combined clean reads were assembled into 147,295 unigenes. The Canna unigenes showed maximum homology with Populus trichocarpa (26.79%). A total of 2702 unigenes expressed differentially between the two cultivars of which 1972 were up-regulated and 730 were down-regulated in RP. Phenylpropanoid and flavonoid biosynthetic processes were the significant processes in RP. Expression of a vast number of transcription factors including MYB, bHLH, ARF, and WRKY were higher in RP than TS. The expression analysis of RNA sequencing data was validated by qRT-PCR analysis. Further, concentration of measured anthocyanidins and flavonols were very low or absent in TS, corroborating largely with our transcriptome data. These findings may help in understanding flower color development in Canna and in future crop breeding program.

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.

Small RNA Transcriptome of Hibiscus Syriacus Provides Insights into the Potential Influence of microRNAs in Flower Development and Terpene Synthesis

MicroRNAs (miRNAs) are essential small RNA molecules that regulate the expression of target mRNAs in plants and animals. Here, we aimed to identify miRNAs and their putative targets in Hibiscus syriacus, the national flower of South Korea. We employed high-throughput sequencing of small RNAs obtained from four different tissues (i.e., leaf, root, flower, and ovary) and identified 33 conserved and 30 novel miRNA families, many of which showed differential tissue-specific expressions. In addition, we computationally predicted novel targets of miRNAs and validated some of them using 5' rapid amplification of cDNA ends analysis. One of the validated novel targets of miR477 was a terpene synthase, the primary gene involved in the formation of disease-resistant terpene metabolites such as sterols and phytoalexins. In addition, a predicted target of conserved miRNAs, miR396, is SHORT VEGETATIVE PHASE, which is involved in flower initiation and is duplicated in H. syriacus. Collectively, this study provides the first reliable draft of the H. syriacus miRNA transcriptome that should constitute a basis for understanding the biological roles of miRNAs in H. syriacus.

Comparative transcriptomic analysis of the evolution and development of flower size in Saltugilia (Polemoniaceae)

BACKGROUND

Flower size varies dramatically across angiosperms, representing innovations over the course of >130 million years of evolution and contributing substantially to relationships with pollinators. However, the genetic underpinning of flower size is not well understood. Saltugilia (Polemoniaceae) provides an excellent non-model system for extending the genetic study of flower size to interspecific differences that coincide with variation in pollinators.

RESULTS

Using targeted gene capture methods, we infer phylogenetic relationships among all members of Saltugilia to provide a framework for investigating the genetic control of flower size differences via RNA-Seq de novo assembly. Nuclear concatenation and species tree inference methods provide congruent topologies. The inferred evolutionary trajectory of flower size is from small flowers to larger flowers. We identified 4 to 10,368 transcripts that are differentially expressed during flower development, with many unigenes associated with cell wall modification and components of the auxin and gibberellin pathways.

CONCLUSIONS

Saltugilia is an excellent model for investigating covarying floral and pollinator evolution. Four candidate genes from model systems (BIG BROTHER, BIG PETAL, GASA, and LONGIFOLIA) show differential expression during development of flowers in Saltugilia, and four other genes (FLOWERING-PROMOTING FACTOR 1, PECTINESTERASE, POLYGALACTURONASE, and SUCROSE SYNTHASE) fit into hypothesized organ size pathways. Together, these gene sets provide a strong foundation for future functional studies to determine their roles in specifying interspecific differences in flower size.

New technologies accelerate the exploration of non-coding RNAs in horticultural plants

Non-coding RNAs (ncRNAs), that is, RNAs not translated into proteins, are crucial regulators of a variety of biological processes in plants. While protein-encoding genes have been relatively well-annotated in sequenced genomes, accounting for a small portion of the genome space in plants, the universe of plant ncRNAs is rapidly expanding. Recent advances in experimental and computational technologies have generated a great momentum for discovery and functional characterization of ncRNAs. Here we summarize the classification and known biological functions of plant ncRNAs, review the application of next-generation sequencing (NGS) technology and ribosome profiling technology to ncRNA discovery in horticultural plants and discuss the application of new technologies, especially the new genome-editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) systems, to functional characterization of plant ncRNAs.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Towards development of new ornamental plants: status and progress in wide hybridization

The present review provides insights into the key findings of the hybridization process, crucial factors affecting the adaptation of new technologies within wide hybridization of ornamental plants and presents perspectives of further development of this strategy. Wide hybridization is one of the oldest breeding techniques that contributed enormously to the development of modern plant cultivars. Within ornamental breeding, it represents the main source of genetic variation. During the long history of wide hybridization, a number of methods were implemented allowing the evolution from a conventional breeding tool into a modern methodology. Nowadays, the research on model plants and crop species increases our understanding of reproductive isolation among distant species and partly explains the background of the traditional approaches previously used for overcoming hybridization barriers. Characterization of parental plants and hybrids is performed using molecular and cytological techniques that strongly facilitate breeding processes. Molecular markers and sequencing technologies are used for the assessment of genetic relationships among plants, as the genetic distance is typically depicted as one of the most important factors influencing cross-compatibility in hybridization processes. Furthermore, molecular marker systems are frequently applied for verification of hybrid state of the progeny. The flow cytometry and genomic in situ hybridization are used in the assessment of hybridization partners and characterization of hybrid progeny in relation to genome stabilization as well as genome recombination and introgression. In the future, new research and technologies are likely to provide more detailed information about genes and pathways responsible for interspecific reproductive isolation. Ultimately, this knowledge will enable development of strategies for obtaining compatible lines for hybrid production. Recent development in sequencing technologies and availability of sequence data will also facilitate creation of new molecular markers that will advance marker-assisted selection in hybridization process.

In silico search and biological validation of microRNAs related to drought response in peach and almond

Plant responses to drought stress are regulated at the transcriptional and post-transcriptional levels through noncoding endogenous microRNAs. These microRNAs play key roles in gene expression, mainly by down-regulating target mRNAs. In this work, an in silico search and validation for microRNAs related to drought response in peach ('G.H. Hill'), almond ('Sefied') and an interspecific peach-almond hybrid ('GN 15') has been performed. We used qPCR to analyse the gene expression of several miRNAs described as being related to drought response in peach, including miR156, miR159, miR160, miR167, miR171, miR172, miR398, miR403, miR408, miR842 and miR2275 under mild and severe water deficit. These miRNAs were in silico selected on the basis of previous works, their conservation in plants and their drought response. qPCR analysis confirmed the implication of these miRNAs in the dehydration stress response in the three assayed genotypes. Comparison of miRNA expression patterns in the three evaluated genotypes indicated that the hybrid GN 15 showed higher expression levels of specific miRNAs which should be related to the observed drought tolerance. mRNA target transcripts of the miRNAs studied were predicted using the Rose database, which includes transcription factors that regulate plant growth and development. In addition, results showed that the promoter region contains responsive elements to hormone-mediated regulatory elements. Network analysis not only unravelled the interaction between miRNAs and their predicted gene targets but also highlighted the roles of miRNAs in response to drought stress.

Identification and Expression Analyses of miRNAs from Two Contrasting Flower Color Cultivars of Canna by Deep Sequencing

miRNAs are endogenous small RNA (sRNA) that play critical roles in plant development processes. Canna is an ornamental plant belonging to family Cannaceae. Here, we report for the first time the identification and differential expression of miRNAs in two contrasting flower color cultivars of Canna, Tropical sunrise and Red president. A total of 313 known miRNAs belonging to 78 miRNA families were identified from both the cultivars. Thirty one miRNAs (17 miRNA families) were specific to Tropical sunrise and 43 miRNAs (10 miRNA families) were specific to Red president. Thirty two and 18 putative new miRNAs were identified from Tropical sunrise and Red president, respectively. One hundred and nine miRNAs were differentially expressed in the two cultivars targeting 1343 genes. Among these, 16 miRNAs families targeting60 genes were involved in flower development related traits and five miRNA families targeting five genes were involved in phenyl propanoid and pigment metabolic processes. We further validated the expression analysis of a few miRNA and their target genes by qRT-PCR. Transcription factors were the major miRNA targets identified. Target validation of a few randomly selected miRNAs by RLM-RACE was performed but was successful with only miR162. These findings will help in understanding flower development processes, particularly the color development in Canna.

Combined Effects of Retinoic Acid and Hydro-Alcoholic Extract of Rosa Damascena Mill on Wound in Diabetic Rats

Retinoic acid and Rosa damascena are compounds that have considerable effects in the cellular proliferation and synthesis of extracellular matrix. The present study was designed to assess the combined effects of retinoic acid and Rosa damascena mill on wound in diabetic rats. Seventy-two rats were used in this study. Diabetes was induced by a single intraperitoneal injection of streptozotocin (60 mg. Kg(-1)). Then, a full thickness wound was created on dorsal surface of all animals. After that, rats were divided, into three groups; control (normal saline), positive control (Phenytoin), and  case (combined of 0.1% Tretinoein lotion and hydro-alcoholic extract of Rosa damascena mill). Afterward, wounds were evaluated macroscopically and microscopically on days 5, 10 and 15. Macroscopic and microscopic evaluations showed a significant improvement (p<0.05) of wounds in case group on 5(th) and 10(th) days when compared to positive control and control groups. The combination of Retinoic acid and hydro-alcholic extract of Rosa damascena mill can accelerate wound healing in diabetic rats.

Transcriptomic Analysis of Paeonia delavayi Wild Population Flowers to Identify Differentially Expressed Genes Involved in Purple-Red and Yellow Petal Pigmentation

Tree peony (Paeonia suffruticosa Andrews) is a very famous traditional ornamental plant in China. P. delavayi is a species endemic to Southwest China that has aroused great interest from researchers as a precious genetic resource for flower color breeding. However, the current understanding of the molecular mechanisms of flower pigmentation in this plant is limited, hindering the genetic engineering of novel flower color in tree peonies. In this study, we conducted a large-scale transcriptome analysis based on Illumina HiSeq sequencing of cDNA libraries generated from yellow and purple-red P. delavayi petals. A total of 90,202 unigenes were obtained by de novo assembly, with an average length of 721 nt. Using Blastx, 44,811 unigenes (49.68%) were found to have significant similarity to accessions in the NR, NT, and Swiss-Prot databases. We also examined COG, GO and KEGG annotations to better understand the functions of these unigenes. Further analysis of the two digital transcriptomes revealed that 6,855 unigenes were differentially expressed between yellow and purple-red flower petals, with 3,430 up-regulated and 3,425 down-regulated. According to the RNA-Seq data and qRT-PCR analysis, we proposed that four up-regulated key structural genes, including F3H, DFR, ANS and 3GT, might play an important role in purple-red petal pigmentation, while high co-expression of THC2'GT, CHI and FNS II ensures the accumulation of pigments contributing to the yellow color. We also found 50 differentially expressed transcription factors that might be involved in flavonoid biosynthesis. This study is the first to report genetic information for P. delavayi. The large number of gene sequences produced by transcriptome sequencing and the candidate genes identified using pathway mapping and expression profiles will provide a valuable resource for future association studies aimed at better understanding the molecular mechanisms underlying flower pigmentation in tree peonies.

Using RNA-Seq to assemble a rose transcriptome with more than 13,000 full-length expressed genes and to develop the WagRhSNP 68k Axiom SNP array for rose (Rosa L.)

In order to develop a versatile and large SNP array for rose, we set out to mine ESTs from diverse sets of rose germplasm. For this RNA-Seq libraries containing about 700 million reads were generated from tetraploid cut and garden roses using Illumina paired-end sequencing, and from diploid Rosa multiflora using 454 sequencing. Separate de novo assemblies were performed in order to identify single nucleotide polymorphisms (SNPs) within and between rose varieties. SNPs among tetraploid roses were selected for constructing a genotyping array that can be employed for genetic mapping and marker-trait association discovery in breeding programs based on tetraploid germplasm, both from cut roses and from garden roses. In total 68,893 SNPs were included on the WagRhSNP Axiom array. Next, an orthology-guided assembly was performed for the construction of a non-redundant rose transcriptome database. A total of 21,740 transcripts had significant hits with orthologous genes in the strawberry (Fragaria vesca L.) genome. Of these 13,390 appeared to contain the full-length coding regions. This newly established transcriptome resource adds considerably to the currently available sequence resources for the Rosaceae family in general and the genus Rosa in particular.

Disease resistance breeding in rose: current status and potential of biotechnological tools

The cultivated rose is a multispecies complex for which a high level of disease protection is needed due to the low tolerance of blemishes in ornamental plants. The most important fungal diseases are black spot, powdery mildew, botrytis and downy mildew. Rose rosette, a lethal viral pathogen, is emerging as a devastating disease in North America. Currently rose breeders use a recurrent phenotypic selection approach and perform selection for disease resistance for most pathogen issues in a 2-3 year field trial. Marker assisted selection could accelerate this breeding process. Thus far markers have been identified for resistance to black spot (Rdrs) and powdery mildew and with the ability of genotyping by sequencing to generate 1000s of markers our ability to identify markers useful in plant improvement should increase exponentially. Transgenic rose lines with various fungal resistance genes inserted have shown limited success and RNAi technology has potential to provide virus resistance. Roses, as do other plants, have sequences homologous to characterized R-genes in their genomes, some which have been related to specific disease resistance. With improving next generation sequencing technology, our ability to do genomic and transcriptomic studies of the resistance related genes in both the rose and the pathogens to reveal novel gene targets to develop resistant roses will accelerate. Finally, the development of designer nucleases opens up a potentially non-GMO approach to directly modify a rose's DNA to create a disease resistant rose. Although there is much potential, at present rose breeders are not using marker assisted breeding primarily because a good suite of marker/trait associations (MTA) that would ensure a path to stable disease resistance is not available. As our genomic analytical tools improve, so will our ability to identify useful genes and linked markers. Once these MTAs are available, it will be the cost savings, both in time and money, that will convince the breeders to use the technology.

Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

The transcript pool of a plant part, under any given condition, is a collection of mRNAs that will pave the way for a biochemical reaction of the plant to stimuli. Over the past decades, transcriptome study has advanced from Northern blotting to RNA sequencing (RNA-seq), through other techniques, of which real-time quantitative polymerase chain reaction (PCR) and microarray are the most significant ones. The questions being addressed by such studies have also matured from a solitary process to expression atlas and marker-assisted genetic enhancement. Not only genes and their networks involved in various developmental processes of plant parts have been elucidated, but also stress tolerant genes have been highlighted. The transcriptome of a plant with altered expression of a target gene has given information about the downstream genes. Marker information has been used for breeding improved varieties. Fortunately, the data generated by transcriptome analysis has been made freely available for ample utilization and comparison. The review discusses this wide variety of transcriptome data being generated in plants, which includes developmental stages, abiotic and biotic stress, effect of altered gene expression, as well as comparative transcriptomics, with a special emphasis on microarray and RNA-seq. Such data can be used to determine the regulatory gene networks, which can subsequently be utilized for generating improved plant varieties.

Transcriptome analysis by Illumina high-throughout paired-end sequencing reveals the complexity of differential gene expression during in vitro plantlet growth and flowering in Amaranthus tricolor L

Amaranthus tricolor L. is a C4 plant, which is consumed as a major leafy vegetable in some tropical countries. Under conditions of high temperature and short daylight, Am. tricolor readily bolts and blooms, degrading leaf quality. A preliminary in vitro flowering study demonstrated that the flowering control pathway in Am. tricolor may differ from that of Arabidopsis. Nevertheless, no transcriptome analysis of the flowering process in Amaranthus has been conducted. To study Am. tricolor floral regulatory mechanisms, we conducted a large-scale transcriptome analysis--based on Illumina HiSeq sequencing of cDNA libraries generated from Am. tricolor at young seedling (YSS), adult seedling (ASS), flower bud (FBS), and flowering (FS) stages. A total of 99,312 unigenes were obtained. Using BLASTX, 43,088 unigenes (43.39%) were found to have significant similarity with accessions in Nr, Nt, and Swiss-Prot databases. Of these unigenes, 11,291 were mapped to 266 KEGG pathways. Further analysis of the four digital transcriptomes revealed that 735, 17,184, 274, and 206 unigenes were specifically expressed during YSS, ASS, FBS, and FS, respectively, with 59,517 unigenes expressed throughout the four stages. These unigenes were involved in many metabolic pathways related to in vitro flowering. Among these pathways, 259 unigenes were associated with ubiquitin-mediated proteolysis, indicating its importance for in vitro flowering in Am. tricolor. Other pathways, such as circadian rhythm and cell cycle, also had important roles. Finally, 26 unigenes were validated by qRT-PCR in samples from Am. tricolor at YSS, ASS, FBS, and FS; their differential expressions at the various stages indicate their possible roles in Am. tricolor growth and development, but the results were somewhat similar to Arabidopsis. Because unigenes involved in many metabolic pathways or of unknown function were revealed to regulate in vitro plantlet growth and flowering in Am. tricolor, the process appears to be highly complex in this species.

TRV-GFP: a modified Tobacco rattle virus vector for efficient and visualizable analysis of gene function

Virus-induced gene silencing (VIGS) is a useful tool for functional characterization of genes in plants. Unfortunately, the efficiency of infection by Tobacco rattle virus (TRV) is relatively low for some non-Solanaceae plants, which are economically important, such as rose (Rosa sp.). Here, to generate an easy traceable TRV vector, a green fluorescent protein (GFP) gene was tagged to the 3' terminus of the coat protein gene in the original TRV2 vector, and the silencing efficiency of the modified TRV-GFP vector was tested in several plants, including Nicotiana benthamiana, Arabidopsis thaliana, rose, strawberry (Fragaria ananassa), and chrysanthemum (Dendranthema grandiflorum). The results showed that the efficiency of infection by TRV-GFP was equal to that of the original TRV vector in each tested plant. Spread of the modified TRV virus was easy to monitor by using fluorescent microscopy and a hand-held UV lamp. When TRV-GFP was used to silence the endogenous phytoene desaturase (PDS) gene in rose cuttings and seedlings, the typical photobleached phenotype was observed in 75-80% plants which were identified as GFP positive by UV lamp. In addition, the abundance of GFP protein, which represented the concentration of TRV virus, was proved to correlate negatively with the level of the PDS gene, suggesting that GFP could be used as an indicator of the degree of silencing of a target gene. Taken together, this work provides a visualizable and efficient tool to predict positive gene silencing plants, which is valuable for research into gene function in plants, especially for non-Solanaceae plants.

Filling gaps with construction of a genetic linkage map in tetraploid roses

Rose (Rosa sp.) is one of the most economically important ornamental crops worldwide. The present work contains a genetic linkage map for tetraploid roses that was constructed from an F1 segregation population using AFLPs and SSRs on 189 individuals. The preliminary 'Yunzheng Xiawei' and 'Sun City' maps consisted of 298 and 255 markers arranged into 26 and 32 linkage groups, respectively. The recombined parental maps covered 737 and 752 cM of the genome, respectively. The integrated linkage map was composed of 295 polymorphic markers that spanned 874 cM, and it had a mean intermarker distance of 2.9 cM. In addition, a set of newly developed EST-SSRs that are distributed evenly throughout the mapping population were released. The work identified 67 anchoring points that came from 43 common SSRs. The results that were produced from a large number of individuals (189) and polymorphic SSRs (242) will enhance the ability to construct higher density consensus maps with the available diploid level rose maps, and they will definitely serve as a tool for accurate QTL detection and marker assisted selection.

Molecular genetics and genomics of the Rosoideae: state of the art and future perspectives

The Rosoideae is a subfamily of the Rosaceae that contains a number of species of economic importance, including the soft fruit species strawberry (Fragaria ×ananassa), red (Rubus idaeus) and black (Rubus occidentalis) raspberries, blackberries (Rubus spp.) and one of the most economically important cut flower genera, the roses (Rosa spp.). Molecular genetics and genomics resources for the Rosoideae have developed rapidly over the past two decades, beginning with the development and application of a number of molecular marker types including restriction fragment length polymorphisms, amplified fragment length polymorphisms and microsatellites, and culminating in the recent publication of the genome sequence of the woodland strawberry, Fragaria vesca, and the development of high throughput single nucleotide polymorphism (SNP)-genotyping resources for Fragaria, Rosa and Rubus. These tools have been used to identify genes and other functional elements that control traits of economic importance, to study the evolution of plant genome structure within the subfamily, and are beginning to facilitate genomic-assisted breeding through the development and deployment of markers linked to traits such as aspects of fruit quality, disease resistance and the timing of flowering. In this review, we report on the developments that have been made over the last 20 years in the field of molecular genetics and structural genomics within the Rosoideae, comment on how the knowledge gained will improve the efficiency of cultivar development and discuss how these advances will enhance our understanding of the biological processes determining agronomically important traits in all Rosoideae species.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

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.

The Antioxidants Changes in Ornamental Flowers during Development and Senescence

The concentration of antioxidant compounds is constitutive and variable from species to species and is also variable considering the development of the plant tissue. In this review, we take into consideration the antioxidant changes and the physiological, biochemical and molecular factors that are able to modulate the accumulation of antioxidant compounds in ornamental flowers during the whole development process until the senescence. Many ornamental flowers are natural sources of very important bioactive compounds with benefit to the human health and their possible role as dietary components has been reported. The most part of antioxidants are flower pigments such as carotenoids and polyphenols, often present in higher concentration compared with the most common fruits and vegetables. The antioxidants content changes during development and during senescence many biochemical systems and molecular mechanisms are activated to counteract the increase of reactive oxygen species and free radicals. There is a tight correlation between antioxidants and senescence processes and this aspect is detailed and appropriately discussed.

Transcriptome/degradome-wide identification of R. glutinosa miRNAs and their targets: the role of miRNA activity in the replanting disease

Rehmannia glutinosa, a traditional Chinese medicine herb, is unable to grow normally in a soil where the same species has recently been cultivated. The biological basis of this so called "replanting disease" is unknown, but it may involve the action of microRNAs (miRNAs), which are known to be important regulators of plant growth and development. High throughput Solexa/Illumina sequencing was used to generate a transcript library of the R. glutinosa transcriptome and degradome in order to identify possible miRNAs and their targets implicated in the replanting disease. A total of 87,665 unigenes and 589 miRNA families (17 of which have not been identified in plants to date) was identified from the libraries made from a first year (FP) and a second year (SP) crop. A comparison between the FP and SP miRNAs showed that the abundance of eight of the novel and 295 of the known miRNA families differed between the FP and SP plants. Sequencing of the degradome sampled from FP and SP plants led to the identification of 165 transcript targets of 85 of the differentially abundant miRNA families. The interaction of some of these miRNAs with their target(s) is likely to form an important part of the molecular basis of the replanting disease of R. glutinosa.

The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: a foundation for understanding MicroRNA functional roles in hot pepper

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5′ rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.

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.

MiRAuto: an automated user-friendly microRNA prediction tool utilizing plant small RNA sequencing data

MicroRNAs (miRNAs) are a class of small RNAs that post-transcriptionally regulate gene expression in animals and plants. The recent rapid advancement in miRNA biology, including high-throughput sequencing of small RNA libraries, inspired the development of a bioinformatics software, miRAuto, which predicts putative miRNAs in model plant genomes computationally. Furthermore, miRAuto enables users to identify miRNAs in non-model plant species whose genomes have yet to be fully sequenced. miRAuto analyzes the expression of the 5'-end position of mapped small RNAs in reference sequences to prevent the possibility of mRNA fragments being included as candidate miRNAs. We validated the utility of miRAuto on a small RNA dataset, and the results were compared to other publicly available miRNA prediction programs. In conclusion, miRAuto is a fully automated user-friendly tool for predicting miRNAs from small RNA sequencing data in both model and non-model plant species. miRAuto is available at http://nature.snu.ac.kr/software/miRAuto.htm.