PRISEs (progesterone 5β-reductase and/or iridoid synthase-like 1,4-enone reductases): Catalytic and substrate promiscuity allows for realization of multiple pathways in plant metabolism

Four enzymes control natural variation in the steroid core of Erysimum cardenolides

Plants commonly produce families of structurally related metabolites with similar defensive functions. This apparent redundancy raises the question of underlying molecular mechanisms and adaptive benefits of such chemical variation. Cardenolides, a class defensive compounds found in the wallflower genus Erysimum (L., Brassicaceae) and scattered across other plant families, show substantial structural variation, with glycosylation and hydroxylation being common modifications of a steroid core, which itself may vary in terms of stereochemistry and saturation. Through a combination of chemical mutagenesis and analysis of gene coexpression networks, we identified four enzymes involved in cardenolide biosynthesis in Erysimum that work together to determine stereochemistry at carbon 5 of the steroid core: Ec3βHSD, a 3β-hydroxysteroid dehydrogenase, Ec3KSI, a ketosteroid isomerase, EcP5βR2, a progesterone 5β-reductase, and EcDET2, a steroid 5α-reductase. We biochemically characterized the activity of these enzymes in vitro and generated CRISPR/Cas9 knockout lines to confirm activity in vivo. Cardenolide biosynthesis was not eliminated in any of the knockouts. Instead, mutant plants accumulated cardenolides with altered saturation and stereochemistry of the steroid core. Furthermore, we found variation in carbon 5 configuration among the cardenolides of 44 species of Erysimum, where the occurrence of some 5β-cardenolides is associated with the expression and sequence of P5βR2. This may have allowed Erysimum species to fine-tune their defensive profiles to target specific herbivore populations over the course of evolution.

SIGNIFICANCE STATEMENT

Plants use an array of toxic compounds to defend themselves from attack against insects and other herbivores. One mechanism through which plants may evolve more toxic compounds is through modifications to the structure of compounds they already produce. In this study, we show how plants in the wallflower genus Erysimum use four enzymes to fine-tune the structure of toxic metabolites called cardenolides. Natural variation in the sequence and expression of a single enzyme called progesterone 5β-reductase 2 partly explains the variation in cardenolides observed across the Erysimum genus. These alterations to cardenolide structure over the course of evolution suggests that there may be context-dependent benefits to Erysimum to invest in one cardenolide variant over another.

A fused hybrid enzyme of 8-hydroxygeraniol oxidoreductase (8HGO) from Gardenia jasminoides and iridoid synthase (ISY) from Catharanthus roseus significantly enhances nepetalactol and iridoid production

MAIN CONCLUSION

The operation of 8HGO-ISY fusion enzymes can increase nepetalactol flux to iridoid biosynthesis, and the Gj8HGO-CrISY expression in Gardenia jasminoides indicates that seco-iridoids and closed-ring iridoids share a nepetalactol pool. Nepetalactol is a common precursor of (seco)iridoids and their derivatives, which are a group of noncanonical monoterpenes. Functional characterization of an 8HGO (8-hydroxygeraniol oxidoreductase) from Catharanthus roseus, a seco-iridoids producing plant, has been reported; however, the 8HGO from G. jasminoides with plenty of closed-ring iridoids remains uninvestigated. In this work, a Gj8HGO was cloned and biochemically characterized. In addition, the relatively low production of nepetalactol in plants and engineered microbial host is likely to be attributed to the fact that Cr8HGO and CrISY (iridoid synthase) are substrate-promiscuous enzymes catalyzing unexpected substrates to the undesired products. Herein, a bifunctional enzyme consisting of an 8HGO fused to an ISY was designed for the proximity to the substrate and recycling of NADP+ and NADPH cofactor to reduce the undesired intermediate in the synthesis of nepetalactol. Of four fusion enzymes (i.e., Gj8HGO-GjISY, Gj8HGO-GjISY2, Gj8HGO-GjISY4, and Gj8HGO-CrISY), interestingly, only the last one can enable cascade reaction to form cis-trans-nepetalactol. Furthermore, we establish a reliable Agrobacterium-mediated transformation system. The expression of Gj8HGO-CrISY in G. jasminoides led to a significant enhancement of nepetalactol production, about 19-fold higher than that in wild-type plants, which further resulted in the twofold to fivefold increase of total iridoids and representative iridoid such as geniposide, indicating that seco-iridoids in C. roseus and closed-ring iridoids in G. jasminoides share a nepetalactol pool. All results suggest that 8HGO and ISY can be manipulated to maximize metabolic flux for nepetalactol and iridoid production.

Citronellol biosynthesis in pelargonium is a multistep pathway involving progesterone 5β-reductase and/or iridoid synthase-like enzymes

Citronellol is a pleasant-smelling compound produced in rose (Rosa spp.) flowers and in the leaves of many aromatic plants, including pelargoniums (Pelargonium spp.). Although geraniol production has been well studied in several plants, citronellol biosynthesis has been documented only in crab-lipped spider orchid (Caladenia plicata) and its mechanism remains open to question in other species. We therefore profiled 10 pelargonium accessions using RNA sequencing and gas chromatography-MS analysis. Three enzymes from the progesterone 5β-reductase and/or iridoid synthase-like enzymes (PRISE) family were characterized in vitroand subsequently identified as citral reductases (named PhCIRs). Transgenic RNAi lines supported a role for PhCIRs in the biosynthesis of citronellol as well as in the production of mint-scented terpenes. Despite their high amino acid sequence identity, the 3 enzymes showed contrasting stereoselectivity, either producing mainly (S)-citronellal or a racemate of both (R)- and (S)-citronellal. Using site-directed mutagenesis, we identified a single amino acid substitution as being primarily responsible for the enzyme's enantioselectivity. Phylogenetic analysis of pelargonium PRISEs revealed 3 clades and 7 groups of orthologs. PRISEs from different groups exhibited differential affinities toward substrates (citral and progesterone) and cofactors (NADH/NADPH), but most were able to reduce both substrates, prompting hypotheses regarding the evolutionary history of PhCIRs. Our results demonstrate that pelargoniums evolved citronellol biosynthesis independently through a 3-step pathway involving PRISE homologs and both citral and citronellal as intermediates. In addition, these enzymes control the enantiomeric ratio of citronellol thanks to small alterations of the catalytic site.

A novel steroid hydroxylase from Nigrospora sphaerica with various hydroxylation capabilities to different steroid substrates

Fungal hydroxylation of steroids is a key step in the industrial production of various steroid drugs. The main enzymes that enable these reactions are Cytochrome P450s (CYP), though very few industrially important CYPs have been identified and characterized. In this study, we identified a CYP enzyme (CYP-N2) and a cytochrome P450 reductase (CPRns) from Nigrospora sphaerica 722 by a combination of transcriptome sequencing and heterologous expression in Pichia pastoris. Gene CYP-N2 co-expressed with CPRns in Pichia pastoris GS115 showed 6β- and 15α-hydroxylation activities on progesterone. Different hydroxylation specificity of CYP-N2 was observed on different steroid substrates. CYP-N2 showed 1α-hydroxylation on cortisone and 1α-hydroxylation and 6β-hydroxylation activities on androstenedione (AD). With dehydroepiandrosterone (DHEA) as a substrate, the hydroxylated products of CYP-N2 included 7α-hydroxy-DHEA and 7α,15α-dihydroxy-DHEA. In order to precisely elucidate CYP-N2 biological function and find out the key amino acids influencing its hydroxylation capabilities in the binding pocket, new generation artificial intelligence technology AlphaFold 2 was used to predict the function-structure of CYP-N2 with high reliability. Through molecular docking, it was concluded that the residues almost binding all substrates were located in the same substrate binding pocket and the various hydroxylation abilities might be due to the different binding conformations of different substrates in the binding pocket. Alanine scanning mutagenesis was used to verify key amino acids identified by the molecular docking with steroid substrates. The 128 THR mutation resulted in conversion rate increase for substrates AD and cortisone by 2.6-fold and 2.1-fold respectively. The information obtained in this study is beneficial to facilitating the engineering of more efficient steroid hydroxylases for industrial applications.

Cardenolide Increase in Foxglove after 2,1,3-Benzothiadiazole Treatment Reveals a Potential Link between Cardenolide and Phytosterol Biosynthesis

Cardenolides are steroidal metabolites in Digitalis lanata with potent cardioactive effects on animals. In plants, cardenolides are likely involved in various stress responses. However, the molecular mechanism of cardenolide increase during stresses is mostly unknown. Additionally, cardenolides are proposed to arise from cholesterol, but indirect results show that phytosterols may also be substrates for cardenolide biosynthesis. Here, we show that cardenolides increased after methyl jasmonate (MJ), sorbitol, potassium chloride (KCl) and salicylic acid analog [2,1,3-benzothiadiazole (BTH)] treatments. However, the expression of three known genes for cardenolide biosynthesis did not correlate well with these increases. Specifically, the expression of progesterone-5β-reductases (P5βR and P5βR2) did not correlate with the cardenolide increase. The expression of 3β-hydroxysteroid dehydrogenase (3βHSD) correlated with changes in cardenolide levels only during the BTH treatment. Mining the D. lanata transcriptome identified genes involved in cholesterol and phytosterol biosynthesis: C24 sterol sidechain reductase 1 (SSR1), C4 sterol methyl oxidase 1, and 3 (SMO1 and SMO3). Surprisingly, the expression of all three genes correlated well with the cardenolide increase after the BTH treatment. Phylogenetic analysis showed that SSR1 is likely involved in both cholesterol and phytosterol biosynthesis. In addition, SMO1 is likely specific to phytosterol biosynthesis, and SMO3 is specific to cholesterol biosynthesis. These results suggest that stress-induced increase of cardenolides in foxglove may correlate with cholesterol and phytosterol biosynthesis. In summary, this work shows that cardenolides are important for stress responses in D. lanata and reveals a potential link between phytosterol and cardenolide biosynthesis.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

Identification and functional characterization of three iridoid synthases in Gardenia jasminoides

The discovery of three iridoid synthases (GjISY, GjISY2 and GjISY4) from Gardenia jasminoides and their functional characterization increase the understanding of iridoid scaffold/iridoid glycoside biosynthesis in iridoid-producing plants. Iridoids are a class of noncanonical monoterpenes that are found naturally in the plant kingdom mostly as glycosides. Over 40 iridoid glycosides (e.g., geniposide, gardenoside and shanzhiside) have been isolated from Gardenia jasminoides. They have multiple pharmacological properties and health-promoting effects. However, their biosynthetic pathway is poorly understood, and the iridoid synthase (ISY) responsible for the cyclization of the core scaffold remains unclear. In this study, three homologs of ISYs from G. jasminoides (GjISY, GjISY2 and GjISY4) were identified on the basis of transcriptomic data and functionally characterized. The genomic structure and intron-exon arrangement revealed that all three ISYs contained an intron. Biochemical assays indicated that all three recombinant enzymes reduced 8-oxogeranial to nepetalactol and its open forms (iridodials) as the products of the classical CrISY (Catharanthus roseus). In addition, all three enzymes reduced progesterone to 5-β-prognane-3,20-dione. However, only GjISY2 and GjISY4 reduced 2-cyclohexen-1-one to cyclohexanone. Overall, the GjISY2 expression levels in the flowers and fruits were similar to the GjISY and GjISY4 expression levels. By contrast, the GjISY2 expression levels in the upper and lower leaves were substantially higher than the GjISY and GjISY4 expression levels. Among the three, GjISY2 exhibited the highest catalytic efficiency for 8-oxogeranial. GjISY2 might be the major contributor to iridoid biosynthesis in G. jasminoides. Collectively, our results advance the understanding of iridoid scaffold/iridoid glycoside biosynthesis in G. jasminoides and provide a potential target for metabolic engineering and breeding.

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.

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).

The Progesterone 5β-Reductase/Iridoid Synthase Family: A Catalytic Reservoir for Specialized Metabolism across Land Plants

Iridoids are plant-derived terpenoids with a rich array of bioactivities. The key step in iridoid skeleton formation is the reduction of 8-oxogeranial by certain members of the progesterone 5β-reductase/iridoid synthase (PRISE) family of short-chain alcohol dehydrogenases. Other members of the PRISE family have previously been implicated in the biosynthesis of the triterpenoid class of cardenolides, which requires the reduction of progesterone. Here, we explore the occurrence and activity of PRISE across major lineages of plants. We observed trace activities toward either 8-oxogeranial or progesterone in all PRISEs, including those from nonseed plants and green algae. Phylogenetic analysis, coupled with enzymatic assays, show that these activities appear to have become specialized in specific angiosperm lineages. This broad analysis of the PRISE family provides insight into how these enzymes evolved in plants and also suggests that iridoid synthase activity is an ancestral trait in all land plants, which might have contributed to the rise of iridoid metabolites.

Unveiling the genetic basis of Sclerotinia head rot resistance in sunflower

BACKGROUND

Sclerotinia sclerotiorum is a necrotrophic fungus that causes Sclerotinia head rot (SHR) in sunflower, with epidemics leading to severe yield losses. In this work, we present an association mapping (AM) approach to investigate the genetic basis of natural resistance to SHR in cultivated sunflower, the fourth most widely grown oilseed crop in the world.

RESULTS

Our association mapping population (AMP), which comprises 135 inbred breeding lines (ILs), was genotyped using 27 candidate genes, a panel of 9 Simple Sequence Repeat (SSR) markers previously associated with SHR resistance via bi-parental mapping, and a set of 384 SNPs located in genes with molecular functions related to stress responses. Moreover, given the complexity of the trait, we evaluated four disease descriptors (i.e, disease incidence, disease severity, area under the disease progress curve for disease incidence, and incubation period). As a result, this work constitutes the most exhaustive AM study of disease resistance in sunflower performed to date. Mixed linear models accounting for population structure and kinship relatedness were used for the statistical analysis of phenotype-genotype associations, allowing the identification of 13 markers associated with disease reduction. The number of favourable alleles was negatively correlated to disease incidence, disease severity and area under the disease progress curve for disease incidence, whereas it was positevily correlated to the incubation period.

CONCLUSIONS

Four of the markers identified here as associated with SHR resistance (HA1848, HaCOI_1, G33 and G34) validate previous research, while other four novel markers (SNP117, SNP136, SNP44, SNP128) were consistently associated with SHR resistance, emerging as promising candidates for marker-assisted breeding. From the germplasm point of view, the five ILs carrying the largest combination of resistance alleles provide a valuable resource for sunflower breeding programs worldwide.

In vitro prototyping of limonene biosynthesis using cell-free protein synthesis

Metabolic engineering of microorganisms to produce sustainable chemicals has emerged as an important part of the global bioeconomy. Unfortunately, efforts to design and engineer microbial cell factories are challenging because design-build-test cycles, iterations of re-engineering organisms to test and optimize new sets of enzymes, are slow. To alleviate this challenge, we demonstrate a cell-free approach termed in vitro Prototyping and Rapid Optimization of Biosynthetic Enzymes (or iPROBE). In iPROBE, a large number of pathway combinations can be rapidly built and optimized. The key idea is to use cell-free protein synthesis (CFPS) to manufacture pathway enzymes in separate reactions that are then mixed to modularly assemble multiple, distinct biosynthetic pathways. As a model, we apply our approach to the 9-step heterologous enzyme pathway to limonene in extracts from Escherichia coli. In iterative cycles of design, we studied the impact of 54 enzyme homologs, multiple enzyme levels, and cofactor concentrations on pathway performance. In total, we screened over 150 unique sets of enzymes in 580 unique pathway conditions to increase limonene production in 24 h from 0.2 to 4.5 mM (23-610 mg/L). Finally, to demonstrate the modularity of this pathway, we also synthesized the biofuel precursors pinene and bisabolene. We anticipate that iPROBE will accelerate design-build-test cycles for metabolic engineering, enabling data-driven multiplexed cell-free methods for testing large combinations of biosynthetic enzymes to inform cellular design.

Identification of key sites determining the cofactor specificity and improvement of catalytic activity of a steroid 5β-reductase from Capsella rubella

Progesterone 5β-reductases (P5βRs) are involved in 5β-cardenolide formation by stereo-specific reduction of the △^4,5 double bond of steroid precursors. In this study a steroid 5β-reductase was identified in Capsella rubella (CrSt5βR1) and its function in steroid 5β-reduction was validated experimentally. CrSt5βR1 is capable of enantioselectively reducing the activated CC bond of broad substrates such as steroids and enones by using NADPH as a cofactor and therefore has the potential as a biocatalyst in organic synthesis. However, for industrial purposes the cheaper NADH is the preferred cofactor. By applying rational design based on literature and complementary mutagenesis strategies, we successfully identified two key amino acid residues determining the cofactor specificity of the enzyme. The R63 K mutation enables the enzyme to convert progesterone to 5β-pregnane-3,20-dione with NADH as cofactor, whereas the wild-type CrSt5βR1 is strictly NADPH-dependent. By further introducing the R64H mutation, the double mutant R63K_R64H of CrSt5βR1 was shown to increase enzymatic activity by13.8-fold with NADH as a cofactor and to increase the NADH/NADPH conversion ratio by 10.9-fold over the R63 K single mutant. This finding was successfully applied to change the cofactor specificity and to improve activity of other members of the same enzyme family, AtP5βR and DlP5βR. CrSt5βR1 mutants are expected to have the potential for biotechnological applications in combination with the well-established NADH regeneration systems.

Biosynthetic approach to combine the first steps of cardenolide formation in Saccharomyces cerevisiae

A yeast expression plasmid was constructed containing a cardenolide biosynthetic module, referred to as CARD II, using the AssemblX toolkit, which enables the assembly of large DNA constructs. The genes cloned into the vector were (a) a Δ⁵‐3β‐hydroxysteroid dehydrogenase gene from Digitalis lanata, (b) a steroid Δ⁵‐isomerase gene from Comamonas testosteronii, (c) a mutated steroid‐5β‐reductase gene from Arabidopsis thaliana, and (d) a steroid 21‐hydroxylase gene from Mus musculus. A second plasmid bearing an ADR/ADX fusion gene from Bos taurus was also constructed. A Saccharomyces cerevisiae strain bearing these two plasmids was generated. This strain, termed “CARD II yeast”, was capable of producing 5β‐pregnane‐3β,21‐diol‐20‐one, a central intermediate in 5β‐cardenolide biosynthesis, starting from pregnenolone which was added to the culture medium. Using this approach, five consecutive steps in cardenolide biosynthesis were realized in baker's yeast.

Exploiting enzyme promiscuity to shape plant specialized metabolism

The amazing variability of plant metabolism and its rapid divergence during evolution pose fundamental questions as to the driving forces, mechanisms, and players in metabolic differentiation. This review examines concepts that help us understand adaptive pathway evolution, with a particular emphasis on plant specialized metabolism, previously often termed secondary metabolism. Following a general introduction to pathway and metabolite evolution, the focus is directed to enzyme promiscuity and its classification. Promiscuous enzymes (or substrates), 'silent' elements of the metabolome, and the 'underground metabolism' may be used and combined to evolve 'new' metabolic pathways. It appears that new pathways rarely appear from scratch, but instead emerge from 'floppy' enzymes and elements of a 'messy' metabolism, and in this way a range of metabolites is generated, some of which may provide benefits to the plant.