Phytoalexins of the crucifer Barbarea vulgaris: Structural profile and correlation with glucosinolate turnover.

Phytoalexins are antimicrobial plant metabolites elicited by microbial attack or abiotic stress. We investigated phytoalexin profiles after foliar abiotic elicitation in the crucifer Barbarea vulgaris and interactions with the glucosinolate-myrosinase system. The treatment for abiotic elicitation was a foliar spray with CuCl2 solution, a usual eliciting agent, and three independent experiments were carried out. Two genotypes of B. vulgaris (G-type and P-type) accumulated the same three major phytoalexins in rosette leaves after treatment: phenyl-containing nasturlexin D and indole-containing cyclonasturlexin and cyclobrassinin. Phytoalexin levels were investigated daily by UHPLC-QToF MS and tended to differ among plant types and individual phytoalexins. In roots, phytoalexins were low or not detected. In treated leaves, typical total phytoalexin levels were in the range 1-10 nmol/g fresh wt. during three days after treatment while typical total glucosinolate (GSL) levels were three orders of magnitude higher. Levels of some minor GSLs responded to the treatment: phenethylGSL (PE) and 4-substituted indole GSLs. Levels of PE, a suggested nasturlexin D precursor, were lower in treated plants than controls. Another suggested precursor GSL, 3-hydroxyPE, was not detected, suggesting PE hydrolysis to be a key biosynthetic step. Levels of 4-substituted indole GSLs differed markedly between treated and control plants in most experiments, but not in a consistent way. The dominant GSLs, glucobarbarins, are not believed to be phytoalexin precursors. We observed statistically significant linear correlations between total major phytoalexins and the glucobarbarin products barbarin and resedine, suggesting that GSL turnover for phytoalexin biosynthesis was unspecific. In contrast, we did not find correlations between total major phytoalexins and raphanusamic acid or total glucobarbarins and barbarin. In conclusion, two groups of phytoalexins were detected in B. vulgaris, apparently derived from the GSLs PE and indol-3-ylmethylGSL. Phytoalexin biosynthesis was accompanied by depletion of the precursor PE and by turnover of major non-precursor GSLs to resedine. This work paves the way for identifying and characterizing genes and enzymes in the biosyntheses of phytoalexins and resedine.

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.

The saponin bomb: a nucleolar-localized ß-glucosidase hydrolyzes triterpene saponins in Medicago truncatula.

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a ß-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned ß-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Isolate-Specific Effect of Entomopathogenic Endophytic Fungi on Population Growth of Two-Spotted Spider Mite (Tetranychus urticae Koch) and Levels of Steroidal Glycoalkaloids in Tomato.

Entomopathogenic fungi (EPF) can be experimentally established in several plant species as endophytes. Ecological effects of EPF inoculations on plant growth and plant-herbivore interactions have been demonstrated, potentially by altering plant physiological responses. However, the role of these responses in plant-fungus-herbivore tripartite interactions has not been well elucidated. Steroidal glycoalkaloids (SGAs) are plant specialized metabolites with bioactive properties against arthropod herbivores. Here, the effects of seed treatments by three EPF isolates, representing Beauveria bassiana, Metarhizium brunneum, and M. robertsii, on population growth of two-spotted spider mites (Tetranychus urticae Koch) were evaluated on tomato (Solanum lycopersicum). The levels of two SGAs, α-tomatine and dehydrotomatine, were determined in tomato leaves by LC-MS with and without T. urticae infestations after EPF inoculations. Interestingly, the population growth of T. urticae was significantly highest with M. brunneum and lowest with M. robertsii and B. bassiana at 15 days after infestation. Overall there was a significant negative correlation between SGAs content and the number of T. urticae. The levels of SGAs were significantly induced by T. urticae presence in all treatments, while only M. robertsii showed significantly higher levels of SGAs than M. brunneum and control in one of two experiments. Contrastingly, the effects on SGAs accumulation and population growth of T. urticae did not directly correlate with EPF endophytic colonization patterns of the inoculated plants. This study suggests a link between ecological effects and physiological responses mediated by EPF inoculations and T. urticae infestation with potential implications for plant protection.

Analysis of wild tomato introgression lines elucidates the genetic basis of transcriptome and metabolome variation underlying fruit traits and pathogen response.

Wild tomato species represent a rich gene pool for numerous desirable traits lost during domestication. Here, we exploited an introgression population representing wild desert-adapted species and a domesticated cultivar to establish the genetic basis of gene expression and chemical variation accompanying the transfer of wild-species-associated fruit traits. Transcriptome and metabolome analysis of 580 lines coupled to pathogen sensitivity assays resulted in the identification of genomic loci associated with levels of hundreds of transcripts and metabolites. These associations occurred in hotspots representing coordinated perturbation of metabolic pathways and ripening-related processes. Here, we identify components of the Solanum alkaloid pathway, as well as genes and metabolites involved in pathogen defense and linking fungal resistance with changes in the fruit ripening regulatory network. Our results outline a framework for understanding metabolism and pathogen resistance during tomato fruit ripening and provide insights into key fruit quality traits.

Pathways to defense metabolites and evading fruit bitterness in genus Solanum evolved through 2-oxoglutarate-dependent dioxygenases.

The genus Solanum comprises three food crops (potato, tomato, and eggplant), which are consumed on daily basis worldwide and also producers of notorious anti-nutritional steroidal glycoalkaloids (SGAs). Hydroxylated SGAs (i.e. leptinines) serve as precursors for leptines that act as defenses against Colorado Potato Beetle (Leptinotarsa decemlineata Say), an important pest of potato worldwide. However, SGA hydroxylating enzymes remain unknown. Here, we discover that 2-OXOGLUTARATE-DEPENDENT-DIOXYGENASE (2-ODD) enzymes catalyze SGA-hydroxylation across various Solanum species. In contrast to cultivated potato, Solanum chacoense, a widespread wild potato species, has evolved a 2-ODD enzyme leading to the formation of leptinines. Furthermore, we find a related 2-ODD in tomato that catalyzes the hydroxylation of the bitter α-tomatine to hydroxytomatine, the first committed step in the chemical shift towards downstream ripening-associated non-bitter SGAs (e.g. esculeoside A). This 2-ODD enzyme prevents bitterness in ripe tomato fruit consumed today which otherwise would remain unpleasant in taste and more toxic.

The cytochrome P450 CYP72A552 is key to production of hederagenin-based saponins that mediate plant defense against herbivores.

Plants continuously evolve new defense compounds. One class of such compounds is triterpenoid saponins. A few species in the Barbarea genus produce saponins as the only ones in the large crucifer family. However, the molecular mechanism behind saponin biosynthesis and their role in plant defense remains unclear. We used pathway reconstitution in planta, enzymatic production of saponins in vitro, insect feeding assays, and bioinformatics to identify a missing gene involved in saponin biosynthesis and saponin-based herbivore defense. A tandem repeat of eight CYP72A cytochromes P450 colocalise with a quantitative trait locus (QTL) for saponin accumulation and flea beetle resistance in Barbarea vulgaris. We found that CYP72A552 oxidises oleanolic acid at position C-23 to hederagenin. In vitro-produced hederagenin monoglucosides reduced larval feeding by up to 90% and caused 75% larval mortality of the major crucifer pest diamondback moth and the tobacco hornworm. Sequence analysis indicated that CYP72A552 evolved through gene duplication and has been under strong selection pressure. In conclusion, CYP72A552 has evolved to catalyse the formation of hederagenin-based saponins that mediate plant defense against herbivores. Our study highlights the evolution of chemical novelties by gene duplication and selection for enzyme innovations, and the importance of chemical modification in plant defense evolution.

Evolution of Structural Diversity of Triterpenoids.

Plants have evolved to produce a blend of specialized metabolites that serve functional roles in plant adaptation. Among them, triterpenoids are one of the largest subclasses of such specialized metabolites, with more than 14,000 known structures. They play a role in plant defense and development and have potential applications within food and pharma. Triterpenoids are cyclized from oxidized squalene precursors by oxidosqualene cyclases, creating more than 100 different cyclical triterpene scaffolds. This limited number of scaffolds is the first step towards creating the vast structural diversity of triterpenoids followed by extensive diversification, in particular, by oxygenation and glycosylation. Gene duplication, divergence, and selection are major forces that drive triterpenoid structural diversification. The triterpenoid biosynthetic genes can be organized in non-homologous gene clusters, such as in Avena spp., Cucurbitaceae and Solanum spp., or scattered along plant chromosomes as in Barbarea vulgaris. Paralogous genes organized as tandem repeats reflect the extended gene duplication activities in the evolutionary history of the triterpenoid saponin pathways, as seen in B. vulgaris. We review and discuss examples of convergent and divergent evolution in triterpenoid biosynthesis, and the apparent mechanisms occurring in plants that drive their increasing structural diversity within and across species. Using B. vulgaris' saponins as examples, we discuss the impact a single structural modification can have on the structure of a triterpenoid and how this affect its biological properties. These examples provide insight into how plants continuously evolve their specialized metabolome, opening the way to study uncharacterized triterpenoid biosynthetic pathways.

Corrigendum: Plant cholesterol biosynthetic pathway overlaps with phytosterol metabolism.

This corrects the article DOI: 10.1038/nplants.2016.205.

Plant cholesterol biosynthetic pathway overlaps with phytosterol metabolism.

The amount of cholesterol made by many plants is not negligible. Whereas cholesterogenesis in animals was elucidated decades ago, the plant pathway has remained enigmatic. Among other roles, cholesterol is a key precursor for thousands of bioactive plant metabolites, including the well-known Solanum steroidal glycoalkaloids. Integrating tomato transcript and protein co-expression data revealed candidate genes putatively associated with cholesterol biosynthesis. A combination of functional assays including gene silencing, examination of recombinant enzyme activity and yeast mutant complementation suggests the cholesterol pathway comprises 12 enzymes acting in 10 steps. It appears that half of the cholesterogenesis-specific enzymes evolved through gene duplication and divergence from phytosterol biosynthetic enzymes, whereas others act reciprocally in both cholesterol and phytosterol metabolism. Our findings provide a unique example of nature's capacity to exploit existing protein folds and catalytic machineries from primary metabolism to assemble a new, multi-step metabolic pathway. Finally, the engineering of a 'high-cholesterol' model plant underscores the future value of our gene toolbox to produce high-value steroidal compounds via synthetic biology.

GAME9 regulates the biosynthesis of steroidal alkaloids and upstream isoprenoids in the plant mevalonate pathway.

Steroidal glycoalkaloids (SGAs) are cholesterol-derived molecules produced by solanaceous species. They contribute to pathogen defence but are toxic to humans and considered as anti-nutritional compounds. Here we show that GLYCOALKALOID METABOLISM 9 (GAME9), an APETALA2/Ethylene Response Factor, related to regulators of alkaloid production in tobacco and Catharanthus roseus, controls SGA biosynthesis. GAME9 knockdown and overexpression in tomato and potato alters expression of SGAs and upstream mevalonate pathway genes including the cholesterol biosynthesis gene STEROL SIDE CHAIN REDUCTASE 2 (SSR2). Levels of SGAs, C24-alkylsterols and the upstream mevalonate and cholesterol pathways intermediates are modified in these plants. Δ(7)-STEROL-C5(6)-DESATURASE (C5-SD) in the hitherto unresolved cholesterol pathway is a direct target of GAME9. Transactivation and promoter-binding assays show that GAME9 exerts its activity either directly or cooperatively with the SlMYC2 transcription factor as in the case of the C5-SD gene promoter. Our findings provide insight into the regulation of SGA biosynthesis and means for manipulating these metabolites in crops.

Retention of triplicated phytoene synthase (PSY) genes in Brassica napus L. and its diploid progenitors during the evolution of the Brassiceae.

The extent of genome redundancy exhibited by Brassica species provides a model to study the evolutionary fate of multi-copy genes and the effects of polyploidy in economically important crops. Phytoene synthase (PSY) catalyzes the first committed reaction of the carotenoid biosynthetic pathway, which has been shown to be rate-limiting in Brassica napus seeds. In Arabidopsis thaliana, a single PSY gene (AtPSY) regulates phytoene synthesis in all tissues. Considering that diploid Brassica genomes contain three Arabidopsis-like subgenomes, the objectives of the present work were to determine whether PSY gene families exist in B. napus (AACC) and its diploid progenitor species, Brassica rapa (AA) and Brassica oleracea (CC); to establish the level of retention of Brassica PSY genes; to map PSY gene family members in the A and C genomes and to compare Brassica PSY gene expression patterns. A total of 12 PSY homologues were identified, 6 in B. napus (BnaX.PSY.a-f) and 3 in B. rapa (BraA.PSY.a-c) and B. oleracea (BolC.PSY.a-c). Indeed, with six members, B. napus has the largest PSY gene family described to date. Sequence comparison between AtPSY and Brassica PSY genes revealed a highly conserved gene structure and identity percentages above 85% at the coding sequence (CDS) level. Altogether, our data indicate that PSY gene family expansion preceded the speciation of B. rapa and B. oleracea, dating back to the paralogous subgenome triplication event. In these three Brassica species, all PSY homologues are expressed, exhibiting overlapping redundancy and signs of subfunctionalization among photosynthetic and non-photosynthetic tissues. This evidence supports the hypothesis that functional divergence of PSY gene expression facilitates the accumulation of high levels of carotenoids in chromoplast-rich tissues. Thus, functional retention of triplicated Brassica PSY genes could be at least partially explained by the selective advantage provided by increased levels of gene product in floral organs. A better understanding of carotenogenesis in Brassica will aid in the future development of transgenic and conventional cultivars with carotenoid-enriched oil.