Transcriptome Analysis of Syringa oblata Lindl. Inflorescence Identifies Genes Associated with Pigment Biosynthesis and Scent Metabolism

Volatile secondary metabolome and transcriptome analysis reveals distinct regulation mechanism of aroma biosynthesis in Syringa oblata and S. vulgaris

Lilacs have high ornamental value due to their strong aroma. However, the molecular regulatory mechanisms of aroma biosynthesis and metabolism in lilac were largely unclear. In this study, two varieties with distinct aroma, Syringa oblata 'Zi Kui' (faint aroma) and Syringa vulgaris 'Li Fei' (strong aroma), were used for exploring the regulation mechanism of aroma difference. Via GC-MS analysis, a total of 43 volatile components were identified. Terpene volatiles was the most abundant volatiles constituting the aroma of two varieties. Notably, 3 volatile secondary metabolites were unique in 'Zi Kui' and 30 volatile secondary metabolites were unique in 'Li Fei'. Then, a transcriptome analysis was performed to clarify the regulation mechanism of aroma metabolism difference between these two varieties, and identified 6411 differentially expressed genes (DEGs). Interestingly, ubiquinone and other terpenoid-quinone biosynthesis genes were significantly enriched in DEGs. We further conducted a correlation analysis between the volatile metabolome and transcriptome and found that TPS, GGPPS, and HMGS genes might be the key contributors to the differences in floral fragrance composition between the two lilac varieties. Our study improves the understanding in the regulation mechanism of Lilac aroma and would help improve the aroma of ornamental crops by metabolic engineering.

Syringa oblata genome provides new insights into molecular mechanism of flower color differences among individuals and biosynthesis of its flower volatiles

Syringa oblata is a high ornamental value tree owing to its elegant colors, unique aromas and wide adaptability, however, studies on the molecular mechanism underlying the formation of its ornamental traits are still lacking. Here, we presented a chromosome-scale genome assembly of S. oblata and the final genome size was 1.11 Gb with a contig N50 of 4.75 Mb, anchored on 23 chromosomes and was a better reference for S. oblata transcriptome assembly. Further by integrating transcriptomic and metabolic data, it was concluded that F3H, F3'H, 4CL and PAL, especially the F3'H, were important candidates involved in the formation of floral color differences among S. oblata individuals. Genome-wide identification and analysis revealed that the TPS-b subfamily was the most abundant subfamily of TPS family in S. oblata, which together with the CYP76 family genes determined the formation of the major floral volatiles of S. oblata. Overall, our results provide an important reference for mechanistic studies on the main ornamental traits and molecular breeding in S. oblata.

Lilac (Syringa oblata) genome provides insights into its evolution and molecular mechanism of petal color change

Color change during flower opening is common; however, little is understood on the biochemical and molecular basis related. Lilac (Syringa oblata), a well-known woody ornamental plant with obvious petal color changes, is an ideal model. Here, we presented chromosome-scale genome assembly for lilac, resolved the flavonoids metabolism, and identified key genes and potential regulatory networks related to petal color change. The genome assembly is 1.05 Gb anchored onto 23 chromosomes, with a BUSCO score of 96.6%. Whole-genome duplication (WGD) event shared within Oleaceae was revealed. Metabolome quantification identified delphinidin-3-O-rutinoside (Dp3Ru) and cyanidin-3-O-rutinoside (Cy3Ru) as the major pigments; gene co-expression networks indicated WRKY an essential regulation factor at the early flowering stage, ERF more important in the color transition period (from violet to light nearly white), while the MBW complex participated in the entire process. Our results provide a foundation for functional study and molecular breeding in lilac.

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.

Integration of Transcriptome and Methylome Analyses Provides Insight Into the Pathway of Floral Scent Biosynthesis in Prunus mume

DNA methylation is a common epigenetic modification involved in regulating many biological processes. However, the epigenetic mechanisms involved in the formation of floral scent have rarely been reported within a famous traditional ornamental plant Prunus mume emitting pleasant fragrance in China. By combining whole-genome bisulfite sequencing and RNA-seq, we determined the global change in DNA methylation and expression levels of genes involved in the biosynthesis of floral scent in four different flowering stages of P. mume. During flowering, the methylation status in the “CHH” sequence context (with H representing A, T, or C) in the promoter regions of genes showed the most significant change. Enrichment analysis showed that the differentially methylated genes (DMGs) were widely involved in eight pathways known to be related to floral scent biosynthesis. As the key biosynthesis pathway of the dominant volatile fragrance of P. mume, the phenylpropane biosynthesis pathway contained the most differentially expressed genes (DEGs) and DMGs. We detected 97 DMGs participated in the most biosynthetic steps of the phenylpropane biosynthesis pathway. Furthermore, among the previously identified genes encoding key enzymes in the biosynthesis of the floral scent of P. mume, 47 candidate genes showed an expression pattern matching the release of floral fragrances and 22 of them were differentially methylated during flowering. Some of these DMGs may or have already been proven to play an important role in biosynthesis of the key floral scent components of P. mume, such as PmCFAT1a/1c, PmBEAT36/37, PmPAL2, PmPAAS3, PmBAR8/9/10, and PmCNL1/3/5/6/14/17/20. In conclusion, our results for the first time revealed that DNA methylation is widely involved in the biosynthesis of floral scent and may play critical roles in regulating the floral scent biosynthesis of P. mume. This study provided insights into floral scent metabolism for molecular breeding.

Transcriptome sequencing reveals terpene biosynthesis pathway genes accounting for volatile terpene of tree peony

Transcriptomic and volatile component analyses showed that high expression levels of genes from the terpenoid backbone biosynthesis pathway and the monoterpene metabolic pathway can strengthen the floral fragrance of tree peony. Floral fragrance is a crucial ornamental trait whose improvement is one of the primary objectives of tree peony breeding. So far, exploration of the floral fragrance of tree peony has focused on the identification of its volatile components, but the molecular mechanisms responsible for their formation remain unclear. Here, we identified 128 volatile components from the petals of tree peony and found that they consisted primarily of terpenes, alcohols, and esters. Based on the distribution pattern of these major fragrance components, 24 tree peony cultivars were classified into 4 types: grassy scent (ocimene), woody scent (longifolene), lily of the valley scent (linalool), and fruity scent (2-ethyl hexanol). We used RNA-seq to explore the mechanistic basis of terpenoid metabolism in tree peony petals with various scents. The expression levels of AACT, HMGR, PMK, DXS, DXR, HDS, HDR, and GGPS, which encode key enzymes of terpenoid backbone biosynthesis, were upregulated in 'Huangguan' (strong fragrance) compared to 'Fengdan' (faint fragrance). Moreover, the transcript abundance of LIS and MYS, two monoterpene synthase genes, was also enhanced in petals of 'Huangguan' compared to those of 'Fengdan'. Together, these results demonstrate that differences in the expression of genes from the monoterpene synthesis and terpenoid backbone pathways are associated with differences in the fragrance of tree peony. This research provides crucial genetic resources for fragrance improvement and also lays a foundation for further clarification of the mechanisms that underlie tree peony fragrance.

Identification of suitable reference genes for studies of Syringa pinnatifolia Hemsl

Syringa pinnatifolia Hemsl. (Oleaceae) is a species of shrub with a limited distribution in China. Several chemical compounds with pharmacological effects have been isolated from S. pinnatifolia, including new lignans and sesquiterpenes. Studies of gene expression in this species require the identification of suitable reference genes that are stably expressed under different conditions and in different tissues. To identify candidate reference genes, here we used the geNorm, NormFinder, and BestKeeper algorithms to analyze the stability of 12 candidate genes. The geometric mean of the rankings generated with these algorithms was used to obtain a comprehensive ranking. TBP and PP2A were found to be appropriate reference genes for calli treated with different external stimuli, and TIP41 and TBP were found to be appropriate reference genes in differentiated tissues. When calli and differentiated tissues were considered together, TBP and TIP41 were found to be the most reliable reference genes. The selected genes were validated by analysis of HMGR expression in calli and differentiated tissues. This study is the first to screen candidate reference genes in the genus Syringa and could help guide future molecular studies in this genus.

The spatio-temporal biosynthesis of floral flavonols is controlled by differential phylogenetic MYB regulators in Freesia hybrida

Floral flavonols play specific pivotal roles in pollinator attraction, pollen germination and fertility, in addition to other functions in vegetative organs. For many plants, the process of flavonol biosynthesis in late flower development stages and in mature flower tissues is poorly understood, in contrast to early flower development stages. It is thought that this process may be regulated independently of subgroup 7 R2R3 MYB (SG7 MYB) transcription factors. In this study, two FLS genes were shown to be expressed synchronously with the flower development-specific and tissue-specific biosynthesis of flavonols in Freesia hybrida. FhFLS1 contributed to flavonol biosynthesis in early flower buds, toruses and calyxes, and was regulated by four well-known SG7 MYB proteins, designated as FhMYBFs, with at least partial regulatory redundancy. FhFLS2 accounted for flavonols in late developed flowers and in the petals, stamens and pistils, and was targeted directly by non SG7 MYB protein FhMYB21L2. In parallel, AtMYB21 and AtMYB24 also activated AtFLS1, a gene highly expressed in Arabidopsis anthers and pollen, indicating the conserved regulatory roles of MYB21 against FLS genes in these two evolutionarily divergent angiosperm plants. Our results reveal a novel regulatory and synthetic mechanism underlying flavonol biosynthesis in floral organs and tissues which may be exploited to investigate supplementary roles of flavonols in flowers.

Analysis of Relative Scent Intensity, Volatile Compounds and Gene Expression in Freesia "Shiny Gold"

Scent is one of the most important economic traits in Freesia hybrida. "Shiny Gold", a popular cultivar in South Korea, is widely cultivated for its scent. The relative scent intensity of "Shiny Gold" was approximately 16% higher in full-bloomed flower when compared to the yellow bud stage, while tissue-specifically, tepals showed higher intensity in electronic-nose (e-nose) analysis. E-nose analysis also showed that the scent intensity of "Shiny Gold" was higher and lower than "10C3-424" and "10C3-894", respectively, and was similar to "Yvonne". These results correlated to those of the olfactory tests. In total, 19 volatile compounds, including linalool, β-ocimene, D-limonene, trans-β-ionone were detected in gas chromatography-mass spectrometry analysis. Among these, linalool was the major volatile compound, accounting for 38.7% in "Shiny Gold". Linalool synthase and TPS gene expression corresponded to the scent intensity of the four cultivars, with the lowest expression in the "10C3-424". TPS 2, TPS 3, TPS 5, TPS 6 and TPS 8 were highly expressed in both bud and flower in "Shiny Gold", while the expression of TPS 4 was lower, relative to other TPS genes in both the flowering stages. These results may aid in enhancing scent composition in Freesia cultivars using marker-assisted selection.

Comparative transcriptome analysis to identify candidate genes involved in 2-methoxy-1,4-naphthoquinone (MNQ) biosynthesis in Impatiens balsamina L

Impatiens balsamina L. is a tropical ornamental and traditional medicinal herb rich in natural compounds, especially 2-methoxy-1,4-naphthoquinone (MNQ) which is a bioactive compound with tested anticancer activities. Characterization of key genes involved in the shikimate and 1,4-dihydroxy-2-naphthoate (DHNA) pathways responsible for MNQ biosynthesis and their expression profiles in I. balsamina will facilitate adoption of genetic/metabolic engineering or synthetic biology approaches to further increase production for pre-commercialization. In this study, HPLC analysis showed that MNQ was present in significantly higher quantities in the capsule pericarps throughout three developmental stages (early-, mature- and postbreaker stages) whilst its immediate precursor, 2-hydroxy-1,4-naphthoquinone (lawsone) was mainly detected in mature leaves. Transcriptomes of I. balsamina derived from leaf, flower, and three capsule developmental stages were generated, totalling 59.643 Gb of raw reads that were assembled into 94,659 unigenes (595,828 transcripts). A total of 73.96% of unigenes were functionally annotated against seven public databases and 50,786 differentially expressed genes (DEGs) were identified. Expression profiles of 20 selected genes from four major secondary metabolism pathways were studied and validated using qRT-PCR method. Majority of the DHNA pathway genes were found to be significantly upregulated in early stage capsule compared to flower and leaf, suggesting tissue-specific synthesis of MNQ. Correlation analysis identified 11 candidate unigenes related to three enzymes (NADH-quinone oxidoreductase, UDP-glycosyltransferases and S-adenosylmethionine-dependent O-methyltransferase) important in the final steps of MNQ biosynthesis based on genes expression profiles consistent with MNQ content. This study provides the first molecular insight into the dynamics of MNQ biosynthesis and accumulation across different tissues of I. balsamina and serves as a valuable resource to facilitate further manipulation to increase production of MNQ.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

MYB21 interacts with MYC2 to control the expression of terpene synthase genes in flowers of Freesia hybrida and Arabidopsis thaliana

Previously, linalool was found to be the most abundant component among the cocktail of volatiles released from flowers of Freesia hybrida. Linalool formation is catalysed by monoterpene synthase TPS1. However, the regulatory network developmentally modulating the expression of the TPS1 gene in Freesia hybrida remains unexplored. In this study, three regulatory genes, FhMYB21L1, FhMYB21L2, and FhMYC2, were screened from 52 candidates. Two MYB transcription factor genes were synchronously expressed with FhTPS1 and could activate its expression significantly when overexpressed, and the binding of FhMYB21L2 to the MYBCORE sites in the FhTPS1 promoter was further confirmed, indicating a direct role in activation. FhMYC2 showed an inverse expression pattern compared with FhTPS1; its expression led to a decreased binding of FhMYB21 to the FhTPS1 promoter to reduce its activation capacity when co-expressed, suggesting a role for an MYB-bHLH complex in the regulation of the FhTPS1 gene. In Arabidopsis, both MYB21 and MYC2 regulators were shown to activate the expression of sesquiterpene synthase genes, and the regulatory roles of AtMYB21 and AtMYC2 in the expression of the linalool synthase gene were also confirmed, implying conserved functions of the MYB-bHLH complex in these two evolutionarily divergent plants. Moreover, the expression ratio between MYB21 and MYC2 orthologues might be a determinant factor in floral linalool emission.

Volatile Organic Compounds from Orchids: From Synthesis and Function to Gene Regulation

Orchids are one of the most significant plants that have ecologically adapted to every habitat on earth. Orchids show a high level of variation in their floral morphologies, which makes them popular as ornamental plants in the global market. Floral scent and color are key traits for many floricultural crops. Volatile organic compounds (VOCs) play vital roles in pollinator attraction, defense, and interaction with the environment. Recent progress in omics technology has led to the isolation of genes encoding candidate enzymes responsible for the biosynthesis and regulatory circuits of plant VOCs. Uncovering the biosynthetic pathways and regulatory mechanisms underlying the production of floral scents is necessary not only for a better understanding of the function of relevant genes but also for the generation of new cultivars with desirable traits through molecular breeding approaches. However, little is known about the pathways responsible for floral scents in orchids because of their long life cycle as well as the complex and large genome; only partial terpenoid pathways have been reported in orchids. Here, we review the biosynthesis and regulation of floral volatile compounds in orchids. In particular, we focused on the genes responsible for volatile compounds in various tissues and developmental stages in Cymbidium orchids. We also described the emission of orchid floral volatiles and their function in pollination ecology. Taken together, this review will provide a broad scope for the study of orchid floral scents.

Development of SSR markers in Paeonia based on De Novo transcriptomic assemblies

Peony is a famous ornamental and medicinal plant in China, and peony hybrid breeding is an important means of germplasm innovation. However, research on the genome of this species is limited, thereby hindering the genetic and breeding research on peony. In the present study, simple sequence repeat (SSR) locus analysis was performed on expressed sequence tags obtained by the transcriptome sequencing of Paeonia using Microsatellite software. Primers with polymorphism were obtained via polymerase chain reaction amplification and electrophoresis. As a result, a total of 86,195 unigenes were obtained by assembling the transcriptome data of Paeonia. Functional annotations were obtained in seven functional databases including 49,172 (Non-Redundant Protein Sequence Database: 57.05%), 38,352 (Nucleotide Sequence Database: 44.49%), 36,477 (Swiss Prot: 42.32%), 38,905 (Clusters of Orthologous Groups for Eukaryotic Complete Genomes: 45.14%), 37,993 (Kyoto Encyclopedia of Genes and Genomes: 44.08%), 26,832 (Gene Ontology: 31.13%) and 37,758 (Pfam: 43.81%) unigenes. Meanwhile, 21,998 SSR loci were distributed in 17,567 unigenes containing SSR sequences, and the SSR distribution frequency was 25.52%, with an average of one SSR sequence per 4.66 kb. Mononucleotide, dinucleotide, and trinucleotide were the main repeat types, accounting for 55.74%, 25.58%, and 13.21% of the total repeat times, respectively. Forty-five pairs of the 100 pairs of primers selected randomly could amplify clear polymorphic bands. The polymorphic primers of these 45 pairs were used to cluster and analyze 16 species of peony. The new SSR molecular markers can be useful for the study of genetic diversity and marker-assisted breeding of peony.

Transcriptomic analysis reveals flavonoid biosynthesis of Syringa oblata Lindl. in response to different light intensity

BACKGROUND

Hazy weather significantly increase air pollution and affect light intensity which may also affect medicinal plants growth. Syringa oblata Lindl. (S. oblata), an effective anti-biofilm medicinal plants, is also vulnerable to changes in plant photoperiods and other abiotic stress responses. Rutin, one of the flavonoids, is the main bioactive ingredient in S. oblata that inhibits Streptococcus suis biofilm formation. Thus, the present study aims to explore the biosynthesis and molecular basis of flavonoids in S. oblata in response to different light intensity.

RESULTS

In this study, it was shown that compared with natural (Z0) and 25% ~ 35% (Z2) light intensities, the rutin content of S. oblata under 50% ~ 60% (Z1) light intensity increased significantly. In addition, an integrated analysis of metabolome and transcriptome was performed using light intensity stress conditions from two kinds of light intensities which S. oblata was subjected to: Z0 and Z1. The results revealed that differential metabolites and genes were mainly related to the flavonoid biosynthetic pathway. We found out that 13 putative structural genes and a transcription factor bHLH were significantly up-regulated in Z1. Among them, integration analysis showed that 3 putative structural genes including 4CL1, CYP73A and CYP75B1 significantly up-regulated the rutin biosynthesis, suggesting that these putative genes may be involved in regulating the flavonoid biosynthetic pathway, thereby making them key target genes in the whole metabolic process.

CONCLUSIONS

The present study provided helpful information to search for the novel putative genes that are potential targets for S. oblata in response to light intensity.

RNA sequencing analysis of Cymbidium goeringii identifies floral scent biosynthesis related genes

BACKGROUND

Cymbidium goeringii belongs to the Orchidaceae, which is one of the most abundant angiosperm families. Cymbidium goeringii consist with high economic value and characteristics include fragrance and multiple flower colors. Floral scent is one of the important strategies for ensuring fertilization. However, limited genetic data is available in this non-model plant, and little known about the molecular mechanism responsible for floral scent in this orchid. Transcriptome and expression profiling data are needed to identify genes and better understand the biological mechanisms of floral scents in this species. Present transcriptomic data provides basic information on the genes and enzymes related to and pathways involved in flower secondary metabolism in this plant.

RESULTS

In this study, RNA sequencing analyses were performed to identify changes in gene expression and biological pathways related scent metabolism. Three cDNA libraries were obtained from three developmental floral stages: closed bud, half flowering stage and full flowering stage. Using Illumina technique 159,616,374 clean reads were obtained and were assembled into 85,868 final unigenes (average length 1194 nt), 33.85% of which were annotated in the NCBI non redundant protein database. Among this unigenes 36,082 were assigned to gene ontology and 23,164 were combined with COG groups. Total 33,417 unigenes were assigned in 127 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. According these transcriptomic data we identified number of candidates genes which differentially expressed in different developmental stages of flower related to fragrance biosynthesis. In q-RT-PCR most of the fragrance related genes highly expressed in half flowering stage.

CONCLUSIONS

RNA-seq and DEG data provided comprehensive gene expression information at the transcriptional level that could be facilitate the molecular mechanisms of floral biosynthesis pathways in three developmental phase's flowers in Cymbidium goeringii, moreover providing useful information for further analysis on C. goeringii, and other plants of genus Cymbidium.

CYP79D73 Participates in Biosynthesis of Floral Scent Compound 2-Phenylethanol in Plumeria rubra

Plumeria (Plumeria rubra), well known for its brightly colored and fragrant flowers, emits a number of floral volatile organic compounds (VOCs). Plumeria flowers emit a total of 43 VOCs including nine phenylpropanoids/benzenoids, such as 2-phenylethanol (2PE), benzaldehyde, 2-phenylacetaldehyde (PAld), (E/Z)-phenylacetaldoxime (PAOx), benzyl nitrile (BN), and 2-phenylnitroethane (PN). To identify genes and pathways involved in the production of the major compound 2PE, we analyzed the plumeria floral transcriptome and found a highly expressed, flower-specific gene encoding a cytochrome P450 family 79D protein (PrCYP79D73), which catalyzed the formation of (E/Z)-PAOx. Feeding experiments with deuterated phenylalanine or deuterated (E/Z)-PAOx showed that (E/Z)-PAOx is an intermediate in the biosynthesis of 2PE, as are two nitrogen-containing volatiles, BN and PN, in plumeria flowers. Crude enzyme extracts from plumeria flowers converted l-phenylalanine to (E/Z)-PAOx, PAld, 2PE, BN, and PN. The biosynthesis of these compounds increased with addition of PrCYP79D73-enriched microsomes but was blocked by pretreatment with 4-phenylimidazole, an inhibitor of cytochrome P450 enzymes. Moreover, overexpression of PrCYP79D73 in Nicotiana benthamiana resulted in the emission of (E/Z)-PAOx as well as PAld, 2PE, BN, and PN, all of which were also found among plumeria floral VOCs. Taken together, our results demonstrate that PrCYP79D73 is a crucial player in the biosynthesis of the major floral VOC 2PE and other nitrogen-containing volatiles. These volatiles may be required for plant defense as well as to attract pollinators for the successful reproduction of plumeria.

De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using Illumina platform

Cornus officinalis is one of the most widely used medicinal plants in China and other East Asian countries to cure diseases such as liver, kidney, cardiovascular diseases and frequent urination for thousands of years. It is a Level 3 protected species, and is one of the 42 national key protected wild species of animals and plants in China. However, the genetics and molecular biology of C. officinalis are poorly understood, which has hindered research on the molecular mechanism of its metabolism and utilization. Hence, enriching its genomic data and information is very important. In recent years, the fast-growing technology of next generation sequencing has provided an effective path to gain genomic information from nonmodel species. This study is the first to explore the leaf and fruit tissue transcriptome of C. officinalis using the Illumina HiSeq 4000 platform. A total of 57,954,134 and 60,971,652 clean reads from leaf and fruit were acquired, respectively (GenBank number SRP115440). The pooled reads from all two libraries were assembled into 56,392 unigenes with an average length 856 bp. Among these, 41,146 unigenes matched with sequences in the NCBI nonredundant protein database. The Gene Ontology database assigned 24,336 unigenes with biological process (83.26%), cellular components (53.58%), and molecular function (83.93%). In addition, 10,808 unigenes were assigned a KOG functional classification by the KOG database. Searching against the KEGG pathway database indicated that 18,435 unigenes were mapped to 371 KEGG pathways. Moreover, the edgeR database identified 4,585 significant differentially expressed genes (DEGs), of which 1,392 were up-regulated and 3,193 were down-regulated in fruit tissue compared with leaf tissue. Finally, we explored 581 transcription factors with 50 transcription factor gene families. Most DEGs and transcription factors were related to terpene biosynthesis and secondary metabolic regulation. This study not only represented the first de novo transcriptomic analysis of C. officinalis but also provided fundamental information on its genes and biosynthetic pathway. These findings will help us explore the molecular metabolism mechanism of terpene biosynthesis in C. officinalis.

Molecular characterization and functional analysis of chalcone synthase from Syringa oblata Lindl. in the flavonoid biosynthetic pathway

The flower color of Syringa oblata Lindl., which is often modulated by the flavonoid content, varies and is an important ornamental feature. Chalcone synthase (CHS) catalyzes the first key step in the flavonoid biosynthetic pathway. However, little is known about the role of S. oblata CHS (SoCHS) in flavonoid biosynthesis in this species. Here, we isolate and analyze the cDNA (SoCHS1) that encodes CHS in S. oblata. We also sought to analyzed the molecular characteristics and function of flavonoid metabolism by SoCHS1. We successfully isolated the CHS-encoding genomic DNA (gDNA) in S. oblata (SoCHS1), and the gene structural analysis indicated it had no intron. The opening reading frame (ORF) sequence of SoCHS1 was 1170bp long and encoded a 389-amino acid polypeptide. Multiple sequence alignment revealed that both the conserved CHS active site residues and CHS signature sequence were in the deduced amino acid sequence of SoCHS1. Crystallographic analysis revealed that the protein structure of SoCHS1 is highly similar to that of FnCHS1 in Freesia hybrida. The quantitative real-time polymerase chain reaction (PCR) performed to detect the SoCHS1 transcript expression levels in flowers, and other tissues revealed the expression was significantly correlated with anthocyanin accumulation during flower development. The ectopic expression results of Nicotiana tabacum showed that SoCHS1 overexpression in transgenic tobacco changed the flower color from pale pink to pink. In conclusion, these results suggest that SoCHS1 plays an essential role in flavonoid biosynthesis in S. oblata, and could be used to modify flavonoid components in other plant species.

De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl

Sorbus pohuashanensis is a native tree species of northern China that is used for a variety of ecological purposes. The species is often grown as an ornamental landscape tree because of its beautiful form, silver flowers in early summer, attractive pinnate leaves in summer, and red leaves and fruits in autumn. However, development and further utilization of the species are hindered by the lack of comprehensive genetic information, which impedes research into its genetics and molecular biology. Recent advances in de novo transcriptome sequencing (RNA-seq) technology have provided an effective means to obtain genomic information from non-model species. Here, we applied RNA-seq for sequencing S. pohuashanensis leaves and obtained a total of 137,506 clean reads. After assembly, 96,213 unigenes with an average length of 770 bp were obtained. We found that 64.5% of the unigenes could be annotated using bioinformatics tools to analyze gene function and alignment with the NCBI database. Overall, 59,089 unigenes were annotated using the Nr database(non-redundant protein database), 35,225 unigenes were annotated using the GO (Gene Ontology categories) database, and 33,168 unigenes were annotated using COG (Cluster of Orthologous Groups). Analysis of the unigenes using the KEGG (Kyoto Encyclopedia of Genes and Genomes) database indicated that 13,953 unigenes were involved in 322 metabolic pathways. Finally, simple sequence repeat (SSR) site detection identified 6,604 unigenes that included EST-SSRs and a total of 7,473 EST-SSRs in the unigene sequences. Fifteen polymorphic SSRs were screened and found to be of use for future genetic research. These unigene sequences will provide important genetic resources for genetic improvement and investigation of biochemical processes in S. pohuashanensis.

Transcriptome Sequencing Analysis Reveals a Difference in Monoterpene Biosynthesis between Scented Lilium 'Siberia' and Unscented Lilium 'Novano'

Lilium is a world famous fragrant bulb flower with high ornamental and economic values, and significant differences in fragrance are found among different Lilium genotypes. In order to explore the mechanism underlying the different fragrances, the floral scents of Lilium 'Sibeia', with a strong fragrance, and Lilium 'Novano', with a very faint fragrance, were collected in vivo using a dynamic headspace technique. These scents were identified using automated thermal desorption-gas chromatography/mass spectrometry (ATD-GC/MS) at different flowering stages. We used RNA-Seq technique to determine the petal transcriptome at the full-bloom stage and analyzed differentially expressed genes (DEGs) to investigate the molecular mechanism of floral scent biosynthesis. The results showed that a significantly higher amount of Lilium 'Siberia' floral scent was released compared with Lilium 'Novano'. Moreover, monoterpenes played a dominant role in the floral scent of Lilium 'Siberia'; therefore, it is believed that the different emissions of monoterpenes mainly contributed to the difference in the floral scent between the two Lilium genotypes. Transcriptome sequencing analysis indicated that ~29.24 Gb of raw data were generated and assembled into 124,233 unigenes, of which 35,749 unigenes were annotated. Through a comparison of gene expression between these two Lilium genotypes, 6,496 DEGs were identified. The genes in the terpenoid backbone biosynthesis pathway showed significantly different expression levels. The gene expressions of 1-deoxy-D-xylulose 5-phosphate synthase (DXS), 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR), 4-hydroxy-3-methylbut-2-enyl diphosphate synthase (HDS), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (HDR), isopentenyl diphosphate isomerase (IDI), and geranyl diphosphate synthase (GPS/GGPS), were upregulated in Lilium 'Siberia' compared to Lilium 'Novano', and two monoterpene synthase genes, ocimene synthase gene (OCS) and myrcene synthase gene (MYS), were also expressed at higher levels in the tepals of Lilium 'Siberia', which was consistent with the monoterpene release amounts. We demonstrated that the high activation levels of the pathways contributed to monoterpene biosynthesis in Lilium 'Siberia' resulting in high accumulations and emissions of monoterpenes, which led to the difference in fragrance between these two Lilium genotypes.

Transcriptomic Analysis Reveals Mechanisms of Sterile and Fertile Flower Differentiation and Development in Viburnum macrocephalum f. keteleeri

Sterile and fertile flowers are an important evolutionary developmental (evo-devo) phenotype in angiosperm flowers, playing important roles in pollinator attraction and sexual reproductive success. However, the gene regulatory mechanisms underlying fertile and sterile flower differentiation and development remain largely unknown. Viburnum macrocephalum f. keteleeri, which possesses fertile and sterile flowers in a single inflorescence, is a useful candidate species for investigating the regulatory networks in differentiation and development. We developed a de novo-assembled flower reference transcriptome. Using RNA sequencing (RNA-seq), we compared the expression patterns of fertile and sterile flowers isolated from the same inflorescence over its rapid developmental stages. The flower reference transcriptome consisted of 105,683 non-redundant transcripts, of which 5,675 transcripts showed significant differential expression between fertile and sterile flowers. Combined with morphological and cytological changes between fertile and sterile flowers, we identified expression changes of many genes potentially involved in reproductive processes, phytohormone signaling, and cell proliferation and expansion using RNA-seq and qRT-PCR. In particular, many transcription factors (TFs), including MADS-box family members and ABCDE-class genes, were identified, and expression changes in TFs involved in multiple functions were analyzed and highlighted to determine their roles in regulating fertile and sterile flower differentiation and development. Our large-scale transcriptional analysis of fertile and sterile flowers revealed the dynamics of transcriptional networks and potentially key components in regulating differentiation and development of fertile and sterile flowers in Viburnum macrocephalum f. keteleeri. Our data provide a useful resource for Viburnum transcriptional research and offer insights into gene regulation of differentiation of diverse evo-devo processes in flowers.

Comparative Transcriptome Reveals Benzenoid Biosynthesis Regulation as Inducer of Floral Scent in the Woody Plant Prunus mume

Mei (Prunus mume) is a peculiar woody ornamental plant famous for its inviting fragrance in winter. However, in this valuable plant, the mechanism behind floral volatile development remains poorly defined. Therefore, to explore the floral scent formation, a comparative transcriptome was conducted in order to identify the global transcripts specifying flower buds and blooming flowers of P. mume. Differentially expressed genes were identified between the two different stages showing great discrepancy in floral volatile production. Moreover, according to the expression specificity among the organs (stem, root, fruit, leaf), we summarized one gene cluster regulating the benzenoid floral scent. Significant gene changes were observed in accordance with the formation of benzenoid, thus pointing the pivotal roles of genes as well as cytochrome-P450s and short chain dehydrogenases in the benzenoid biosynthetic process. Further, transcription factors like EMISSION OF BENZENOID I and ODORANT I performed the same expression pattern suggesting key roles in the management of the downstream genes. Taken together, these data provide potential novel anchors for the benzenoid pathway, and the insight for the floral scent induction and regulation mechanism in woody plants.

Analysis of transcriptome in hickory (Carya cathayensis), and uncover the dynamics in the hormonal signaling pathway during graft process

BACKGROUND

Hickory (Carya cathayensis), a woody plant with high nutritional and economic value, is widely planted in China. Due to its long juvenile phase, grafting is a useful technique for large-scale cultivation of hickory. To reveal the molecular mechanism during the graft process, we sequenced the transcriptomes of graft union in hickory.

RESULTS

In our study, six RNA-seq libraries yielded a total of 83,676,860 clean short reads comprising 4.19 Gb of sequence data. A large number of differentially expressed genes (DEGs) at three time points during the graft process were identified. In detail, 777 DEGs in the 7 d vs 0 d (day after grafting) comparison were classified into 11 enriched Gene Ontology (GO) categories, and 262 DEGs in the 14 d vs 0 d comparison were classified into 15 enriched GO categories. Furthermore, an overview of the PPI network was constructed by these DEGs. In addition, 20 genes related to the auxin-and cytokinin-signaling pathways were identified, and some were validated by qRT-PCR analysis.

CONCLUSIONS

Our comprehensive analysis provides basic information on the candidate genes and hormone signaling pathways involved in the graft process in hickory and other woody plants.