Lessons learned from metabolic engineering in hairy roots: Transcriptome and metabolic profile changes caused by Rhizobium-mediated plant transformation in Cucurbitaceae species.

Cucurbitaceae species are used in traditional medicine around the world. Cucurbitacins are highly oxygenated triterpenoids found in Cucurbitaceae species and exhibit potent anticancer activity alone and in combination with other existing chemotherapeutic drugs. Therefore, increasing production of these specialized metabolites is of great relevance. We recently showed that hairy roots of Cucurbita pepo can be used as a platform for metabolic engineering of cucurbitacins to modify their structure and increase their production. To study the changes in cucurbitacin accumulation upon formation of hairy roots, an empty vector (EV) control and Cucurbitacin inducing bHLH transcription factor 1 (CpCUCbH1)-overexpressing hairy roots of C. pepo were compared to untransformed (WT) roots. Whilst CpCUCbH1-overexpression increased production of cucurbitacins I and B by 5-fold, and cucurbitacin E by 3-fold when compared to EV lines, this increase was not significantly different when compared to WT roots. This indicated that Rhizobium rhizogenes transformation lowered the cucurbitacins levels in hairy roots, but that increasing expression of cucurbitacin biosynthetic genes by CpCUCbH1-overexpression restored cucurbitacin production to WT levels. Subsequent metabolomic and RNA-seq analysis indicated that the metabolic profile and transcriptome of hairy roots was significantly changed when compared to WT roots. Interestingly, it was observed that 11% of the differentially expressed genes were transcription factors. It was noteworthy that the majority of transcripts showing highest Pearson correlation coefficients to the Rhizobium rhizogenes genes rolB, rolC and ORF13a, were predicted to be transcription factors. In summary, hairy roots are an excellent platform for metabolic engineering of plant specialized metabolites, but these extensive transcriptome and metabolic profile changes should be considered in subsequent studies.

A Combination of Metabolomics and Machine Learning Results in the Identification of a New Cyst Nematode Hatching Factor.

Potato Cyst Nematodes (PCNs) are an economically important pest for potato growers. A crucial event in the life cycle of the nematode is hatching, after which the juvenile will move toward the host root and infect it. The hatching of PCNs is induced by known and unknown compounds in the root exudates of host plant species, called hatching factors (HFs, induce hatching independently), such as solanoeclepin A (solA), or hatching stimulants (HSs, enhance hatching activity of HFs). Unraveling the identity of unknown HSs and HFs and their natural variation is important for the selection of cultivars that produce low amounts of HFs and HSs, thus contributing to more sustainable agriculture. In this study, we used a new approach aimed at the identification of new HFs and HSs for PCNs in potato. Hereto, root exudates of a series of different potato cultivars were analyzed for their PCN hatch-inducing activity and their solA content. The exudates were also analyzed using untargeted metabolomics, and subsequently the data were integrated using machine learning, specifically random forest feature selection, and Pearson's correlation testing. As expected, solA highly correlates with hatching. Furthermore, this resulted in the discovery of a number of metabolite features present in the root exudate that correlate with hatching and solA content, and one of these is a compound of m/z 526.18 that predicts hatching even better than solA with both data methods. This compound's involvement in hatch stimulation was confirmed by the fractionation of three representative root exudates and hatching assays with the resulting fractions. Moreover, the compound shares mass fragmentation similarity with solA, and we therefore assume it has a similar structure. With this work, we show that potato likely produces a solA analogue, and we contribute to unraveling the hatch-inducing cocktail exuded by plant roots.

The tomato cytochrome P450 CYP712G1 catalyses the double oxidation of orobanchol en route to the rhizosphere signalling strigolactone, solanacol.

Strigolactones (SLs) are rhizosphere signalling molecules and phytohormones. The biosynthetic pathway of SLs in tomato has been partially elucidated, but the structural diversity in tomato SLs predicts that additional biosynthetic steps are required. Here, root RNA-seq data and co-expression analysis were used for SL biosynthetic gene discovery. This strategy resulted in a candidate gene list containing several cytochrome P450s. Heterologous expression in Nicotiana benthamiana and yeast showed that one of these, CYP712G1, can catalyse the double oxidation of orobanchol, resulting in the formation of three didehydro-orobanchol (DDH) isomers. Virus-induced gene silencing and heterologous expression in yeast showed that one of these DDH isomers is converted to solanacol, one of the most abundant SLs in tomato root exudate. Protein modelling and substrate docking analysis suggest that hydroxy-orbanchol is the likely intermediate in the conversion from orobanchol to the DDH isomers. Phylogenetic analysis demonstrated the occurrence of CYP712G1 homologues in the Eudicots only, which fits with the reports on DDH isomers in that clade. Protein modelling and orobanchol docking of the putative tobacco CYP712G1 homologue suggest that it can convert orobanchol to similar DDH isomers as tomato.

UPLC-MS/MS analysis and biological activity of the potato cyst nematode hatching stimulant, solanoeclepin A, in the root exudate of Solanum spp.

MAIN CONCLUSION: Solanoeclepin A is a hatching stimulant for potato cyst nematode in very low (pM) concentrations. We report a highly sensitive method for the analysis of SolA in plant root exudates using UHPLC-MS/MS and show that there is considerable natural variation in SolA production in Solanum spp. corresponding with their hatching inducing activity. Potato cyst nematode (PCN) is a plant root sedentary endoparasite, specialized in the infection of solanaceous species such as potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Earlier reports (Mulder et al. in Hatching agent for the potato cyst nematode, Patent application No. PCT/NL92/00126, 1996; Schenk et al. in Croat Chem Acta 72:593-606, 1999) showed that solanoeclepin A (SolA), a triterpenoid metabolite that was isolated from the root exudate of potato, induces the hatching of PCN. Its low concentration in potato root exudate has hindered progress in fully understanding its hatching inducing activity and exploitation in the control of PCN. To further investigate the role of SolA in hatching of PCN, the establishment of a highly sensitive analytical method is a prerequisite. Here we present the efficient single-step extraction and UHPLC-MS/MS based analysis for rapid determination of SolA in sub-nanomolar concentrations in tomato root exudate. This method was used to analyze SolA production in different tomato cultivars and related solanaceous species, including the trap crop Solanum sisymbriifolium. Hatching assays with PCN, Globodera pallida, with root exudates of tomato genotypes revealed a significant positive correlation between SolA concentration and hatching activity. Our results demonstrate that there is natural variation in SolA production within solanaceous species and that this has an effect on PCN hatching. The analytical method we have developed can potentially be used to support breeding for crop genotypes that induce less hatching and may therefore display reduced infection by PCN.

Plant triterpenoids with bond-missing skeletons: biogenesis, distribution and bioactivity.

Covering: up to December 2018 The polycyclic ABCD(E) framework of triterpenoids can miss a single endocyclic C-C bond as a result of a modification of the cyclization cascade that triggers their formation (interrupted- or diverted cascades), or can be the result of post-cyclization ring cleavage by late-stage oxidative modifications (seco-triterpenoids). Because of mechanistic and biogenetic differences, ring opening associated with loss of a skeletal fragment, as in nor-seco-triterpenoids (limonoids, quassinoids), will not be covered, nor will compounds where ring opening is part of a fragmentation cascade or of a multiple diversion from it. Even with these limitations, 342 bond-missing triterpenoids could be retrieved from the literature, with transversal distribution in the plant kingdom. Their structural diversity translates into a variety of biological targets, with dominance of potential applications in the realm of cancer, neuroprotection, and anti-infective therapy. In addition to the bioactivity and chemotaxonomic relevance of bond-missing triterpenoids, current knowledge on the genetic basis of interrupted- and diverted oxidosqualene cyclases will be summarized. This untapped source of enzymes could be useful to selectively modify triterpenoids by metabolic engineering, circumventing the bottlenecks of their isolation (poor yield or inadequate supply chain) to explore new areas of their chemical space.

Co-expression of squalene epoxidases with triterpene cyclases boosts production of triterpenoids in plants and yeast.

Triterpene cyclases catalyze the first committed step in triterpene biosynthesis, by forming mono- to pentacyclic backbone structures from oxygenated C30 isoprenoid precursors. Squalene epoxidase precedes this cyclization by providing the oxygenated and activated substrate for triterpene biosynthesis. Three squalene epoxidases from Cucurbita pepo (CpSEs) were isolated and shown to have evolved under purifying selection with signs of sites under positive selection in their N- and C-termini. They all localize to the Endoplasmic Reticulum (ER) and produce 2,3-oxidosqualene and 2,3:22,23-dioxidosqualene when expressed in a yeast erg1 (squalene epoxidase) erg7 (lanosterol synthase) double mutant. Co-expression of the CpSEs with four different triterpene cyclases, either transiently in Nicotiana benthamiana or constitutively in yeast, showed that CpSEs boost triterpene production. CpSE2 was the best performing in this regard, which could reflect either increased substrate production or superior channeling of the substrate to the triterpene cyclases. Fluorescence Lifetime Imaging Microscopy (FLIM) analysis with C. pepo cucurbitadienol synthase (CpCPQ) revealed a specific interaction with CpSE2 but not with the other CpSEs. When CpSE2 was transformed into C. pepo hairy root lines, cucurbitacin E production was increased two folds compared to empty vector control lines. This study provides new insight into the importance of SEs in triterpene biosynthesis, suggesting that they may facilitate substrate channeling, and demonstrates that SE overexpression is a new tool for increasing triterpene production in plants and yeast.

A Single Oxidosqualene Cyclase Produces the Seco-Triterpenoid α-Onocerin.

8,14-seco-Triterpenoids are characterized by their unusual open C-ring. Their distribution in nature is rare and scattered in taxonomically unrelated plants. The 8,14-seco-triterpenoid α-onocerin is only known from the evolutionarily distant clubmoss genus Lycopodium and the leguminous genus Ononis, which makes the biosynthesis of this seco-triterpenoid intriguing from an evolutionary standpoint. In our experiments with Ononis spinosa, α-onocerin was detected only in the roots. Through transcriptome analysis of the roots, an oxidosqualene cyclase, OsONS1, was identified that produces α-onocerin from squalene-2,3;22,23-dioxide when transiently expressed in Nicotiana bethamiana In contrast, in Lycopodium clavatum, two sequential cyclases, LcLCC and LcLCD, are required to produce α-onocerin in the N. benthamiana transient expression system. Expression of OsONS1 in the lanosterol synthase knockout yeast strain GIL77, which accumulates squalene-2,3;22,23-dioxide, verified the α-onocerin production. A phylogenetic analysis predicts that OsONS1 branches off from specific lupeol synthases and does not group with the known L. clavatum α-onocerin cyclases. Both the biochemical and phylogenetic analyses of OsONS1 suggest convergent evolution of the α-onocerin pathways. When OsONS1 was coexpressed in N. benthamiana leaves with either of the two O. spinosa squalene epoxidases, OsSQE1 or OsSQE2, α-onocerin production was boosted, most likely because the epoxidases produce higher amounts of squalene-2,3;22,23-dioxide. Fluorescence lifetime imaging microscopy analysis demonstrated specific protein-protein interactions between OsONS1 and both O. spinosa squalene epoxidases. Coexpression of OsONS1 with the two OsSQEs suggests that OsSQE2 is the preferred partner of OsONS1 in planta. Our results provide an example of the convergent evolution of plant specialized metabolism.

Monoterpene biosynthesis potential of plant subcellular compartments.

Subcellular monoterpene biosynthesis capacity based on local geranyl diphosphate (GDP) availability or locally boosted GDP production was determined for plastids, cytosol and mitochondria. A geraniol synthase (GES) was targeted to plastids, cytosol, or mitochondria. Transient expression in Nicotiana benthamiana indicated local GDP availability for each compartment but resulted in different product levels. A GDP synthase from Picea abies (PaGDPS1) was shown to boost GDP production. PaGDPS1 was also targeted to plastids, cytosol or mitochondria and PaGDPS1 and GES were coexpressed in all possible combinations. Geraniol and geraniol-derived products were analyzed by GC-MS and LC-MS, respectively. GES product levels were highest for plastid-targeted GES, followed by mitochondrial- and then cytosolic-targeted GES. For each compartment local boosting of GDP biosynthesis increased GES product levels. GDP exchange between compartments is not equal: while no GDP is exchanged from the cytosol to the plastids, 100% of GDP in mitochondria can be exchanged to plastids, while only 7% of GDP from plastids is available for mitochondria. This suggests a direct exchange mechanism for GDP between plastids and mitochondria. Cytosolic PaGDPS1 competes with plastidial GES activity, suggesting an effective drain of isopentenyl diphosphate from the plastids to the cytosol.

The seco-iridoid pathway from Catharanthus roseus.

The (seco)iridoids and their derivatives, the monoterpenoid indole alkaloids (MIAs), form two large families of plant-derived bioactive compounds with a wide spectrum of high-value pharmacological and insect-repellent activities. Vinblastine and vincristine, MIAs used as anticancer drugs, are produced by Catharanthus roseus in extremely low levels, leading to high market prices and poor availability. Their biotechnological production is hampered by the fragmentary knowledge of their biosynthesis. Here we report the discovery of the last four missing steps of the (seco)iridoid biosynthesis pathway. Expression of the eight genes encoding this pathway, together with two genes boosting precursor formation and two downstream alkaloid biosynthesis genes, in an alternative plant host, allows the heterologous production of the complex MIA strictosidine. This confirms the functionality of all enzymes of the pathway and highlights their utility for synthetic biology programmes towards a sustainable biotechnological production of valuable (seco)iridoids and alkaloids with pharmaceutical and agricultural applications.

Evaluation of tobacco (Nicotiana tabacum L. cv. Petit Havana SR1) hairy roots for the production of geraniol, the first committed step in terpenoid indole alkaloid pathway.

The terpenoid indole alkaloids are one of the major classes of plant-derived natural products and are well known for their many applications in the pharmaceutical, fragrance and cosmetics industries. Hairy root cultures are useful for the production of plant secondary metabolites because of their genetic and biochemical stability and their rapid growth in hormone-free media. Tobacco (Nicotiana tabacum L. cv. Petit Havana SR1) hairy roots, which do not produce geraniol naturally, were engineered to express a plastid-targeted geraniol synthase gene originally isolated from Valeriana officinalis L. (VoGES). A SPME-GC-MS screening tool was developed for the rapid evaluation of production clones. The GC-MS analysis revealed that the free geraniol content in 20 hairy root clones expressing VoGES was an average of 13.7 µg/g dry weight (DW) and a maximum of 31.3 µg/g DW. More detailed metabolic analysis revealed that geraniol derivatives were present in six major glycoside forms, namely the hexose and/or pentose conjugates of geraniol and hydroxygeraniol, resulting in total geraniol levels of up to 204.3 µg/g DW following deglycosylation. A benchtop-scale process was developed in a 20-L wave-mixed bioreactor eventually yielding hundreds of grams of biomass and milligram quantities of geraniol per cultivation bag.

Metabolic flux phenotype of tobacco hairy roots engineered for increased geraniol production.

The goal of this study was to characterise the metabolic flux phenotype of transgenic tobacco (Nicotiana tabacum) hairy roots engineered for increased biosynthesis of geraniol, an intermediate of the terpenoid indole alkaloid pathway. Steady state, stable isotope labelling was used to determine flux maps of central carbon metabolism for transgenic lines over-expressing (i) plastid-targeted geraniol synthase (pGES) from Valeriana officinalis, and (ii) pGES in combination with plastid-targeted geranyl pyrophosphate synthase from Arabidopsis thaliana (pGES+pGPPS), as well as for wild type and control-vector-transformed roots. Fluxes were constrained by the redistribution of label from [1-¹³C]-, [2-¹³C]- or [¹³C6]glucose into amino acids, sugars and organic acids at isotopic steady state, and by biomass output fluxes determined from the fractionation of [U-¹4C]glucose into insoluble polymers. No significant differences in growth and biomass composition were observed between the lines. The pGES line accumulated significant amounts of geraniol/geraniol glycosides (151±24 ng/mg dry weight) and the de novo synthesis of geraniol in pGES was confirmed by ¹³C labelling analysis. The pGES+pGPPS also accumulated geraniol and geraniol glycosides, but to lower levels than the pGES line. Although there was a distinct impact of the transgenes at the level of geraniol synthesis, other network fluxes were unaffected, reflecting the capacity of central metabolism to meet the relatively modest demand for increased precursors in the transgenic lines. It is concluded that re-engineering of the terpenoid indole alkaloid pathway will only require simultaneous manipulation of the steps producing the pathway precursors that originate in central metabolism in tissues engineered to produce at least an order of magnitude more geraniol than has been achieved so far.

Characterization of two geraniol synthases from Valeriana officinalis and Lippia dulcis: similar activity but difference in subcellular localization.

Two geraniol synthases (GES), from Valeriana officinalis (VoGES) and Lippia dulcis (LdGES), were isolated and were shown to have geraniol biosynthetic activity with Km values of 32 µM and 51 µM for GPP, respectively, upon expression in Escherichia coli. The in planta enzymatic activity and sub-cellular localization of VoGES and LdGES were characterized in stable transformed tobacco and using transient expression in Nicotiana benthamiana. Transgenic tobacco expressing VoGES or LdGES accumulate geraniol, oxidized geraniol compounds like geranial, geranic acid and hexose conjugates of these compounds to similar levels. Geraniol emission of leaves was lower than that of flowers, which could be related to higher levels of competing geraniol-conjugating activities in leaves. GFP-fusions of the two GES proteins show that VoGES resides (as expected) predominantly in the plastids, while LdGES import into to the plastid is clearly impaired compared to that of VoGES, resulting in both cytosolic and plastidic localization. Geraniol production by VoGES and LdGES in N. benthamiana was nonetheless very similar. Expression of a truncated version of VoGES or LdGES (cytosolic targeting) resulted in the accumulation of 30% less geraniol glycosides than with the plastid targeted VoGES and LdGES, suggesting that the substrate geranyl diphosphate is readily available, both in the plastids as well as in the cytosol. The potential role of GES in the engineering of the TIA pathway in heterologous hosts is discussed.

Geraniol hydroxylase and hydroxygeraniol oxidase activities of the CYP76 family of cytochrome P450 enzymes and potential for engineering the early steps of the (seco)iridoid pathway.

The geraniol-derived (seco)iridoid skeleton is a precursor for a large group of bioactive compounds with diverse therapeutic applications, including the widely used anticancer molecule vinblastine. Despite of this economic prospect, the pathway leading to iridoid biosynthesis from geraniol is still unclear. The first geraniol hydroxylation step has been reported to be catalyzed by cytochrome P450 enzymes such as CYP76B6 from Catharanthus roseus and CYP76C1 from Arabidopsis thaliana. In the present study, an extended functional analysis of CYP76 family members was carried-out to identify the most effective enzyme to be used for pathway reconstruction. This disproved CYP76C1 activity and led to the characterization of CYP76C4 from A. thaliana as a geraniol 9- or 8-hydroxylase. CYP76B6 emerged as a highly specialized multifunctional enzyme catalyzing two sequential oxidation steps leading to the formation of 8-oxogeraniol from geraniol. This dual function was confirmed in planta using a leaf-disc assay. The first step, geraniol hydroxylation, was very efficient and fast enough to outcompete geraniol conjugation in plant tissues. When the enzyme was expressed in leaf tissues, 8-oxogeraniol was converted into further oxidized and/or reduced compounds in the absence of the next enzyme of the iridoid pathway.