Expression of 3beta-HSD and P5betaR, genes respectively coding for Delta5-3beta-hydroxysteroid dehydrogenase and progesterone 5beta-reductase, in leaves and cell cultures of Digitalis lanata EHRH

Integrated transcriptomic and metabolomic analysis reveals the molecular basis of tissue-specific accumulation of bioactive steroidal alkaloids in Fritillaria unibracteata

Fritillaria unibracteata is an endangered medicinal plant whose bulb has been used as a Chinese herb to suppress cough, asthma and excessive phlegm for centuries. Steroidal alkaloids, which are synthesized via the steroid synthesis pathways, are their significant bioactive constituents. However, few studies on genes involved in steroidal alkaloid biosynthesis in F. unibracteata have been reported, mainly due to the lack of the F. unibracteata genome. In this paper, comparative transcriptomic and metabolomic analyses of four different tissues of F. unibracteata (leaves, flowers, stems, and bulbs) were performed. Imperialine, peiminine, and peimisine were among the significant bioactive compounds that were considerably abundant in bulb tissue, according to the metabolomic findings. Then, 83.60 Gb transcriptome sequencing of four different tissues was performed, of which one gene encoding phosphomevalonate kinase was directly functionally characterized to verify the accuracy of sequences obtained from the transcriptome. A total of 9217 differentially expressed unigenes (DEGs) were identified in four different tissues of F. unibracteata. GO and KEGG enrichments revealed that phenylpropanoid biosynthesis, MVA-mediated terpenoid backbone biosynthesis, and steroid biosynthesis were enriched in bulb tissue, whereas enrichment of MEP-mediated terpenoid backbone biosynthesis, photosynthesis, photosynthesis-antenna protein and carotenoid biosynthesis was observed in aerial tissues. Moreover, clustering analysis indicated that the downstream steroid biosynthesis pathway was more important in steroidal alkaloid biosynthesis compared to the upstream terpenoid backbone biosynthesis pathway. Hence, MVA-mediated biosynthesis of steroidal alkaloids was proposed, in which 15 bulb-clustered DEGs were positively correlated with a high accumulation of bioactive steroid alkaloids, further validating our proposal. In addition, 36 CYP450s showing a positive correlation with bioactive steroidal alkaloids provided candidate enzymes to catalyze the subsequent steps of steroidal alkaloid biosynthesis. In addition, the transcription factors and ABC transporters clustered in bulb tissue might be responsible for the regulation and transportation of steroidal alkaloid biosynthesis. Protein-protein interaction analysis implied a highly complex steroid alkaloid biosynthesis network in which delta (24)-sterol reductase might be one of the central catalysts. Based on the integrated transcriptome and metabolome, this current study is a first step in understanding the tissue-specific biosynthesis of steroidal alkaloids in F. unibracteata. Furthermore, key genes and regulators identified herein could facilitate metabolic engineering to improve steroidal alkaloids in F. unibracteata.

Overexpression and RNAi-mediated Knockdown of Two 3β-hydroxy-Δ5-steroid dehydrogenase Genes in Digitalis lanata Shoot Cultures Reveal Their Role in Cardenolide Biosynthesis

3β-hydroxy-Δ5-steroid dehydrogenases (3βHSDs) are supposed to be involved in 5β-cardenolide biosynthesis. Here, a novel 3βHSD (Dl3βHSD2) was isolated from Digitalis lanata shoot cultures and expressed in E. coli. Recombinant Dl3βHSD1 and Dl3βHSD2 shared 70% amino acid identity, reduced various 3-oxopregnanes and oxidised 3-hydroxypregnanes, but only rDl3βHSD2 converted small ketones and secondary alcohols efficiently. To explain these differences in substrate specificity, we established homology models using borneol dehydrogenase of Salvia rosmarinus (6zyz) as the template. Hydrophobicity and amino acid residues in the binding pocket may explain the difference in enzyme activities and substrate preferences. Compared to Dl3βHSD1, Dl3βHSD2 is weakly expressed in D. lanata shoots. High constitutive expression of Dl3βHSDs was realised by Agrobacterium-mediated transfer of Dl3βHSD genes fused to the CaMV-35S promotor into the genome of D. lanata wild type shoot cultures. Transformed shoots (35S:Dl3βHSD1 and 35S:Dl3βHSD2) accumulated less cardenolides than controls. The levels of reduced glutathione (GSH), which is known to inhibit cardenolide formation, were higher in the 35S:Dl3βHSD1 lines than in the controls. In the 35S:Dl3βHSD1 lines cardenolide levels were restored after adding of the substrate pregnane-3,20-dione in combination with buthionine-sulfoximine (BSO), an inhibitor of GSH formation. RNAi-mediated knockdown of the Dl3βHSD1 yielded several shoot culture lines with strongly reduced cardenolide levels. In these lines, cardenolide biosynthesis was fully restored after addition of the downstream precursor pregnan-3β-ol-20-one, whereas upstream precursors such as progesterone had no effect, indicating that no shunt pathway could overcome the Dl3βHSD1 knockdown. These results can be taken as the first direct proof that Dl3βHSD1 is indeed involved in 5β-cardenolide biosynthesis.

Knockout of Arabidopsis thaliana VEP1, Encoding a PRISE (Progesterone 5β-Reductase/Iridoid Synthase-Like Enzyme), Leads to Metabolic Changes in Response to Exogenous Methyl Vinyl Ketone (MVK)

Small or specialized natural products (SNAPs) produced by plants vary greatly in structure and function, leading to selective advantages during evolution. With a limited number of genes available, a high promiscuity of the enzymes involved allows the generation of a broad range of SNAPs in complex metabolic networks. Comparative metabolic studies may help to understand why—or why not—certain SNAPs are produced in plants. Here, we used the wound-induced, vein patterning regulating VEP1 (AtStR1, At4g24220) and its paralogue gene on locus At5g58750 (AtStR2) from Arabidopsis to study this issue. The enzymes encoded by VEP1-like genes were clustered under the term PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes) as it was previously demonstrated that they are involved in cardenolide and/or iridoid biosynthesis in other plants. In order to further understand the general role of PRISEs and to detect additional more “accidental” roles we herein characterized A. thaliana steroid reductase 1 (AtStR1) and compared it to A. thaliana steroid reductase 2 (AtStR2). We used A. thaliana Col-0 wildtype plants as well as VEP1 knockout mutants and VEP1 knockout mutants overexpressing either AtStR1 or AtStR2 to investigate the effects on vein patterning and on the stress response after treatment with methyl vinyl ketone (MVK). Our results added evidence to the assumption that AtStR1 and AtStR2, as well as PRISEs in general, play specific roles in stress and defense situations and may be responsible for sudden metabolic shifts.

Plastidial Expression of 3β-Hydroxysteroid Dehydrogenase and Progesterone 5β-Reductase Genes Confer Enhanced Salt Tolerance in Tobacco

The short-chain dehydrogenase/reductase (SDR) gene family is widely distributed in all kingdoms of life. The SDR genes, 3β-hydroxysteroid dehydrogenase (3β-HSD) and progesterone 5-β-reductases (P5βR1, P5βR2) play a crucial role in cardenolide biosynthesis pathway in the Digitalis species. However, their role in plant stress, especially in salinity stress management, remains unexplored. In the present study, transplastomic tobacco plants were developed by inserting the 3β-HSD, P5βR1 and P5βR2 genes. The integration of transgenes in plastomes, copy number and transgene expression at transcript and protein level in transplastomic plants were confirmed by PCR, end-to-end PCR, qRT-PCR and Western blot analysis, respectively. Subcellular localization analysis showed that 3β-HSD and P5βR1 are cytoplasmic, and P5βR2 is tonoplast-localized. Transplastomic lines showed enhanced growth in terms of biomass and chlorophyll content compared to wild type (WT) under 300 mM salt stress. Under salt stress, transplastomic lines remained greener without negative impact on shoot or root growth compared to the WT. The salt-tolerant transplastomic lines exhibited enhanced levels of a series of metabolites (sucrose, glutamate, glutamine and proline) under control and NaCl stress. Furthermore, a lower Na+/K+ ratio in transplastomic lines was also observed. The salt tolerance, mediated by plastidial expression of the 3β-HSD, P5βR1 and P5βR2 genes, could be due to the involvement in the upregulation of nitrogen assimilation, osmolytes as well as lower Na+/K+ ratio. Taken together, the plastid-based expression of the SDR genes leading to enhanced salt tolerance, which opens a window for developing saline-tolerant plants via plastid genetic engineering.

RNAi-mediated gene knockdown of progesterone 5β-reductases in Digitalis lanata reduces 5β-cardenolide content

Studying RNAi-mediated DlP5βR1 and DlP5βR2 knockdown shoot culture lines of Digitalis lanata, we here provide direct evidence for the participation of PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes) in 5β-cardenolide formation. Progesterone 5β-reductases (P5βR) are assumed to catalyze the reduction of progesterone to 5β-pregnane-3,20-dione, which is a crucial step in the biosynthesis of the 5β-cardenolides. P5βRs are encoded by VEP1-like genes occurring ubiquitously in embryophytes. P5βRs are substrate-promiscuous enone-1,4-reductases recently termed PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes). Two PRISE genes, termed DlP5βR1 (AY585867.1) and DlP5βR2 (HM210089.1) were isolated from Digitalis lanata. To give experimental evidence for the participation of PRISEs in 5β-cardenolide formation, we here established several RNAi-mediated DlP5βR1 and DlP5βR2 knockdown shoot culture lines of D. lanata. Cardenolide contents were lower in D. lanata P5βR-RNAi lines than in wild-type shoots. We considered that the gene knockdowns may have had pleiotropic effects such as an increase in glutathione (GSH) which is known to inhibit cardenolide formation. GSH levels and expression of glutathione reductase (GR) were measured. Both were higher in the Dl P5βR-RNAi lines than in the wild-type shoots. Cardenolide biosynthesis was restored by buthionine sulfoximine (BSO) treatment in Dl P5βR2-RNAi lines but not in Dl P5βR1-RNAi lines. Since progesterone is a precursor of cardenolides but can also act as a reactive electrophile species (RES), we here discriminated between these by comparing the effects of progesterone and methyl vinyl ketone, a small RES but not a precursor of cardenolides. To the best of our knowledge, we here demonstrated for the first time that P5βR1 is involved in cardenolide formation. We also provide further evidence that PRISEs are also important for plants dealing with stress by detoxifying reactive electrophile species (RES).

De novo transcriptome analysis deciphered polyoxypregnane glycoside biosynthesis pathway in Gymnema sylvestre

Gymnema sylvestre is an important medicinal plant containing antidiabetic activity. Through de novo transcriptomic study, the pathways of polyoxypregnane glycosides were explored and candidate genes of these pathways were identified in G. sylvestre. High-quality raw reads were assembled into transcripts which resulted in 193,615 unigenes. These unigenes further decoded 58,274 coding DNA sequences (CDSs). Functional annotation of predicted CDSs was carried out using the protein databases, i.e., NCBI's non-redundant, Uniprot and Pfam. Eukaryotic orthologous group (KOG) classification and transcription factor analysis has revealed most CDS-enriched categories as "Signal transduction mechanism" and "Basic Helix loop helix" (bHLH) transcription factor family, respectively. A total of 16,569 CDSs were assigned minimum one Gene Ontology (GO) term. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis disclosed 235 CDSs which represented total 27 genes of pregnane glycoside pathways and 19 CDSs represented 10 important enzymes of polyoxypregnane glycoside biosynthesis, i.e., sterol 24-C-methyltransferase, cycloeucalenol cycloisomerase, Δ¹⁴-sterol reductase, C-8,7 sterol isomerase, sterol methyltransferase 2, C-5 sterol desaturase, sterol Δ⁷ reductase, Δ²⁴ sterol reductase, 3β-hydroxysteroid dehydrogenase and progesterone 5β reductase (5βPOR). This transcriptome analysis provided an important resource for future functional genomic studies in G. sylvestre.

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.

Biotransformation of 21-O-acetyl-deoxycorticosterone by cell suspension cultures of Digitalis lanata (strain W.1.4)

Cell cultures of Digitalis species are known to accept exogenous substrates for biotransformation reactions. We here report the biotransformation of 21-O-acetyl-deoxycorticosterone (1) by cell suspension cultures of Digitalis lanata strain W.1.4. Nine derivatives of 1 were obtained and their chemical structures determined by spectroscopic methods. 2β-Hydroxylation and C-21-glucosylation of the steroidal nucleus were described for the first time in suspension-cultured plant cells. Steroid 5α- and 5β-reduction products were also observed. Among the compounds isolated and structures elucidated were 2β,3β,21-trihydroxy-4-pregnen-20-one, 2β,3α,21-trihydroxy-4-pregnen-20-one and 3β,21-dihydroxy-5α-pregnan-20-one-3β-O-β-glucoside.