De novo transcriptome assembly and the putative biosynthetic pathway of steroidal sapogenins of Dioscorea composita

Endophyte-inoculated rhizomes of Paris polyphylla improve polyphyllin biosynthesis and yield: a transcriptomic analysis of the underlying mechanism

Introduction

Polyphyllin from Paris polyphylla var. yunnanensis exhibits anti-inflammatory, analgesic, antibacterial, and antiviral properties. However, the current production of polyphyllin can barely meet market demand. To improve the content of polyphyllin produced by P. polyphylla, two endophyte strains, Bacillus cereus LgD2 and Fusarium oxysporum TPB, were isolated from Paris fargesii Franch. and inoculated in the roots of P. polyphylla. Both symbiotic strains significantly promoted the accumulation of saponins in P. polyphylla.

Methods

The content of polyphyllin in rhizomes of P. polyphylla treated with TPB with LgD2 strain was determined using High Performance Liquid Chromatography and the expressed genes were analyzed by RNA-seq. Gene Ontology and Kyoto Encyclopedia of Genes annotations were performed on the differentially expressed genes, a clustering tree of UDP-glycosyltransferase (UGT) and cytochrome P450 (CYP450) gene families was constructed, and UGT and CYP450 involved in the biosynthesis of polyphyllin were predicted using weighted correlation network analysis (WGCNA).

Results

RNA-seq and qRT-PCR analyses showed that endophytic inoculation did not promote polyphyllin accumulation by enhancing the upstream terpene biosynthesis pathway, but probably by up-regulating the downstream CYP450 and UGT genes associated with polyphyllin biosynthesis. Genomes enrichment analyses of differentially expressed genes indicated that inoculation with LgD2 and TPB played a positive role in promoting the defense against pathogenic bacteria, enhancing the biosynthesis of carbohydrates, attenuating the process of nitrogen metabolism, and maintaining the equilibrium of the redox reaction homeostasis, potentially indirectly enhancing the polyphyllin yield of P. polyphylla. By combining differentially expressed genes screening, WGCNA, and phylogenetic tree analyses, 17 CYP450 and 2 UGT candidate genes involved in the biosynthesis of polyphyllin I, polyphyllin II, polyphyllin VII, polyphyllin D, and polyphyllin H were identified. These results suggest that endophytes probably effectively promote the accumulation of polyphyllin by regulating key downstream genes in biosynthetic pathways.

Discussion

This study provides a new approach for investigating the regulatory mechanisms of endophytes that promote the production and accumulation of polyphyllin in P. polyphylla, providing a basis for further elucidating the mechanisms of plant-endophyte interactions.

Comparative transcriptome analysis of different tissues of Solanum khasianum reveals candidate genes involved in steroidal glycoalkaloid biosynthesis

Fruits and leaves of Solanum khasianum C. B. Clarke have long been used as a common Chinese herbal medicine. Steroidal glycoalkaloids (SGAs), the main active ingredient in S. khasianum, exhibit various pharmacological effects. However, genes involved in the SGA biosynthetic pathway in S. khasianum have not yet been identified. Genes encoding potential key SGA biosynthesis enzymes were identified through comprehensive RNA sequencing analysis (RNA-seq) of S. khasianum leaves, stems, and fruits. A total of 123,704 unigenes were obtained, of which 109,775 (88.74%) were annotated in seven public databases. Among these, 54 unigenes potentially involved in SGA biosynthesis were identified. Additionally, 23,636 differentially expressed genes were identified by comparing gene expression levels among the fruits, stems, and leaves of S. khasianum. The structural characteristics and phylogenetic relationship of cycloartenol synthase involved in SGA biosynthesis were further analyzed. Solasodine constituent was detected by high-performance liquid chromatography. This is the first study to report the comparative transcriptome analysis of different tissues of S. khasianum that identifies valuable genes potentially involved in SGA biosynthesis in this species.

Dioscorea composita WRKY5 positively regulates AtSOD1 and AtABF2 to enhance drought and salt tolerances

KEY MESSAGE

DcWRKY5 increases the antioxidant enzyme activity and proline accumulation, oppositely, reduces the accumulation of ROS and MDA, through directly activating the genes expression, finally enhances the salt and drought tolerance. Drought and salinity are two main environmental factors that limit the large-scale cultivation of the medicinal plant Dioscorea composita (D. composita). WRKY transcription factors (TFs) play vital roles in regulating drought and salt tolerance in plants. Nevertheless, the molecular mechanism of WRKY TF mediates drought and salt resistance of D. composita remains largely unknown. Here, we isolated and characterized a WRKY TF from D. composita, namely DcWRKY5, which was localized to the nucleus and bound to the W-box cis-acting elements. Expression pattern analysis showed that it was highly expressed in root and significantly up-regulated in the presence of salt, polyethylene glycol-6000 (PEG-6000) and abscisic acid (ABA). Heterologous expression of DcWRKY5 increased salt and drought tolerance in Arabidopsis, but was insensitive to ABA. In addition, compared with the wild type, the DcWRKY5 overexpressing transgenic lines had more proline, higher antioxidant enzyme (POD, SOD, and CAT) activities, less reactive oxygen species (ROS) and malondialdehyde (MDA). Correspondingly, the overexpression of DcWRKY5 modulated the expression of genes related to salt and drought stresses, such as AtSS1, AtP5CS1, AtCAT, AtSOD1, AtRD22, and AtABF2. Dual luciferase assay and Y1H were further confirmed that DcWRKY5 activate the promoter of AtSOD1 and AtABF2 through directly binding to the enrichment region of the W-box cis-acting elements. These results suggest that DcWRKY5 is a positive regulator of the drought and salt tolerance in D. composita and has potential applications in transgenic breeding.

The Identification of SQS/SQE/OSC Gene Families in Regulating the Biosynthesis of Triterpenes in Potentilla anserina

The tuberous roots of Potentilla anserina (Pan) are an edible and medicinal resource in Qinghai-Tibetan Plateau, China. The triterpenoids from tuberous roots have shown promising anti-cancer, hepatoprotective, and anti-inflammatory properties. In this study, we carried out phylogenetic analysis of squalene synthases (SQSs), squalene epoxidases (SQEs), and oxidosqualene cyclases (OSCs) in the pathway of triterpenes. In total, 6, 26, and 20 genes of SQSs, SQEs, and OSCs were retrieved from the genome of Pan, respectively. Moreover, 6 SQSs and 25 SQEs genes expressed in two sub-genomes (A and B) of Pan. SQSs were not expanded after whole-genome duplication (WGD), and the duplicated genes were detected in SQEs. Twenty OSCs were divided into two clades of cycloartenol synthases (CASs) and β-amyrin synthases (β-ASs) by a phylogenetic tree, characterized with gene duplication and evolutionary divergence. We speculated that β-ASs and CASs may participate in triterpenes synthesis. The data presented act as valuable references for future studies on the triterpene synthetic pathway of Pan.

The genome of Dioscorea zingiberensis sheds light on the biosynthesis, origin and evolution of the medicinally important diosgenin saponins

Diosgenin saponins isolated from Dioscorea species such as D. zingiberensis exhibit a broad spectrum of pharmacological activities. Diosgenin, the aglycone of diosgenin saponins, is an important starting material for the production of steroidal drugs. However, how plants produce diosgenin saponins and the origin and evolution of the diosgenin saponin biosynthetic pathway remain a mystery. Here we report a high-quality, 629-Mb genome of D. zingiberensis anchored on 10 chromosomes with 30 322 protein-coding genes. We reveal that diosgenin is synthesized in leaves ('source'), then converted into diosgenin saponins, and finally transported to rhizomes ('sink') for storage in plants. By evaluating the distribution and evolutionary patterns of diosgenin saponins in Dioscorea species, we find that diosgenin saponin-containing may be an ancestral trait in Dioscorea and is selectively retained. The results of comparative genomic analysis indicate that tandem duplication coupled with a whole-genome duplication event provided key evolutionary resources for the diosgenin saponin biosynthetic pathway in the D. zingiberensis genome. Furthermore, comparative transcriptome and metabolite analysis among 13 Dioscorea species suggests that specific gene expression patterns of pathway genes promote the differential evolution of the diosgenin saponin biosynthetic pathway in Dioscorea species. Our study provides important insights and valuable resources for further understanding the biosynthesis, evolution, and utilization of plant specialized metabolites such as diosgenin saponins.

Illuminating the Plant Rhabdovirus Landscape through Metatranscriptomics Data

Rhabdoviruses infect a large number of plant species and cause significant crop diseases. They have a negative-sense, single-stranded unsegmented or bisegmented RNA genome. The number of plant-associated rhabdovirid sequences has grown in the last few years in concert with the extensive use of high-throughput sequencing platforms. Here, we report the discovery of 27 novel rhabdovirus genomes associated with 25 different host plant species and one insect, which were hidden in public databases. These viral sequences were identified through homology searches in more than 3000 plant and insect transcriptomes from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) using known plant rhabdovirus sequences as the query. The identification, assembly and curation of raw SRA reads resulted in sixteen viral genome sequences with full-length coding regions and ten partial genomes. Highlights of the obtained sequences include viruses with unique and novel genome organizations among known plant rhabdoviruses. Phylogenetic analysis showed that thirteen of the novel viruses were related to cytorhabdoviruses, one to alphanucleorhabdoviruses, five to betanucleorhabdoviruses, one to dichorhaviruses and seven to varicosaviruses. These findings resulted in the most complete phylogeny of plant rhabdoviruses to date and shed new light on the phylogenetic relationships and evolutionary landscape of this group of plant viruses. Furthermore, this study provided additional evidence for the complexity and diversity of plant rhabdovirus genomes and demonstrated that analyzing SRA public data provides an invaluable tool to accelerate virus discovery, gain evolutionary insights and refine virus taxonomy.

Cloning and Functional Characterization of a Squalene Synthase from Paris polyphylla var. yunnanensis

Paris polyphylla Smith var. yunnanensis (Franch.) Hand. - Mazz. is a precious traditional Chinese medicine, and steroidal saponins are its major bioactive constituents possessing extensive biological activities. Squalene synthase (SQS) catalyzes the first dedicated step converting two molecular of farnesyl diphosphate (FDP) into squalene, a key intermediate in the biosynthetic pathway of steroidal saponins. In this study, a squalene synthase gene (PpSQS1) was cloned and functionally characterized from P. polyphylla var. yunnanensis, representing the first identified SQS from the genus Paris. The open reading frame of PpSQS1 is 1239 bp, which encodes a protein of 412 amino acids showing high similarity to those of other plant SQSs. Expression of PpSQS1 in Escherichia coli resulted in production of soluble recombinant proteins. Gas chromatography-mass spectrometry analysis showed that the purified recombinant PpSQS1 protein could produce squalene using FDP as a substrate in the in vitro enzymatic assay. qRT-PCR analysis indicated that PpSQS1 was highly expressed in rhizomes, consistent with the dominant accumulation of steroidal saponins there, suggesting that PpSQS1 is likely involved in the biosynthesis of steroidal saponins in the plant. The findings lay a foundation for further investigation on the biosynthesis and regulation of steroidal saponins, and also provide an alternative gene for manipulation of steroid production using synthetic biology.

A customized nuclear target enrichment approach for developing a phylogenomic baseline for Dioscorea yams (Dioscoreaceae)

PREMISE

We developed a target enrichment panel for phylogenomic studies of Dioscorea, an economically important genus with incompletely resolved relationships.

METHODS

Our bait panel comprises 260 low- to single-copy nuclear genes targeted to work in Dioscorea, assessed here using a preliminary taxon sampling that includes both distantly and closely related taxa, including several yam crops and potential crop wild relatives. We applied coalescent-based and maximum likelihood phylogenomic inference approaches to the pilot taxon set, incorporating new and published transcriptome data from additional species.

RESULTS

The custom panel retrieved ~94% of targets and >80% of full gene length from 88% and 68% of samples, respectively. It has minimal gene overlap with existing panels designed for angiosperm-wide studies and generally recovers longer and more variable targets. Pilot phylogenomic analyses consistently resolve most deep and recent relationships with strong support across analyses and point to revised relationships between the crop species D. alata and candidate crop wild relatives.

DISCUSSION

Our customized panel reliably retrieves targeted loci from Dioscorea, is informative for resolving relationships in denser samplings, and is suitable for refining our understanding of the independent origins of cultivated yam species; the panel likely has broader promise for phylogenomic studies across Dioscoreales.

Recent advances in steroidal saponins biosynthesis and in vitro production

Steroidal saponins exhibited numerous pharmacological activities due to the modification of their backbone by different cytochrome P450s (P450) and UDP glycosyltransferases (UGTs). Plant-derived steroidal saponins are not sufficient for utilizing them for commercial purpose so in vitro production of saponin by tissue culture, root culture, embryo culture, etc, is necessary for its large-scale production. Saponin glycosides are the important class of plant secondary metabolites, which consists of either steroidal or terpenoidal backbone. Due to the existence of a wide range of medicinal properties, saponin glycosides are pharmacologically very important. This review is focused on important medicinal properties of steroidal saponin, its occurrence, and biosynthesis. In addition to this, some recently identified plants containing steroidal saponins in different parts were summarized. The high throughput transcriptome sequencing approach elaborates our understanding related to the secondary metabolic pathway and its regulation even in the absence of adequate genomic information of non-model plants. The aim of this review is to encapsulate the information related to applications of steroidal saponin and its biosynthetic enzymes specially P450s and UGTs that are involved at later stage modifications of saponin backbone. Lastly, we discussed the in vitro production of steroidal saponin as the plant-based production of saponin is time-consuming and yield a limited amount of saponins. A large amount of plant material has been used to increase the production of steroidal saponin by employing in vitro culture technique, which has received a lot of attention in past two decades and provides a way to conserve medicinal plants as well as to escape them for being endangered.

Ethanol production from mixtures of sugarcane bagasse and Dioscorea composita extracted residue with high solid loading

Various mixing ratios of alkali pretreated sugarcane bagasse and starch-rich waste Dioscorea composita hemls extracted residue (DER) were evaluated via simultaneous saccharification and fermentation (SSF) with 12% (w/w) solid loading, and the mixture ratio of 1:1 achieved the highest ethanol concentration and yield. When the solid loading was increased from 12% to 32%, the ethanol concentration was increased to 72.04 g/L, whereas the ethanol yield was reduced from 84.40% to 73.71%. With batch feeding and the addition of 0.1% (w/v) Tween 80, the final ethanol concentration and yield of SSF at 34% loading were 82.83 g/L and 77.22%, respectively. Due to the integration with existing starch-based ethanol industry, the co-fermentation is expected to be a competitive alternative form for cellulosic ethanol production.

Identification, characterization and expression analysis of genes involved in steroidal saponin biosynthesis in Dracaena cambodiana

Dracaena cambodiana is a traditional medicinal plant used for producing dragon's blood. The plants and dragon's blood of D. cambodiana contain a rich variety of steroidal saponins. However, little is known about steroidal saponin biosynthesis and its regulation in D. cambodiana. Here, 122 genes encoding enzymes involved in steroidal saponin biosynthesis were identified based on transcriptome data, with 29 of them containing complete open reading frames (ORF). Transcript expression analysis revealed that several genes related to steroidal saponin biosynthesis showed distinct tissue-specific expression patterns; the expression levels of genes encoding the key enzymes involved in the biosynthesis and early modification of steroidal saponins were significantly down-regulated in the stems in response to the inducer of dragon's blood, exhibiting positive correlations with the content of steroidal saponins. These results provide insights on the steroidal saponins biosynthetic pathway and mechanisms underlying induced formation of dragon's blood in D. cambodiana.

De novo assembly and annotation of the Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) transcriptome in different growth stages

Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) is a traditional medicinal herb in the family Liliaceae that has significant pharmacological effects on immunity and cardiovascular disease. In this study, three different growth stages of Zhe-Maidong were investigated using RNA-seq, and a total of 16.4 Gb of raw data was obtained. After filtering and assembling, 96,738 unigenes with an average length of 605.3 bp were ultimately generated. A total of 77,300 unigenes were annotated using information from five databases, including the NT, NR, SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases. Additionally, the mechanisms of flavonoid, saponin and polysaccharide biosynthesis and of accumulation at different stages of tuber development were also characterized. From the first to third years, the contents of flavonoids, saponins and polysaccharides all increased, whereas the expression levels of related genes decreased in the flavonoid and saponin pathways and first increased and then decreased in the polysaccharide pathway. The results of this study provide the most comprehensive expressed sequence resource for Zhe-Maidong and will expand the available O. japonicus gene library and facilitate further genome-wide research and analyses of this species.

Next-generation sequencing of representational difference analysis products for identification of genes involved in diosgenin biosynthesis in fenugreek (Trigonella foenum-graecum)

Representational difference analysis of cDNA was performed and differential products were sequenced and annotated. Candidate genes involved in biosynthesis of diosgenin in fenugreek were identified. Detailed mechanism of diosgenin synthesis was proposed. Fenugreek (Trigonella foenum-graecum L.) is a valuable medicinal and crop plant. It belongs to Fabaceae family and has a unique potential to synthesize valuable steroidal saponins, e.g., diosgenin. Elicitation (methyl jasmonate) and precursor feeding (cholesterol and squalene) were used to enhance the content of sterols and steroidal sapogenins in in vitro grown plants for representational difference analysis of cDNA (cDNA-RDA). To identify candidate genes involved in diosgenin biosynthesis, differential, factor-specific libraries were subject to the next-generation sequencing. Approximately 9.9 million reads were obtained, trimmed, and assembled into 31,491 unigenes with an average length of 291 bp. Then, functional annotation and gene ontogeny enrichment analysis was performed by aligning all-unigenes with public databases. Within the transcripts related to sterol and steroidal saponin biosynthesis, we discovered novel candidate genes of diosgenin biosynthesis and validated their expression using quantitative RT-PCR analysis. Based on these findings, we supported the idea that diosgenin is biosynthesized from cycloartenol via cholesterol. This is the first report on the next-generation sequencing of cDNA-RDA products. Analysis of the transcriptomes enriched in low copy sequences contributed substantially to our understanding of the biochemical pathways of steroid synthesis in fenugreek.

Spatial transcriptome analysis provides insights of key gene(s) involved in steroidal saponin biosynthesis in medicinally important herb Trillium govanianum

Trillium govanianum, an endangered medicinal herb native to the Himalaya, is less studied at the molecular level due to the non-availability of genomic resources. To facilitate the basic understanding of the key genes and regulatory mechanism of pharmaceutically important biosynthesis pathways, first spatial transcriptome sequencing of T. govanianum was performed. 151,622,376 (~11.5 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 69,174 transcripts. Functional annotation with multiple public databases identified array of genes involved in steroidal saponin biosynthesis and other secondary metabolite pathways including brassinosteroid, carotenoid, diterpenoid, flavonoid, phenylpropanoid, steroid and terpenoid backbone biosynthesis, and important TF families (bHLH, MYB related, NAC, FAR1, bZIP, B3 and WRKY). Differentially expressed large number of transcripts, together with CYPs and UGTs suggests involvement of these candidates in tissue specific expression. Combined transcriptome and expression analysis revealed that leaf and fruit tissues are the main site of steroidal saponin biosynthesis. In conclusion, comprehensive genomic dataset created in the current study will serve as a resource for identification of potential candidates for genetic manipulation of targeted bioactive metabolites and also contribute for development of functionally relevant molecular marker resource to expedite molecular breeding and conservation efforts in T. govanianum.