Identification of a second 16-hydroxytabersonine-O-methyltransferase suggests an evolutionary relationship between alkaloid and flavonoid metabolisms in Catharanthus roseus.

The medicinal plant Catharanthus roseus biosynthesizes many important drugs for human health, including the anticancer monoterpene indole alkaloids (MIAs) vinblastine and vincristine. Over the past decades, the continuous increase in pharmaceutical demand has prompted several research groups to characterize MIA biosynthetic pathways for considering future metabolic engineering processes of supply. In line with previous work suggesting that diversification can potentially occur at various steps along the vindoline branch, we were here interested in investigating the involvement of distinct isoforms of tabersonine-16-O-methyltransferase (16OMT) which plays a pivotal role in the MIA biosynthetic pathway. By combining homology searches based on the previously characterized 16OMT1, phylogenetic analyses, functional assays in yeast, and biochemical and in planta characterizations, we identified a second isoform of 16OMT, referred to as 16OMT2. 16OMT2 appears to be a multifunctional enzyme working on both MIA and flavonoid substrates, suggesting that a constrained evolution of the enzyme for accommodating the MIA substrate has probably occurred to favor the apparition of 16OMT2 from an ancestral specific flavonoid-O-methyltransferase. Since 16OMT1 and 16OMT2 displays a high sequence identity and similar kinetic parameters for 16-hydroxytabersonine, we postulate that 16OMT1 may result from a later 16OMT2 gene duplication accompanied by a continuous neofunctionalization leading to an almost complete loss of flavonoid O-methyltransferase activity. Overall, these results participate in increasing our knowledge on the evolutionary processes that have likely led to enzyme co-optation for MIA synthesis.

Optimization of Tabersonine Methoxylation to Increase Vindoline Precursor Synthesis in Yeast Cell Factories.

Plant specialized metabolites are widely used in the pharmaceutical industry, including the monoterpene indole alkaloids (MIAs) vinblastine and vincristine, which both display anticancer activity. Both compounds can be obtained through the chemical condensation of their precursors vindoline and catharanthine extracted from leaves of the Madagascar periwinkle. However, the extensive use of these molecules in chemotherapy increases precursor demand and results in recurrent shortages, explaining why the development of alternative production approaches, such microbial cell factories, is mandatory. In this context, the precursor-directed biosynthesis of vindoline from tabersonine in yeast-expressing heterologous biosynthetic genes is of particular interest but has not reached high production scales to date. To circumvent production bottlenecks, the metabolic flux was channeled towards the MIA of interest by modulating the copy number of the first two genes of the vindoline biosynthetic pathway, namely tabersonine 16-hydroxylase and tabersonine-16-O-methyltransferase. Increasing gene copies resulted in an optimized methoxylation of tabersonine and overcame the competition for tabersonine access with the third enzyme of the pathway, tabersonine 3-oxygenase, which exhibits a high substrate promiscuity. Through this approach, we successfully created a yeast strain that produces the fourth biosynthetic intermediate of vindoline without accumulation of other intermediates or undesired side-products. This optimization will probably pave the way towards the future development of yeast cell factories to produce vindoline at an industrial scale.

Enhanced bioproduction of anticancer precursor vindoline by yeast cell factories.

The pharmaceutical industry faces a growing demand and recurrent shortages in many anticancer plant drugs given their extensive use in human chemotherapy. Efficient alternative strategies of supply of these natural products such as bioproduction by microorganisms are needed to ensure stable and massive manufacturing. Here, we developed and optimized yeast cell factories efficiently converting tabersonine to vindoline, a precursor of the major anticancer alkaloids vinblastine and vincristine. First, fine-tuning of heterologous gene copies restrained side metabolites synthesis towards vindoline production. Tabersonine to vindoline bioconversion was further enhanced through a rational medium optimization (pH, composition) and a sequential feeding strategy. Finally, a vindoline titre of 266 mg l-1 (88% yield) was reached in an optimized fed-batch bioreactor. This precursor-directed synthesis of vindoline thus paves the way towards future industrial bioproduction through the valorization of abundant tabersonine resources.

Alternative splicing creates a pseudo-strictosidine ß-d-glucosidase modulating alkaloid synthesis in Catharanthus roseus.

Deglycosylation is a key step in the activation of specialized metabolites involved in plant defense mechanisms. This reaction is notably catalyzed by ß-glucosidases of the glycosyl hydrolase 1 (GH1) family such as strictosidine ß-d-glucosidase (SGD) from Catharanthus roseus. SGD catalyzes the deglycosylation of strictosidine, forming a highly reactive aglycone involved in the synthesis of cytotoxic monoterpene indole alkaloids (MIAs) and in the crosslinking of aggressor proteins. By exploring C. roseus transcriptomic resources, we identified an alternative splicing event of the SGD gene leading to the formation of a shorter isoform of this enzyme (shSGD) that lacks the last 71-residues and whose transcript ratio with SGD ranges from 1.7% up to 42.8%, depending on organs and conditions. Whereas it completely lacks ß-glucosidase activity, shSGD interacts with SGD and causes the disruption of SGD multimers. Such disorganization drastically inhibits SGD activity and impacts downstream MIA synthesis. In addition, shSGD disrupts the metabolic channeling of downstream biosynthetic steps by hampering the recruitment of tetrahydroalstonine synthase in cell nuclei. shSGD thus corresponds to a pseudo-enzyme acting as a regulator of MIA biosynthesis. These data shed light on a peculiar control mechanism of ß-glucosidase multimerization, an organization common to many defensive GH1 members.

Exploiting Spermidine N-Hydroxycinnamoyltransferase Diversity and Substrate Promiscuity to Produce Various Trihydroxycinnamoyl Spermidines and Analogues in Engineered Yeast.

Trihydroxycinnamoyl spermidines (THCSpd) are plant specialized metabolites with promising pharmacological activities as antifungals, antibacterial, antiviral, and antidepressant drugs. However, their characterization and potential pharmaceutical exploitation are greatly impaired by the sourcing of these compounds, restricted to the pollen of core Eudicot plant species. In this work, we developed a precursor-directed biosynthesis of THCSpd in yeast using a dual enzymatic system based on 4-coumarate-CoA ligases (4CL) and spermidine N-hydroxycinnamoyltransferases (SHT). The system relies on the yeast endogenous spermidine pool and only requires hydroxycinnamic acids as exogenous precursors. By exploring 4CL isoforms and SHT diversity among plants, we have driven the production of 8 natural THCSpd, using single or mixed hydroxycinnamic acid precursors. Substrate promiscuities of 4CL and SHT were genuinely exploited to produce 8 new-to-nature THCSpd from exotic hydroxycinnamic and dihydrohydroxycinnamic acids, together with 3 new-to-nature THCSpd containing halogenated hydroxycinnamoyl moieties. In this work, we established a versatile and modular biotechnological production platform allowing the tailor-made THCSpd synthesis, constituting pioneer metabolic engineering for access to these valuable natural products.

UPLC-HRMS Analysis Revealed the Differential Accumulation of Antioxidant and Anti-Aging Lignans and Neolignans in In Vitro Cultures of Linum usitatissimum L.

Over the last few decades, methods relating to plant tissue culture have become prevalent within the cosmetic industry. Forecasts predict the cosmetic industry to grow to an annual turnover of around a few hundred billion US dollars. Here we focused on Linum usitatissimum L., a plant that is well-known for its potent cosmetic properties. Following the a) establishment of cell cultures from three distinct initial explant origins (root, hypocotyl, and cotyledon) and b) selection of optimal hormonal concentrations, two in vitro systems (callus vs cell suspensions) were subjected to different light conditions. Phytochemical analysis by UPLC-HRMS not only confirmed high (neo)lignan accumulation capacity of this species with high concentrations of seven newly described (neo)lignans. Evaluation over 30 days revealed strong variations between the two different in vitro systems cultivated under light or dark, in terms of their growth kinetics and phytochemical composition. Additionally, antioxidant (i.e. four different in vitro assays based on hydrogen-atom transfer or electron transfer mechanism) and anti-aging (i.e. four in vitro inhibition potential of the skin remodeling enzymes: elastase, hyaluronidase, collagenase and tyrosinase) properties were evaluated for the two different in vitro systems cultivated under light or dark. A prominent hydrogen-atom transfer antioxidant mechanism was illustrated by the DPPH and ABTS assays. Potent tyrosinase and elastase inhibitory activities were also observed, which was strongly influenced by the in vitro system and light conditions. Statistical treatments of the data showed relationship of some (neo)lignans with these biological activities. These results confirmed the accumulation of flax (neo)lignans in different in vitro systems that were subjected to distinct light conditions. Furthermore, we showed the importance of optimizing these parameters for specific applications within the cosmetic industry.

A Biolistic-Mediated Virus-Induced Gene Silencing in Apocynaceae to Map Biosynthetic Pathways of Alkaloids.

Monoterpene indole alkaloids (MIAs) are specialized metabolites synthesized in many plants of the Apocynaceae family including Catharanthus roseus and Rauvolfia sp. MIAs are part of the chemical arsenal that plants evolved to face pet and herbivore attacks, and their high biological activities also confer pharmaceutical properties exploited in human pharmacopeia. Developing robust and straightforward tools to elucidate each step of MIA biosynthetic pathways thus constitutes a prerequisite to the understanding of Apocynaceae defense mechanisms and to the exploitation of MIA cytotoxicity through their production by metabolic engineering. While protocols of virus-induced gene silencing (VIGS) based on Agrobacterium-based transformation have emerged, the recalcitrance of Apocynaceae to this type of transformation prompted us to develop an universal procedure of VIGS vector inoculation. Such procedure relies on the delivery of the transforming plasmids through a particle bombardment performed using a biolistic device and offers the possibility to overcome host specificity to silence genes in any plant species. Using silencing of geissoschizine oxidase as an example, we described the main steps of this biolistic mediated VIGS in C. roseus and R. tetraphylla.

Identifying Missing Biosynthesis Enzymes of Plant Natural Products.

Elucidating plant-specialized biosynthetic pathways has always constituted a laborious task, notably for natural products with high pharmaceutical values. Here, we discuss emerging omics-based strategies that facilitate the identification of genes from these complex metabolic pathways, paving the way to engineered supplies of these compounds through synthetic biology approaches.

Improved gene co-expression network quality through expression dataset down-sampling and network aggregation.

Large-scale gene co-expression networks are an effective methodology to analyze sets of co-expressed genes and discover new gene functions or associations. Distances between genes are estimated according to their expression profiles and are visualized in networks that may be further partitioned to reveal communities of co-expressed genes. Creating expression profiles is now eased by the large amounts of publicly available expression data (microarrays and RNA-seq). Although many distance calculation methods have been intensively compared and reviewed in the past, it is unclear how to proceed when many samples reflecting a wide range of different conditions are available. Should as many samples as possible be integrated into network construction or be partitioned into smaller sets of more related samples? Previous studies have indicated a saturation in network performances to capture known associations once a certain number of samples is included in distance calculations. Here, we examined the influence of sample size on co-expression network construction using microarray and RNA-seq expression data from three plant species. We tested different down-sampling methods and compared network performances in recovering known gene associations to networks obtained from full datasets. We further examined how aggregating networks may help increase this performance by testing six aggregation methods.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes.

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

Genome-wide identification and biochemical characterization of the UGT88F subfamily in Malus x domestica Borkh.

The UDP-glycosyltransferase UGT88F subfamily has been described first in Malus x domestica with the characterization of UGT88F1. Up to now UGT88F1 was one of the most active UGT glycosylating dihydrochalcones in vitro. The involvement of UGT88F1 in phloridzin (phloretin 2'-O-glucoside) synthesis, the main apple tree dihydrochalcone, was further confirmed in planta. Since the characterization of UGT88F1, this new UGT subfamily has been poorly studied probably because it seemed restricted to Maloideae. In the present study, we investigate the apple tree genome to identify and biochemically characterize the whole UGT88F subfamily. The apple tree genome contains five full-length UGT88F genes out of which three newly identified members (UGT88F6, UGT88F7 and UGT88F8) and a pseudogene. These genes are organized into two genomic clusters resulting from the recent global genomic duplication event in the apple tree. We show that recombinant UGT88F8 protein specifically glycosylates phloretin in the 2'OH position to synthetize phloridzin in vitro and was therefore named UDP-glucose: phloretin 2'-O-glycosyltransferase. The Km values of UGT88F8 are 7.72 µM and 10.84 µM for phloretin and UDP-glucose respectively and are in the same range as UGT88F1 catalytic parameters thus constituting two isoforms. Co-expression patterns of both UGT88F1 and UGT88F8 argue for a redundant function in phloridzin biosynthesis in planta. Contrastingly, recombinant UGT88F6 protein is able to glycosylate in vitro a wide range of flavonoids including flavonols, flavones, flavanones, chalcones and dihydrochalcones, although flavonols are the preferred substrates, e.g. Km value for kaempferol is 2.1 µM. Depending on the flavonoid, glycosylation occurs at least on the 3-OH and 7-OH positions. Therefore UGT88F6 corresponds to an UDP-glucose: flavonoid 3/7-O-glycosyltransferase. Finally, a molecular modeling study highlights a very high substitution rate of residues in the acceptor binding pocket between UGT88F8 and UGT88F6 which is responsible for the enzymes divergence in substrate and regiospecificity, despite an overall high protein homology.

Setting-up a fast and reliable cytokinin biosensor based on a plant histidine kinase receptor expressed in Saccharomyces cerevisiae.

Cytokinins (CK) have been extensively studied for their roles in plant development. Recently, they also appeared to ensure crucial functions in the pathogenicity of some bacterial and fungal plant pathogens. Thus, identifying cytokinin-producing pathogens is a prerequisite to gain a better understanding of their role in pathogenicity. Taking advantage of the cytokinin perception properties of Malus domestica CHASE Histidine Kinase receptor 2 (MdCHK2), we thereby developed a selective and highly sensitive yeast biosensor for the application of cytokinin detection in bacterial samples. The biosensor is based on the mutated sln1Δ Saccharomyces cerevisiae strain expressing MdCHK2. The biosensor does not require any extraction or purification steps of biological samples, enabling cytokinin analysis directly from crude bacterial supernatants. For the first time, the production of cytokinin was shown in the well-known plant pathogenic bacteria Erwinia amylovora and was also revealed in human pathogens Staphylococcus aureus and Streptococcus agalactiae. Importantly, this biosensor was shown to be an efficient tool for unraveling certain steps in cytokinin biosynthesis by micro-organisms since this it was successfully used to unveil the role of ygdH22, a LOG-like gene, that is probably involved in cytokinin biosynthesis pathway in Escherichia coli. Overall, we demonstrated that our biosensor displays several advantages including time- and cost-effectiveness by allowing a rapid and specific detection of cytokinins in bacterial supernatants These results also support its scalability to high-throughput formats.

Field-Based Metabolomics of Vitis vinifera L. Stems Provides New Insights for Genotype Discrimination and Polyphenol Metabolism Structuring.

Grape accumulates numerous polyphenols with abundant health benefit and organoleptic properties that in planta act as key components of the plant defense system against diseases. Considerable advances have been made in the chemical characterization of wine metabolites particularly volatile and polyphenolic compounds. However, the metabotyping (metabolite-phenotype characterization) of grape varieties, from polyphenolic-rich vineyard by-product is unprecedented. As this composition might result from the complex interaction between genotype, environment and viticultural practices, a field experiment was setting up with uniform pedo-climatic factors and viticultural practices of growing vines to favor the genetic determinism of polyphenol expression. As a result, UPLC-MS-based targeted metabolomic analyses of grape stems from 8 Vitis vinifera L. cultivars allowed the determination of 42 polyphenols related to phenolic acids, flavonoids, procyanidins, and stilbenoids as resveratrol oligomers (degree of oligomerization 1-4). Using a partial least-square discriminant analysis approach, grape stem chemical profiles were discriminated according to their genotypic origin showing that polyphenol profile express a varietal signature. Furthermore, hierarchical clustering highlights various degree of polyphenol similarity between grape varieties that were in agreement with the genetic distance using clustering analyses of 22 microsatellite DNA markers. Metabolite correlation network suggested that several polyphenol subclasses were differently controlled. The present polyphenol metabotyping approach coupled to multivariate statistical analyses might assist grape selection programs to improve metabolites with both health-benefit potential and plant defense traits.

Ranking genome-wide correlation measurements improves microarray and RNA-seq based global and targeted co-expression networks.

Co-expression networks are essential tools to infer biological associations between gene products and predict gene annotation. Global networks can be analyzed at the transcriptome-wide scale or after querying them with a set of guide genes to capture the transcriptional landscape of a given pathway in a process named Pathway Level Coexpression (PLC). A critical step in network construction remains the definition of gene co-expression. In the present work, we compared how Pearson Correlation Coefficient (PCC), Spearman Correlation Coefficient (SCC), their respective ranked values (Highest Reciprocal Rank (HRR)), Mutual Information (MI) and Partial Correlations (PC) performed on global networks and PLCs. This evaluation was conducted on the model plant Arabidopsis thaliana using microarray and differently pre-processed RNA-seq datasets. We particularly evaluated how dataset × distance measurement combinations performed in 5 PLCs corresponding to 4 well described plant metabolic pathways (phenylpropanoid, carbohydrate, fatty acid and terpene metabolisms) and the cytokinin signaling pathway. Our present work highlights how PCC ranked with HRR is better suited for global network construction and PLC with microarray and RNA-seq data than other distance methods, especially to cluster genes in partitions similar to biological subpathways.

Two Tabersonine 6,7-Epoxidases Initiate Lochnericine-Derived Alkaloid Biosynthesis in Catharanthus roseus.

Lochnericine is a major monoterpene indole alkaloid (MIA) in the roots of Madagascar periwinkle (Catharanthus roseus). Lochnericine is derived from the stereoselective C6,C7-epoxidation of tabersonine and can be metabolized further to generate other complex MIAs. While the enzymes responsible for its downstream modifications have been characterized, those involved in lochnericine biosynthesis remain unknown. By combining gene correlation studies, functional assays, and transient gene inactivation, we identified two highly conserved P450s that efficiently catalyze the epoxidation of tabersonine: tabersonine 6,7-epoxidase isoforms 1 and 2 (TEX1 and TEX2). Both proteins are quite divergent from the previously characterized tabersonine 2,3-epoxidase and are more closely related to tabersonine 16-hydroxylase, involved in vindoline biosynthesis in leaves. Biochemical characterization of TEX1/2 revealed their strict substrate specificity for tabersonine and their inability to epoxidize 19-hydroxytabersonine, indicating that they catalyze the first step in the pathway leading to hörhammericine production. TEX1 and TEX2 displayed complementary expression profiles, with TEX1 expressed mainly in roots and TEX2 in aerial organs. Our results suggest that TEX1 and TEX2 originated from a gene duplication event and later acquired divergent, organ-specific regulatory elements for lochnericine biosynthesis throughout the plant, as supported by the presence of lochnericine in flowers. Finally, through the sequential expression of TEX1 and up to four other MIA biosynthetic genes in yeast, we reconstituted the 19-acetylhörhammericine biosynthetic pathway and produced tailor-made MIAs by mixing enzymatic modules that are naturally spatially separated in the plant. These results lay the groundwork for the metabolic engineering of tabersonine/lochnericine derivatives of pharmaceutical interest.

Vacuole-Targeted Proteins: Ins and Outs of Subcellular Localization Studies.

Accurate and efficient demonstrations of protein localizations to the vacuole or tonoplast remain strict prerequisites to decipher the role of vacuoles in the whole plant cell biology and notably in defence processes. In this chapter, we describe a reliable procedure of protein subcellular localization study through transient transformations of Catharanthus roseus or onion cells and expression of fusions with fluorescent proteins allowing minimizing artefacts of targeting.

A BAHD acyltransferase catalyzing 19-O-acetylation of tabersonine derivatives in roots of Catharanthus roseus enables combinatorial synthesis of monoterpene indole alkaloids.

While the characterization of the biosynthetic pathway of monoterpene indole alkaloids (MIAs) in leaves of Catharanthus roseus is now reaching completion, only two enzymes from the root counterpart dedicated to tabersonine metabolism have been identified to date, namely tabersonine 19-hydroxylase (T19H) and minovincine 19-O-acetyltransferase (MAT). Albeit the recombinant MAT catalyzes MIA acetylation at low efficiency in vitro, we demonstrated that MAT was inactive when expressed in yeast and in planta, suggesting an alternative function for this enzyme. Therefore, through transcriptomic analysis of periwinkle adventitious roots, several other BAHD acyltransferase candidates were identified based on the correlation of their expression profile with T19H and found to localize in small genomic clusters. Only one, named tabersonine derivative 19-O-acetyltransferase (TAT) was able to acetylate the 19-hydroxytabersonine derivatives from roots, such as minovincinine and hörhammericine, following expression in yeast. Kinetic studies also showed that the recombinant TAT was specific for root MIAs and displayed an up to 200-fold higher catalytic efficiency than MAT. In addition, gene expression analysis, protein subcellular localization and heterologous expression in Nicotiana benthamiana were in agreement with the prominent role of TAT in acetylation of root-specific MIAs, thereby redefining the molecular determinants of the root MIA biosynthetic pathway. Finally, identification of TAT provided a convenient tool for metabolic engineering of MIAs in yeast enabling efficiently mixing different biosynthetic modules spatially separated in the whole plant. This combinatorial synthesis associating several enzymes from Catharanthus roseus resulted in the conversion of tabersonine in tailor-made MIAs bearing both leaf and root-type decorations.

Mechanical stress rapidly induces E-resveratrol and E-piceatannol biosynthesis in grape canes stored as a freshly-pruned byproduct.

Grape canes represent a promising source of bioactive phytochemicals. However the stabilization of the raw material after pruning remains challenging. We recently reported the induction of stilbenoid metabolism after winter pruning including a strong accumulation of E-resveratrol and E-piceatannol during the first six weeks of storage. In the present study, the effect of mechanical wounding on freshly-pruned canes was tested to increase the induction of stilbenoid metabolism. Cutting the grape canes in short segments immediately after pruning triggered a transient expression of phenylalanine ammonia-lyase (PAL) and stilbene synthase (STS) genes, followed by a rapid accumulation of E-resveratrol and E-piceatannol. The degree of stilbenoid induction was related to the intensity of mechanical wounding. Data suggest that a global defense response is triggered involving jasmonate signaling, PR proteins and stilbenoid metabolism. Mechanical wounding of freshly-pruned canes drastically shortens the time required to reach maximal stilbenoid accumulation from 6 to 2weeks.