Background: Siraitia grosvenorii (Luohanguo) is an herbaceous perennial plant native to southern China and most prevalent in Guilin city. Its fruit contains a sweet, fleshy, edible pulp that is widely used in traditional Chinese medicine. The major bioactive constituents in the fruit extract are the cucurbitane-type triterpene saponins known as mogrosides. Among them, mogroside V is nearly 300 times sweeter than sucrose. However, little is known about mogrosides biosynthesis in S. grosvenorii, especially the late steps of the pathway. Results: In this study, a cDNA library generated from of equal amount of RNA taken from S. grosvenorii fruit at 50 days after flowering (DAF) and 70 DAF were sequenced using Illumina/Solexa platform. More than 48,755,516 high-quality reads from a cDNA library were generated that was assembled into 43,891 unigenes. De novo assembly and gap-filling generated 43,891 unigenes with an average sequence length of 668 base pairs. A total of 26,308 (59.9%) unique sequences were annotated and 11,476 of the unique sequences were assigned to specific metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes. cDNA sequences for all of the known enzymes involved in mogrosides backbone synthesis were identified from our library. Additionally, a total of eighty-five cytochrome P450 (CYP450) and ninety UDP-glucosyltransferase (UDPG) unigenes were identified, some of which appear to encode enzymes responsible for the conversion of the mogroside backbone into the various mogrosides. Digital gene expression profile (DGE) analysis using Solexa sequencing was performed on three important stages of fruit development, and based on their expression pattern, seven CYP450 s and five UDPG s were selected as the candidates most likely to be involved in mogrosides biosynthesis. Conclusion: A combination of RNA-seq and DGE analysis based on the next generation sequencing technology was shown to be a powerful method for identifying candidate genes encoding enzymes responsible for the biosynthesis of novel secondary metabolites in a non-model plant. Seven CYP450 s and five UDPG s were selected as potential candidates involved in mogrosides biosynthesis. The transcriptome data from this study provides an important resource for understanding the formation of major bioactive constituents in the fruit extract from S. grosvenorii.

MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing. P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing. Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages. This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle.

WRKY transcription factors, as the largest gene family in higher plants, play an important role in various biological processes including growth and development, regulation of secondary metabolites, and stress response. In this study, we performed genome-wide identification and analysis of WRKY transcription factors in S. siamensis. A total of 59 SsWRKY genes were identified that were distributed on all 14 chromosomes, and these were classified into three major groups based on phylogenetic relationships. Each of these groups had similar conserved motifs and gene structures. We compared all the S. siamensis SsWRKY genes with WRKY genes identified from three diverse plant species, and the results implied that segmental duplication and tandem duplication play an important roles in the evolution processes of the WRKY gene family. Promoter region analysis revealed that SsWRKY genes included many cis-acting elements related to plant growth and development, phytohormone response, and both abiotic and biotic stress. Expression profiles originating from the transcriptome database showed expression patterns of these SsWRKY genes in four different tissues and revealed that most genes are expressed in plant roots. Fifteen SsWRKY genes with low-temperature response motifs were surveyed for their gene expression under cold stress, showing that most genes displayed continuous up-regulation during cold treatment. Our study provides a foundation for further study on the function and regulatory mechanism of the SsWRKY gene family.

Siraitia siamensis is a traditional Chinese medicinal herb. In this study, using S. siamensis cultivated in vitro, twelve candidate reference genes under various treatments were analyzed for their expression stability by using algorithms such as GeNorm, NormFinder, BestKeeper, Delta CT, and RefFinder. The selected reference genes were then used to characterize the gene expression of cucurbitadienol synthase, which is a rate-limiting enzyme for mogroside biosynthesis. The results showed that CDC6 and NCBP2 expression was the most stable across all treatments and are the best reference genes under the tested conditions. Utilizing the validated reference genes, we analyzed the expression profiles of genes related to the synthesis pathway of mogroside in S. siamensis in response to a range of abiotic stresses. The findings of this study provide clear standards for gene expression normalization in Siraitia plants and exploring the rationale behind differential gene expression related to mogroside synthesis pathways.