Transcriptome analysis of carnation (Dianthus caryophyllus L.) based on next-generation sequencing technology

Distinguish Dianthus species or varieties based on chloroplast genomes

Most plants belonging to the widely distributed genus Dianthus are used for gardening. Interspecific hybridization of different Dianthus species leads to blurred genetic backgrounds. To obtain more genomic resources and understand the phylogenetic relationships among Dianthus species, the chloroplast genomes of 12 Dianthus species, including nine Dianthus gratianopolitanus varieties, were analyzed. The chloroplast genomes of these 12 species exhibited similar sizes (149,474-149,735 bp), with Dianthus caryophyllus having a chloroplast genome size of 149,604 bp marked by a significant contraction in inverted repeats. In the chloroplast genome of Dianthus, we identified 124-126 annotated genes, including 83-84 protein-coding genes. Notably, D. caryophyllus had 83 protein-coding genes but lacked rpl2. The repeat sequences of the chloroplast genome were consistent among species, and variations in the sequence were limited and not prominent. However, notable gene replacements were observed in the boundary region. Phylogenetic analysis of Dianthus indicated that D. caryophyllus and D. gratianopolitanus were most closely related, suggesting that the degree of variation within nine Dianthus varieties was no less than the variation observed between species. These differences provide a theoretical foundation for a more comprehensive understanding of the diversity within Dianthus species.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

Low-Oxygen Responses of Cut Carnation Flowers Associated with Modified Atmosphere Packaging

Gaseous factors affect post-harvest physiological processes in horticultural crops, including ornamental flowers. However, the molecular responses of cut flowers to the low-oxygen conditions associated with modified atmosphere packaging (MAP) have not yet been elucidated. Here, we show that storage of cut carnation flowers in a sealed polypropylene bag decreased the oxygen concentration in the bag to 3-5% and slowed flower opening. The vase life of carnation flowers after storage for seven days under MAP conditions was comparable to that without storage and was improved by the application of a commercial-quality preservative. The adenylate energy charge (AEC) was maintained at high levels in petals from florets stored under MAP conditions. This was accompanied by the upregulation of four hypoxia-related genes, among which the HYPOXIA-RESPONSIVE ETHYLENE RESPONSE FACTOR and PHYTOGLOBIN genes (DcERF19 and DcPGB1) were newly identified. These results suggest that hypoxia-responsive genes contribute to the maintenance of the energy status in carnation flowers stored under MAP conditions, making this gas-controlling technique potentially effective for maintaining cut flower quality without cooling.

Floral Scents and Fruit Aromas: Functions, Compositions, Biosynthesis, and Regulation

Floral scents and fruit aromas are crucial volatile organic compounds (VOCs) in plants. They are used in defense mechanisms, along with mechanisms to attract pollinators and seed dispersers. In addition, they are economically important for the quality of crops, as well as quality in the perfume, cosmetics, food, drink, and pharmaceutical industries. Floral scents and fruit aromas share many volatile organic compounds in flowers and fruits. Volatile compounds are classified as terpenoids, phenylpropanoids/benzenoids, fatty acid derivatives, and amino acid derivatives. Many genes and transcription factors regulating the synthesis of volatiles have been discovered. In this review, we summarize recent progress in volatile function, composition, biosynthetic pathway, and metabolism regulation. We also discuss unresolved issues and research perspectives, providing insight into improvements and applications of plant VOCs.

DcWRKY75 promotes ethylene induced petal senescence in carnation (Dianthus caryophyllus L.)

Carnation (Dianthus caryophyllus L.) is one of the most important and typical ethylene sensitive cut flowers worldwide, although how ethylene influences the petal senescence process in carnation remains largely unknown. Here, we screened out one of the key transcription factors, DcWRKY75, using a constructed ethylene induced petal senescence transcriptome in carnation and found that it shows quick induction by ethylene treatment. Silencing of DcWRKY75 delays ethylene induced petal senescence in carnation. Molecular evidence confirms that DcWRKY75 can bind to the promoter regions of two main ethylene biosynthetic genes (DcACS1 and DcACO1) and a couple of senescence associated genes (DcSAG12 and DcSAG29) to activate their expression. Furthermore, we show that DcWRKY75 is a direct target gene of DcEIL3-1, which is a homolog of the ethylene signaling core transcription factor EIN3 in Arabidopsis. DcEIL3-1 can physically interact with DcWRKY75 and silencing of DcEIL3-1 also delays ethylene induced petal senescence in carnation and inhibits the ethylene induced expression of DcWRKY75 and its target genes. The present study demonstrates that the transcriptional regulation network is vitally important for ethylene induced petal senescence process in carnation and potentially in other ethylene sensitive cut flowers.

Transcriptome of Cenostigma pyramidale roots, a woody legume, under different salt stress times

Salinity stress has a significant impact on the gain of plant biomass. Our study provides the first root transcriptome of Cenostigma pyramidale, a tolerant woody legume from a tropical dry forest, under three different salt stress times (30 min, 2 h, and 11 days). The transcriptome was assembled using the RNA sequencing (RNA-Seq) de novo pipeline from GenPipes. We observed 932, 804, and 3157 upregulated differentially expressed genes (DEGs) and 164, 273, and 1332 downregulated DEGs for salt over 30 min, 2 h, and 11 days, respectively. For DEGs annotated with the Viridiplantae clade in the early stress periods, the response to salt stress was mainly achieved by stabilizing homeostasis of such ions like Na⁺ and K⁺ , signaling by Ca²⁺ , transcription factor modulation, water transport, and oxidative stress. For salt stress at 11 days, we observed a higher modulation of transcription factors including the WRKY, MYB, bHLH, NAC, HSF, and AP2-EREBP families, as well as DEGs involved in hormonal responses, water transport, sugar metabolism, proline, and reactive oxygen scavenging mechanisms. Five selected DEGs (K⁺ transporter, aquaporin, glutathione S-transferase, cyclic nucleotide-gated channel, and superoxide dismutase) were validated by qPCR. Our results indicated that C. pyramidale had an early perception of salt stress modulating ionic channels and transporters, and as the stress progressed, the focus turned to the antioxidant system, aquaporins, and complex hormone responses. The results of this first root transcriptome provide clues on how this native species modulate gene expression to achieve salt stress tolerance.

Recent Progress in Enhancing Fungal Disease Resistance in Ornamental Plants

Fungal diseases pose a major threat to ornamental plants, with an increasing percentage of pathogen-driven host losses. In ornamental plants, management of the majority of fungal diseases primarily depends upon chemical control methods that are often non-specific. Host basal resistance, which is deficient in many ornamental plants, plays a key role in combating diseases. Despite their economic importance, conventional and molecular breeding approaches in ornamental plants to facilitate disease resistance are lagging, and this is predominantly due to their complex genomes, limited availability of gene pools, and degree of heterozygosity. Although genetic engineering in ornamental plants offers feasible methods to overcome the intrinsic barriers of classical breeding, achievements have mainly been reported only in regard to the modification of floral attributes in ornamentals. The unavailability of transformation protocols and candidate gene resources for several ornamental crops presents an obstacle for tackling the functional studies on disease resistance. Recently, multiomics technologies, in combination with genome editing tools, have provided shortcuts to examine the molecular and genetic regulatory mechanisms underlying fungal disease resistance, ultimately leading to the subsequent advances in the development of novel cultivars with desired fungal disease-resistant traits, in ornamental crops. Although fungal diseases constitute the majority of ornamental plant diseases, a comprehensive overview of this highly important fungal disease resistance seems to be insufficient in the field of ornamental horticulture. Hence, in this review, we highlight the representative mechanisms of the fungal infection-related resistance to pathogens in plants, with a focus on ornamental crops. Recent progress in molecular breeding, genetic engineering strategies, and RNAi technologies, such as HIGS and SIGS for the enhancement of fungal disease resistance in various important ornamental crops, is also described.

Comprehensive analysis of sucrolytic enzyme gene families in carnation (Dianthus caryophyllus L.)

Sucrose plays crucial roles in growth and responses of plants to the environment, including those in ornamental species. During post-harvest handling of cut flowers, sucrose degradation is an essential process of inter- and intra-cellular carbon partitioning affecting flower opening and senescence and, subsequently, flower quality. However, complete information about the molecular basis of sucrose degradation in ornamental flowers, which can be catalyzed by two kinds of sucrolytic enzymes, invertase (INV), and sucrose synthase (SUS), is not available from past reports. The present study shows that sucrose treatment of carnation (Dianthus caryophyllus L.) florets increased starch content in petals, accompanied by decreased vacuolar INV (VIN) activity and increased SUS activity. However, hypoxic treatment of carnation florets decreased sucrose content and cell-wall INV (CWIN) activity in petals. In silico analysis using the carnation genome database identified six CWIN, three VIN, eight cytoplasmic INV (CIN), and five SUS genes. Real-time RT-PCR analysis confirmed that these genes are differentially expressed in carnation petals in response to sucrose and hypoxic treatments, partially corresponding to the changes in enzyme activities. In contrast to DcSUS1 (Dca4507.1), a SUS gene already reported in carnation, which showed preferential expression under aerated conditions, the expression of DcSUS2 (Dca22218.1), an undescribed carnation SUS gene, was enhanced under hypoxia similarly to an alcohol dehydrogenase gene DcADH1 (Dca18671.1). These results suggest that sugar metabolism in carnation petals is regulated in response to environmental cues, accompanied by modulated activities and gene expression of a set of sucrolytic enzymes.

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.

Integrating physiological and metabolites analysis to identify ethylene involvement in petal senescence in Tulipa gesneriana

Flower senescence is classified into ethylene-dependent and ethylene-independent manners and determines the flower longevity which is valuable for ornamental plants. However, the manner of petal senescence in tulip is still less defined. In this study, we characterized the physiological indexes in the process of petal senescence, as well as metabolic and ethylene responses in tulip cultivar 'American Dream', and further identified the role of ethylene biosynthesis genes TgACS by transgenic and transient assays. Primary metabolites profiling revealed that sugars, amino acids and organic acids preferentially accumulated in senescent petals. Additionally, senescence-associated genes were identified and significantly up-regulated, coupled with increased ROS contents, rapid water loss and accelerated cell membrane breakdown. Moreover, ethylene production was stimulated as evidenced by increasing in ACS activity and ethylene biosynthesis-related genes expression. Exogenous treatment of cutting flowers with 1-MCP or ethephon resulted in delayed or enhanced petal senescence, respectively. Transient down-regulation of TgACS by VIGS assay in tulip petals delayed senescence, while over-expressed TgACS1 in tobacco promoted leaf senescence. Taken together, this study provides evidences to certify ethylene roles and TgACS functions during flower senescence in tulip.

Transcriptome assembly and abiotic related gene expression analysis of date palm reveal candidate genes involved in response to cadmium stress

Date palm in Tunisia is of major economic importance but are also factors of social, environmental and economic stability. An annotated assembly of the transcriptome of cultivar Deglet Nour was reported. RNA was isolated from plant Cd-contaminated leaves, and 37,049 unique Illumina RNA-seq reads were used in a transcriptome assembly. The draft transcriptome assembly consists of 6789 contigs and 17.285 singletons with a means length of 858 bp and 1.042 bp, respectively. The final assembly was functionally annotated using Blast2GO software, allowing the identification of putative genes controlling important agronomic traits. The annotated transcriptome data sets were used to query all known Kyoto Encyclopedia of Genes and Genomes pathways. The most represented molecular functions and biological processes were nucleotide binding and transcription, transport and response to stress and abiotic and biotic stimuli. A prediction of the genes interaction network was proposed by selecting corresponding functionally similar genes from Arabidopsis datasets, downloaded by GeneMANIA version 2.1. Several Cd-responsive genes expression was monitored in in vitro isolated explant of Cd stressed Deglet Nour. Some chelators encoding genes were upregulated confirming in silico findings. Genes encoding HMs transporters in date palm showed expression enhancement more pronounced after 20 days of exposure. P. dactylifera transcriptome provides a valuable resource for future functional analysis of candidate genes involved in metal stress response.

The study of transcriptome sequencing for flower coloration in different anthesis stages of alpine ornamental herb (Meconopsis 'Lingholm')

Meconopsis (Papaveraceae) is an interesting alpine herb, mainly distributed in the mountainous area of southwest China and high altitude zone in Tibetan-Himalaya. Different Meconopsis species showed a flower color alteration in different anthesis stages, Meconopsis 'Lingholm' is one of the localized species whose petal color changes from purple to blue during the flowering process. In general, the blue color flower is a rare kind, and usually hard to cultivate artificially. The molecular mechanism of flower color formation and color alteration of alpine flowers were reported by many research workers. To find critical genes that regulate Meconopsis 'Lingholm' color alteration and the mechanism of environmental adaptation, the current study performed transcriptome sequencing by using Meconopsis 'Lingholm' petals from different anthesis stages. There were totally 91,615 unigenes obtained from 31.4 Gb sequencing data, and differentially expressed genes between two consecutive flowering stages were obtained. Bioinformatics studies showed genes regulating petal color alteration were activated. Moreover, the functional analysis showed that Meconopsis 'Lingholm' showed a stress response to mechanical damage, non-biological stimulation and water deficiency in the bud stage, as well as showed a stress response to the cold from cracking stage to blooming stage. Furthermore, RNA-Seq results were verified using nine randomly selected genes by qPCR, which showed same expression trend with sequencing results. During this study, 20 candidate genes identified for further studies, which included five petal color related genes and 15 environmental response genes.

Cysmotif Searcher Pipeline for Antimicrobial Peptide Identification in Plant Transcriptomes

In this paper, we present the new Cysmotif searcher pipeline for identification of various antimicrobial peptides (AMPs), the most important components of innate immunity, in plant transcriptomes. Cysmotif searcher reveals and classifies short cysteine-rich amino acid sequences containing an open reading frame and a signal peptide cleavage site. Due to the combination of various search methods, Cysmotif searcher allows to obtain the most complete repertoire of AMPs for one or more transcriptomes in a short amount of time. The pipeline performance is estimated on the model plant Arabidopsis thaliana and nine other plants, including cultivated and wild species. The obtained results are compared to the existing annotation (A. thaliana) and results of conventional homology search (other plants). The comparison is carried out for known families of plant AMPs and newly discovered peptides that could not be assigned to existing families. The applicability of Cysmotif searcher in detecting new AMPs is discussed, and some practical recommendations on the pipeline usage for end users are given. The Cysmotif searcher pipeline is free for academic use and can be downloaded from Github (http://github.com/fallandar/cysmotifsearcher).

Genome-Wide Identification and Characterization of Aquaporins and Their Role in the Flower Opening Processes in Carnation (Dianthus caryophyllus)

Aquaporins (AQPs) are associated with the transport of water and other small solutes across biological membranes. Genome-wide identification and characterization will pave the way for further insights into the AQPs’ roles in the commercial carnation (Dianthus caryophyllus). This study focuses on the analysis of AQPs in carnation (DcaAQPs) involved in flower opening processes. Thirty DcaAQPs were identified and grouped to five subfamilies: nine PIPs, 11 TIPs, six NIPs, three SIPs, and one XIP. Subsequently, gene structure, protein motifs, and co-expression network of DcaAQPs were analyzed and substrate specificity of DcaAQPs was predicted. qRT-PCR, RNA-seq, and semi-qRTRCR were used for DcaAQP genes expression analysis. The analysis results indicated that DcaAQPs were relatively conserved in gene structure and protein motifs, that DcaAQPs had significant differences in substrate specificity among different subfamilies, and that DcaAQP genes’ expressions were significantly different in roots, stems, leaves and flowers. Five DcaAQP genes (DcaPIP1;3, DcaPIP2;2, DcaPIP2;5, DcaTIP1;4, and DcaTIP2;2) might play important roles in flower opening process. However, the roles they play are different in flower organs, namely, sepals, petals, stamens, and pistils. Overall, this study provides a theoretical basis for further functional analysis of DcaAQPs.

Characterization of increased cuticular wax mutant and analysis of genes involved in wax biosynthesis in Dianthus spiculifolius

Cuticular wax formation on the surface of plant leaves is associated with drought-stress tolerance. The identification of wax biosynthesis-related genes will contribute to the genetic improvement of drought resistance in plants. In this study, we characterize a novel Dianthus spiculifolius mutant with increased cuticular wax. The mutant exhibited stronger drought resistance as indicated by less leaf wilting and death, higher leaf relative water content and water retention capacity, and slower water loss and chlorophyll extraction than did the wild type during drought treatment. In the mutant leaves, 2 730 upregulated and 2 151 downregulated differentially expressed genes (DEGs) were identified by transcriptome sequencing. A wax biosynthesis pathway of the identified DEGs was significantly enriched. Finally, three key genes (DsCER1, DsMAH1, and DsWSD1) involved in wax biosynthesis were identified and verified by qPCR. These results suggest that differential expression of DEGs involved in wax biosynthesis may be associated with the increase in cuticular wax in the mutant. Taken together, our results help elucidate wax formation patterns in D. spiculifolius. Furthermore, the DEGs involved in wax biosynthesis identified here may be valuable genetic resources for improving plant stress tolerance through increased accumulation of cuticular wax.

Recent progress in whole genome sequencing, high-density linkage maps, and genomic databases of ornamental plants

Genome information is useful for functional analysis of genes, comparative genomic analysis, breeding of new varieties by marker-assisted selection, and map-based gene isolation. Genome-related research in ornamentals plants has been relatively slow to develop because of their heterozygosity or polyploidy. Advances in analytical instruments, such as next-generation sequencers and information processing technologies have revolutionized biology, and have been applied in a large number and variety of species, including ornamental plants. Recently, high-quality whole genome sequences have been reported in plant genetics and physiology studies of model ornamentals, such as those in genus Petunia and Japanese morning glory (Ipomoea nil). In this review, whole genome sequencing and construction of high-density genetic linkage maps based on SNP markers of ornamentals will be discussed. The databases that store this information for ornamentals are also described.

Molecular aspects of flower senescence and strategies to improve flower longevity

Flower longevity is one of the most important traits for ornamental plants. Ethylene plays a crucial role in flower senescence in some plant species. In several species that show ethylene-dependent flower senescence, genetic modification targeting genes for ethylene biosynthesis or signaling has improved flower longevity. Although little is known about regulatory mechanisms of petal senescence in flowers that show ethylene-independent senescence, a recent study of Japanese morning glory revealed that a NAC transcription factor, EPHEMERAL1 (EPH1), is a key regulator in ethylene-independent petal senescence. EPH1 is induced in an age-dependent manner irrespective of ethylene signal, and suppression of EPH1 expression dramatically delays petal senescence. In ethylene-dependent petal senescence, comprehensive transcriptome analyses revealed the involvement of transcription factors, a basic helix-loop-helix protein and a homeodomain-leucine zipper protein, in the transcriptional regulation of the ethylene biosynthesis enzymes. This review summarizes molecular aspects of flower senescence and discusses strategies to improve flower longevity by molecular breeding.

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.

Defense peptide repertoire of Stellaria media predicted by high throughput next generation sequencing

Being perfectly adapted to diverse environments, chickweed (Stellaria media (L.) Vill), a ubiquitous garden weed, grows widely in Europe and North America. As opposed to the model plants, many weeds, and S. media in particular, have been poorly studied, although they are likely to contain promising components of immunity and novel resistance genes. In this study, for the first time RNA-seq analysis of healthy and infected with Fusarium oxysporum chickweed seedlings, as well as de novo transcriptome assembly and annotation, are presented. Note, this research is focused on antimicrobial peptides (AMPs), the major components of plant immune system. Using custom software developed earlier, 145 unique putative AMPs (pAMPs) including defensins, thionins, hevein-like peptides, snakins, alpha-hairpinins, LTPs, and cysteine-rich peptides with novel cysteine motifs were predicted. Furthermore, changes in AMP expression profile in response to fungal infection were traced. In addition, the comparison of chickweed AMP repertoire with those of other Caryophyllaceae plants whose transcriptomes are presently available is made. As a result, alpha-hairpinins and hevein-like peptides which display characteristic modular structure appear to be specific AMPs distinguishing S. media from Dianthus caryophyllus, Silene vulgaris, and Silene latifolia. Finally, revealing several AMPs with proven antimicrobial activity gives opportunity to conclude that the presented method of AMP repertoire analysis reveals highly active AMPs playing vital role in plant immunity.

Transcriptome and Metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera

Calotropis procera is a medicinal plant of immense importance due to its pharmaceutical active components, especially cardiac glycosides (CG). As genomic resources for this plant are limited, the genes involved in CG biosynthetic pathway remain largely unknown till date. Our study on stage and tissue specific metabolite accumulation showed that CG's were maximally accumulated in stems of 3 month old seedlings. De novo transcriptome sequencing of same was done using high throughput Illumina HiSeq platform generating 44074 unigenes with average mean length of 1785 base pair. Around 66.6% of unigenes were annotated by using various public databases and 5324 unigenes showed significant match in the KEGG database involved in 133 different pathways of plant metabolism. Further KEGG analysis resulted in identification of 336 unigenes involved in cardenolide biosynthesis. Tissue specific expression analysis of 30 putative transcripts involved in terpenoid, steroid and cardenolide pathways showed a positive correlation between metabolite and transcript accumulation. Wound stress elevated CG levels as well the levels of the putative transcripts involved in its biosynthetic pathways. This result further validated the involvement of identified transcripts in CGs biosynthesis. The identified transcripts will lay a substantial foundation for further research on metabolic engineering and regulation of cardiac glycosides biosynthesis pathway genes.

Current status and biotechnological advances in genetic engineering of ornamental plants

Cut flower markets are developing in many countries as the international demand for cut flowers is rapidly growing. Developing new varieties with modified characteristics is an important aim in floriculture. Production of transgenic ornamental plants can shorten the time required in the conventional breeding of a cultivar. Biotechnology tools in combination with conventional breeding methods have been used by cut flower breeders to change flower color, plant architecture, post-harvest traits, and disease resistance. In this review, we describe advances in genetic engineering that have led to the development of new cut flower varieties.

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.

Transcriptome Analysis of Two Vicia sativa Subspecies: Mining Molecular Markers to Enhance Genomic Resources for Vetch Improvement

The vetch (Vicia sativa) is one of the most important annual forage legumes globally due to its multiple uses and high nutritional content. Despite these agronomical benefits, many drawbacks, including cyano-alanine toxin, has reduced the agronomic value of vetch varieties. Here, we used 454 technology to sequence the two V. sativa subspecies (ssp. sativa and ssp. nigra) to enrich functional information and genetic marker resources for the vetch research community. A total of 86,532 and 47,103 reads produced 35,202 and 18,808 unigenes with average lengths of 735 and 601 bp for V. sativa sativa and V. sativa nigra, respectively. Gene Ontology annotations and the cluster of orthologous gene classes were used to annotate the function of the Vicia transcriptomes. The Vicia transcriptome sequences were then mined for simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. About 13% and 3% of the Vicia unigenes contained the putative SSR and SNP sequences, respectively. Among those SSRs, 100 were chosen for the validation and the polymorphism test using the Vicia germplasm set. Thus, our approach takes advantage of the utility of transcriptomic data to expedite a vetch breeding program.

Identification of heat-responsive genes in carnation (Dianthus caryophyllus L.) by RNA-seq

Carnation (Dianthus caryophyllus L.) is an important flower crop, having substantial commercial value as a cut-flower due to the long vase-life and wide array of flower colors and forms. Standard carnation varieties perform well under cool climates but are very susceptible to high temperatures which adversely affect the yield and the quality of the cut-flowers. Despite several studies of carnation contributing to the number of expressed sequence tags (ESTs), transcriptomic information of this species remains very limited, particularly regarding abiotic stress-related genes. Here, transcriptome analysis was performed to generate expression profiles of heat stress (HS)-responsive genes in carnation. We sequenced a cDNA library constructed with mixed RNA from carnation leaves subjected to 42°C HS (0, 0.5, 1, and 2 h) and 46°C HS (0.5, 1, and 2 h), and obtained 45,604,882 high quality paired-end reads. After de novo assembly and quantitative assessment 99,255 contigs were generated with an average length of 1053 bp. We then obtained functional annotations by aligning contigs with public protein databases including NR, SwissProt, KEGG, and COG. Using the above carnation transcriptome as the reference, we compared the effects of high temperature treatments (42°C: duration 0.5, 2, or 12 h) delivered to aseptic carnation seedlings, relative to untreated controls, using the FPKM metric. Overall, 11,471 genes were identified which showed a significant response to one or more of the three HS treatment times. In addition, based on GO and metabolic pathway enrichment analyses, a series of candidate genes involved in thermo-tolerance responses were selected and characterized. This study represents the first expression profiling analysis of D. caryophyllus under heat stress treatments. Numerous genes were found to be induced in response to HS, the study of which may advance our understanding of heat response of carnation.

Identification of genes associated with chlorophyll accumulation in flower petals

Plants have an ability to prevent chlorophyll accumulation, which would mask the bright flower color, in their petals. In contrast, leaves contain substantial amounts of chlorophyll, as it is essential for photosynthesis. The mechanisms of organ-specific chlorophyll accumulation are unknown. To identify factors that determine the chlorophyll content in petals, we compared the expression of genes related to chlorophyll metabolism in different stages of non-green (red and white) petals (very low chlorophyll content), pale-green petals (low chlorophyll content), and leaves (high chlorophyll content) of carnation (Dianthus caryophyllus L.). The expression of many genes encoding chlorophyll biosynthesis enzymes, in particular Mg-chelatase, was lower in non-green petals than in leaves. Non-green petals also showed higher expression of genes involved in chlorophyll degradation, including STAY-GREEN gene and pheophytinase. These data suggest that the absence of chlorophylls in carnation petals may be caused by the low rate of chlorophyll biosynthesis and high rate of degradation. Similar results were obtained by the analysis of Arabidopsis microarray data. In carnation, most genes related to chlorophyll biosynthesis were expressed at similar levels in pale-green petals and leaves, whereas the expression of chlorophyll catabolic genes was higher in pale-green petals than in leaves. Therefore, we hypothesize that the difference in chlorophyll content between non-green and pale-green petals is due to different levels of chlorophyll biosynthesis. Our study provides a basis for future molecular and genetic studies on organ-specific chlorophyll accumulation.

Characterization of flower-bud transcriptome and development of genic SSR markers in Asian lotus (Nelumbo nucifera Gaertn.)

Background

Asian lotus (Nelumbo nucifera Gaertn.) is the national flower of India, Vietnam, and one of the top ten traditional Chinese flowers. Although lotus is highly valued for its ornamental, economic and cultural uses, genomic information, particularly the expressed sequence based (genic) markers is limited. High-throughput transcriptome sequencing provides large amounts of transcriptome data for promoting gene discovery and development of molecular markers.

Results

In this study, 68,593 unigenes were assembled from 1.34 million 454 GS-FLX sequence reads of a mixed flower-bud cDNA pool derived from three accessions of N. nucifera. A total of 5,226 SSR loci were identified, and 3,059 primer pairs were designed for marker development. Di-nucleotide repeat motifs were the most abundant type identified with a frequency of 65.2%, followed by tri- (31.7%), tetra- (2.1%), penta- (0.5%) and hexa-nucleotide repeats (0.5%). A total of 575 primer pairs were synthesized, of which 514 (89.4%) yielded PCR amplification products. In eight Nelumbo accessions, 109 markers were polymorphic. They were used to genotype a sample of 44 accessions representing diverse wild and cultivated genotypes of Nelumbo. The number of alleles per locus varied from 2 to 9 alleles and the polymorphism information content values ranged from 0.6 to 0.9. We performed genetic diversity analysis using 109 polymorphic markers. A UPGMA dendrogram was constructed based on Jaccard’s similarity coefficients revealing distinct clusters among the 44 accessions.

Conclusions

Deep transcriptome sequencing of lotus flower buds developed 3,059 genic SSRs, making a significant addition to the existing SSR markers in lotus. Among them, 109 polymorphic markers were successfully validated in 44 accessions of Nelumbo. This comprehensive set of genic SSR markers developed in our study will facilitate analyses of genetic diversity, construction of linkage maps, gene mapping, and marker-assisted selection breeding for lotus.

De novo transcriptome analysis of petal senescence in Gardenia jasminoides Ellis

BACKGROUND

The petal senescence of ethylene insensitive species has not been investigated thoroughly while little is known about the temporal and tissue specific expression patterns of transcription factors (TFs) in this developmental process. Even less is known on flower senescence of the ornamental pot plant Gardenia jasminoides, a non climacteric flower with significant commercial value.

RESULTS

We initiated a de novo transcriptome study to investigate the petal senescence in four developmental stages of cut gardenia flowers considering that the visible symptoms of senescence appear within 4 days of flower opening. De novo assembly of transcriptome sequencing resulted in 102,263 contigs with mean length of 360 nucleotides that generated 57,503 unigenes. These were further clustered into 20,970 clusters and 36,533 singletons. The comparison of the consecutive developmental stages resulted in 180 common, differentially expressed unigenes. A large number of Simple Sequence Repeats were also identified comprising a large number of dinucleotides and trinucleotides. The prevailing families of differentially expressed TFs comprise the AP2/EREBP, WRKY and the bHLH. There are 81 differentially expressed TFs when the symptoms of flower senescence become visible with the most prevailing being the WRKY family with 19 unigenes. No other WRKY TFs had been identified up to now in petal senescence of ethylene insensitive species. A large number of differentially expressed genes were identified at the initiation of visible symptoms of senescence compared to the open flower stage indicating a significant shift in the expression profiles which might be coordinated by up-regulated and/or down-regulated TFs. The expression of 16 genes that belong to the TF families of WRKY, bHLH and the ethylene sensing pathway was validated using qRT--PCR.

CONCLUSION

This de novo transcriptome analysis resulted in the identification of TFs with specific temporal expression patterns such as two WRKYs and one bHLH, which might play the role of senescence progression regulators. Further research is required to investigate their role in gardenia flowers in order to develop tools to delay petal senescence.

Gleditsia sinensis: transcriptome sequencing, construction, and application of its protein-protein interaction network

Gleditsia sinensis is a genus of deciduous tree in the family Caesalpinioideae, native to China, and is of great economic importance. However, despite its economic value, gene sequence information is strongly lacking. In the present study, transcriptome sequencing of G. sinensis was performed resulting in approximately 75.5 million clean reads assembled into 142155 unique transcripts generating 58583 unigenes. The average length of the unigenes was 900 bp, with an N50 of 549 bp. The obtained unigene sequences were then compared to four protein databases to include NCBI nonredundant protein (NRDB), Swiss-prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Cluster of Orthologous Groups (COG). Using BLAST procedure, 31385 unigenes (53.6%) were generated to have functional annotations. Additionally, sequence homologies between identified unigenes and genes of known species in a protein-protein interaction (PPI) network facilitated G. sinensis PPI network construction. Based on this network construction, new stress resistance genes (including cold, drought, and high salinity) were predicted. The present study is the first investigation of genome-wide gene expression in G. sinensis with the results providing a basis for future functional genomic studies relating to this species.

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.

Polyploidy and the petal transcriptome of Gossypium

Background

Genes duplicated by polyploidy (homoeologs) may be differentially expressed in plant tissues. Recent research using DNA microarrays and RNAseq data have described a cacophony of complex expression patterns during development of cotton fibers, petals, and leaves. Because of its highly canalized development, petal tissue has been used as a model tissue for gene expression in cotton. Recent advances in cotton genome annotation and assembly now permit an enhanced analysis of duplicate gene deployment in petals from allopolyploid cotton.

Results

Homoeologous gene expression levels were quantified in diploid and tetraploid flower petals of Gossypium using the Gossypium raimondii genome sequence as a reference. In the polyploid, most homoeologous genes were expressed at equal levels, though a subset had an expression bias of A_T and D_T copies. The direction of gene expression bias was conserved in natural and recent polyploids of cotton. Conservation of direction of bias and additional comparisons between the diploids and tetraploids suggested different regulation mechanisms of gene expression. We described three phases in the evolution of cotton genomes that contribute to gene expression in the polyploid nucleus.

Conclusions

Compared to previous studies, a surprising level of expression homeostasis was observed in the expression patterns of polyploid genomes. Conserved expression bias in polyploid petals may have resulted from cis-acting modifications that occurred prior to polyploidization. Some duplicated genes were intriguing exceptions to general trends. Mechanisms of gene regulation for these and other genes in the cotton genome warrants further investigation.

Transcriptomic analysis of cut tree peony with glucose supply using the RNA-Seq technique

Several unigenes encoding ACS and ERF involved in ethylene biosynthesis and signal transduction were greatly down-regulated in the petal transcriptome of cut tree peony 'Luoyang Hong' with glucose treatment. Glucose also repressed stress-related transcription factor genes DREB, CBF, NAC, WRKY and bHLH. Tree peony (Paeonia suffruticosa Andrews) is a famous traditional flower in China. Glucose supply prolonging vase life of cut tree peony flowers is associated with its role in the suppression of sensitivity to ethylene and ethylene production, but the regulation mechanism of sugar on ethylene biosynthesis and signaling is unclear. In the present work, a normalized cDNA pool was constructed as the reference transcriptome from mixed petals of different developmental cut tree peony 'Luoyang Hong' and sequenced using the Illumina HiSeq™ 2000 platform. We obtained 33,117 unigenes annotated with public protein databases. In addition, the transcriptome change in petals of cut tree peony with glucose supply and the control treatment was investigated. With non-redundant annotation, 173 differentially expressed genes were identified, with 41 up-regulated genes and 132 down-regulated genes. According to RNA-Seq data and real-time quantitative polymerase chain reaction validation, one unigene encoding ACS, a key ethylene synthetic enzyme, and four unigenes encoding ERF, which is involved in ethylene signal transduction was greatly down-regulated with glucose treatment. Furthermore, stress-related transcription factor genes DREB, CBF, NAC, WRKY and bHLH were also repressed with glucose supply, as well as several other stress-responsive and stress-tolerance genes, indicating that glucose supply probably releases the effects induced by various environmental stress. All the results and analysis are valuable resources for better understanding of the beneficial influence of exogenous sugars on cut tree peony.

Sequence analysis of the genome of carnation (Dianthus caryophyllus L.)

The whole-genome sequence of carnation (Dianthus caryophyllus L.) cv. ‘Francesco’ was determined using a combination of different new-generation multiplex sequencing platforms. The total length of the non-redundant sequences was 568 887 315 bp, consisting of 45 088 scaffolds, which covered 91% of the 622 Mb carnation genome estimated by k-mer analysis. The N50 values of contigs and scaffolds were 16 644 bp and 60 737 bp, respectively, and the longest scaffold was 1 287 144 bp. The average GC content of the contig sequences was 36%. A total of 1050, 13, 92 and 143 genes for tRNAs, rRNAs, snoRNA and miRNA, respectively, were identified in the assembled genomic sequences. For protein-encoding genes, 43 266 complete and partial gene structures excluding those in transposable elements were deduced. Gene coverage was ∼98%, as deduced from the coverage of the core eukaryotic genes. Intensive characterization of the assigned carnation genes and comparison with those of other plant species revealed characteristic features of the carnation genome. The results of this study will serve as a valuable resource for fundamental and applied research of carnation, especially for breeding new carnation varieties. Further information on the genomic sequences is available at http://carnation.kazusa.or.jp.

Construction of a reference genetic linkage map for carnation (Dianthus caryophyllus L.)

Background

Genetic linkage maps are important tools for many genetic applications including mapping of quantitative trait loci (QTLs), identifying DNA markers for fingerprinting, and map-based gene cloning. Carnation (Dianthus caryophyllus L.) is an important ornamental flower worldwide. We previously reported a random amplified polymorphic DNA (RAPD)-based genetic linkage map derived from Dianthus capitatus ssp. andrezejowskianus and a simple sequence repeat (SSR)-based genetic linkage map constructed using data from intraspecific F₂ populations; however, the number of markers was insufficient, and so the number of linkage groups (LGs) did not coincide with the number of chromosomes (x = 15). Therefore, we aimed to produce a high-density genetic map to improve its usefulness for breeding purposes and genetic research.

Results

We improved the SSR-based genetic linkage map using SSR markers derived from a genomic library, expression sequence tags, and RNA-seq data. Linkage analysis revealed that 412 SSR loci (including 234 newly developed SSR loci) could be mapped to 17 linkage groups (LGs) covering 969.6 cM. Comparison of five minor LGs covering less than 50 cM with LGs in our previous RAPD-based genetic map suggested that four LGs could be integrated into two LGs by anchoring common SSR loci. Consequently, the number of LGs corresponded to the number of chromosomes (x = 15). We added 192 new SSRs, eight RAPD, and two sequence-tagged site loci to refine the RAPD-based genetic linkage map, which comprised 15 LGs consisting of 348 loci covering 978.3 cM. The two maps had 125 SSR loci in common, and most of the positions of markers were conserved between them. We identified 635 loci in carnation using the two linkage maps. We also mapped QTLs for two traits (bacterial wilt resistance and anthocyanin pigmentation in the flower) and a phenotypic locus for flower-type by analyzing previously reported genotype and phenotype data.

Conclusions

The improved genetic linkage maps and SSR markers developed in this study will serve as reference genetic linkage maps for members of the genus Dianthus, including carnation, and will be useful for mapping QTLs associated with various traits, and for improving carnation breeding programs.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.

De novo assembly and characterization of the transcriptome, and development of SSR markers in wax gourd (Benicasa hispida)

BACKGROUND

Wax gourd is a widely used vegetable of Cucuribtaceae, and also has important medicinal and health values. However, the genomic resources of wax gourd were scarcity, and only a few nucleotide sequences could be obtained in public databases.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we examined transcriptome in wax gourd. More than 44 million of high quality reads were generated from five different tissues of wax gourd using Illumina paired-end sequencing technology. Approximately 4 Gbp data were generated, and de novo assembled into 65,059 unigenes, with an N50 of 1,132 bp. Based on sequence similarity search with known protein database, 36,070 (55.4%) showed significant similarity to known proteins in Nr database, and 24,969 (38.4%) had BLAST hits in Swiss-Prot database. Among the annotated unigenes, 14,994 of wax gourd unigenes were assigned to GO term annotation, and 23,977 were found to have COG classifications. In addition, a total of 18,713 unigenes were assigned to 281 KEGG pathways. Furthermore, 6,242 microsatellites (simple sequence repeats) were detected as potential molecular markers in wax gourd. Two hundred primer pairs for SSRs were designed for validation of the amplification and polymorphism. The result showed that 170 of the 200 primer pairs were successfully amplified and 49 (28.8%) of them exhibited polymorphisms.

CONCLUSION/SIGNIFICANCE

Our study enriches the genomic resources of wax gourd and provides powerful information for future studies. The availability of this ample amount of information about the transcriptome and SSRs in wax gourd could serve as valuable basis for studies on the physiology, biochemistry, molecular genetics and molecular breeding of this important vegetable crop.

Development of 65 novel polymorphic cDNA-SSR markers in common vetch (Vicia sativa subsp. sativa) using next generation sequencing

Vetch (Vicia sativa L.) is one of the most important annual forage legumes in the World due to its multiple uses (i.e., hay, grain, silage and green manure) and high nutritional value. However, detrimental cyanoalanine toxins in its plant parts including seeds and its vulnerability to hard winter conditions are currently reducing the agronomic values of vetch varieties. Moreover, the existence in the public domain of very few genomic resources, especially molecular markers, has further hampered breeding efforts. Polymorphic simple sequence repeat markers from transcript sequences (cDNA; simple sequence repeat [SSR]) were developed for Vicia sativa subsp. sativa. We found 3,811 SSR loci from 31,504 individual sequence reads, and 300 primer pairs were designed and synthesized. In total, 65 primer pairs were found to be consistently scorable when 32 accessions were tested. The numbers of alleles ranged from 2 to 19, frequency of major alleles per locus were 0.27-0.87, the genotype number was 2-19, the overall polymorphism information content (PIC) values were 0.20-0.86, and the observed and expected heterozygosity values were 0.00-0.41 and 0.264-0.852, respectively. These markers provide a useful tool for assessing genetic diversity, population structure, and positional cloning, facilitating vetch breeding programs.